1 |
|
|
!OPTIONS XOPT(NOEVAL) |
2 |
|
1 |
SUBROUTINE SUPHY1(KULOUT) |
3 |
|
|
|
4 |
|
|
!**** *SUPHY1* - Initialize common YOMPHY1 physics land surface |
5 |
|
|
! constants |
6 |
|
|
|
7 |
|
|
! Purpose. |
8 |
|
|
! -------- |
9 |
|
|
! Initialize YOMPHY1, the common that contains the parameters |
10 |
|
|
! for the land surface part of the physics of the model. |
11 |
|
|
|
12 |
|
|
!** Interface. |
13 |
|
|
! ---------- |
14 |
|
|
! *CALL* *SUPHY1(KULOUT) |
15 |
|
|
|
16 |
|
|
! Explicit arguments : |
17 |
|
|
! -------------------- |
18 |
|
|
! KULOUT : Logical unit for the output |
19 |
|
|
|
20 |
|
|
! Implicit arguments : |
21 |
|
|
! -------------------- |
22 |
|
|
! COMMON YOMPHY1 |
23 |
|
|
|
24 |
|
|
! Method. |
25 |
|
|
! ------- |
26 |
|
|
! See documentation |
27 |
|
|
|
28 |
|
|
! Externals. |
29 |
|
|
! ---------- |
30 |
|
|
|
31 |
|
|
! Reference. |
32 |
|
|
! ---------- |
33 |
|
|
! Documentation ARPEGE |
34 |
|
|
|
35 |
|
|
! Author. |
36 |
|
|
! ------- |
37 |
|
|
! J.-F. Geleyn . |
38 |
|
|
|
39 |
|
|
! Modifications. |
40 |
|
|
! -------------- |
41 |
|
|
! Original : 90-9-1 |
42 |
|
|
! Modified 91-02-28 by Michel Deque (Relaxation of deep soil values) |
43 |
|
|
! Modified by R. EL Khatib : 93-04-02 Set-up defaults controled by LECMWF |
44 |
|
|
! Modified by M. Deque : 94-10-19 4-layer soil temperature |
45 |
|
|
! Modified by D. Giard : 94-10-07 Constants for ACSOL |
46 |
|
|
! Control keys : LIMC, LIMW |
47 |
|
|
! Modified by H. Douville : 95-01-13 Snow parameterization |
48 |
|
|
! Modified by D. Giard : 95-03-09 Loop for SODELX modified |
49 |
|
|
! 95-09-08 Constants for LSNV |
50 |
|
|
! Modified by P. Mercier : 97-03-24 Ozone difusion + deposition |
51 |
|
|
! Modified by E. Bazile : 97-05-05 C1 option vapour phase |
52 |
|
|
! Modified by J.M. Piriou : 97-02-26 soil inertia. |
53 |
|
|
! Modified by M. Deque : 97-04-10 Sea-ice parameters |
54 |
|
|
! Modified by J.M. Piriou : 97-04-17 soil inertia default values. |
55 |
|
|
! Modified by D. Giard : 97-11-13 defaults for veget. features |
56 |
|
|
! Modified by E. Bazile : 97-12-08 Soil freezing. |
57 |
|
|
! Modified by E. Bazile : 00-12-12 Default value for soil freezing. |
58 |
|
|
! Modified by E. Bazile : 02-10-29 Defaults for the snow scheme LVGSN. |
59 |
|
|
! Modified by E. Bazile : 04-02-24 Introduce NCHSP. |
60 |
|
|
! M.Hamrud 01-Oct-2003 CY28 Cleaning |
61 |
|
|
! ------------------------------------------------------------------ |
62 |
|
|
|
63 |
|
|
USE PARKIND1 ,ONLY : JPIM ,JPRB |
64 |
|
|
USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
65 |
|
|
|
66 |
|
|
! Ce qui concerne NULNAM commente par MPL le 15.04.09 |
67 |
|
|
!USE YOMLUN , ONLY : NULNAM |
68 |
|
|
USE YOMCST , ONLY : RPI ,RDAY ,RG ,RCS |
69 |
|
|
USE SURFACE_FIELDS , ONLY : YSD_VVD, YSP_SBD |
70 |
|
|
USE YOMCT0B , ONLY : LECMWF |
71 |
|
|
USE YOMPHY1 , ONLY : GF3 ,GF4 ,TREF4 ,RCTVEG ,& |
72 |
|
|
& RGL ,SODELX ,GCZ0H ,ALBGLA ,ALBMAX ,& |
73 |
|
|
& ALBMER ,ALBMED ,ALBMIN ,ALCRIN ,ALRCN1 ,ALRCN2 ,& |
74 |
|
|
& EA ,EC2REF ,EMCRIN ,EMMGLA ,EMMMER ,& |
75 |
|
|
& EWFC ,EWWILT ,GA ,GC1 ,GC1S1 ,& |
76 |
|
|
& GC1S2 ,GC1S3 ,GC1S4 ,GC1Y1 ,GTSVAP ,& |
77 |
|
|
& GVEGMX ,GLAIMX ,GNEIMX ,GWPIMX ,GCGEL ,& |
78 |
|
|
& GC2 ,GC2REF ,GC3 ,GC31 ,GC32 ,& |
79 |
|
|
& GCONV ,GF1 ,GWFC ,GWLEX ,GWLMX ,& |
80 |
|
|
& GWWILT ,G1B ,G1CGSAT ,G1C1SAT ,G1P ,& |
81 |
|
|
& G1WSAT ,G2B ,G2CGSAT ,G2C1SAT ,G2P ,& |
82 |
|
|
& G2WSAT ,G3CGSAT ,GSNC1 ,GSNC2 ,HSOL ,& |
83 |
|
|
& HSOLIWR ,HSOLIT0 ,OMTPRO ,OMWPRO ,RC1MAX ,& |
84 |
|
|
& RCTGLA ,RCGMAX ,RD1 ,RD2GLA ,RD2MER ,& |
85 |
|
|
& RHOMAX ,RHOMIN ,RSMAX ,RTINER ,RZ0GLA ,& |
86 |
|
|
& RZ0MER ,RZHZ0G ,RZHZ0M ,RZHGLA ,RZHMER ,& |
87 |
|
|
& TMERGL ,TOEXP ,& |
88 |
|
|
& TOLIN ,WCRIN ,WCRINC ,WCRING ,WNEW ,& |
89 |
|
|
& WPMX ,WSMX ,XCRINR ,XCRINV ,LIMC ,& |
90 |
|
|
& LIMW ,LC1VAP ,NTVGLA ,NTVMER ,GCGELS ,& |
91 |
|
|
& GVEGMXS ,GLAIMXS ,GNEIMXS ,ALB1 ,ALB2 ,& |
92 |
|
|
& RLAIMX ,RLAI ,NCHSP |
93 |
|
|
USE YOMVDOZ , ONLY : VDHJS ,VDHJH ,VDHNS ,VDHNH ,& |
94 |
|
|
& VDPJS ,VDPJH ,VDPNS ,VDPNH ,VDEJS ,& |
95 |
|
|
& VDEJH ,VDENS ,VDENH ,VDAJS ,VDAJH ,& |
96 |
|
|
& VDANS ,VDANH ,VDNJS ,VDNJH ,VDNNS ,& |
97 |
|
|
& VDNNH ,VOZNJ ,VOZHS ,LRDIFOZ ,LRDEPOZ |
98 |
|
|
|
99 |
|
|
IMPLICIT NONE |
100 |
|
|
|
101 |
|
|
INTEGER(KIND=JPIM),INTENT(IN) :: KULOUT |
102 |
|
|
INTEGER(KIND=JPIM) :: IVEG, J |
103 |
|
|
REAL(KIND=JPRB) :: ZHOOK_HANDLE |
104 |
|
|
#include "abor1.intfb.h" |
105 |
|
|
#include "posnam.intfb.h" |
106 |
|
|
|
107 |
|
|
#include "namphy1.h" |
108 |
|
|
#include "namvdoz.h" |
109 |
|
|
! ------------------------------------------------------------------ |
110 |
|
|
|
111 |
|
|
!* 1. Set default values. |
112 |
|
|
! ------------------- |
113 |
|
|
|
114 |
|
|
! 1.1 Set implicit default values |
115 |
|
|
|
116 |
|
|
! Snow , Sea-Ice |
117 |
✓✗ |
1 |
IF (LHOOK) CALL DR_HOOK('SUPHY1',0,ZHOOK_HANDLE) |
118 |
|
1 |
ALCRIN=0.7_JPRB |
119 |
|
1 |
ALRCN1=1.E-03_JPRB |
120 |
|
1 |
ALRCN2=2.5E-03_JPRB |
121 |
|
1 |
EMCRIN=0.98_JPRB |
122 |
|
1 |
WCRIN=10._JPRB |
123 |
|
1 |
TMERGL=271.23_JPRB |
124 |
|
1 |
NCHSP=0 |
125 |
|
|
! Characteristics of ice and sea |
126 |
|
1 |
NTVGLA=2 |
127 |
|
1 |
NTVMER=1 |
128 |
|
1 |
RD2GLA=8._JPRB |
129 |
|
1 |
RD2MER=8._JPRB |
130 |
|
1 |
ALBGLA=.65_JPRB |
131 |
|
1 |
ALBMER=.07_JPRB |
132 |
|
1 |
ALBMED=.06_JPRB |
133 |
|
1 |
EMMGLA=.97_JPRB |
134 |
|
1 |
EMMMER=.96_JPRB |
135 |
|
1 |
RZ0GLA=.001_JPRB |
136 |
|
1 |
RZ0MER=.001_JPRB |
137 |
|
1 |
RZHGLA=.001_JPRB |
138 |
|
1 |
RZHMER=.001_JPRB |
139 |
|
1 |
RZHZ0G=1.0_JPRB |
140 |
|
1 |
RZHZ0M=1.0_JPRB |
141 |
|
|
! Usual - soil |
142 |
|
1 |
HSOL=1.0E-05_JPRB |
143 |
|
1 |
HSOLIT0=0.35_JPRB |
144 |
|
1 |
HSOLIWR=6._JPRB |
145 |
|
1 |
RTINER=5._JPRB |
146 |
|
1 |
WPMX=100._JPRB |
147 |
|
1 |
WSMX=20._JPRB |
148 |
|
1 |
SODELX(0)=1.0_JPRB/SQRT(1.0_JPRB+2.0_JPRB*RPI) |
149 |
✓✓ |
10 |
DO J=1,9 |
150 |
|
10 |
SODELX(J)=SODELX(J-1)*2.0_JPRB*RPI |
151 |
|
|
ENDDO |
152 |
|
|
! Relaxation |
153 |
|
1 |
OMTPRO=0._JPRB |
154 |
|
1 |
OMWPRO=0._JPRB |
155 |
|
|
! ISBA - soil |
156 |
|
1 |
EA=-0.54_JPRB |
157 |
|
1 |
GA=732.42E-3_JPRB |
158 |
|
1 |
G1B=0.137_JPRB |
159 |
|
1 |
G2B=3.5_JPRB |
160 |
|
1 |
G1P=0.134_JPRB |
161 |
|
1 |
G2P=3.4_JPRB |
162 |
|
1 |
GC1=0.5_JPRB |
163 |
|
1 |
GC2=10._JPRB |
164 |
|
1 |
GC3=8._JPRB |
165 |
|
1 |
GCONV=1.E3_JPRB |
166 |
|
|
|
167 |
|
1 |
G1WSAT=-1.08E-3_JPRB |
168 |
|
1 |
G2WSAT=494.31E-3_JPRB |
169 |
|
1 |
EWFC=0.35_JPRB |
170 |
|
1 |
GWFC=89.0467E-3_JPRB |
171 |
|
1 |
EWWILT=0.5_JPRB |
172 |
|
1 |
GWWILT=37.1342E-3_JPRB |
173 |
|
1 |
EC2REF=-0.95_JPRB |
174 |
|
1 |
GC2REF=13.82_JPRB |
175 |
|
1 |
G1CGSAT=-1.5571E-8_JPRB |
176 |
|
1 |
G2CGSAT=-1.441E-8_JPRB |
177 |
|
1 |
G3CGSAT=4.70217E-6_JPRB |
178 |
|
1 |
G1C1SAT=5.58E-3_JPRB |
179 |
|
1 |
G2C1SAT=84.88E-3_JPRB |
180 |
|
1 |
GC31=5.3275_JPRB |
181 |
|
1 |
GC32=-1.043_JPRB |
182 |
|
|
|
183 |
|
1 |
RD1=1.E-2_JPRB |
184 |
|
1 |
RC1MAX=500._JPRB |
185 |
|
1 |
RCTGLA=5.5E-6_JPRB |
186 |
|
1 |
RCGMAX=0.8E-5_JPRB |
187 |
|
1 |
LIMC=.TRUE. |
188 |
|
1 |
LIMW=.TRUE. |
189 |
|
|
! ISBA - vegetation |
190 |
|
1 |
GF1=0.55_JPRB |
191 |
|
1 |
GWLEX=2.0_JPRB/3._JPRB |
192 |
|
1 |
GWLMX=0.2_JPRB |
193 |
|
1 |
RSMAX=5000._JPRB |
194 |
✓✓ |
19 |
DO J=1,18 |
195 |
|
18 |
GF3(J)=0.0_JPRB |
196 |
|
18 |
GF4(J)=0.0016_JPRB |
197 |
|
18 |
RCTVEG(J)=0.8E-5_JPRB |
198 |
|
18 |
RGL(J)=100._JPRB |
199 |
|
19 |
TREF4(J)=298._JPRB |
200 |
|
|
ENDDO |
201 |
|
1 |
GF3(4)=40._JPRB |
202 |
|
1 |
RGL(4)=30._JPRB |
203 |
|
|
! ISBA - roughness length |
204 |
|
1 |
GCZ0H(0,1)=7.5_JPRB |
205 |
|
1 |
GCZ0H(1,1)=2.39037_JPRB |
206 |
|
1 |
GCZ0H(2,1)=-.28583_JPRB |
207 |
|
1 |
GCZ0H(3,1)=.01074_JPRB |
208 |
|
1 |
GCZ0H(0,2)=0.5_JPRB |
209 |
|
1 |
GCZ0H(1,2)=-.07028_JPRB |
210 |
|
1 |
GCZ0H(2,2)=.01023_JPRB |
211 |
|
1 |
GCZ0H(3,2)=-.00067_JPRB |
212 |
|
1 |
GCZ0H(0,3)=5.0_JPRB |
213 |
|
1 |
GCZ0H(1,3)=4.51268_JPRB |
214 |
|
1 |
GCZ0H(2,3)=.34012_JPRB |
215 |
|
1 |
GCZ0H(3,3)=-.05330_JPRB |
216 |
|
1 |
GCZ0H(0,4)=0.5_JPRB |
217 |
|
1 |
GCZ0H(1,4)=-.09421_JPRB |
218 |
|
1 |
GCZ0H(2,4)=.01463_JPRB |
219 |
|
1 |
GCZ0H(3,4)=-.00099_JPRB |
220 |
|
|
! ISBA - snow |
221 |
|
1 |
ALBMAX=0.85_JPRB |
222 |
|
1 |
ALBMIN=0.50_JPRB |
223 |
|
1 |
RHOMAX=0.3_JPRB |
224 |
|
1 |
RHOMIN=0.1_JPRB |
225 |
|
1 |
TOEXP=0.24_JPRB/86400._JPRB |
226 |
|
1 |
TOLIN=0.008_JPRB/86400._JPRB |
227 |
|
1 |
WCRINC=70._JPRB |
228 |
|
1 |
WCRING=10._JPRB |
229 |
|
1 |
WNEW=10._JPRB |
230 |
|
1 |
XCRINR=1.0_JPRB/RG |
231 |
|
1 |
XCRINV=10000._JPRB |
232 |
|
1 |
GSNC1=RPI/(2.22_JPRB*RCS*RDAY*1000._JPRB) |
233 |
|
1 |
GSNC2=2.885_JPRB |
234 |
|
|
! ISBA - C1 vapour phase |
235 |
|
1 |
LC1VAP=.TRUE. |
236 |
|
1 |
GTSVAP=0._JPRB |
237 |
|
1 |
GC1S1= 1.19_JPRB |
238 |
|
1 |
GC1S2=-5.09_JPRB |
239 |
|
1 |
GC1S3=-1.464E+2_JPRB |
240 |
|
1 |
GC1S4= 17.86E+2_JPRB |
241 |
|
1 |
GC1Y1=10._JPRB |
242 |
|
|
! Deep Soil freezing |
243 |
|
1 |
GVEGMX=5._JPRB |
244 |
|
1 |
GLAIMX=30._JPRB |
245 |
|
1 |
GNEIMX=1.8_JPRB |
246 |
|
1 |
GWPIMX=150._JPRB |
247 |
|
1 |
GCGEL=3.E-5_JPRB |
248 |
|
|
! Surface soil freezing |
249 |
|
1 |
GCGELS=5.E-5_JPRB |
250 |
|
1 |
GVEGMXS=5._JPRB |
251 |
|
1 |
GLAIMXS=30._JPRB |
252 |
|
1 |
GNEIMXS=1.8_JPRB |
253 |
|
|
! OZONE DIFFUSION AND DEPOSITION |
254 |
|
1 |
LRDIFOZ=.FALSE. |
255 |
|
1 |
LRDEPOZ=.FALSE. |
256 |
|
1 |
VOZNJ=1._JPRB |
257 |
|
1 |
VOZHS=1.0_JPRB/86400._JPRB |
258 |
✓✓ |
100 |
DO J=1,99 |
259 |
|
99 |
VDHJS(J)=-999._JPRB |
260 |
|
99 |
VDHJH(J)=-999._JPRB |
261 |
|
99 |
VDHNS(J)=-999._JPRB |
262 |
|
99 |
VDHNH(J)=-999._JPRB |
263 |
|
99 |
VDPJS(J)=-999._JPRB |
264 |
|
99 |
VDPJH(J)=-999._JPRB |
265 |
|
99 |
VDPNS(J)=-999._JPRB |
266 |
|
99 |
VDPNH(J)=-999._JPRB |
267 |
|
99 |
VDEJS(J)=-999._JPRB |
268 |
|
99 |
VDEJH(J)=-999._JPRB |
269 |
|
99 |
VDENS(J)=-999._JPRB |
270 |
|
99 |
VDENH(J)=-999._JPRB |
271 |
|
99 |
VDAJS(J)=-999._JPRB |
272 |
|
99 |
VDAJH(J)=-999._JPRB |
273 |
|
99 |
VDANS(J)=-999._JPRB |
274 |
|
99 |
VDANH(J)=-999._JPRB |
275 |
|
99 |
VDNJS(J)=-999._JPRB |
276 |
|
99 |
VDNJH(J)=-999._JPRB |
277 |
|
99 |
VDNNS(J)=-999._JPRB |
278 |
|
100 |
VDNNH(J)=-999._JPRB |
279 |
|
|
ENDDO |
280 |
|
|
! New snow scheme (LVGSN) |
281 |
|
1 |
ALB1=0.87_JPRB |
282 |
|
1 |
ALB2=0.84_JPRB |
283 |
|
1 |
RLAIMX=7._JPRB |
284 |
|
1 |
RLAI=3._JPRB |
285 |
|
|
|
286 |
|
|
! 1.2 Modify default values according to LECMWF |
287 |
|
|
|
288 |
|
|
IF (LECMWF) THEN |
289 |
|
|
ELSE |
290 |
|
|
ENDIF |
291 |
|
|
|
292 |
|
|
! ------------------------------------------------------------------ |
293 |
|
|
|
294 |
|
|
!* 2. Modify default values. |
295 |
|
|
! ---------------------- |
296 |
|
|
|
297 |
|
|
! Ce qui concerne NAMPHY1 et NAMVDOZ commente par MPL le 15.04.09 |
298 |
|
|
!CALL POSNAM(NULNAM,'NAMPHY1') |
299 |
|
|
!READ(NULNAM,NAMPHY1) |
300 |
|
|
!CALL POSNAM(NULNAM,'NAMVDOZ') |
301 |
|
|
!READ(NULNAM,NAMVDOZ) |
302 |
|
|
|
303 |
|
|
! 2.5 Check consistency |
304 |
|
|
! ----------------- |
305 |
|
|
IF (GC1Y1 > 60._JPRB) THEN |
306 |
|
|
CALL ABOR1 ('GC1Y1 FOR C1-VAPOUR PHASE IS BIGGER THAN 60.') |
307 |
|
|
ENDIF |
308 |
|
|
|
309 |
|
|
IF ((YSD_VVD%NUMFLDS < 8).AND.(RZHZ0M /= 1.0_JPRB)) THEN |
310 |
|
|
CALL ABOR1('YSD_VVD%NUMFLDS<8 IMPLIES RZHZ0M=1.0_JPRB !...') |
311 |
|
|
ENDIF |
312 |
|
|
|
313 |
|
|
! ------------------------------------------------------------------ |
314 |
|
|
|
315 |
|
|
!* 3. Print final values. |
316 |
|
|
! ------------------- |
317 |
|
|
|
318 |
|
1 |
WRITE(UNIT=KULOUT,FMT='('' COMMON YOMPHY1 '')') |
319 |
|
|
|
320 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
321 |
|
|
& '' ALCRIN ='',E10.4,'' ALRCN1 ='',E10.4,'' ALRCN2 ='',E10.4 & |
322 |
|
|
& ,'' EMCRIN ='',E10.4,'' WCRIN ='',E10.4,'' TMERGL ='',E10.4)')& |
323 |
|
1 |
& ALCRIN,ALRCN1,ALRCN2,EMCRIN,WCRIN,TMERGL |
324 |
|
|
|
325 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
326 |
|
|
& '' NCHSP ='',I3,& |
327 |
|
|
& '' NTVMER ='',I3,'' NTVGLA ='',I3 & |
328 |
|
|
& ,'' RD2MER ='',E10.4,'' ALBMER='',E10.4,'' EMMMER ='',E10.4,/& |
329 |
|
|
& ,'' RD2GLA ='',E10.4,'' ALBGLA='',E10.4,'' EMMGLA ='',E10.4,/& |
330 |
|
|
& ,'' RZ0MER ='',E10.4,'' RZHMER='',E10.4 & |
331 |
|
|
& ,'' RZ0GLA ='',E10.4,'' RZHGLA='',E10.4)')& |
332 |
|
1 |
& NCHSP,NTVMER,NTVGLA,RD2MER,ALBMER,EMMMER & |
333 |
|
2 |
& ,RD2GLA,ALBGLA,EMMGLA,RZ0MER,RZHMER,RZ0GLA,RZHGLA |
334 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
335 |
|
|
& '' RTINER ='',E10.4,'' HSOL ='',E10.4,'' HSOLIT0 ='',E10.4 & |
336 |
|
|
& ,'' HSOLIWR ='',E10.4,'' WPMX ='',E10.4 & |
337 |
|
|
& ,'' WSMX ='',E10.4,'' OMTPRO ='',E10.4,'' OMWPRO ='',E10.4)')& |
338 |
|
1 |
& RTINER,HSOL,HSOLIT0,HSOLIWR,WPMX,WSMX,OMTPRO,OMWPRO |
339 |
|
|
|
340 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
341 |
|
|
& '' EA ='',E10.4,'' GA ='',E10.4,'' G1B ='',E10.4 & |
342 |
|
|
& ,'' G2B ='',E10.4,'' G1P ='',E10.4,'' G2P = '',E10.4,/& |
343 |
|
|
& ,'' GC1 ='',E10.4,'' GC2 ='',E10.4,'' GC3 ='',E10.4 & |
344 |
|
|
& ,'' GCONV ='',E10.4)')& |
345 |
|
1 |
& EA,GA,G1B,G2B,G1P,G2P,GC1,GC2,GC3,GCONV |
346 |
|
|
|
347 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
348 |
|
|
& '' G1WSAT ='',E10.4,'' G2WSAT ='',E10.4,'' EWFC ='',E10.4 & |
349 |
|
|
& ,'' GWFC ='',E10.4,/,'' EWWILT ='',E10.4,'' GWWILT ='',E10.4 & |
350 |
|
|
& ,'' EC2REF ='',E10.4,'' GC2REF ='',E10.4,/& |
351 |
|
|
& ,'' G1CGSAT ='',E10.4,'' G2CGSAT ='',E10.4,'' G3CGSAT =''& |
352 |
|
|
& ,E10.4,'' G1C1SAT ='',E10.4,'' G2C1SAT ='',E10.4)')& |
353 |
|
1 |
& G1WSAT,G2WSAT,EWFC,GWFC,EWWILT,GWWILT,EC2REF,GC2REF,& |
354 |
|
2 |
& G1CGSAT,G2CGSAT,G3CGSAT,G1C1SAT,G2C1SAT |
355 |
|
|
|
356 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
357 |
|
|
& '' RD1 ='',E10.4,'' RC1MAX ='',E10.4,'' RCTGLA ='',E10.4 & |
358 |
|
|
& ,'' RCGMAX ='',E10.4,'' LIMC ='',L2,'' LIMW ='',L2)')& |
359 |
|
1 |
& RD1,RC1MAX,RCTGLA,RCGMAX,LIMC,LIMW |
360 |
|
|
|
361 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
362 |
|
|
& '' GC1S1 ='',E10.4,'' GC1S2 ='',E10.4,'' GC1S3 ='',E10.4 & |
363 |
|
|
& ,'' GC1S4 ='',E10.4,'' GC1Y1 ='',E10.4,'' LC1VAP ='',L2 & |
364 |
|
|
& ,'' GTSVAP ='',E10.4)')& |
365 |
|
1 |
& GC1S1,GC1S2,GC1S3,GC1S4,GC1Y1,LC1VAP,GTSVAP |
366 |
|
|
|
367 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
368 |
|
|
& '' GCGEL ='',E10.4,'' GVEGMX ='',E10.4,'' GLAIMX ='',E10.4 & |
369 |
|
|
& ,'' GWPIMX ='',E10.4,'' GNEIMX ='',E10.4)')& |
370 |
|
1 |
& GCGEL,GVEGMX,GLAIMX,GWPIMX,GNEIMX |
371 |
|
|
|
372 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
373 |
|
|
& '' GCGELS ='',E10.4,'' GVEGMXS ='',E10.4 & |
374 |
|
|
& ,'' GLAIMXS ='',E10.4,'' GNEIMXS ='',E10.4)')& |
375 |
|
1 |
& GCGELS,GVEGMXS,GLAIMXS,GNEIMXS |
376 |
|
|
|
377 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
378 |
|
|
& '' ALB1 ='',E10.4,'' ALB2 ='',E10.4 & |
379 |
|
|
& ,'' RLAIMX ='',E10.4,'' RLAI ='',E10.4)')& |
380 |
|
1 |
& ALB1,ALB2,RLAIMX,RLAI |
381 |
|
|
|
382 |
|
1 |
WRITE(UNIT=KULOUT,FMT='('' GCZ0H ='',/,4(1X,4E11.4,/))')GCZ0H |
383 |
|
|
|
384 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
385 |
|
|
& '' GF1 ='',E10.4,'' GWLEX ='',E10.4,'' GWLMX ='',E10.4 & |
386 |
|
|
& ,'' RSMAX ='',E10.4,/& |
387 |
|
|
& ,'' GF3 ='',/,2(1X,9E11.4,/),'' GF4 ='',/,2(1X,9E11.4,/)& |
388 |
|
|
& ,'' RCTVEG ='',/,2(1X,9E11.4,/),'' RGL ='',/,2(1X,9E11.4,/)& |
389 |
|
|
& ,'' TREF4 ='',/,2(1X,9E11.4,/))')& |
390 |
|
1 |
& GF1,GWLEX,GWLMX,RSMAX,GF3,GF4,RCTVEG,RGL,TREF4 |
391 |
|
|
|
392 |
|
|
WRITE(UNIT=KULOUT,FMT='(& |
393 |
|
|
& '' ALBMAX = '',E10.4,'' ALBMIN = '',E10.4 & |
394 |
|
|
& ,'' RHOMAX = '',E10.4,'' RHOMIN = '',E10.4,/& |
395 |
|
|
& ,'' TOEXP = '',E10.4,'' TOLIN = '',E10.4 & |
396 |
|
|
& ,'' WCRINC = '',E10.4,'' WCRING = '',E10.4,/& |
397 |
|
|
& ,'' WNEW = '',E10.4 & |
398 |
|
|
& ,'' XCRINR = '',E10.4,'' XCRINV = '',E10.4)')& |
399 |
|
1 |
& ALBMAX,ALBMIN,RHOMAX,RHOMIN,TOEXP,TOLIN,WCRINC,WCRING,WNEW,& |
400 |
|
2 |
& XCRINR,XCRINV |
401 |
|
|
|
402 |
|
1 |
WRITE(UNIT=KULOUT,FMT='('' SODELX = ''/5E11.4/5E11.4)') SODELX |
403 |
✗✓ |
1 |
IF(YSP_SBD%NLEVS > 9) CALL ABOR1(' TOO MANY SOIL LAYERS !') |
404 |
|
|
|
405 |
|
1 |
WRITE(UNIT=KULOUT,FMT='('' COMMON YOMVDOZ '')') |
406 |
|
|
|
407 |
|
1 |
WRITE(UNIT=KULOUT,FMT='('' LRDIFOZ ='',L2,'' LRDEPOZ ='',L2)')LRDIFOZ,LRDEPOZ |
408 |
|
|
|
409 |
✗✓✗✗
|
1 |
IF (LRDIFOZ.AND.LRDEPOZ) THEN |
410 |
|
|
WRITE(UNIT=KULOUT,FMT='('' VOZNJ ='',F10.5,'' VOZHS ='',F10.5)')VOZNJ,VOZHS |
411 |
|
|
DO IVEG=1,99 |
412 |
|
|
IF((VDHJS(IVEG) /= -999._JPRB).OR.(VDHJH(IVEG) /= -999._JPRB).OR.& |
413 |
|
|
& (VDHNS(IVEG) /= -999._JPRB).OR.(VDHNH(IVEG) /= -999._JPRB).OR.& |
414 |
|
|
& (VDPJS(IVEG) /= -999._JPRB).OR.(VDPJH(IVEG) /= -999._JPRB).OR.& |
415 |
|
|
& (VDPNS(IVEG) /= -999._JPRB).OR.(VDPNH(IVEG) /= -999._JPRB).OR.& |
416 |
|
|
& (VDEJS(IVEG) /= -999._JPRB).OR.(VDEJH(IVEG) /= -999._JPRB).OR.& |
417 |
|
|
& (VDENS(IVEG) /= -999._JPRB).OR.(VDENH(IVEG) /= -999._JPRB).OR.& |
418 |
|
|
& (VDAJS(IVEG) /= -999._JPRB).OR.(VDAJH(IVEG) /= -999._JPRB).OR.& |
419 |
|
|
& (VDANS(IVEG) /= -999._JPRB).OR.(VDANH(IVEG) /= -999._JPRB).OR.& |
420 |
|
|
& (VDNJS(IVEG) /= -999._JPRB).OR.(VDNJH(IVEG) /= -999._JPRB).OR.& |
421 |
|
|
& (VDNNS(IVEG) /= -999._JPRB).OR.(VDNNH(IVEG) /= -999._JPRB))THEN |
422 |
|
|
WRITE(UNIT=KULOUT,FMT='('' IVEJ ='',I2 & |
423 |
|
|
& ,'' VDHJS ='',E9.3,'' VDHJH ='',E9.3,'' VDHNS ='',E9.3 & |
424 |
|
|
& ,'' VDHNH ='',E9.3)')& |
425 |
|
|
& IVEG,VDHJS(IVEG),VDHJH(IVEG),VDHNS(IVEG),VDHNH(IVEG) |
426 |
|
|
WRITE(UNIT=KULOUT,FMT='('' IVEJ ='',I2 & |
427 |
|
|
& ,'' VDPJS ='',E9.3,'' VDPJH ='',E9.3,'' VDPNS ='',E9.3 & |
428 |
|
|
& ,'' VDPNH ='',E9.3)')& |
429 |
|
|
& IVEG,VDPJS(IVEG),VDPJH(IVEG),VDPNS(IVEG),VDPNH(IVEG) |
430 |
|
|
WRITE(UNIT=KULOUT,FMT='('' IVEJ ='',I2 & |
431 |
|
|
& ,'' VDEJS ='',E9.3,'' VDEJH ='',E9.3,'' VDENS ='',E9.3 & |
432 |
|
|
& ,'' VDENH ='',E9.3)')& |
433 |
|
|
& IVEG,VDEJS(IVEG),VDEJH(IVEG),VDENS(IVEG),VDENH(IVEG) |
434 |
|
|
WRITE(UNIT=KULOUT,FMT='('' IVEJ ='',I2 & |
435 |
|
|
& ,'' VDAJS ='',E9.3,'' VDAJH ='',E9.3,'' VDANS ='',E9.3 & |
436 |
|
|
& ,'' VDANH ='',E9.3)')& |
437 |
|
|
& IVEG,VDAJS(IVEG),VDAJH(IVEG),VDANS(IVEG),VDANH(IVEG) |
438 |
|
|
WRITE(UNIT=KULOUT,FMT='('' IVEJ ='',I2 & |
439 |
|
|
& ,'' VDNJS ='',E9.3,'' VDNJH ='',E9.3,'' VDNNS ='',E9.3 & |
440 |
|
|
& ,'' VDNNH ='',E9.3)')& |
441 |
|
|
& IVEG,VDNJS(IVEG),VDNJH(IVEG),VDNNS(IVEG),VDNNH(IVEG) |
442 |
|
|
ENDIF |
443 |
|
|
ENDDO |
444 |
|
|
ENDIF |
445 |
|
|
!* |
446 |
|
|
! ------------------------------------------------------------------ |
447 |
|
|
|
448 |
✓✗ |
1 |
IF (LHOOK) CALL DR_HOOK('SUPHY1',1,ZHOOK_HANDLE) |
449 |
|
1 |
END SUBROUTINE SUPHY1 |