LMDZ
suphy1.F90
Go to the documentation of this file.
1 !OPTIONS XOPT(NOEVAL)
2 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 IF (lhook) CALL dr_hook('SUPHY1',0,zhook_handle)
118 alcrin=0.7_jprb
119 alrcn1=1.e-03_jprb
120 alrcn2=2.5e-03_jprb
121 emcrin=0.98_jprb
122 wcrin=10._jprb
123 tmergl=271.23_jprb
124 nchsp=0
125 ! Characteristics of ice and sea
126 ntvgla=2
127 ntvmer=1
128 rd2gla=8._jprb
129 rd2mer=8._jprb
130 albgla=.65_jprb
131 albmer=.07_jprb
132 albmed=.06_jprb
133 emmgla=.97_jprb
134 emmmer=.96_jprb
135 rz0gla=.001_jprb
136 rz0mer=.001_jprb
137 rzhgla=.001_jprb
138 rzhmer=.001_jprb
139 rzhz0g=1.0_jprb
140 rzhz0m=1.0_jprb
141 ! Usual - soil
142 hsol=1.0e-05_jprb
143 hsolit0=0.35_jprb
144 hsoliwr=6._jprb
145 rtiner=5._jprb
146 wpmx=100._jprb
147 wsmx=20._jprb
148 sodelx(0)=1.0_jprb/sqrt(1.0_jprb+2.0_jprb*rpi)
149 DO j=1,9
150  sodelx(j)=sodelx(j-1)*2.0_jprb*rpi
151 ENDDO
152 ! Relaxation
153 omtpro=0._jprb
154 omwpro=0._jprb
155 ! ISBA - soil
156 ea=-0.54_jprb
157 ga=732.42e-3_jprb
158 g1b=0.137_jprb
159 g2b=3.5_jprb
160 g1p=0.134_jprb
161 g2p=3.4_jprb
162 gc1=0.5_jprb
163 gc2=10._jprb
164 gc3=8._jprb
165 gconv=1.e3_jprb
166 
167 g1wsat=-1.08e-3_jprb
168 g2wsat=494.31e-3_jprb
169 ewfc=0.35_jprb
170 gwfc=89.0467e-3_jprb
171 ewwilt=0.5_jprb
172 gwwilt=37.1342e-3_jprb
173 ec2ref=-0.95_jprb
174 gc2ref=13.82_jprb
175 g1cgsat=-1.5571e-8_jprb
176 g2cgsat=-1.441e-8_jprb
177 g3cgsat=4.70217e-6_jprb
178 g1c1sat=5.58e-3_jprb
179 g2c1sat=84.88e-3_jprb
180 gc31=5.3275_jprb
181 gc32=-1.043_jprb
182 
183 rd1=1.e-2_jprb
184 rc1max=500._jprb
185 rctgla=5.5e-6_jprb
186 rcgmax=0.8e-5_jprb
187 limc=.true.
188 limw=.true.
189 ! ISBA - vegetation
190 gf1=0.55_jprb
191 gwlex=2.0_jprb/3._jprb
192 gwlmx=0.2_jprb
193 rsmax=5000._jprb
194 DO j=1,18
195  gf3(j)=0.0_jprb
196  gf4(j)=0.0016_jprb
197  rctveg(j)=0.8e-5_jprb
198  rgl(j)=100._jprb
199  tref4(j)=298._jprb
200 ENDDO
201 gf3(4)=40._jprb
202 rgl(4)=30._jprb
203 ! ISBA - roughness length
204 gcz0h(0,1)=7.5_jprb
205 gcz0h(1,1)=2.39037_jprb
206 gcz0h(2,1)=-.28583_jprb
207 gcz0h(3,1)=.01074_jprb
208 gcz0h(0,2)=0.5_jprb
209 gcz0h(1,2)=-.07028_jprb
210 gcz0h(2,2)=.01023_jprb
211 gcz0h(3,2)=-.00067_jprb
212 gcz0h(0,3)=5.0_jprb
213 gcz0h(1,3)=4.51268_jprb
214 gcz0h(2,3)=.34012_jprb
215 gcz0h(3,3)=-.05330_jprb
216 gcz0h(0,4)=0.5_jprb
217 gcz0h(1,4)=-.09421_jprb
218 gcz0h(2,4)=.01463_jprb
219 gcz0h(3,4)=-.00099_jprb
220 ! ISBA - snow
221 albmax=0.85_jprb
222 albmin=0.50_jprb
223 rhomax=0.3_jprb
224 rhomin=0.1_jprb
225 toexp=0.24_jprb/86400._jprb
226 tolin=0.008_jprb/86400._jprb
227 wcrinc=70._jprb
228 wcring=10._jprb
229 wnew=10._jprb
230 xcrinr=1.0_jprb/rg
231 xcrinv=10000._jprb
232 gsnc1=rpi/(2.22_jprb*rcs*rday*1000._jprb)
233 gsnc2=2.885_jprb
234 ! ISBA - C1 vapour phase
235 lc1vap=.true.
236 gtsvap=0._jprb
237 gc1s1= 1.19_jprb
238 gc1s2=-5.09_jprb
239 gc1s3=-1.464e+2_jprb
240 gc1s4= 17.86e+2_jprb
241 gc1y1=10._jprb
242 ! Deep Soil freezing
243 gvegmx=5._jprb
244 glaimx=30._jprb
245 gneimx=1.8_jprb
246 gwpimx=150._jprb
247 gcgel=3.e-5_jprb
248 ! Surface soil freezing
249 gcgels=5.e-5_jprb
250 gvegmxs=5._jprb
251 glaimxs=30._jprb
252 gneimxs=1.8_jprb
253 ! OZONE DIFFUSION AND DEPOSITION
254 lrdifoz=.false.
255 lrdepoz=.false.
256 voznj=1._jprb
257 vozhs=1.0_jprb/86400._jprb
258 DO j=1,99
259  vdhjs(j)=-999._jprb
260  vdhjh(j)=-999._jprb
261  vdhns(j)=-999._jprb
262  vdhnh(j)=-999._jprb
263  vdpjs(j)=-999._jprb
264  vdpjh(j)=-999._jprb
265  vdpns(j)=-999._jprb
266  vdpnh(j)=-999._jprb
267  vdejs(j)=-999._jprb
268  vdejh(j)=-999._jprb
269  vdens(j)=-999._jprb
270  vdenh(j)=-999._jprb
271  vdajs(j)=-999._jprb
272  vdajh(j)=-999._jprb
273  vdans(j)=-999._jprb
274  vdanh(j)=-999._jprb
275  vdnjs(j)=-999._jprb
276  vdnjh(j)=-999._jprb
277  vdnns(j)=-999._jprb
278  vdnnh(j)=-999._jprb
279 ENDDO
280 ! New snow scheme (LVGSN)
281 alb1=0.87_jprb
282 alb2=0.84_jprb
283 rlaimx=7._jprb
284 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 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  & 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  & nchsp,ntvmer,ntvgla,rd2mer,albmer,emmmer &
333  & ,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  & 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  & 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  & g1wsat,g2wsat,ewfc,gwfc,ewwilt,gwwilt,ec2ref,gc2ref,&
354  & 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  & 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  & 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  & 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  & 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  & alb1,alb2,rlaimx,rlai
381 
382 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  & 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  & albmax,albmin,rhomax,rhomin,toexp,tolin,wcrinc,wcring,wnew,&
400  & xcrinr,xcrinv
401 
402 WRITE(unit=kulout,fmt='('' SODELX = ''/5E11.4/5E11.4)') sodelx
403 IF(ysp_sbd%NLEVS > 9) CALL abor1(' TOO MANY SOIL LAYERS !')
404 
405 WRITE(unit=kulout,fmt='('' COMMON YOMVDOZ '')')
406 
407 WRITE(unit=kulout,fmt='('' LRDIFOZ ='',L2,'' LRDEPOZ ='',L2)')lrdifoz,lrdepoz
408 
409 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 IF (lhook) CALL dr_hook('SUPHY1',1,zhook_handle)
449 END SUBROUTINE suphy1
yomvdoz::vdhjs
real(kind=jprb), dimension(99) vdhjs
Definition: yomvdoz.F90:51
yomphy1::g2wsat
real(kind=jprb) g2wsat
Definition: yomphy1.F90:255
yomphy1::rhomax
real(kind=jprb) rhomax
Definition: yomphy1.F90:270
yomcst::rpi
real(kind=jprb) rpi
Definition: yomcst.F90:15
yomphy1::rlai
real(kind=jprb) rlai
Definition: yomphy1.F90:303
yomphy1::g1p
real(kind=jprb) g1p
Definition: yomphy1.F90:249
yomphy1::rd2gla
real(kind=jprb) rd2gla
Definition: yomphy1.F90:268
yomphy1::rz0gla
real(kind=jprb) rz0gla
Definition: yomphy1.F90:274
yomphy1::tolin
real(kind=jprb) tolin
Definition: yomphy1.F90:282
yomphy1::glaimxs
real(kind=jprb) glaimxs
Definition: yomphy1.F90:298
yomphy1::hsolit0
real(kind=jprb) hsolit0
Definition: yomphy1.F90:261
yomvdoz::vozhs
real(kind=jprb) vozhs
Definition: yomvdoz.F90:72
yomphy1::gf3
real(kind=jprb), dimension(18) gf3
Definition: yomphy1.F90:200
yomcst::rday
real(kind=jprb) rday
Definition: yomcst.F90:21
yomphy1::gwlmx
real(kind=jprb) gwlmx
Definition: yomphy1.F90:244
yomphy1::rlaimx
real(kind=jprb) rlaimx
Definition: yomphy1.F90:302
yomphy1::hsoliwr
real(kind=jprb) hsoliwr
Definition: yomphy1.F90:260
yomvdoz::vdnjs
real(kind=jprb), dimension(99) vdnjs
Definition: yomvdoz.F90:67
yomvdoz::vdpns
real(kind=jprb), dimension(99) vdpns
Definition: yomvdoz.F90:57
yomphy1::emmmer
real(kind=jprb) emmmer
Definition: yomphy1.F90:219
yomphy1::gc2
real(kind=jprb) gc2
Definition: yomphy1.F90:235
yomphy1::gc32
real(kind=jprb) gc32
Definition: yomphy1.F90:239
yomphy1::gc1s1
real(kind=jprb) gc1s1
Definition: yomphy1.F90:224
yomphy1::wpmx
real(kind=jprb) wpmx
Definition: yomphy1.F90:287
yomphy1::tref4
real(kind=jprb), dimension(18) tref4
Definition: yomphy1.F90:202
yomphy1::gcgels
real(kind=jprb) gcgels
Definition: yomphy1.F90:296
surface_fields::ysd_vvd
type(type_surf_gen) ysd_vvd
Definition: surface_fields.F90:457
yomphy1::g2cgsat
real(kind=jprb) g2cgsat
Definition: yomphy1.F90:252
yomphy1::gsnc2
real(kind=jprb) gsnc2
Definition: yomphy1.F90:258
yomphy1::rd2mer
real(kind=jprb) rd2mer
Definition: yomphy1.F90:269
yomphy1::rctgla
real(kind=jprb) rctgla
Definition: yomphy1.F90:265
yomphy1::gc1y1
real(kind=jprb) gc1y1
Definition: yomphy1.F90:228
yomphy1::ntvgla
integer(kind=jpim) ntvgla
Definition: yomphy1.F90:294
yomphy1::g1cgsat
real(kind=jprb) g1cgsat
Definition: yomphy1.F90:247
surface_fields
Definition: surface_fields.F90:1
yomphy1::g2c1sat
real(kind=jprb) g2c1sat
Definition: yomphy1.F90:253
yomphy1::gvegmx
real(kind=jprb) gvegmx
Definition: yomphy1.F90:230
yomphy1::sodelx
real(kind=jprb), dimension(0:9) sodelx
Definition: yomphy1.F90:205
yomvdoz::voznj
real(kind=jprb) voznj
Definition: yomvdoz.F90:71
yomphy1::gc3
real(kind=jprb) gc3
Definition: yomphy1.F90:237
yomphy1::albmin
real(kind=jprb) albmin
Definition: yomphy1.F90:211
yomphy1::alb1
real(kind=jprb) alb1
Definition: yomphy1.F90:300
abor1
subroutine abor1(CDTEXT)
Definition: abor1.F90:2
yomphy1::wcring
real(kind=jprb) wcring
Definition: yomphy1.F90:285
yomphy1::rzhz0g
real(kind=jprb) rzhz0g
Definition: yomphy1.F90:276
yomphy1::gsnc1
real(kind=jprb) gsnc1
Definition: yomphy1.F90:257
yomphy1::albmax
real(kind=jprb) albmax
Definition: yomphy1.F90:208
yomphy1::xcrinr
real(kind=jprb) xcrinr
Definition: yomphy1.F90:289
yomvdoz::vdanh
real(kind=jprb), dimension(99) vdanh
Definition: yomvdoz.F90:66
yomphy1::rc1max
real(kind=jprb) rc1max
Definition: yomphy1.F90:264
yomcst::rg
real(kind=jprb) rg
Definition: yomcst.F90:29
yomphy1::limc
logical limc
Definition: yomphy1.F90:291
yomphy1::wcrin
real(kind=jprb) wcrin
Definition: yomphy1.F90:283
yomphy1::gcgel
real(kind=jprb) gcgel
Definition: yomphy1.F90:234
yomphy1::rtiner
real(kind=jprb) rtiner
Definition: yomphy1.F90:273
yomphy1::hsol
real(kind=jprb) hsol
Definition: yomphy1.F90:259
yomphy1::g2p
real(kind=jprb) g2p
Definition: yomphy1.F90:254
yomvdoz::vdejh
real(kind=jprb), dimension(99) vdejh
Definition: yomvdoz.F90:60
false
!$Id itapm1 ENDIF!IM on interpole les champs sur les niveaux STD de pression!IM a chaque pas de temps de la physique c!positionnement de l argument logique a false c!pour ne pas recalculer deux fois la meme chose!c!a cet effet un appel a plevel_new a ete deplace c!a la fin de la serie d appels c!la boucle DO nlevSTD a ete internalisee c!dans d ou la creation de cette routine c c!CALL false
Definition: calcul_STDlev.h:26
yomphy1::glaimx
real(kind=jprb) glaimx
Definition: yomphy1.F90:231
yomphy1::limw
logical limw
Definition: yomphy1.F90:292
yomvdoz::vdhns
real(kind=jprb), dimension(99) vdhns
Definition: yomvdoz.F90:53
yomphy1::rhomin
real(kind=jprb) rhomin
Definition: yomphy1.F90:271
yomct0b::lecmwf
logical lecmwf
Definition: yomct0b.F90:15
yomvdoz::vdpjs
real(kind=jprb), dimension(99) vdpjs
Definition: yomvdoz.F90:55
yomphy1::gc1s4
real(kind=jprb) gc1s4
Definition: yomphy1.F90:227
yomphy1::rzhgla
real(kind=jprb) rzhgla
Definition: yomphy1.F90:278
yomphy1::gc31
real(kind=jprb) gc31
Definition: yomphy1.F90:238
yomcst::rcs
real(kind=jprb) rcs
Definition: yomcst.F90:58
yomvdoz::vdans
real(kind=jprb), dimension(99) vdans
Definition: yomvdoz.F90:65
yomphy1::rz0mer
real(kind=jprb) rz0mer
Definition: yomphy1.F90:275
yomvdoz::vdajh
real(kind=jprb), dimension(99) vdajh
Definition: yomvdoz.F90:64
parkind1::jprb
integer, parameter jprb
Definition: parkind1.F90:31
yomvdoz::vdenh
real(kind=jprb), dimension(99) vdenh
Definition: yomvdoz.F90:62
yomphy1::g2b
real(kind=jprb) g2b
Definition: yomphy1.F90:251
yomphy1::albmed
real(kind=jprb) albmed
Definition: yomphy1.F90:210
yomphy1::gc2ref
real(kind=jprb) gc2ref
Definition: yomphy1.F90:236
yomphy1::gtsvap
real(kind=jprb) gtsvap
Definition: yomphy1.F90:229
yomphy1::emmgla
real(kind=jprb) emmgla
Definition: yomphy1.F90:218
yomphy1::gf1
real(kind=jprb) gf1
Definition: yomphy1.F90:241
yomphy1::xcrinv
real(kind=jprb) xcrinv
Definition: yomphy1.F90:290
yomphy1::gconv
real(kind=jprb) gconv
Definition: yomphy1.F90:240
yomphy1::gc1s2
real(kind=jprb) gc1s2
Definition: yomphy1.F90:225
yomphy1::nchsp
integer(kind=jpim) nchsp
Definition: yomphy1.F90:304
yomphy1::rctveg
real(kind=jprb), dimension(18) rctveg
Definition: yomphy1.F90:203
yomvdoz::vdejs
real(kind=jprb), dimension(99) vdejs
Definition: yomvdoz.F90:59
yomphy1::ewwilt
real(kind=jprb) ewwilt
Definition: yomphy1.F90:221
yomphy1::albmer
real(kind=jprb) albmer
Definition: yomphy1.F90:209
yomvdoz::vdajs
real(kind=jprb), dimension(99) vdajs
Definition: yomvdoz.F90:63
yomvdoz::vdnnh
real(kind=jprb), dimension(99) vdnnh
Definition: yomvdoz.F90:70
yomvdoz::vdnns
real(kind=jprb), dimension(99) vdnns
Definition: yomvdoz.F90:69
yomphy1::lc1vap
logical lc1vap
Definition: yomphy1.F90:293
yomphy1::gneimx
real(kind=jprb) gneimx
Definition: yomphy1.F90:232
yomphy1::wsmx
real(kind=jprb) wsmx
Definition: yomphy1.F90:288
yomphy1::omtpro
real(kind=jprb) omtpro
Definition: yomphy1.F90:262
yomphy1::alrcn2
real(kind=jprb) alrcn2
Definition: yomphy1.F90:214
yomphy1::albgla
real(kind=jprb) albgla
Definition: yomphy1.F90:207
yomphy1::alb2
real(kind=jprb) alb2
Definition: yomphy1.F90:301
yomvdoz::vdpjh
real(kind=jprb), dimension(99) vdpjh
Definition: yomvdoz.F90:56
true
!$Id itapm1 ENDIF!IM on interpole les champs sur les niveaux STD de pression!IM a chaque pas de temps de la physique c!positionnement de l argument logique a false c!pour ne pas recalculer deux fois la meme chose!c!a cet effet un appel a plevel_new a ete deplace c!a la fin de la serie d appels c!la boucle DO nlevSTD a ete internalisee c!dans d ou la creation de cette routine c c!CALL ulevSTD CALL &zphi philevSTD CALL &zx_rh rhlevSTD!DO klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon du jour ou toutes les read_climoz CALL true
Definition: calcul_STDlev.h:139
yomvdoz::lrdepoz
logical lrdepoz
Definition: yomvdoz.F90:74
yomphy1::rsmax
real(kind=jprb) rsmax
Definition: yomphy1.F90:272
yomphy1::wnew
real(kind=jprb) wnew
Definition: yomphy1.F90:286
yomvdoz::vdpnh
real(kind=jprb), dimension(99) vdpnh
Definition: yomvdoz.F90:58
yomphy1::rd1
real(kind=jprb) rd1
Definition: yomphy1.F90:267
yomphy1::g1c1sat
real(kind=jprb) g1c1sat
Definition: yomphy1.F90:248
yomphy1::gwwilt
real(kind=jprb) gwwilt
Definition: yomphy1.F90:245
yomhook::lhook
logical lhook
Definition: yomhook.F90:12
yomvdoz::vdnjh
real(kind=jprb), dimension(99) vdnjh
Definition: yomvdoz.F90:68
yomphy1::tmergl
real(kind=jprb) tmergl
Definition: yomphy1.F90:280
yomphy1::gneimxs
real(kind=jprb) gneimxs
Definition: yomphy1.F90:299
yomphy1::g1wsat
real(kind=jprb) g1wsat
Definition: yomphy1.F90:250
yomphy1::omwpro
real(kind=jprb) omwpro
Definition: yomphy1.F90:263
yomphy1::gcz0h
real(kind=jprb), dimension(0:3, 4) gcz0h
Definition: yomphy1.F90:206
yomvdoz::vdhnh
real(kind=jprb), dimension(99) vdhnh
Definition: yomvdoz.F90:54
yomphy1::gc1
real(kind=jprb) gc1
Definition: yomphy1.F90:223
yomhook::dr_hook
subroutine dr_hook(CDNAME, KSWITCH, PKEY)
Definition: yomhook.F90:17
yomphy1::gwlex
real(kind=jprb) gwlex
Definition: yomphy1.F90:243
yomphy1::gf4
real(kind=jprb), dimension(18) gf4
Definition: yomphy1.F90:201
yomphy1::ntvmer
integer(kind=jpim) ntvmer
Definition: yomphy1.F90:295
surface_fields::ysp_sbd
type(type_surf_gen) ysp_sbd
Definition: surface_fields.F90:414
yomphy1::wcrinc
real(kind=jprb) wcrinc
Definition: yomphy1.F90:284
yomvdoz::vdens
real(kind=jprb), dimension(99) vdens
Definition: yomvdoz.F90:61
yomphy1::alcrin
real(kind=jprb) alcrin
Definition: yomphy1.F90:212
yomphy1::alrcn1
real(kind=jprb) alrcn1
Definition: yomphy1.F90:213
yomphy1::ea
real(kind=jprb) ea
Definition: yomphy1.F90:215
parkind1::jpim
integer, parameter jpim
Definition: parkind1.F90:13
yomvdoz::vdhjh
real(kind=jprb), dimension(99) vdhjh
Definition: yomvdoz.F90:52
yomphy1::ga
real(kind=jprb) ga
Definition: yomphy1.F90:222
yomphy1::gwfc
real(kind=jprb) gwfc
Definition: yomphy1.F90:242
yomphy1::toexp
real(kind=jprb) toexp
Definition: yomphy1.F90:281
suphy1
subroutine suphy1(KULOUT)
Definition: suphy1.F90:3
yomphy1::ec2ref
real(kind=jprb) ec2ref
Definition: yomphy1.F90:216
yomcst
Definition: yomcst.F90:1
yomphy1::g1b
real(kind=jprb) g1b
Definition: yomphy1.F90:246
yomphy1::gvegmxs
real(kind=jprb) gvegmxs
Definition: yomphy1.F90:297
yomphy1::g3cgsat
real(kind=jprb) g3cgsat
Definition: yomphy1.F90:256
yomphy1::gwpimx
real(kind=jprb) gwpimx
Definition: yomphy1.F90:233
unit
!$Header!integer nvarmx s s unit
Definition: gradsdef.h:20
yomphy1::gc1s3
real(kind=jprb) gc1s3
Definition: yomphy1.F90:226
yomphy1::rgl
real(kind=jprb), dimension(18) rgl
Definition: yomphy1.F90:204
yomphy1::emcrin
real(kind=jprb) emcrin
Definition: yomphy1.F90:217
yomphy1::rcgmax
real(kind=jprb) rcgmax
Definition: yomphy1.F90:266
yomphy1::rzhz0m
real(kind=jprb) rzhz0m
Definition: yomphy1.F90:277
yomphy1::rzhmer
real(kind=jprb) rzhmer
Definition: yomphy1.F90:279
yomvdoz::lrdifoz
logical lrdifoz
Definition: yomvdoz.F90:73
yomphy1::ewfc
real(kind=jprb) ewfc
Definition: yomphy1.F90:220