| Directory: | ./ |
|---|---|
| File: | rad/suphy3.f90 |
| Date: | 2022-01-11 19:19:34 |
| Exec | Total | Coverage | |
|---|---|---|---|
| Lines: | 218 | 267 | 81.6% |
| Branches: | 3 | 6 | 50.0% |
| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | !OPTIONS XOPT(NOEVAL) | ||
| 2 | 1 | SUBROUTINE SUPHY3(KULOUT) | |
| 3 | |||
| 4 | !**** *SUPHY3* - Initialize common YOMPHY3 physics radiative | ||
| 5 | ! constants | ||
| 6 | |||
| 7 | ! Purpose. | ||
| 8 | ! -------- | ||
| 9 | ! Initialize YOMPHY3, the common that contains the parameters | ||
| 10 | ! for the radiation part of the physics of the model. | ||
| 11 | |||
| 12 | !** Interface. | ||
| 13 | ! ---------- | ||
| 14 | ! *CALL* *SUPHY3(KULOUT) | ||
| 15 | |||
| 16 | ! Explicit arguments : | ||
| 17 | ! -------------------- | ||
| 18 | ! KULOUT : Logical unit for the output | ||
| 19 | |||
| 20 | ! Implicit arguments : | ||
| 21 | ! -------------------- | ||
| 22 | ! COMMON YOMPHY3 | ||
| 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 | ! Ajout de GCE4 et | ||
| 43 | ! modification de toutes les constantes optiques : 91-2-1 | ||
| 44 | ! Ajout des nuages de glace : 92-04, L. Labbe. | ||
| 45 | ! Modified by R. EL Khatib : 93-04-02 Set-up defaults controled by LECMWF | ||
| 46 | ! Modified 97-04-17 J.M. Piriou: radiative coefficients default values, and new QSUSX default. | ||
| 47 | ! Modified 2000-08 R. Randriamampianina, J.F. Geleyn et J.M. Piriou: change infra-red interaction between layers (GIREC*). | ||
| 48 | ! Modified 2000-12-12 by E. Bazile: CYCORA's default value. | ||
| 49 | ! M.Hamrud 01-Oct-2003 CY28 Cleaning | ||
| 50 | ! R.Brozkova 24-Sep-2006 ALARO-0 (new cloud model) | ||
| 51 | ! ------------------------------------------------------------------ | ||
| 52 | |||
| 53 | USE PARKIND1 ,ONLY : JPIM ,JPRB | ||
| 54 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK | ||
| 55 | |||
| 56 | ! Ce qui concerne NULNAM commente par MPL le 15.04.09 | ||
| 57 | !USE YOMLUN , ONLY : NULNAM | ||
| 58 | USE YOMCT0B , ONLY : LECMWF | ||
| 59 | USE YOMPHY , ONLY : LCLSATUR | ||
| 60 | USE YOMPHY3 , ONLY : GCA ,GCB ,GCC ,VDP ,& | ||
| 61 | & VNP ,BSFSA ,BSFSI ,BSFSN ,BSFTA ,& | ||
| 62 | & BSFTI ,BSFTN ,EARRT ,EOASA ,EOASI ,& | ||
| 63 | & EOASN ,EOATA ,EOATI ,EOATN ,EODSA ,& | ||
| 64 | & EODSI ,EODSN ,EODTA ,EODTI ,EODTN ,& | ||
| 65 | & EORAY ,GCD4 ,GCE4 ,QCO2 ,QLIMI ,& | ||
| 66 | & QLIP0 ,RII0 ,USAA ,USAI ,& | ||
| 67 | & USAN ,USBA ,USBI ,USBN ,& | ||
| 68 | & GIREC1 ,GIREC2 ,GIREC3 ,GIREC4 ,& | ||
| 69 | & FCM_DEL_A,FCM_DEL_D,FCM_MU_A ,FCM_MU_D ,FCM_N_I ,& | ||
| 70 | & FCM_N_L ,FCM_P_AI ,FCM_P_AL ,FCM_P_DI ,FCM_P_DL ,& | ||
| 71 | & FCM_P_GI ,FCM_P_GL ,FCM_Q_AI ,FCM_Q_AL ,FCM_Q_DI ,& | ||
| 72 | & FCM_Q_DL ,FCM_Q_GI ,FCM_Q_GL ,N_SPBAND ,REXP_NEB | ||
| 73 | |||
| 74 | IMPLICIT NONE | ||
| 75 | |||
| 76 | INTEGER(KIND=JPIM),INTENT(IN) :: KULOUT | ||
| 77 | REAL(KIND=JPRB) :: ZHOOK_HANDLE | ||
| 78 | |||
| 79 | INTERFACE | ||
| 80 | SUBROUTINE POSNAM(KULNAM,CDNAML) | ||
| 81 | USE PARKIND1 ,ONLY : JPIM ,JPRB | ||
| 82 | INTEGER(KIND=JPIM),INTENT(IN) :: KULNAM | ||
| 83 | CHARACTER(LEN=*) ,INTENT(IN) :: CDNAML | ||
| 84 | END SUBROUTINE POSNAM | ||
| 85 | END INTERFACE | ||
| 86 | |||
| 87 | ! ------------------------------------------------------------------ | ||
| 88 | NAMELIST/NAMPHY3/& | ||
| 89 | &BSFSA,BSFSI,BSFSN,BSFTA,BSFTI,BSFTN,EARRT,EOASA,EOASI,EOASN & | ||
| 90 | &,EOATA,EOATI,EOATN,EODSA,EODSI,EODSN,EODTA,EODTI,EODTN,EORAY & | ||
| 91 | &,GCA,GCB,GCC,GCD4,GCE4,GIREC1,GIREC2,GIREC3,GIREC4 & | ||
| 92 | &,QCO2,QLIMI,QLIP0,RII0,USAA,USAI,USAN,USBA,USBI,USBN,VDP,VNP & | ||
| 93 | &,FCM_DEL_A,FCM_DEL_D,FCM_MU_A,FCM_MU_D,FCM_N_I,FCM_N_L & | ||
| 94 | &,FCM_P_AI,FCM_P_AL,FCM_P_DI,FCM_P_DL,FCM_P_GI,FCM_P_GL & | ||
| 95 | &,FCM_Q_AI,FCM_Q_AL,FCM_Q_DI,FCM_Q_DL,FCM_Q_GI,FCM_Q_GL & | ||
| 96 | &,REXP_NEB | ||
| 97 | ! ------------------------------------------------------------------ | ||
| 98 | |||
| 99 | ! ------------------------------------------------------------------ | ||
| 100 | |||
| 101 | !* 1. Set default values. | ||
| 102 | ! ------------------- | ||
| 103 | |||
| 104 | ! 1.1 Set implicit default values | ||
| 105 | |||
| 106 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | IF (LHOOK) CALL DR_HOOK('SUPHY3',0,ZHOOK_HANDLE) |
| 107 | 1 | BSFSA=.3490_JPRB | |
| 108 | 1 | BSFSI=.3238_JPRB | |
| 109 | 1 | BSFSN=.3300_JPRB | |
| 110 | 1 | BSFTA=.3471_JPRB | |
| 111 | 1 | BSFTI=.3231_JPRB | |
| 112 | 1 | BSFTN=.3584_JPRB | |
| 113 | 1 | EARRT=0.001324_JPRB | |
| 114 | 1 | EOASA=.6727E-01_JPRB | |
| 115 | 1 | EOASI=.1692E-01_JPRB | |
| 116 | 1 | EOASN=.3353E-01_JPRB | |
| 117 | 1 | EOATA=.5037E-01_JPRB | |
| 118 | 1 | EOATI=.3455E+01_JPRB | |
| 119 | 1 | EOATN=.9077E+01_JPRB | |
| 120 | 1 | EODSA=.3665E+00_JPRB | |
| 121 | 1 | EODSI=.1406E+01_JPRB | |
| 122 | 1 | EODSN=.6558E+01_JPRB | |
| 123 | 1 | EODTA=.1766E-01_JPRB | |
| 124 | 1 | EODTI=.8435E+00_JPRB | |
| 125 | 1 | EODTN=.4319E+01_JPRB | |
| 126 | 1 | EORAY=.8606E-06_JPRB | |
| 127 | 1 | GCA(1)=.8041E-01_JPRB | |
| 128 | 1 | GCA(2)=.1456E+00_JPRB | |
| 129 | 1 | GCA(3)=.4787E+01_JPRB | |
| 130 | 1 | GCA(4)=.2102E+04_JPRB | |
| 131 | 1 | GCA(5)=.1334E+01_JPRB | |
| 132 | 1 | GCA(6)=.1551E-04_JPRB | |
| 133 | 1 | GCB(1)=.8968E+07_JPRB | |
| 134 | 1 | GCB(2)=.2413E+10_JPRB | |
| 135 | 1 | GCB(3)=.3548E+05_JPRB | |
| 136 | 1 | GCB(4)=.6370E+10_JPRB | |
| 137 | 1 | GCB(5)=.8499E+11_JPRB | |
| 138 | 1 | GCB(6)=.1012E+06_JPRB | |
| 139 | 1 | GCC(1)=.5925E-10_JPRB | |
| 140 | 1 | GCC(2)=.1842E-10_JPRB | |
| 141 | 1 | GCC(3)=.2532E-07_JPRB | |
| 142 | 1 | GCC(4)=.1953E+07_JPRB | |
| 143 | 1 | GCC(5)=.1734E-11_JPRB | |
| 144 | 1 | GCC(6)=.1225E-16_JPRB | |
| 145 | 1 | GCD4=.3608E-69_JPRB | |
| 146 | 1 | GCE4=.7563E+04_JPRB | |
| 147 | 1 | GIREC1=0.0_JPRB | |
| 148 | 1 | GIREC2=0.0_JPRB | |
| 149 | 1 | GIREC3=0.0_JPRB | |
| 150 | 1 | GIREC4=0.0_JPRB | |
| 151 | 1 | QCO2=.5366E-03_JPRB | |
| 152 | 1 | QLIMI=15000._JPRB | |
| 153 | 1 | QLIP0=8.E+06_JPRB | |
| 154 | 1 | USAA=-.3020_JPRB | |
| 155 | 1 | USAI=-.3524_JPRB | |
| 156 | 1 | USAN=-.3400_JPRB | |
| 157 | 1 | USBA=0._JPRB | |
| 158 | 1 | USBI=0._JPRB | |
| 159 | 1 | USBN=0._JPRB | |
| 160 | 1 | VDP(1,1)=.21868E+02_JPRB | |
| 161 | 1 | VDP(2,1)=.17453E+02_JPRB | |
| 162 | 1 | VDP(3,1)=.68918E+00_JPRB | |
| 163 | 1 | VDP(4,1)=.23456E-03_JPRB | |
| 164 | 1 | VDP(5,1)=.22317E-09_JPRB | |
| 165 | 1 | VDP(1,2)=.48401E+02_JPRB | |
| 166 | 1 | VDP(2,2)=.18648E+02_JPRB | |
| 167 | 1 | VDP(3,2)=.35199E-01_JPRB | |
| 168 | 1 | VDP(4,2)=.63691E-06_JPRB | |
| 169 | 1 | VDP(5,2)=.86395E-13_JPRB | |
| 170 | 1 | VDP(1,3)=.76948E+02_JPRB | |
| 171 | 1 | VDP(2,3)=.41056E+01_JPRB | |
| 172 | 1 | VDP(3,3)=.42667E-03_JPRB | |
| 173 | 1 | VDP(4,3)=0._JPRB | |
| 174 | 1 | VDP(5,3)=0._JPRB | |
| 175 | 1 | VDP(1,4)=.27241E+03_JPRB | |
| 176 | 1 | VDP(2,4)=.57091E+04_JPRB | |
| 177 | 1 | VDP(3,4)=.14393E+05_JPRB | |
| 178 | 1 | VDP(4,4)=.29879E+04_JPRB | |
| 179 | 1 | VDP(5,4)=.25382E+02_JPRB | |
| 180 | 1 | VDP(1,5)=.86942E+02_JPRB | |
| 181 | 1 | VDP(2,5)=.32186E+03_JPRB | |
| 182 | 1 | VDP(3,5)=.10775E+03_JPRB | |
| 183 | 1 | VDP(4,5)=.21261E+01_JPRB | |
| 184 | 1 | VDP(5,5)=.40003E-02_JPRB | |
| 185 | 1 | VDP(1,6)=.24408E+02_JPRB | |
| 186 | 1 | VDP(2,6)=.81919E+01_JPRB | |
| 187 | 1 | VDP(3,6)=.72193E+01_JPRB | |
| 188 | 1 | VDP(4,6)=.56230E+00_JPRB | |
| 189 | 1 | VDP(5,6)=.25384E-02_JPRB | |
| 190 | 1 | VNP(1,1)=.69926E+01_JPRB | |
| 191 | 1 | VNP(2,1)=.63915E+00_JPRB | |
| 192 | 1 | VNP(3,1)=.28896E-03_JPRB | |
| 193 | 1 | VNP(4,1)=.10050E-08_JPRB | |
| 194 | 1 | VNP(5,1)=.99037E-16_JPRB | |
| 195 | 1 | VNP(1,2)=.31105E+01_JPRB | |
| 196 | 1 | VNP(2,2)=.14225E-01_JPRB | |
| 197 | 1 | VNP(3,2)=.69355E-06_JPRB | |
| 198 | 1 | VNP(4,2)=.36087E-12_JPRB | |
| 199 | 1 | VNP(5,2)=.44113E-20_JPRB | |
| 200 | 1 | VNP(1,3)=.23659E+01_JPRB | |
| 201 | 1 | VNP(2,3)=.11139E-02_JPRB | |
| 202 | 1 | VNP(3,3)=.10618E-05_JPRB | |
| 203 | 1 | VNP(4,3)=0._JPRB | |
| 204 | 1 | VNP(5,3)=0._JPRB | |
| 205 | 1 | VNP(1,4)=.19810E+03_JPRB | |
| 206 | 1 | VNP(2,4)=.46954E+04_JPRB | |
| 207 | 1 | VNP(3,4)=.22512E+04_JPRB | |
| 208 | 1 | VNP(4,4)=.52461E+02_JPRB | |
| 209 | 1 | VNP(5,4)=.11645E+00_JPRB | |
| 210 | 1 | VNP(1,5)=.46348E+02_JPRB | |
| 211 | 1 | VNP(2,5)=.35630E+02_JPRB | |
| 212 | 1 | VNP(3,5)=.33005E+01_JPRB | |
| 213 | 1 | VNP(4,5)=.18045E-01_JPRB | |
| 214 | 1 | VNP(5,5)=.88667E-05_JPRB | |
| 215 | 1 | VNP(1,6)=.47413E+01_JPRB | |
| 216 | 1 | VNP(2,6)=.16334E+01_JPRB | |
| 217 | 1 | VNP(3,6)=.48164E+00_JPRB | |
| 218 | 1 | VNP(4,6)=.56140E-02_JPRB | |
| 219 | 1 | VNP(5,6)=.67790E-04_JPRB | |
| 220 | |||
| 221 | ! Default values for cloud model - first index denotes spectral band: | ||
| 222 | ! 1 - solar | ||
| 223 | ! 2 - thermal | ||
| 224 | 1 | FCM_DEL_A(1) = 9.55496e-01_JPRB | |
| 225 | 1 | FCM_DEL_A(2) = 1.79152e+02_JPRB | |
| 226 | 1 | FCM_DEL_D(1) = 5.72959e+01_JPRB | |
| 227 | 1 | FCM_DEL_D(2) = 2.00159e+02_JPRB | |
| 228 | 1 | FCM_MU_A(1) = 8.66200e-01_JPRB | |
| 229 | 1 | FCM_MU_A(2) = 6.65532e-01_JPRB | |
| 230 | 1 | FCM_MU_D(1) = 1.00000e+00_JPRB | |
| 231 | 1 | FCM_MU_D(2) = 6.70517e-01_JPRB | |
| 232 | 1 | FCM_N_I = 1.00000e-01_JPRB | |
| 233 | 1 | FCM_N_L = 1.00000e-01_JPRB | |
| 234 | 1 | FCM_P_AI(1,0) = 1.48920e+01_JPRB | |
| 235 | 1 | FCM_P_AI(1,1) = -3.71379e+01_JPRB | |
| 236 | 1 | FCM_P_AI(1,2) = 0.00000e+00_JPRB | |
| 237 | 1 | FCM_P_AI(1,3) = 0.00000e+00_JPRB | |
| 238 | 1 | FCM_P_AI(2,0) = 9.76920e+00_JPRB | |
| 239 | 1 | FCM_P_AI(2,1) = -7.54840e+00_JPRB | |
| 240 | 1 | FCM_P_AI(2,2) = 0.00000e+00_JPRB | |
| 241 | 1 | FCM_P_AI(2,3) = 0.00000e+00_JPRB | |
| 242 | 1 | FCM_P_AL(1,0) = 1.42062e+00_JPRB | |
| 243 | 1 | FCM_P_AL(1,1) = -3.54037e+00_JPRB | |
| 244 | 1 | FCM_P_AL(1,2) = 1.73535e+00_JPRB | |
| 245 | 1 | FCM_P_AL(1,3) = 0.00000e+00_JPRB | |
| 246 | 1 | FCM_P_AL(2,0) = 5.08522e+00_JPRB | |
| 247 | 1 | FCM_P_AL(2,1) = -1.08144e+01_JPRB | |
| 248 | 1 | FCM_P_AL(2,2) = 5.16331e+00_JPRB | |
| 249 | 1 | FCM_P_AL(2,3) = 0.00000e+00_JPRB | |
| 250 | 1 | FCM_P_DI(1,0) = 1.30487e+01_JPRB | |
| 251 | 1 | FCM_P_DI(1,1) = 1.98670e+01_JPRB | |
| 252 | 1 | FCM_P_DI(1,2) = -5.08578e+01_JPRB | |
| 253 | 1 | FCM_P_DI(1,3) = 0.00000e+00_JPRB | |
| 254 | 1 | FCM_P_DI(2,0) = 1.23173e+01_JPRB | |
| 255 | 1 | FCM_P_DI(2,1) = -3.32808e+01_JPRB | |
| 256 | 1 | FCM_P_DI(2,2) = 5.77295e+01_JPRB | |
| 257 | 1 | FCM_P_DI(2,3) = -7.35728e+01_JPRB | |
| 258 | 1 | FCM_P_DL(1,0) = 4.68646e+00_JPRB | |
| 259 | 1 | FCM_P_DL(1,1) = -8.77652e+00_JPRB | |
| 260 | 1 | FCM_P_DL(1,2) = 2.12789e+00_JPRB | |
| 261 | 1 | FCM_P_DL(1,3) = 0.00000e+00_JPRB | |
| 262 | 1 | FCM_P_DL(2,0) = 4.52292e+00_JPRB | |
| 263 | 1 | FCM_P_DL(2,1) = -1.43327e+01_JPRB | |
| 264 | 1 | FCM_P_DL(2,2) = 1.16650e+01_JPRB | |
| 265 | 1 | FCM_P_DL(2,3) = 0.00000e+00_JPRB | |
| 266 | 1 | FCM_P_GI(1,0) = -2.01472e+00_JPRB | |
| 267 | 1 | FCM_P_GI(1,1) = 9.16284e+01_JPRB | |
| 268 | 1 | FCM_P_GI(1,2) = 7.82275e+01_JPRB | |
| 269 | 1 | FCM_P_GI(1,3) = 0.00000e+00_JPRB | |
| 270 | 1 | FCM_P_GI(2,0) = -7.55124e-01_JPRB | |
| 271 | 1 | FCM_P_GI(2,1) = -4.68837e+00_JPRB | |
| 272 | 1 | FCM_P_GI(2,2) = 6.72877e+01_JPRB | |
| 273 | 1 | FCM_P_GI(2,3) = 0.00000e+00_JPRB | |
| 274 | 1 | FCM_P_GL(1,0) = 8.82686e-01_JPRB | |
| 275 | 1 | FCM_P_GL(1,1) = -2.02445e+00_JPRB | |
| 276 | 1 | FCM_P_GL(1,2) = 1.45259e+00_JPRB | |
| 277 | 1 | FCM_P_GL(1,3) = 0.00000e+00_JPRB | |
| 278 | 1 | FCM_P_GL(2,0) = 1.93784e-01_JPRB | |
| 279 | 1 | FCM_P_GL(2,1) = -6.07676e-01_JPRB | |
| 280 | 1 | FCM_P_GL(2,2) = 5.37734e-01_JPRB | |
| 281 | 1 | FCM_P_GL(2,3) = 0.00000e+00_JPRB | |
| 282 | 1 | FCM_Q_AI(1,1) = 1.86588e+01_JPRB | |
| 283 | 1 | FCM_Q_AI(1,2) = 0.00000e+00_JPRB | |
| 284 | 1 | FCM_Q_AI(1,3) = 0.00000e+00_JPRB | |
| 285 | 1 | FCM_Q_AI(2,1) = 2.78063e+00_JPRB | |
| 286 | 1 | FCM_Q_AI(2,2) = 0.00000e+00_JPRB | |
| 287 | 1 | FCM_Q_AI(2,3) = 0.00000e+00_JPRB | |
| 288 | 1 | FCM_Q_AL(1,1) = -2.81210e+00_JPRB | |
| 289 | 1 | FCM_Q_AL(1,2) = 2.35758e+00_JPRB | |
| 290 | 1 | FCM_Q_AL(1,3) = 0.00000e+00_JPRB | |
| 291 | 1 | FCM_Q_AL(2,1) = -2.31252e+00_JPRB | |
| 292 | 1 | FCM_Q_AL(2,2) = 1.56848e+00_JPRB | |
| 293 | 1 | FCM_Q_AL(2,3) = 0.00000e+00_JPRB | |
| 294 | 1 | FCM_Q_DI(1,1) = 1.04582e+01_JPRB | |
| 295 | 1 | FCM_Q_DI(1,2) = 6.15641e+00_JPRB | |
| 296 | 1 | FCM_Q_DI(1,3) = 0.00000e+00_JPRB | |
| 297 | 1 | FCM_Q_DI(2,1) = 1.18370e+00_JPRB | |
| 298 | 1 | FCM_Q_DI(2,2) = 2.23187e+00_JPRB | |
| 299 | 1 | FCM_Q_DI(2,3) = 7.00424e+00_JPRB | |
| 300 | 1 | FCM_Q_DL(1,1) = -2.18205e+00_JPRB | |
| 301 | 1 | FCM_Q_DL(1,2) = 1.34831e+00_JPRB | |
| 302 | 1 | FCM_Q_DL(1,3) = 0.00000e+00_JPRB | |
| 303 | 1 | FCM_Q_DL(2,1) = -3.12310e+00_JPRB | |
| 304 | 1 | FCM_Q_DL(2,2) = 2.57959e+00_JPRB | |
| 305 | 1 | FCM_Q_DL(2,3) = 0.00000e+00_JPRB | |
| 306 | 1 | FCM_Q_GI(1,1) = 7.89999e+01_JPRB | |
| 307 | 1 | FCM_Q_GI(1,2) = 5.24929e+01_JPRB | |
| 308 | 1 | FCM_Q_GI(1,3) = 0.00000e+00_JPRB | |
| 309 | 1 | FCM_Q_GI(2,1) = 4.93801e+00_JPRB | |
| 310 | 1 | FCM_Q_GI(2,2) = 2.38096e+01_JPRB | |
| 311 | 1 | FCM_Q_GI(2,3) = 0.00000e+00_JPRB | |
| 312 | 1 | FCM_Q_GL(1,1) = -2.20830e+00_JPRB | |
| 313 | 1 | FCM_Q_GL(1,2) = 1.39832e+00_JPRB | |
| 314 | 1 | FCM_Q_GL(1,3) = 0.00000e+00_JPRB | |
| 315 | 1 | FCM_Q_GL(2,1) = -3.27770e+00_JPRB | |
| 316 | 1 | FCM_Q_GL(2,2) = 2.78446e+00_JPRB | |
| 317 | 1 | FCM_Q_GL(2,3) = 0.00000e+00_JPRB | |
| 318 | 1 | REXP_NEB = 8.0_JPRB | |
| 319 | |||
| 320 | ! 1.2 Modify default values according to LECMWF | ||
| 321 | |||
| 322 | IF (LECMWF) THEN | ||
| 323 | ELSE | ||
| 324 | ENDIF | ||
| 325 | |||
| 326 | ! Remark : the values or RII0 is calculated and not set up | ||
| 327 | |||
| 328 | ! ------------------------------------------------------------------ | ||
| 329 | |||
| 330 | !* 2. Modify default values. | ||
| 331 | ! ---------------------- | ||
| 332 | |||
| 333 | ! Ce qui concerne NAMPHY3 commente par MPL le 15.04.09 | ||
| 334 | !CALL POSNAM(NULNAM,'NAMPHY3') | ||
| 335 | !READ(NULNAM,NAMPHY3) | ||
| 336 | ! ------------------------------------------------------------------ | ||
| 337 | |||
| 338 | !* 3. Print final values. | ||
| 339 | ! ------------------- | ||
| 340 | |||
| 341 | 1 | WRITE(UNIT=KULOUT,FMT='('' COMMON YOMPHY3 '')') | |
| 342 | WRITE(UNIT=KULOUT,FMT='('' BSFSA = '',E10.4,'' BSFSI = '',E10.4 & | ||
| 343 | & ,'' BSFSN = '',E10.4,'' BSFTA = '',E10.4,'' BSFTI = '',E10.4,/& | ||
| 344 | & ,'' BSFTN = '',E10.4,'' EARRT = '',E10.4,'' EOASA = ''& | ||
| 345 | & ,E10.4,'' EOASI = '',E10.4,'' EOASN = '',E10.4 & | ||
| 346 | & ,'' EOATA = '',E10.4,/,'' EOATI = '',E10.4,'' EOATN = '',E10.4 & | ||
| 347 | & ,'' EODSA = '',E10.4,'' EODSI = '',E10.4,'' EODSN = '',E10.4 & | ||
| 348 | & ,'' EODTA = '',E10.4,/,'' EODTI = '',E10.4 & | ||
| 349 | & ,'' EODTN = '',E10.4,'' EORAY = '',E10.4,/& | ||
| 350 | & ,'' GCA(6) = '',6E11.4,/,'' GCB(6) = '',6E11.4,/,'' GCC(6) = ''& | ||
| 351 | & ,6E11.4,/,'' GCD4 = '',E11.4,'' GCE4 = '',E11.4,/& | ||
| 352 | & ,'' QCO2 = '',E10.4,'' QLIP0 = '',E10.4 & | ||
| 353 | & ,'' USAA = '',E10.4,'' USAI = '',E10.4,'' USAN = '',E10.4 & | ||
| 354 | & ,'' USBA = '',E10.4,/,'' USBI = '',E10.4,'' USBN = '',E10.4,/& | ||
| 355 | & ,'' VDP(5,6) = '',5E12.5,/,5(12X,5E12.5,/)& | ||
| 356 | & ,'' VNP(5,6) = '',5E12.5,/,5(12X,5E12.5,/)& | ||
| 357 | & ,'' QLIMI = '',E10.4 & | ||
| 358 | & ,'' GIREC1 = '',E10.4,'' GIREC2 = '',E10.4,'' GIREC3 = '',E10.4 & | ||
| 359 | & ,'' GIREC4 = '',E10.4 & | ||
| 360 | & )')& | ||
| 361 | 1 | & BSFSA,BSFSI,BSFSN,BSFTA,BSFTI,BSFTN,EARRT,EOASA,EOASI,EOASN,& | |
| 362 | 1 | & EOATA,EOATI,EOATN,EODSA,EODSI,EODSN,EODTA,EODTI,EODTN,EORAY,& | |
| 363 | 1 | & GCA,GCB,GCC,GCD4,GCE4,QCO2,QLIP0,USAA,USAI,USAN,USBA,USBI,USBN,& | |
| 364 | 2 | & VDP,VNP,QLIMI,GIREC1,GIREC2,GIREC3,GIREC4 | |
| 365 | |||
| 366 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (LCLSATUR) THEN |
| 367 | WRITE(KULOUT,'('' N_SPBAND = '',I2,'' REXP_NEB = '',F5.2)') & | ||
| 368 | ✗ | & N_SPBAND,REXP_NEB | |
| 369 | ✗ | WRITE(KULOUT,'('' FCM_DEL_A(:) ='',2(1X,ES12.5))') FCM_DEL_A(:) | |
| 370 | ✗ | WRITE(KULOUT,'('' FCM_DEL_D(:) ='',2(1X,ES12.5))') FCM_DEL_D(:) | |
| 371 | ✗ | WRITE(KULOUT,'('' FCM_MU_A (:) ='',2(1X,ES12.5))') FCM_MU_A (:) | |
| 372 | ✗ | WRITE(KULOUT,'('' FCM_MU_D (:) ='',2(1X,ES12.5))') FCM_MU_D (:) | |
| 373 | ✗ | WRITE(KULOUT,'('' FCM_N_I ='',1X,ES12.5)') FCM_N_I | |
| 374 | ✗ | WRITE(KULOUT,'('' FCM_N_L ='',1X,ES12.5)') FCM_N_L | |
| 375 | ✗ | WRITE(KULOUT,'('' FCM_P_AI (:,0) ='',2(1X,ES12.5))') FCM_P_AI (:,0) | |
| 376 | ✗ | WRITE(KULOUT,'('' FCM_P_AI (:,1) ='',2(1X,ES12.5))') FCM_P_AI (:,1) | |
| 377 | ✗ | WRITE(KULOUT,'('' FCM_P_AI (:,2) ='',2(1X,ES12.5))') FCM_P_AI (:,2) | |
| 378 | ✗ | WRITE(KULOUT,'('' FCM_P_AI (:,3) ='',2(1X,ES12.5))') FCM_P_AI (:,3) | |
| 379 | ✗ | WRITE(KULOUT,'('' FCM_P_AL (:,0) ='',2(1X,ES12.5))') FCM_P_AL (:,0) | |
| 380 | ✗ | WRITE(KULOUT,'('' FCM_P_AL (:,1) ='',2(1X,ES12.5))') FCM_P_AL (:,1) | |
| 381 | ✗ | WRITE(KULOUT,'('' FCM_P_AL (:,2) ='',2(1X,ES12.5))') FCM_P_AL (:,2) | |
| 382 | ✗ | WRITE(KULOUT,'('' FCM_P_AL (:,3) ='',2(1X,ES12.5))') FCM_P_AL (:,3) | |
| 383 | ✗ | WRITE(KULOUT,'('' FCM_P_DI (:,0) ='',2(1X,ES12.5))') FCM_P_DI (:,0) | |
| 384 | ✗ | WRITE(KULOUT,'('' FCM_P_DI (:,1) ='',2(1X,ES12.5))') FCM_P_DI (:,1) | |
| 385 | ✗ | WRITE(KULOUT,'('' FCM_P_DI (:,2) ='',2(1X,ES12.5))') FCM_P_DI (:,2) | |
| 386 | ✗ | WRITE(KULOUT,'('' FCM_P_DI (:,3) ='',2(1X,ES12.5))') FCM_P_DI (:,3) | |
| 387 | ✗ | WRITE(KULOUT,'('' FCM_P_DL (:,0) ='',2(1X,ES12.5))') FCM_P_DL (:,0) | |
| 388 | ✗ | WRITE(KULOUT,'('' FCM_P_DL (:,1) ='',2(1X,ES12.5))') FCM_P_DL (:,1) | |
| 389 | ✗ | WRITE(KULOUT,'('' FCM_P_DL (:,2) ='',2(1X,ES12.5))') FCM_P_DL (:,2) | |
| 390 | ✗ | WRITE(KULOUT,'('' FCM_P_DL (:,3) ='',2(1X,ES12.5))') FCM_P_DL (:,3) | |
| 391 | ✗ | WRITE(KULOUT,'('' FCM_P_GI (:,0) ='',2(1X,ES12.5))') FCM_P_GI (:,0) | |
| 392 | ✗ | WRITE(KULOUT,'('' FCM_P_GI (:,1) ='',2(1X,ES12.5))') FCM_P_GI (:,1) | |
| 393 | ✗ | WRITE(KULOUT,'('' FCM_P_GI (:,2) ='',2(1X,ES12.5))') FCM_P_GI (:,2) | |
| 394 | ✗ | WRITE(KULOUT,'('' FCM_P_GI (:,3) ='',2(1X,ES12.5))') FCM_P_GI (:,3) | |
| 395 | ✗ | WRITE(KULOUT,'('' FCM_P_GL (:,0) ='',2(1X,ES12.5))') FCM_P_GL (:,0) | |
| 396 | ✗ | WRITE(KULOUT,'('' FCM_P_GL (:,1) ='',2(1X,ES12.5))') FCM_P_GL (:,1) | |
| 397 | ✗ | WRITE(KULOUT,'('' FCM_P_GL (:,2) ='',2(1X,ES12.5))') FCM_P_GL (:,2) | |
| 398 | ✗ | WRITE(KULOUT,'('' FCM_P_GL (:,3) ='',2(1X,ES12.5))') FCM_P_GL (:,3) | |
| 399 | ✗ | WRITE(KULOUT,'('' FCM_Q_AI (:,1) ='',2(1X,ES12.5))') FCM_Q_AI (:,1) | |
| 400 | ✗ | WRITE(KULOUT,'('' FCM_Q_AI (:,2) ='',2(1X,ES12.5))') FCM_Q_AI (:,2) | |
| 401 | ✗ | WRITE(KULOUT,'('' FCM_Q_AI (:,3) ='',2(1X,ES12.5))') FCM_Q_AI (:,3) | |
| 402 | ✗ | WRITE(KULOUT,'('' FCM_Q_AL (:,1) ='',2(1X,ES12.5))') FCM_Q_AL (:,1) | |
| 403 | ✗ | WRITE(KULOUT,'('' FCM_Q_AL (:,2) ='',2(1X,ES12.5))') FCM_Q_AL (:,2) | |
| 404 | ✗ | WRITE(KULOUT,'('' FCM_Q_AL (:,3) ='',2(1X,ES12.5))') FCM_Q_AL (:,3) | |
| 405 | ✗ | WRITE(KULOUT,'('' FCM_Q_DI (:,1) ='',2(1X,ES12.5))') FCM_Q_DI (:,1) | |
| 406 | ✗ | WRITE(KULOUT,'('' FCM_Q_DI (:,2) ='',2(1X,ES12.5))') FCM_Q_DI (:,2) | |
| 407 | ✗ | WRITE(KULOUT,'('' FCM_Q_DI (:,3) ='',2(1X,ES12.5))') FCM_Q_DI (:,3) | |
| 408 | ✗ | WRITE(KULOUT,'('' FCM_Q_DL (:,1) ='',2(1X,ES12.5))') FCM_Q_DL (:,1) | |
| 409 | ✗ | WRITE(KULOUT,'('' FCM_Q_DL (:,2) ='',2(1X,ES12.5))') FCM_Q_DL (:,2) | |
| 410 | ✗ | WRITE(KULOUT,'('' FCM_Q_DL (:,3) ='',2(1X,ES12.5))') FCM_Q_DL (:,3) | |
| 411 | ✗ | WRITE(KULOUT,'('' FCM_Q_GI (:,1) ='',2(1X,ES12.5))') FCM_Q_GI (:,1) | |
| 412 | ✗ | WRITE(KULOUT,'('' FCM_Q_GI (:,2) ='',2(1X,ES12.5))') FCM_Q_GI (:,2) | |
| 413 | ✗ | WRITE(KULOUT,'('' FCM_Q_GI (:,3) ='',2(1X,ES12.5))') FCM_Q_GI (:,3) | |
| 414 | ✗ | WRITE(KULOUT,'('' FCM_Q_GL (:,1) ='',2(1X,ES12.5))') FCM_Q_GL (:,1) | |
| 415 | ✗ | WRITE(KULOUT,'('' FCM_Q_GL (:,2) ='',2(1X,ES12.5))') FCM_Q_GL (:,2) | |
| 416 | ✗ | WRITE(KULOUT,'('' FCM_Q_GL (:,3) ='',2(1X,ES12.5))') FCM_Q_GL (:,3) | |
| 417 | ENDIF | ||
| 418 | ! ------------------------------------------------------------------ | ||
| 419 | |||
| 420 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | IF (LHOOK) CALL DR_HOOK('SUPHY3',1,ZHOOK_HANDLE) |
| 421 | 1 | END SUBROUTINE SUPHY3 | |
| 422 |