[3331] | 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 |
---|