| [3331] | 1 | MODULE YOESW |
|---|
| 2 | |
|---|
| 3 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
|---|
| 4 | |
|---|
| 5 | IMPLICIT NONE |
|---|
| 6 | |
|---|
| 7 | SAVE |
|---|
| 8 | |
|---|
| 9 | ! ------------------------------------------------------------------ |
|---|
| 10 | !* ** *YOESW* - COEFFICIENTS FOR SHORTWAVE RADIATION TRANSFER |
|---|
| 11 | ! ------------------------------------------------------------------ |
|---|
| 12 | |
|---|
| 13 | REAL(KIND=JPRB) :: APAD(6,3,7) |
|---|
| 14 | REAL(KIND=JPRB) :: BPAD(6,3,7) |
|---|
| 15 | REAL(KIND=JPRB) :: RRAY(6,6) |
|---|
| 16 | REAL(KIND=JPRB), ALLOCATABLE :: RSUN(:) |
|---|
| 17 | !$OMP THREADPRIVATE(RSUN) |
|---|
| 18 | REAL(KIND=JPRB) :: RPDH1 |
|---|
| 19 | REAL(KIND=JPRB) :: RPDU1 |
|---|
| 20 | REAL(KIND=JPRB) :: RPNH |
|---|
| 21 | REAL(KIND=JPRB) :: RPNU |
|---|
| 22 | REAL(KIND=JPRB) :: RSWCE(6) |
|---|
| 23 | REAL(KIND=JPRB) :: RSWCP(6) |
|---|
| 24 | REAL(KIND=JPRB) :: RTDH2O |
|---|
| 25 | REAL(KIND=JPRB) :: RTDUMG |
|---|
| 26 | REAL(KIND=JPRB) :: RTH2O |
|---|
| 27 | REAL(KIND=JPRB) :: RTUMG |
|---|
| 28 | REAL(KIND=JPRB) :: D(6,3) |
|---|
| 29 | REAL(KIND=JPRB) :: REXPO3(6,2,7) |
|---|
| 30 | INTEGER(KIND=JPIM) :: NEXPO3(6) |
|---|
| 31 | |
|---|
| 32 | REAL(KIND=JPRB) :: RYFWCA(6) |
|---|
| 33 | REAL(KIND=JPRB) :: RYFWCB(6) |
|---|
| 34 | REAL(KIND=JPRB) :: RYFWCC(6) |
|---|
| 35 | REAL(KIND=JPRB) :: RYFWCD(6) |
|---|
| 36 | REAL(KIND=JPRB) :: RYFWCE(6) |
|---|
| 37 | REAL(KIND=JPRB) :: RYFWCF(6) |
|---|
| 38 | |
|---|
| 39 | REAL(KIND=JPRB) :: REBCUA(6) |
|---|
| 40 | REAL(KIND=JPRB) :: REBCUB(6) |
|---|
| 41 | REAL(KIND=JPRB) :: REBCUC(6) |
|---|
| 42 | REAL(KIND=JPRB) :: REBCUD(6) |
|---|
| 43 | REAL(KIND=JPRB) :: REBCUE(6) |
|---|
| 44 | REAL(KIND=JPRB) :: REBCUF(6) |
|---|
| 45 | REAL(KIND=JPRB) :: REBCUG(16) |
|---|
| 46 | REAL(KIND=JPRB) :: REBCUH(16) |
|---|
| 47 | REAL(KIND=JPRB) :: REBCUI(6) |
|---|
| 48 | REAL(KIND=JPRB) :: REBCUJ(6) |
|---|
| 49 | |
|---|
| 50 | REAL(KIND=JPRB) :: RASWCA(6) |
|---|
| 51 | REAL(KIND=JPRB) :: RASWCB(6) |
|---|
| 52 | REAL(KIND=JPRB) :: RASWCC(6) |
|---|
| 53 | REAL(KIND=JPRB) :: RASWCD(6) |
|---|
| 54 | REAL(KIND=JPRB) :: RASWCE(6) |
|---|
| 55 | REAL(KIND=JPRB) :: RASWCF(6) |
|---|
| 56 | |
|---|
| 57 | REAL(KIND=JPRB) :: RFUETA(16,3),RFUETB(16,4), RFUETC(16,4) |
|---|
| 58 | REAL(KIND=JPRB) :: RFULIO(16,3) |
|---|
| 59 | REAL(KIND=JPRB) :: RHSAVI(16,3) |
|---|
| 60 | REAL(KIND=JPRB) :: RLILIA(16,5),RLILIB(16,4) |
|---|
| 61 | |
|---|
| 62 | REAL(KIND=JPRB) :: RFLAA0(6) |
|---|
| 63 | REAL(KIND=JPRB) :: RFLAA1(6) |
|---|
| 64 | REAL(KIND=JPRB) :: RFLBB0(6) |
|---|
| 65 | REAL(KIND=JPRB) :: RFLBB1(6) |
|---|
| 66 | REAL(KIND=JPRB) :: RFLBB2(6) |
|---|
| 67 | REAL(KIND=JPRB) :: RFLBB3(6) |
|---|
| 68 | REAL(KIND=JPRB) :: RFLCC0(6) |
|---|
| 69 | REAL(KIND=JPRB) :: RFLCC1(6) |
|---|
| 70 | REAL(KIND=JPRB) :: RFLCC2(6) |
|---|
| 71 | REAL(KIND=JPRB) :: RFLCC3(6) |
|---|
| 72 | |
|---|
| 73 | REAL(KIND=JPRB) :: RFUAA0(6) |
|---|
| 74 | REAL(KIND=JPRB) :: RFUAA1(6) |
|---|
| 75 | REAL(KIND=JPRB) :: RFUBB0(6) |
|---|
| 76 | REAL(KIND=JPRB) :: RFUBB1(6) |
|---|
| 77 | REAL(KIND=JPRB) :: RFUBB2(6) |
|---|
| 78 | REAL(KIND=JPRB) :: RFUBB3(6) |
|---|
| 79 | REAL(KIND=JPRB) :: RFUCC0(6) |
|---|
| 80 | REAL(KIND=JPRB) :: RFUCC1(6) |
|---|
| 81 | REAL(KIND=JPRB) :: RFUCC2(6) |
|---|
| 82 | REAL(KIND=JPRB) :: RFUCC3(6) |
|---|
| 83 | REAL(KIND=JPRB) :: RFLDD0(6) |
|---|
| 84 | REAL(KIND=JPRB) :: RFLDD1(6) |
|---|
| 85 | REAL(KIND=JPRB) :: RFLDD2(6) |
|---|
| 86 | REAL(KIND=JPRB) :: RFLDD3(6) |
|---|
| 87 | |
|---|
| 88 | REAL(KIND=JPRB) :: RSUSHE(6) |
|---|
| 89 | REAL(KIND=JPRB) :: RSUSHF(6) |
|---|
| 90 | REAL(KIND=JPRB) :: RSUSHH(6) |
|---|
| 91 | REAL(KIND=JPRB) :: RSUSHK(6) |
|---|
| 92 | REAL(KIND=JPRB) :: RSUSHA(6) |
|---|
| 93 | REAL(KIND=JPRB) :: RSUSHG(6) |
|---|
| 94 | REAL(KIND=JPRB) :: RSUSHFA(4) |
|---|
| 95 | REAL(KIND=JPRB) :: RSUSHC |
|---|
| 96 | REAL(KIND=JPRB) :: RSUSHD |
|---|
| 97 | |
|---|
| 98 | REAL(KIND=JPRB) :: REFFIA |
|---|
| 99 | REAL(KIND=JPRB) :: REFFIB |
|---|
| 100 | REAL(KIND=JPRB) :: RTIW |
|---|
| 101 | REAL(KIND=JPRB) :: RRIW |
|---|
| 102 | REAL(KIND=JPRB) :: RROMA(6) |
|---|
| 103 | REAL(KIND=JPRB) :: RROMB(6) |
|---|
| 104 | REAL(KIND=JPRB) :: RRASY(6) |
|---|
| 105 | |
|---|
| 106 | REAL(KIND=JPRB) :: RHSRA(6) |
|---|
| 107 | REAL(KIND=JPRB) :: RHSRB(6) |
|---|
| 108 | REAL(KIND=JPRB) :: RHSRC(6) |
|---|
| 109 | REAL(KIND=JPRB) :: RHSRD(6) |
|---|
| 110 | REAL(KIND=JPRB) :: RHSRE(6) |
|---|
| 111 | REAL(KIND=JPRB) :: RHSRF(6) |
|---|
| 112 | REAL(KIND=JPRB) :: RHSRTA |
|---|
| 113 | REAL(KIND=JPRB) :: RHSRTB |
|---|
| 114 | |
|---|
| 115 | REAL(KIND=JPRB) :: RTAUA(6,6) |
|---|
| 116 | REAL(KIND=JPRB) :: RPIZA(6,6) |
|---|
| 117 | REAL(KIND=JPRB) :: RCGA(6,6) |
|---|
| 118 | REAL(KIND=JPRB) :: RAER(6,6) |
|---|
| 119 | |
|---|
| 120 | REAL(KIND=JPRB) :: RTWEIGHT(8) |
|---|
| 121 | REAL(KIND=JPRB) :: RWEIGS(6) |
|---|
| 122 | REAL(KIND=JPRB) :: RWEIGV(6) |
|---|
| 123 | |
|---|
| 124 | REAL(KIND=JPRB) :: RWEIGHT(6,8) |
|---|
| 125 | INTEGER(KIND=JPIM) :: NMPSRTM(14), NTYPS |
|---|
| 126 | |
|---|
| 127 | REAL(KIND=JPRB) :: RADJUST |
|---|
| 128 | |
|---|
| 129 | ! * E.C.M.W.F. PHYSICS PACKAGE * |
|---|
| 130 | |
|---|
| 131 | ! J.-J. MORCRETTE E.C.M.W.F. 89/07/14 |
|---|
| 132 | |
|---|
| 133 | ! NAME TYPE PURPOSE |
|---|
| 134 | ! ---- : ---- : --------------------------------------------------- |
|---|
| 135 | ! APAD : REAL PADE APPROXIMANTS NUMERATOR |
|---|
| 136 | ! BPAD : REAL PADE APPROXIMANTS DENOMINATOR |
|---|
| 137 | ! D : REAL TRANSMISSION LIMIT FOR INFINITE ABSORBER AMOUNT |
|---|
| 138 | ! RRAY : REAL RAYLEIGH SCATTERING COEFFICIENTS |
|---|
| 139 | ! RSUN : REAL SOLAR FRACTION IN SPECTRAL INTERVALS |
|---|
| 140 | ! RPDH1 : 1 + EXPONENT PRESSURE DEPENDENCE H2O |
|---|
| 141 | ! RPDU1 : 1 + EXPONENT PRESSURE DEPENDENCE UNIFORMLY MIXED GASES |
|---|
| 142 | ! RPNH : REFERENCE PRESSURE FACTOR FOR H2O |
|---|
| 143 | ! RPNU : REFERENCE PRESSURE FACTOR FOR UNIFORMLY MIXED GASES |
|---|
| 144 | ! RSWCE : E-TYPE, H2O CONTINUUM ABSORPTION COEFFICIENT |
|---|
| 145 | ! RSWCP : P-TYPE, H2O CONTINUUM ABSORPTION COEFFICIENT |
|---|
| 146 | ! RTDH2O: EXPONENT TEMPERATURE DEPENDENCE H2O |
|---|
| 147 | ! RTDUMG: EXPONENT TEMPERATURE DEPENDENCE UNIFORMLY MIXED GASES |
|---|
| 148 | ! RTH2O : REFERENCE TEMPERATURE H2O |
|---|
| 149 | ! RTUMG : REFERENCE TEMPERATURE UNIFORMLY MIXED GASES |
|---|
| 150 | ! ----------------------------------------------------------------- |
|---|
| 151 | |
|---|
| 152 | ! * E.C.M.W.F. PHYSICS PACKAGE * |
|---|
| 153 | |
|---|
| 154 | ! J.-J. MORCRETTE E.C.M.W.F. 89/07/14 |
|---|
| 155 | |
|---|
| 156 | ! NAME TYPE PURPOSE |
|---|
| 157 | ! ---- : ---- : --------------------------------------------------- |
|---|
| 158 | !* FOUQUART (1987) WATER CLOUD OPTICAL PROPERTIES |
|---|
| 159 | |
|---|
| 160 | ! RYFWCA : REAL : C1 IN OPTICAL THICKNESS FORMULA |
|---|
| 161 | ! RYFWCB : REAL : C2 IN OPTICAL THICKNESS FORMULA |
|---|
| 162 | ! RYFWCC : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
|---|
| 163 | ! RYFWCD : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
|---|
| 164 | ! RYFWCE : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
|---|
| 165 | ! RYFWCF : REAL : ASSYMETRY FACTOR |
|---|
| 166 | |
|---|
| 167 | !* SLINGO (1989) WATER CLOUD OPTICAL PROPERTIES |
|---|
| 168 | |
|---|
| 169 | ! RASWCA : REAL : C1 IN OPTICAL THICKNESS FORMULA |
|---|
| 170 | ! RASWCB : REAL : C2 IN OPTICAL THICKNESS FORMULA |
|---|
| 171 | ! RASWCC : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
|---|
| 172 | ! RASWCD : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
|---|
| 173 | ! RASWCE : REAL : SINGLE SCATTERING ALBEDO PARAMETER |
|---|
| 174 | ! RASWCF : REAL : ASSYMETRY FACTOR |
|---|
| 175 | |
|---|
| 176 | !* SAVIJARVI (1998) WATER CLOUD OPTICAL PROPERTIES (RRTM) |
|---|
| 177 | |
|---|
| 178 | ! RHSAVI : REAL : MASS ABSORPTION COEFFICIENTS (POLYNOMIAL DEVELOPM) |
|---|
| 179 | |
|---|
| 180 | !* LINDNER,LI (2000) WATER CLOUD OPTICAL PROPERTIES (RRTM) |
|---|
| 181 | |
|---|
| 182 | ! RLILIA : REAL : MASS ABSORPTION COEFFICIENTS (POLYNOMIAL DEVELOPM) |
|---|
| 183 | ! RLILIB : REAL : 1-SSA COEFFICIENTS (POLYNOMIAL DEVELOPM) |
|---|
| 184 | |
|---|
| 185 | !* ICE CLOUD OPTICAL PROPERTIES DERIVED FROM EBERT-CURRY (1992) |
|---|
| 186 | |
|---|
| 187 | ! REBCUA : REAL : C1 IN OPTICAL THICKNESS FORMULA |
|---|
| 188 | ! REBCUB : REAL : C2 IN OPTICAL THICKNESS FORMULA |
|---|
| 189 | ! REBCUC : REAL : 1-C3 IN SINGLE SCATTERING ALBEDO FORMULA |
|---|
| 190 | ! REBCUD : REAL : C4 IN SINGLE SCATTERING ALBEDO FORMULA |
|---|
| 191 | ! REBCUE : REAL : C5 IN ASSYMETRY FACTOR FORMULA |
|---|
| 192 | ! REBCUF : REAL : C6 IN ASSYMETRY FACTOR FORMULA |
|---|
| 193 | ! REBCUG : REAL : C7 IN MASS ABSORPTION COEFFICIENT FORMULA |
|---|
| 194 | ! REBCUH : REAL : C8 IN MASS ABSORPTION COEFFICIENT FORMULA |
|---|
| 195 | ! REBCUI : REAL : C7 IN MASS ABSORPTION COEFFICIENT SPECTRAL FORMULA |
|---|
| 196 | ! REBCUJ : REAL : C8 IN MASS ABSORPTION COEFFICIENT SPECTRAL FORMULA |
|---|
| 197 | |
|---|
| 198 | !* ICE CLOUD OPTICAL PROPERTIES DERIVED FROM SUN-SHINE (1995) |
|---|
| 199 | |
|---|
| 200 | ! RSHSUE : REAL : E IN SINGLE SCATTERING ALBEDO FORMULA |
|---|
| 201 | ! RSHSUF : REAL : F IN SINGLE SCATTERING ALBEDO FORMULA |
|---|
| 202 | ! RSHSUH : REAL : H IN ASSYMETRY FACTOR FORMULA |
|---|
| 203 | ! RSHSUK : REAL : K IN ASSYMETRY FACTOR FORMULA |
|---|
| 204 | ! RSHSUA : REAL : ALPHA IN SSA CORRECTION FACTOR FORMULA |
|---|
| 205 | ! RSHSUG : REAL : GAMMA IN ASSYMETRY CORRECTION FACTOR FORMULA |
|---|
| 206 | ! RSHSUFA: REAL : COEFFICIENTS IN TEMPERATURE CORRECTION FACTOR |
|---|
| 207 | |
|---|
| 208 | ! REFFIA : REAL : C9 IN EFFECTIVE RADIUS FORMULA |
|---|
| 209 | ! REFFIB : REAL : C10 IN EFFECTIVE RADIUS FORMULA |
|---|
| 210 | |
|---|
| 211 | !* ICE CLOUD OPTICAL PROPERTIES DERIVED FROM FU-LIOU (1993) |
|---|
| 212 | |
|---|
| 213 | ! RFULIO : REAL : COEFFICIENTS IN EXPRESSION FOR LW EXTINCTION COEFF. |
|---|
| 214 | ! RFLAA : REAL : COEFFICIENTS IN EXPRESSION FOR SW EXTINCTION COEFF. |
|---|
| 215 | ! RFLBB : REAL : COEFFICIENTS IN EXPRESSION FOR SW SINGLE SCATT.ALB. |
|---|
| 216 | ! RFLCC : REAL : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR |
|---|
| 217 | ! RFLDD : REAL : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR |
|---|
| 218 | |
|---|
| 219 | !* ICE CLOUD OPTICAL PROPERTIES DERIVED FROM FU (1996) & FU ET AL. (1998) |
|---|
| 220 | |
|---|
| 221 | ! RFUETA : REAL : COEFFICIENTS IN EXPRESSION FOR LW EXTINCTION COEFF. |
|---|
| 222 | ! RFUAA : REAL : COEFFICIENTS IN EXPRESSION FOR SW EXTINCTION COEFF. |
|---|
| 223 | ! RFUBB : REAL : COEFFICIENTS IN EXPRESSION FOR SW SINGLE SCATT.ALB. |
|---|
| 224 | ! RFUCC : REAL : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR |
|---|
| 225 | |
|---|
| 226 | !* TRANSITION BETWEEN LIQUID AND SOLID WATER |
|---|
| 227 | |
|---|
| 228 | ! RTIW : REAL : TEMPERATURE THRESHOLD |
|---|
| 229 | ! RRIW : REAL : TRANSITION RANGE |
|---|
| 230 | |
|---|
| 231 | !* RAIN OPTICAL PROPERTIES FROM SAVIJARVI (1996) |
|---|
| 232 | |
|---|
| 233 | ! RROMA : REAL : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO |
|---|
| 234 | ! RROMB : REAL : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO |
|---|
| 235 | ! RRASY : REAL : COEFFICIENTS FOR ASSYMETRY FACTOR |
|---|
| 236 | ! RHSRA : REAL : COEFFICIENTS FOR OPTICAL THICKNESS |
|---|
| 237 | ! RHSRB : REAL : COEFFICIENTS FOR OPTICAL THICKNESS |
|---|
| 238 | ! RHSRC : REAL : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO |
|---|
| 239 | ! RHSRD : REAL : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO |
|---|
| 240 | ! RHSRE : REAL : COEFFICIENTS FOR ASSYMETRY FACTOR |
|---|
| 241 | ! RHSRF : REAL : COEFFICIENTS FOR ASSYMETRY FACTOR |
|---|
| 242 | ! RHSRTA : REAL : COEFFICIENTS FOR OPTICAL THICKNESS |
|---|
| 243 | ! RHSRTB : REAL : COEFFICIENTS FOR OPTICAL THICKNESS |
|---|
| 244 | ! ----------------------------------------------------------------- |
|---|
| 245 | |
|---|
| 246 | ! * E.C.M.W.F. PHYSICS PACKAGE * |
|---|
| 247 | |
|---|
| 248 | ! J.-J. MORCRETTE E.C.M.W.F. 89/07/14 |
|---|
| 249 | |
|---|
| 250 | ! NAME TYPE PURPOSE |
|---|
| 251 | ! ---- : ---- : ------- |
|---|
| 252 | ! RTAUA : REAL S.W. NORMALIZED OPTICAL THICKNESS AT 0.55 MICRON |
|---|
| 253 | ! RPIZA : REAL S.W. SINGLE SCATTERING ALBEDO |
|---|
| 254 | ! RCGA : REAL S.W. ASSYMETRY FACTOR |
|---|
| 255 | ! RAER : REAL L.W. ABSORPTION COEFFICIENTS |
|---|
| 256 | ! ----------------------------------------------------------------- |
|---|
| 257 | |
|---|
| 258 | ! * E.C.M.W.F. PHYSICS PACKAGE * |
|---|
| 259 | |
|---|
| 260 | ! J.-J. MORCRETTE E.C.M.W.F. 89/07/14 |
|---|
| 261 | |
|---|
| 262 | ! NAME TYPE PURPOSE |
|---|
| 263 | ! ---- : ---- : ------- |
|---|
| 264 | !RTWEIGHT: REAL S.W. INTEGRATED WEIGHT |
|---|
| 265 | ! NMPSRTM: INTEGER : Indices for mapping SW[1:6] albedo into SRTM[1:14] |
|---|
| 266 | ! ----------------------------------------------------------------- |
|---|
| 267 | !$OMP THREADPRIVATE(apad,bpad,d,nexpo3,nmpsrtm,ntyps,radjust,raer,raswca) |
|---|
| 268 | !$OMP THREADPRIVATE(raswcb,raswcc,raswcd,raswce,raswcf,rcga,rebcua,rebcub) |
|---|
| 269 | !$OMP THREADPRIVATE(rebcuc,rebcud,rebcue,rebcuf,rebcug,rebcuh,rebcui,rebcuj) |
|---|
| 270 | !$OMP THREADPRIVATE(reffia,reffib,rexpo3,rflaa0,rflaa1,rflbb0,rflbb1,rflbb2) |
|---|
| 271 | !$OMP THREADPRIVATE(rflbb3,rflcc0,rflcc1,rflcc2,rflcc3,rfldd0,rfldd1,rfldd2) |
|---|
| 272 | !$OMP THREADPRIVATE(rfldd3,rfuaa0,rfuaa1,rfubb0,rfubb1,rfubb2,rfubb3,rfucc0) |
|---|
| 273 | !$OMP THREADPRIVATE(rfucc1,rfucc2,rfucc3,rfueta,rfuetb,rfuetc,rfulio,rhsavi) |
|---|
| 274 | !$OMP THREADPRIVATE(rhsra,rhsrb,rhsrc,rhsrd,rhsre,rhsrf,rhsrta,rhsrtb,rlilia) |
|---|
| 275 | !$OMP THREADPRIVATE(rlilib,rpdh1,rpdu1,rpiza,rpnh,rpnu,rrasy,rray,rriw,rroma) |
|---|
| 276 | !$OMP THREADPRIVATE(rromb,rsusha,rsushc,rsushd,rsushe,rsushf,rsushfa,rsushg) |
|---|
| 277 | !$OMP THREADPRIVATE(rsushh,rsushk,rswce,rswcp,rtaua,rtdh2o,rtdumg,rth2o,rtiw) |
|---|
| 278 | !$OMP THREADPRIVATE(rtumg,rtweight,rweight,rweigs,rweigv,ryfwca,ryfwcb,ryfwcc,ryfwcd,ryfwce,ryfwcf) |
|---|
| 279 | END MODULE YOESW |
|---|
