source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/rrtm/suphy1.F90

Last change on this file was 3331, checked in by acozic, 7 years ago

Add modification for isotopes

  • Property svn:executable set to *
File size: 13.7 KB
RevLine 
[3331]1!OPTIONS XOPT(NOEVAL)
2SUBROUTINE 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
63USE PARKIND1  ,ONLY : JPIM     ,JPRB
64USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
65
66! Ce qui concerne NULNAM commente par MPL le 15.04.09
67!USE YOMLUN   , ONLY : NULNAM
68USE YOMCST   , ONLY : RPI      ,RDAY     ,RG       ,RCS
69USE SURFACE_FIELDS   , ONLY : YSD_VVD, YSP_SBD
70USE YOMCT0B  , ONLY : LECMWF
71USE 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 
93USE 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
99IMPLICIT NONE
100
101INTEGER(KIND=JPIM),INTENT(IN)    :: KULOUT
102INTEGER(KIND=JPIM) :: IVEG, J
103REAL(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
117IF (LHOOK) CALL DR_HOOK('SUPHY1',0,ZHOOK_HANDLE)
118ALCRIN=0.7_JPRB
119ALRCN1=1.E-03_JPRB
120ALRCN2=2.5E-03_JPRB
121EMCRIN=0.98_JPRB
122WCRIN=10._JPRB
123TMERGL=271.23_JPRB
124NCHSP=0
125!  Characteristics of ice and sea
126NTVGLA=2
127NTVMER=1
128RD2GLA=8._JPRB
129RD2MER=8._JPRB
130ALBGLA=.65_JPRB
131ALBMER=.07_JPRB
132ALBMED=.06_JPRB
133EMMGLA=.97_JPRB
134EMMMER=.96_JPRB
135RZ0GLA=.001_JPRB
136RZ0MER=.001_JPRB
137RZHGLA=.001_JPRB
138RZHMER=.001_JPRB
139RZHZ0G=1.0_JPRB
140RZHZ0M=1.0_JPRB
141!  Usual - soil
142HSOL=1.0E-05_JPRB
143HSOLIT0=0.35_JPRB
144HSOLIWR=6._JPRB
145RTINER=5._JPRB
146WPMX=100._JPRB
147WSMX=20._JPRB
148SODELX(0)=1.0_JPRB/SQRT(1.0_JPRB+2.0_JPRB*RPI)
149DO J=1,9
150  SODELX(J)=SODELX(J-1)*2.0_JPRB*RPI
151ENDDO
152!  Relaxation
153OMTPRO=0._JPRB
154OMWPRO=0._JPRB
155!  ISBA - soil
156EA=-0.54_JPRB
157GA=732.42E-3_JPRB
158G1B=0.137_JPRB
159G2B=3.5_JPRB
160G1P=0.134_JPRB
161G2P=3.4_JPRB
162GC1=0.5_JPRB
163GC2=10._JPRB
164GC3=8._JPRB
165GCONV=1.E3_JPRB
166
167G1WSAT=-1.08E-3_JPRB
168G2WSAT=494.31E-3_JPRB
169EWFC=0.35_JPRB
170GWFC=89.0467E-3_JPRB
171EWWILT=0.5_JPRB
172GWWILT=37.1342E-3_JPRB
173EC2REF=-0.95_JPRB
174GC2REF=13.82_JPRB
175G1CGSAT=-1.5571E-8_JPRB
176G2CGSAT=-1.441E-8_JPRB
177G3CGSAT=4.70217E-6_JPRB
178G1C1SAT=5.58E-3_JPRB
179G2C1SAT=84.88E-3_JPRB
180GC31=5.3275_JPRB
181GC32=-1.043_JPRB
182
183RD1=1.E-2_JPRB
184RC1MAX=500._JPRB
185RCTGLA=5.5E-6_JPRB
186RCGMAX=0.8E-5_JPRB
187LIMC=.TRUE.
188LIMW=.TRUE.
189!  ISBA - vegetation
190GF1=0.55_JPRB
191GWLEX=2.0_JPRB/3._JPRB
192GWLMX=0.2_JPRB
193RSMAX=5000._JPRB
194DO 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
200ENDDO
201GF3(4)=40._JPRB
202RGL(4)=30._JPRB
203!  ISBA - roughness length
204GCZ0H(0,1)=7.5_JPRB
205GCZ0H(1,1)=2.39037_JPRB
206GCZ0H(2,1)=-.28583_JPRB
207GCZ0H(3,1)=.01074_JPRB
208GCZ0H(0,2)=0.5_JPRB
209GCZ0H(1,2)=-.07028_JPRB
210GCZ0H(2,2)=.01023_JPRB
211GCZ0H(3,2)=-.00067_JPRB
212GCZ0H(0,3)=5.0_JPRB
213GCZ0H(1,3)=4.51268_JPRB
214GCZ0H(2,3)=.34012_JPRB
215GCZ0H(3,3)=-.05330_JPRB
216GCZ0H(0,4)=0.5_JPRB
217GCZ0H(1,4)=-.09421_JPRB
218GCZ0H(2,4)=.01463_JPRB
219GCZ0H(3,4)=-.00099_JPRB
220!  ISBA - snow
221ALBMAX=0.85_JPRB
222ALBMIN=0.50_JPRB
223RHOMAX=0.3_JPRB
224RHOMIN=0.1_JPRB
225TOEXP=0.24_JPRB/86400._JPRB
226TOLIN=0.008_JPRB/86400._JPRB
227WCRINC=70._JPRB
228WCRING=10._JPRB
229WNEW=10._JPRB
230XCRINR=1.0_JPRB/RG
231XCRINV=10000._JPRB
232GSNC1=RPI/(2.22_JPRB*RCS*RDAY*1000._JPRB)
233GSNC2=2.885_JPRB
234!   ISBA - C1 vapour phase
235LC1VAP=.TRUE.
236GTSVAP=0._JPRB
237GC1S1= 1.19_JPRB
238GC1S2=-5.09_JPRB
239GC1S3=-1.464E+2_JPRB
240GC1S4= 17.86E+2_JPRB
241GC1Y1=10._JPRB
242! Deep Soil freezing
243GVEGMX=5._JPRB
244GLAIMX=30._JPRB
245GNEIMX=1.8_JPRB
246GWPIMX=150._JPRB
247GCGEL=3.E-5_JPRB
248! Surface soil freezing
249GCGELS=5.E-5_JPRB
250GVEGMXS=5._JPRB
251GLAIMXS=30._JPRB
252GNEIMXS=1.8_JPRB
253!  OZONE DIFFUSION AND DEPOSITION
254LRDIFOZ=.FALSE.
255LRDEPOZ=.FALSE.
256VOZNJ=1._JPRB
257VOZHS=1.0_JPRB/86400._JPRB
258DO 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
279ENDDO
280! New snow scheme (LVGSN)
281ALB1=0.87_JPRB
282ALB2=0.84_JPRB
283RLAIMX=7._JPRB
284RLAI=3._JPRB
285
286!        1.2 Modify default values according to LECMWF
287
288IF (LECMWF) THEN
289ELSE
290ENDIF
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!              -----------------
305IF (GC1Y1 > 60._JPRB) THEN
306  CALL ABOR1 ('GC1Y1 FOR C1-VAPOUR PHASE IS BIGGER THAN 60.')
307ENDIF
308
309IF ((YSD_VVD%NUMFLDS < 8).AND.(RZHZ0M /= 1.0_JPRB)) THEN
310  CALL ABOR1('YSD_VVD%NUMFLDS<8 IMPLIES RZHZ0M=1.0_JPRB !...')
311ENDIF
312
313!     ------------------------------------------------------------------
314
315!*       3.    Print final values.
316!              -------------------
317
318WRITE(UNIT=KULOUT,FMT='('' COMMON YOMPHY1 '')')
319
320WRITE(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
325WRITE(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 
334WRITE(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
340WRITE(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
347WRITE(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
356WRITE(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
361WRITE(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
367WRITE(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
372WRITE(UNIT=KULOUT,FMT='(&
373 & '' GCGELS  ='',E10.4,'' GVEGMXS ='',E10.4 &
374 & ,'' GLAIMXS ='',E10.4,'' GNEIMXS ='',E10.4)')&
375 & GCGELS,GVEGMXS,GLAIMXS,GNEIMXS 
376
377WRITE(UNIT=KULOUT,FMT='(&
378 & '' ALB1  ='',E10.4,'' ALB2 ='',E10.4 &
379 & ,'' RLAIMX ='',E10.4,'' RLAI ='',E10.4)')&
380 & ALB1,ALB2,RLAIMX,RLAI 
381
382WRITE(UNIT=KULOUT,FMT='('' GCZ0H ='',/,4(1X,4E11.4,/))')GCZ0H
383
384WRITE(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
392WRITE(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
402WRITE(UNIT=KULOUT,FMT='('' SODELX = ''/5E11.4/5E11.4)') SODELX
403IF(YSP_SBD%NLEVS > 9)  CALL ABOR1(' TOO MANY SOIL LAYERS !')
404
405WRITE(UNIT=KULOUT,FMT='('' COMMON YOMVDOZ '')')
406
407WRITE(UNIT=KULOUT,FMT='('' LRDIFOZ ='',L2,'' LRDEPOZ ='',L2)')LRDIFOZ,LRDEPOZ
408
409IF (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
444ENDIF
445!*
446!     ------------------------------------------------------------------
447
448IF (LHOOK) CALL DR_HOOK('SUPHY1',1,ZHOOK_HANDLE)
449END SUBROUTINE SUPHY1
Note: See TracBrowser for help on using the repository browser.