| Directory: | ./ |
|---|---|
| File: | rad/suphy1.f90 |
| Date: | 2022-01-11 19:19:34 |
| Exec | Total | Coverage | |
|---|---|---|---|
| Lines: | 176 | 184 | 95.7% |
| Branches: | 10 | 58 | 17.2% |
| Line | Branch | Exec | Source |
|---|---|---|---|
| 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 | INTERFACE | ||
| 105 | SUBROUTINE ABOR1(CDTEXT) | ||
| 106 | CHARACTER(LEN=*) :: CDTEXT | ||
| 107 | END SUBROUTINE ABOR1 | ||
| 108 | END INTERFACE | ||
| 109 | INTERFACE | ||
| 110 | SUBROUTINE POSNAM(KULNAM,CDNAML) | ||
| 111 | USE PARKIND1 ,ONLY : JPIM ,JPRB | ||
| 112 | INTEGER(KIND=JPIM),INTENT(IN) :: KULNAM | ||
| 113 | CHARACTER(LEN=*) ,INTENT(IN) :: CDNAML | ||
| 114 | END SUBROUTINE POSNAM | ||
| 115 | END INTERFACE | ||
| 116 | |||
| 117 | ! ------------------------------------------------------------------ | ||
| 118 | NAMELIST/NAMPHY1/& | ||
| 119 | &ALBGLA , ALBMAX , ALBMER , ALBMIN , ALCRIN , ALRCN1 ,& | ||
| 120 | &ALRCN2 , EA , EC2REF , EMCRIN , EMMGLA , EMMMER ,& | ||
| 121 | &EWFC , EWWILT , GA , GC1 ,& | ||
| 122 | &GC1S1 , GC1S2 , GC1S3 , GC1S4 , GC1Y1 , GTSVAP ,& | ||
| 123 | &GVEGMX , GLAIMX , GNEIMX , GWPIMX , GCGEL ,& | ||
| 124 | &GC2 , GC2REF , GC3 , GC31 , GC32 ,& | ||
| 125 | &GCZ0H , GF1 , GF3 , GF4 , GWFC , GWLEX ,& | ||
| 126 | &GWLMX , GWWILT , G1B ,& | ||
| 127 | &G1CGSAT , G1C1SAT , G1P , G1WSAT , G2B , G2CGSAT ,& | ||
| 128 | &G2C1SAT , G2P , G2WSAT , G3CGSAT , GSNC1 , GSNC2 ,& | ||
| 129 | &HSOL , HSOLIWR , HSOLIT0 , OMTPRO , OMWPRO , RC1MAX ,& | ||
| 130 | &RCTGLA , RCGMAX , RCTVEG , RD1 , RD2GLA , RD2MER ,& | ||
| 131 | &RGL , RHOMAX , RHOMIN , RSMAX , RTINER , RZ0GLA ,& | ||
| 132 | &RZ0MER , RZHZ0G , RZHZ0M , TMERGL , TOEXP , TOLIN ,& | ||
| 133 | &TREF4 , SODELX , WCRIN , WCRINC , WCRING , WNEW ,& | ||
| 134 | &WPMX , WSMX , XCRINR , XCRINV ,& | ||
| 135 | &LIMC , LIMW , LC1VAP ,& | ||
| 136 | &NTVGLA , NTVMER , GCGELS , GVEGMXS , GLAIMXS , GNEIMXS ,& | ||
| 137 | &ALB1 , ALB2 , RLAIMX , RLAI ,NCHSP | ||
| 138 | ! --------------------------------------------------------------- | ||
| 139 | NAMELIST/NAMVDOZ/& | ||
| 140 | &VDHJS,VDHJH,VDHNS,VDHNH & | ||
| 141 | &,VDPJS,VDPJH,VDPNS,VDPNH & | ||
| 142 | &,VDEJS,VDEJH,VDENS,VDENH & | ||
| 143 | &,VDAJS,VDAJH,VDANS,VDANH & | ||
| 144 | &,VDNJS,VDNJH,VDNNS,VDNNH & | ||
| 145 | &,VOZNJ,VOZHS & | ||
| 146 | &,LRDIFOZ,LRDEPOZ | ||
| 147 | |||
| 148 | |||
| 149 | ! ------------------------------------------------------------------ | ||
| 150 | |||
| 151 | !* 1. Set default values. | ||
| 152 | ! ------------------- | ||
| 153 | |||
| 154 | ! 1.1 Set implicit default values | ||
| 155 | |||
| 156 | ! Snow , Sea-Ice | ||
| 157 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | IF (LHOOK) CALL DR_HOOK('SUPHY1',0,ZHOOK_HANDLE) |
| 158 | 1 | ALCRIN=0.7_JPRB | |
| 159 | 1 | ALRCN1=1.E-03_JPRB | |
| 160 | 1 | ALRCN2=2.5E-03_JPRB | |
| 161 | 1 | EMCRIN=0.98_JPRB | |
| 162 | 1 | WCRIN=10._JPRB | |
| 163 | 1 | TMERGL=271.23_JPRB | |
| 164 | 1 | NCHSP=0 | |
| 165 | ! Characteristics of ice and sea | ||
| 166 | 1 | NTVGLA=2 | |
| 167 | 1 | NTVMER=1 | |
| 168 | 1 | RD2GLA=8._JPRB | |
| 169 | 1 | RD2MER=8._JPRB | |
| 170 | 1 | ALBGLA=.65_JPRB | |
| 171 | 1 | ALBMER=.07_JPRB | |
| 172 | 1 | ALBMED=.06_JPRB | |
| 173 | 1 | EMMGLA=.97_JPRB | |
| 174 | 1 | EMMMER=.96_JPRB | |
| 175 | 1 | RZ0GLA=.001_JPRB | |
| 176 | 1 | RZ0MER=.001_JPRB | |
| 177 | 1 | RZHGLA=.001_JPRB | |
| 178 | 1 | RZHMER=.001_JPRB | |
| 179 | 1 | RZHZ0G=1.0_JPRB | |
| 180 | 1 | RZHZ0M=1.0_JPRB | |
| 181 | ! Usual - soil | ||
| 182 | 1 | HSOL=1.0E-05_JPRB | |
| 183 | 1 | HSOLIT0=0.35_JPRB | |
| 184 | 1 | HSOLIWR=6._JPRB | |
| 185 | 1 | RTINER=5._JPRB | |
| 186 | 1 | WPMX=100._JPRB | |
| 187 | 1 | WSMX=20._JPRB | |
| 188 | 1 | SODELX(0)=1.0_JPRB/SQRT(1.0_JPRB+2.0_JPRB*RPI) | |
| 189 |
2/2✓ Branch 0 taken 9 times.
✓ Branch 1 taken 1 times.
|
10 | DO J=1,9 |
| 190 | 10 | SODELX(J)=SODELX(J-1)*2.0_JPRB*RPI | |
| 191 | ENDDO | ||
| 192 | ! Relaxation | ||
| 193 | 1 | OMTPRO=0._JPRB | |
| 194 | 1 | OMWPRO=0._JPRB | |
| 195 | ! ISBA - soil | ||
| 196 | 1 | EA=-0.54_JPRB | |
| 197 | 1 | GA=732.42E-3_JPRB | |
| 198 | 1 | G1B=0.137_JPRB | |
| 199 | 1 | G2B=3.5_JPRB | |
| 200 | 1 | G1P=0.134_JPRB | |
| 201 | 1 | G2P=3.4_JPRB | |
| 202 | 1 | GC1=0.5_JPRB | |
| 203 | 1 | GC2=10._JPRB | |
| 204 | 1 | GC3=8._JPRB | |
| 205 | 1 | GCONV=1.E3_JPRB | |
| 206 | |||
| 207 | 1 | G1WSAT=-1.08E-3_JPRB | |
| 208 | 1 | G2WSAT=494.31E-3_JPRB | |
| 209 | 1 | EWFC=0.35_JPRB | |
| 210 | 1 | GWFC=89.0467E-3_JPRB | |
| 211 | 1 | EWWILT=0.5_JPRB | |
| 212 | 1 | GWWILT=37.1342E-3_JPRB | |
| 213 | 1 | EC2REF=-0.95_JPRB | |
| 214 | 1 | GC2REF=13.82_JPRB | |
| 215 | 1 | G1CGSAT=-1.5571E-8_JPRB | |
| 216 | 1 | G2CGSAT=-1.441E-8_JPRB | |
| 217 | 1 | G3CGSAT=4.70217E-6_JPRB | |
| 218 | 1 | G1C1SAT=5.58E-3_JPRB | |
| 219 | 1 | G2C1SAT=84.88E-3_JPRB | |
| 220 | 1 | GC31=5.3275_JPRB | |
| 221 | 1 | GC32=-1.043_JPRB | |
| 222 | |||
| 223 | 1 | RD1=1.E-2_JPRB | |
| 224 | 1 | RC1MAX=500._JPRB | |
| 225 | 1 | RCTGLA=5.5E-6_JPRB | |
| 226 | 1 | RCGMAX=0.8E-5_JPRB | |
| 227 | 1 | LIMC=.TRUE. | |
| 228 | 1 | LIMW=.TRUE. | |
| 229 | ! ISBA - vegetation | ||
| 230 | 1 | GF1=0.55_JPRB | |
| 231 | 1 | GWLEX=2.0_JPRB/3._JPRB | |
| 232 | 1 | GWLMX=0.2_JPRB | |
| 233 | 1 | RSMAX=5000._JPRB | |
| 234 |
2/2✓ Branch 0 taken 18 times.
✓ Branch 1 taken 1 times.
|
19 | DO J=1,18 |
| 235 | 18 | GF3(J)=0.0_JPRB | |
| 236 | 18 | GF4(J)=0.0016_JPRB | |
| 237 | 18 | RCTVEG(J)=0.8E-5_JPRB | |
| 238 | 18 | RGL(J)=100._JPRB | |
| 239 | 19 | TREF4(J)=298._JPRB | |
| 240 | ENDDO | ||
| 241 | 1 | GF3(4)=40._JPRB | |
| 242 | 1 | RGL(4)=30._JPRB | |
| 243 | ! ISBA - roughness length | ||
| 244 | 1 | GCZ0H(0,1)=7.5_JPRB | |
| 245 | 1 | GCZ0H(1,1)=2.39037_JPRB | |
| 246 | 1 | GCZ0H(2,1)=-.28583_JPRB | |
| 247 | 1 | GCZ0H(3,1)=.01074_JPRB | |
| 248 | 1 | GCZ0H(0,2)=0.5_JPRB | |
| 249 | 1 | GCZ0H(1,2)=-.07028_JPRB | |
| 250 | 1 | GCZ0H(2,2)=.01023_JPRB | |
| 251 | 1 | GCZ0H(3,2)=-.00067_JPRB | |
| 252 | 1 | GCZ0H(0,3)=5.0_JPRB | |
| 253 | 1 | GCZ0H(1,3)=4.51268_JPRB | |
| 254 | 1 | GCZ0H(2,3)=.34012_JPRB | |
| 255 | 1 | GCZ0H(3,3)=-.05330_JPRB | |
| 256 | 1 | GCZ0H(0,4)=0.5_JPRB | |
| 257 | 1 | GCZ0H(1,4)=-.09421_JPRB | |
| 258 | 1 | GCZ0H(2,4)=.01463_JPRB | |
| 259 | 1 | GCZ0H(3,4)=-.00099_JPRB | |
| 260 | ! ISBA - snow | ||
| 261 | 1 | ALBMAX=0.85_JPRB | |
| 262 | 1 | ALBMIN=0.50_JPRB | |
| 263 | 1 | RHOMAX=0.3_JPRB | |
| 264 | 1 | RHOMIN=0.1_JPRB | |
| 265 | 1 | TOEXP=0.24_JPRB/86400._JPRB | |
| 266 | 1 | TOLIN=0.008_JPRB/86400._JPRB | |
| 267 | 1 | WCRINC=70._JPRB | |
| 268 | 1 | WCRING=10._JPRB | |
| 269 | 1 | WNEW=10._JPRB | |
| 270 | 1 | XCRINR=1.0_JPRB/RG | |
| 271 | 1 | XCRINV=10000._JPRB | |
| 272 | 1 | GSNC1=RPI/(2.22_JPRB*RCS*RDAY*1000._JPRB) | |
| 273 | 1 | GSNC2=2.885_JPRB | |
| 274 | ! ISBA - C1 vapour phase | ||
| 275 | 1 | LC1VAP=.TRUE. | |
| 276 | 1 | GTSVAP=0._JPRB | |
| 277 | 1 | GC1S1= 1.19_JPRB | |
| 278 | 1 | GC1S2=-5.09_JPRB | |
| 279 | 1 | GC1S3=-1.464E+2_JPRB | |
| 280 | 1 | GC1S4= 17.86E+2_JPRB | |
| 281 | 1 | GC1Y1=10._JPRB | |
| 282 | ! Deep Soil freezing | ||
| 283 | 1 | GVEGMX=5._JPRB | |
| 284 | 1 | GLAIMX=30._JPRB | |
| 285 | 1 | GNEIMX=1.8_JPRB | |
| 286 | 1 | GWPIMX=150._JPRB | |
| 287 | 1 | GCGEL=3.E-5_JPRB | |
| 288 | ! Surface soil freezing | ||
| 289 | 1 | GCGELS=5.E-5_JPRB | |
| 290 | 1 | GVEGMXS=5._JPRB | |
| 291 | 1 | GLAIMXS=30._JPRB | |
| 292 | 1 | GNEIMXS=1.8_JPRB | |
| 293 | ! OZONE DIFFUSION AND DEPOSITION | ||
| 294 | 1 | LRDIFOZ=.FALSE. | |
| 295 | 1 | LRDEPOZ=.FALSE. | |
| 296 | 1 | VOZNJ=1._JPRB | |
| 297 | 1 | VOZHS=1.0_JPRB/86400._JPRB | |
| 298 |
2/2✓ Branch 0 taken 99 times.
✓ Branch 1 taken 1 times.
|
100 | DO J=1,99 |
| 299 | 99 | VDHJS(J)=-999._JPRB | |
| 300 | 99 | VDHJH(J)=-999._JPRB | |
| 301 | 99 | VDHNS(J)=-999._JPRB | |
| 302 | 99 | VDHNH(J)=-999._JPRB | |
| 303 | 99 | VDPJS(J)=-999._JPRB | |
| 304 | 99 | VDPJH(J)=-999._JPRB | |
| 305 | 99 | VDPNS(J)=-999._JPRB | |
| 306 | 99 | VDPNH(J)=-999._JPRB | |
| 307 | 99 | VDEJS(J)=-999._JPRB | |
| 308 | 99 | VDEJH(J)=-999._JPRB | |
| 309 | 99 | VDENS(J)=-999._JPRB | |
| 310 | 99 | VDENH(J)=-999._JPRB | |
| 311 | 99 | VDAJS(J)=-999._JPRB | |
| 312 | 99 | VDAJH(J)=-999._JPRB | |
| 313 | 99 | VDANS(J)=-999._JPRB | |
| 314 | 99 | VDANH(J)=-999._JPRB | |
| 315 | 99 | VDNJS(J)=-999._JPRB | |
| 316 | 99 | VDNJH(J)=-999._JPRB | |
| 317 | 99 | VDNNS(J)=-999._JPRB | |
| 318 | 100 | VDNNH(J)=-999._JPRB | |
| 319 | ENDDO | ||
| 320 | ! New snow scheme (LVGSN) | ||
| 321 | 1 | ALB1=0.87_JPRB | |
| 322 | 1 | ALB2=0.84_JPRB | |
| 323 | 1 | RLAIMX=7._JPRB | |
| 324 | 1 | RLAI=3._JPRB | |
| 325 | |||
| 326 | ! 1.2 Modify default values according to LECMWF | ||
| 327 | |||
| 328 | IF (LECMWF) THEN | ||
| 329 | ELSE | ||
| 330 | ENDIF | ||
| 331 | |||
| 332 | ! ------------------------------------------------------------------ | ||
| 333 | |||
| 334 | !* 2. Modify default values. | ||
| 335 | ! ---------------------- | ||
| 336 | |||
| 337 | ! Ce qui concerne NAMPHY1 et NAMVDOZ commente par MPL le 15.04.09 | ||
| 338 | !CALL POSNAM(NULNAM,'NAMPHY1') | ||
| 339 | !READ(NULNAM,NAMPHY1) | ||
| 340 | !CALL POSNAM(NULNAM,'NAMVDOZ') | ||
| 341 | !READ(NULNAM,NAMVDOZ) | ||
| 342 | |||
| 343 | ! 2.5 Check consistency | ||
| 344 | ! ----------------- | ||
| 345 | IF (GC1Y1 > 60._JPRB) THEN | ||
| 346 | CALL ABOR1 ('GC1Y1 FOR C1-VAPOUR PHASE IS BIGGER THAN 60.') | ||
| 347 | ENDIF | ||
| 348 | |||
| 349 | IF ((YSD_VVD%NUMFLDS < 8).AND.(RZHZ0M /= 1.0_JPRB)) THEN | ||
| 350 | CALL ABOR1('YSD_VVD%NUMFLDS<8 IMPLIES RZHZ0M=1.0_JPRB !...') | ||
| 351 | ENDIF | ||
| 352 | |||
| 353 | ! ------------------------------------------------------------------ | ||
| 354 | |||
| 355 | !* 3. Print final values. | ||
| 356 | ! ------------------- | ||
| 357 | |||
| 358 | 1 | WRITE(UNIT=KULOUT,FMT='('' COMMON YOMPHY1 '')') | |
| 359 | |||
| 360 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 361 | & '' ALCRIN ='',E10.4,'' ALRCN1 ='',E10.4,'' ALRCN2 ='',E10.4 & | ||
| 362 | & ,'' EMCRIN ='',E10.4,'' WCRIN ='',E10.4,'' TMERGL ='',E10.4)')& | ||
| 363 | 1 | & ALCRIN,ALRCN1,ALRCN2,EMCRIN,WCRIN,TMERGL | |
| 364 | |||
| 365 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 366 | & '' NCHSP ='',I3,& | ||
| 367 | & '' NTVMER ='',I3,'' NTVGLA ='',I3 & | ||
| 368 | & ,'' RD2MER ='',E10.4,'' ALBMER='',E10.4,'' EMMMER ='',E10.4,/& | ||
| 369 | & ,'' RD2GLA ='',E10.4,'' ALBGLA='',E10.4,'' EMMGLA ='',E10.4,/& | ||
| 370 | & ,'' RZ0MER ='',E10.4,'' RZHMER='',E10.4 & | ||
| 371 | & ,'' RZ0GLA ='',E10.4,'' RZHGLA='',E10.4)')& | ||
| 372 | 1 | & NCHSP,NTVMER,NTVGLA,RD2MER,ALBMER,EMMMER & | |
| 373 | 2 | & ,RD2GLA,ALBGLA,EMMGLA,RZ0MER,RZHMER,RZ0GLA,RZHGLA | |
| 374 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 375 | & '' RTINER ='',E10.4,'' HSOL ='',E10.4,'' HSOLIT0 ='',E10.4 & | ||
| 376 | & ,'' HSOLIWR ='',E10.4,'' WPMX ='',E10.4 & | ||
| 377 | & ,'' WSMX ='',E10.4,'' OMTPRO ='',E10.4,'' OMWPRO ='',E10.4)')& | ||
| 378 | 1 | & RTINER,HSOL,HSOLIT0,HSOLIWR,WPMX,WSMX,OMTPRO,OMWPRO | |
| 379 | |||
| 380 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 381 | & '' EA ='',E10.4,'' GA ='',E10.4,'' G1B ='',E10.4 & | ||
| 382 | & ,'' G2B ='',E10.4,'' G1P ='',E10.4,'' G2P = '',E10.4,/& | ||
| 383 | & ,'' GC1 ='',E10.4,'' GC2 ='',E10.4,'' GC3 ='',E10.4 & | ||
| 384 | & ,'' GCONV ='',E10.4)')& | ||
| 385 | 1 | & EA,GA,G1B,G2B,G1P,G2P,GC1,GC2,GC3,GCONV | |
| 386 | |||
| 387 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 388 | & '' G1WSAT ='',E10.4,'' G2WSAT ='',E10.4,'' EWFC ='',E10.4 & | ||
| 389 | & ,'' GWFC ='',E10.4,/,'' EWWILT ='',E10.4,'' GWWILT ='',E10.4 & | ||
| 390 | & ,'' EC2REF ='',E10.4,'' GC2REF ='',E10.4,/& | ||
| 391 | & ,'' G1CGSAT ='',E10.4,'' G2CGSAT ='',E10.4,'' G3CGSAT =''& | ||
| 392 | & ,E10.4,'' G1C1SAT ='',E10.4,'' G2C1SAT ='',E10.4)')& | ||
| 393 | 1 | & G1WSAT,G2WSAT,EWFC,GWFC,EWWILT,GWWILT,EC2REF,GC2REF,& | |
| 394 | 2 | & G1CGSAT,G2CGSAT,G3CGSAT,G1C1SAT,G2C1SAT | |
| 395 | |||
| 396 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 397 | & '' RD1 ='',E10.4,'' RC1MAX ='',E10.4,'' RCTGLA ='',E10.4 & | ||
| 398 | & ,'' RCGMAX ='',E10.4,'' LIMC ='',L2,'' LIMW ='',L2)')& | ||
| 399 | 1 | & RD1,RC1MAX,RCTGLA,RCGMAX,LIMC,LIMW | |
| 400 | |||
| 401 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 402 | & '' GC1S1 ='',E10.4,'' GC1S2 ='',E10.4,'' GC1S3 ='',E10.4 & | ||
| 403 | & ,'' GC1S4 ='',E10.4,'' GC1Y1 ='',E10.4,'' LC1VAP ='',L2 & | ||
| 404 | & ,'' GTSVAP ='',E10.4)')& | ||
| 405 | 1 | & GC1S1,GC1S2,GC1S3,GC1S4,GC1Y1,LC1VAP,GTSVAP | |
| 406 | |||
| 407 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 408 | & '' GCGEL ='',E10.4,'' GVEGMX ='',E10.4,'' GLAIMX ='',E10.4 & | ||
| 409 | & ,'' GWPIMX ='',E10.4,'' GNEIMX ='',E10.4)')& | ||
| 410 | 1 | & GCGEL,GVEGMX,GLAIMX,GWPIMX,GNEIMX | |
| 411 | |||
| 412 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 413 | & '' GCGELS ='',E10.4,'' GVEGMXS ='',E10.4 & | ||
| 414 | & ,'' GLAIMXS ='',E10.4,'' GNEIMXS ='',E10.4)')& | ||
| 415 | 1 | & GCGELS,GVEGMXS,GLAIMXS,GNEIMXS | |
| 416 | |||
| 417 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 418 | & '' ALB1 ='',E10.4,'' ALB2 ='',E10.4 & | ||
| 419 | & ,'' RLAIMX ='',E10.4,'' RLAI ='',E10.4)')& | ||
| 420 | 1 | & ALB1,ALB2,RLAIMX,RLAI | |
| 421 | |||
| 422 | 1 | WRITE(UNIT=KULOUT,FMT='('' GCZ0H ='',/,4(1X,4E11.4,/))')GCZ0H | |
| 423 | |||
| 424 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 425 | & '' GF1 ='',E10.4,'' GWLEX ='',E10.4,'' GWLMX ='',E10.4 & | ||
| 426 | & ,'' RSMAX ='',E10.4,/& | ||
| 427 | & ,'' GF3 ='',/,2(1X,9E11.4,/),'' GF4 ='',/,2(1X,9E11.4,/)& | ||
| 428 | & ,'' RCTVEG ='',/,2(1X,9E11.4,/),'' RGL ='',/,2(1X,9E11.4,/)& | ||
| 429 | & ,'' TREF4 ='',/,2(1X,9E11.4,/))')& | ||
| 430 | 1 | & GF1,GWLEX,GWLMX,RSMAX,GF3,GF4,RCTVEG,RGL,TREF4 | |
| 431 | |||
| 432 | WRITE(UNIT=KULOUT,FMT='(& | ||
| 433 | & '' ALBMAX = '',E10.4,'' ALBMIN = '',E10.4 & | ||
| 434 | & ,'' RHOMAX = '',E10.4,'' RHOMIN = '',E10.4,/& | ||
| 435 | & ,'' TOEXP = '',E10.4,'' TOLIN = '',E10.4 & | ||
| 436 | & ,'' WCRINC = '',E10.4,'' WCRING = '',E10.4,/& | ||
| 437 | & ,'' WNEW = '',E10.4 & | ||
| 438 | & ,'' XCRINR = '',E10.4,'' XCRINV = '',E10.4)')& | ||
| 439 | 1 | & ALBMAX,ALBMIN,RHOMAX,RHOMIN,TOEXP,TOLIN,WCRINC,WCRING,WNEW,& | |
| 440 | 2 | & XCRINR,XCRINV | |
| 441 | |||
| 442 | 1 | WRITE(UNIT=KULOUT,FMT='('' SODELX = ''/5E11.4/5E11.4)') SODELX | |
| 443 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF(YSP_SBD%NLEVS > 9) CALL ABOR1(' TOO MANY SOIL LAYERS !') |
| 444 | |||
| 445 | 1 | WRITE(UNIT=KULOUT,FMT='('' COMMON YOMVDOZ '')') | |
| 446 | |||
| 447 | 1 | WRITE(UNIT=KULOUT,FMT='('' LRDIFOZ ='',L2,'' LRDEPOZ ='',L2)')LRDIFOZ,LRDEPOZ | |
| 448 | |||
| 449 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
1 | IF (LRDIFOZ.AND.LRDEPOZ) THEN |
| 450 | ✗ | WRITE(UNIT=KULOUT,FMT='('' VOZNJ ='',F10.5,'' VOZHS ='',F10.5)')VOZNJ,VOZHS | |
| 451 | ✗ | DO IVEG=1,99 | |
| 452 | IF((VDHJS(IVEG) /= -999._JPRB).OR.(VDHJH(IVEG) /= -999._JPRB).OR.& | ||
| 453 | & (VDHNS(IVEG) /= -999._JPRB).OR.(VDHNH(IVEG) /= -999._JPRB).OR.& | ||
| 454 | & (VDPJS(IVEG) /= -999._JPRB).OR.(VDPJH(IVEG) /= -999._JPRB).OR.& | ||
| 455 | & (VDPNS(IVEG) /= -999._JPRB).OR.(VDPNH(IVEG) /= -999._JPRB).OR.& | ||
| 456 | & (VDEJS(IVEG) /= -999._JPRB).OR.(VDEJH(IVEG) /= -999._JPRB).OR.& | ||
| 457 | & (VDENS(IVEG) /= -999._JPRB).OR.(VDENH(IVEG) /= -999._JPRB).OR.& | ||
| 458 | & (VDAJS(IVEG) /= -999._JPRB).OR.(VDAJH(IVEG) /= -999._JPRB).OR.& | ||
| 459 | & (VDANS(IVEG) /= -999._JPRB).OR.(VDANH(IVEG) /= -999._JPRB).OR.& | ||
| 460 | & (VDNJS(IVEG) /= -999._JPRB).OR.(VDNJH(IVEG) /= -999._JPRB).OR.& | ||
| 461 | ✗ | & (VDNNS(IVEG) /= -999._JPRB).OR.(VDNNH(IVEG) /= -999._JPRB))THEN | |
| 462 | WRITE(UNIT=KULOUT,FMT='('' IVEJ ='',I2 & | ||
| 463 | & ,'' VDHJS ='',E9.3,'' VDHJH ='',E9.3,'' VDHNS ='',E9.3 & | ||
| 464 | & ,'' VDHNH ='',E9.3)')& | ||
| 465 | ✗ | & IVEG,VDHJS(IVEG),VDHJH(IVEG),VDHNS(IVEG),VDHNH(IVEG) | |
| 466 | WRITE(UNIT=KULOUT,FMT='('' IVEJ ='',I2 & | ||
| 467 | & ,'' VDPJS ='',E9.3,'' VDPJH ='',E9.3,'' VDPNS ='',E9.3 & | ||
| 468 | & ,'' VDPNH ='',E9.3)')& | ||
| 469 | ✗ | & IVEG,VDPJS(IVEG),VDPJH(IVEG),VDPNS(IVEG),VDPNH(IVEG) | |
| 470 | WRITE(UNIT=KULOUT,FMT='('' IVEJ ='',I2 & | ||
| 471 | & ,'' VDEJS ='',E9.3,'' VDEJH ='',E9.3,'' VDENS ='',E9.3 & | ||
| 472 | & ,'' VDENH ='',E9.3)')& | ||
| 473 | ✗ | & IVEG,VDEJS(IVEG),VDEJH(IVEG),VDENS(IVEG),VDENH(IVEG) | |
| 474 | WRITE(UNIT=KULOUT,FMT='('' IVEJ ='',I2 & | ||
| 475 | & ,'' VDAJS ='',E9.3,'' VDAJH ='',E9.3,'' VDANS ='',E9.3 & | ||
| 476 | & ,'' VDANH ='',E9.3)')& | ||
| 477 | ✗ | & IVEG,VDAJS(IVEG),VDAJH(IVEG),VDANS(IVEG),VDANH(IVEG) | |
| 478 | WRITE(UNIT=KULOUT,FMT='('' IVEJ ='',I2 & | ||
| 479 | & ,'' VDNJS ='',E9.3,'' VDNJH ='',E9.3,'' VDNNS ='',E9.3 & | ||
| 480 | & ,'' VDNNH ='',E9.3)')& | ||
| 481 | ✗ | & IVEG,VDNJS(IVEG),VDNJH(IVEG),VDNNS(IVEG),VDNNH(IVEG) | |
| 482 | ENDIF | ||
| 483 | ENDDO | ||
| 484 | ENDIF | ||
| 485 | !* | ||
| 486 | ! ------------------------------------------------------------------ | ||
| 487 | |||
| 488 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | IF (LHOOK) CALL DR_HOOK('SUPHY1',1,ZHOOK_HANDLE) |
| 489 | 1 | END SUBROUTINE SUPHY1 | |
| 490 |