| [1989] | 1 | SUBROUTINE SUSWN (KTSW, KSW) |
|---|
| 2 | |
|---|
| 3 | !**** *SUSW* - INITIALIZE COMMON YOESW |
|---|
| 4 | |
|---|
| 5 | ! PURPOSE. |
|---|
| 6 | ! -------- |
|---|
| 7 | ! INITIALIZE YOESW, THE COMMON THAT CONTAINS COEFFICIENTS |
|---|
| 8 | ! NEEDED TO RUN THE SHORTWAVE RADIATION SUBROUTINES |
|---|
| 9 | |
|---|
| 10 | !** INTERFACE. |
|---|
| 11 | ! ---------- |
|---|
| 12 | ! *CALL* *SUSW |
|---|
| 13 | |
|---|
| 14 | ! EXPLICIT ARGUMENTS : |
|---|
| 15 | ! -------------------- |
|---|
| 16 | ! NONE |
|---|
| 17 | |
|---|
| 18 | ! IMPLICIT ARGUMENTS : |
|---|
| 19 | ! -------------------- |
|---|
| 20 | ! COMMON YOESW |
|---|
| 21 | |
|---|
| 22 | ! METHOD. |
|---|
| 23 | ! ------- |
|---|
| 24 | ! SEE DOCUMENTATION |
|---|
| 25 | |
|---|
| 26 | ! EXTERNALS. |
|---|
| 27 | ! ---------- |
|---|
| 28 | |
|---|
| 29 | ! REFERENCE. |
|---|
| 30 | ! ---------- |
|---|
| 31 | ! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS |
|---|
| 32 | |
|---|
| 33 | ! AUTHOR. |
|---|
| 34 | ! ------- |
|---|
| 35 | ! JEAN-JACQUES MORCRETTE *ECMWF* |
|---|
| 36 | |
|---|
| 37 | ! MODIFICATIONS. |
|---|
| 38 | ! -------------- |
|---|
| 39 | ! ORIGINAL : 88-12-15 |
|---|
| 40 | ! 97-04-16 JJ Morcrette 2 and 4 interval spectral resolution |
|---|
| 41 | ! 00-10-24 JJ Morcrette sea-ice albedo revisited |
|---|
| 42 | ! 00-12-14 JJ Morcrette |
|---|
| 43 | ! and Ph.Dubuisson B.Bonnel 6 spectral interval resolution |
|---|
| 44 | ! 01-04-17 Ph.Dubuisson, B.Bonnel, JJ.Morcrette 6 sp.int.resolu'n |
|---|
| 45 | ! 01-06-28 B.Bonnel, JJ.Morcrette, Ph.Dubuisson Rayleigh (2/4/6) |
|---|
| 46 | ! 01-11-05 Ph.Dubuisson, JJMorcrette (new 2 intervals for TL/AD) |
|---|
| 47 | ! 03-04-24 JJ Morcrette SRTM with 14 spectral intervals |
|---|
| 48 | ! M.Hamrud 01-Oct-2003 CY28 Cleaning |
|---|
| 49 | |
|---|
| 50 | ! ------------------------------------------------------------------ |
|---|
| 51 | |
|---|
| 52 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
|---|
| 53 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
|---|
| 54 | |
|---|
| 55 | USE YOMCST , ONLY : RG |
|---|
| 56 | USE YOERAD , ONLY : NTSW ,LONEWSW ,LRAYL |
|---|
| 57 | USE YOESW , ONLY : APAD ,BPAD ,RRAY ,RSUN ,& |
|---|
| 58 | & RPDH1 ,RPDU1 ,RPNH ,RPNU ,RSWCE ,& |
|---|
| 59 | & RSWCP ,RTDH2O ,RTDUMG ,RTH2O ,RTUMG ,& |
|---|
| 60 | & D ,RROMA ,RROMB ,RRASY ,RHSRA ,& |
|---|
| 61 | & RHSRB ,RHSRC ,RHSRD ,RHSRE ,RHSRF ,& |
|---|
| 62 | & RHSRTA ,RHSRTB ,REXPO3 ,& |
|---|
| 63 | & RADJUST ,NEXPO3 |
|---|
| 64 | |
|---|
| 65 | IMPLICIT NONE |
|---|
| 66 | |
|---|
| 67 | INTEGER(KIND=JPIM),INTENT(IN) :: KTSW |
|---|
| 68 | INTEGER(KIND=JPIM),INTENT(IN) :: KSW |
|---|
| 69 | INTEGER(KIND=JPIM) :: IMPSRTM(14) |
|---|
| 70 | |
|---|
| 71 | ! ---------------------------------------------------------------- |
|---|
| 72 | REAL(KIND=JPRB) :: ZAPAD2(2,3,7) , ZBPAD2(2,3,7) , ZD2(2,3)& |
|---|
| 73 | & , ZRAY2(2,6) , ZSUN2(2) , ZSWCE2(2) , ZSWCP2(2)& |
|---|
| 74 | & , ZROMA2(2) , ZROMB2(2) , ZRASY2(2)& |
|---|
| 75 | & , ZRA2(2) , ZRB2(2) , ZRC2(2)& |
|---|
| 76 | & , ZRD2(2) , ZRE2(2) , ZRF2(2) |
|---|
| 77 | |
|---|
| 78 | REAL(KIND=JPRB) :: ZAPAD4(4,3,7) , ZBPAD4(4,3,7) , ZD4(4,3)& |
|---|
| 79 | & , ZRAY4(4,6) , ZSUN4(4) , ZSWCE4(4) , ZSWCP4(4)& |
|---|
| 80 | & , ZROMA4(4) , ZROMB4(4) , ZRASY4(4)& |
|---|
| 81 | & , ZRA4(4) , ZRB4(4) , ZRC4(4)& |
|---|
| 82 | & , ZRD4(4) , ZRE4(4) , ZRF4(4) |
|---|
| 83 | |
|---|
| 84 | REAL(KIND=JPRB) :: ZAPAD6(6,3,7) , ZBPAD6(6,3,7) , ZD6(6,3)& |
|---|
| 85 | & , ZRAY6(6,6) , ZSUN6(6) , ZSWCE6(6) , ZSWCP6(6)& |
|---|
| 86 | & , ZROMA6(6) , ZROMB6(6) , ZRASY6(6)& |
|---|
| 87 | & , ZRA6(6) , ZRB6(6) , ZRC6(6)& |
|---|
| 88 | & , ZRD6(6) , ZRE6(6) , ZRF6(6)& |
|---|
| 89 | & , ZRBY6(6,6) |
|---|
| 90 | |
|---|
| 91 | |
|---|
| 92 | INTEGER(KIND=JPIM) :: JC3, JC6, JI, JJ, JW, I_K, ISW6 |
|---|
| 93 | |
|---|
| 94 | REAL(KIND=JPRB) :: ZH2O, ZPDHNIS, ZPDH2O, ZPDUNIS,& |
|---|
| 95 | & ZPDUMG, ZPRHNIS, ZPRH2O, & |
|---|
| 96 | & ZPRUNIS, ZPRUMG, ZRTO1, ZRTO2, ZTDHNIS, & |
|---|
| 97 | & ZTDUNIS, ZTHNIS, & |
|---|
| 98 | & ZTUNIS, ZUMG, ZADJUS2, ZADJUS4, ZADJUS6 |
|---|
| 99 | |
|---|
| 100 | REAL(KIND=JPRB) :: ZTDH2IS, ZTDU2IS, ZTH2IS, ZTU2IS, ZPDH2IS, ZPDU2IS, ZPRH2IS, ZPRU2IS |
|---|
| 101 | REAL(KIND=JPRB) :: ZTDH4IS, ZTDU4IS, ZTH4IS, ZTU4IS, ZPDH4IS, ZPDU4IS, ZPRH4IS, ZPRU4IS |
|---|
| 102 | REAL(KIND=JPRB) :: ZTDH6IS, ZTDU6IS, ZTH6IS, ZTU6IS, ZPDH6IS, ZPDU6IS, ZPRH6IS, ZPRU6IS |
|---|
| 103 | REAL(KIND=JPRB) :: ZHOOK_HANDLE |
|---|
| 104 | |
|---|
| 105 | #include "abor1.intfb.h" |
|---|
| 106 | |
|---|
| 107 | ! ---------------------------------------------------------------- |
|---|
| 108 | |
|---|
| 109 | !* 1. CLEAR-SKY ABSORPTION COEFFICIENTS FOR N SPECTRAL INTERVALS |
|---|
| 110 | ! -------------------------------------------------------- |
|---|
| 111 | |
|---|
| 112 | !* spectrally independent |
|---|
| 113 | |
|---|
| 114 | !-- New coefficients |
|---|
| 115 | |
|---|
| 116 | IF (LHOOK) CALL DR_HOOK('SUSWN',0,ZHOOK_HANDLE) |
|---|
| 117 | ZTDHNIS = 0.450_JPRB |
|---|
| 118 | ZTDUNIS = 0.375_JPRB |
|---|
| 119 | ZTHNIS = 273._JPRB |
|---|
| 120 | ZTUNIS = 273._JPRB |
|---|
| 121 | ZPDHNIS = 0.90_JPRB |
|---|
| 122 | ZPDUNIS = 0.75_JPRB |
|---|
| 123 | ZPRHNIS = 101300._JPRB |
|---|
| 124 | ZPRUNIS = 101300._JPRB |
|---|
| 125 | ISW6=6 |
|---|
| 126 | |
|---|
| 127 | ! ---------------------------------------------------------------- |
|---|
| 128 | |
|---|
| 129 | !* 1.1 COEFFICIENTS FOR TWO SPECTRAL INTERVALS |
|---|
| 130 | ! --------------------------------------- |
|---|
| 131 | |
|---|
| 132 | !* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING |
|---|
| 133 | ! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9 |
|---|
| 134 | ! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2) |
|---|
| 135 | ! O3 : unchanged in interval 1, from HITRAN 92 in interval 4 |
|---|
| 136 | |
|---|
| 137 | ZADJUS2 = 1.0_JPRB |
|---|
| 138 | ZTDH2IS = 0.450_JPRB |
|---|
| 139 | ZTDU2IS = 0.375_JPRB |
|---|
| 140 | ZTH2IS = 273._JPRB |
|---|
| 141 | ZTU2IS = 273._JPRB |
|---|
| 142 | ZPDH2IS = 0.90_JPRB |
|---|
| 143 | ZPDU2IS = 0.75_JPRB |
|---|
| 144 | ZPRH2IS = 101300._JPRB |
|---|
| 145 | ZPRU2IS = 101300._JPRB |
|---|
| 146 | |
|---|
| 147 | !* 1st spectral interval: U.V. and Visible (0.25 - 0.69 Micron) |
|---|
| 148 | |
|---|
| 149 | ZSUN2(1) = 0.459760_JPRB |
|---|
| 150 | |
|---|
| 151 | !-- New coefficients |
|---|
| 152 | ZD2(1,:)= (/ 0.800000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB /) |
|---|
| 153 | |
|---|
| 154 | ZAPAD2(1, 1, :) = (/& |
|---|
| 155 | & 0.1095121E+03_JPRB,& |
|---|
| 156 | & 0.8906647E+03_JPRB,& |
|---|
| 157 | & 0.1030601E+03_JPRB,& |
|---|
| 158 | & 0.000000000E-00_JPRB,& |
|---|
| 159 | & 0.000000000E-00_JPRB,& |
|---|
| 160 | & 0.000000000E-00_JPRB,& |
|---|
| 161 | & 0.000000000E-00_JPRB/) |
|---|
| 162 | ZAPAD2(1, 2, :) = (/& |
|---|
| 163 | & 0.5124541E+03_JPRB,& |
|---|
| 164 | & 0.3478800E+04_JPRB,& |
|---|
| 165 | & 0.1244272E+04_JPRB,& |
|---|
| 166 | & 0.000000000E-00_JPRB,& |
|---|
| 167 | & 0.000000000E-00_JPRB,& |
|---|
| 168 | & 0.000000000E-00_JPRB,& |
|---|
| 169 | & 0.000000000E-00_JPRB/) |
|---|
| 170 | ZAPAD2(1, 3, :) = (/& |
|---|
| 171 | & 0.4766857E-04_JPRB,& |
|---|
| 172 | & 0.4019806E-02_JPRB,& |
|---|
| 173 | & 0.1936247E+00_JPRB,& |
|---|
| 174 | & 0.1886522E+02_JPRB,& |
|---|
| 175 | & 0.6498215E+02_JPRB,& |
|---|
| 176 | & 0.1799387E+02_JPRB,& |
|---|
| 177 | & 0.000000000E+00_JPRB/) |
|---|
| 178 | |
|---|
| 179 | ZBPAD2(1, 1, :) = (/& |
|---|
| 180 | & 0.1095121E+03_JPRB,& |
|---|
| 181 | & 0.8919103E+03_JPRB,& |
|---|
| 182 | & 0.1142218E+03_JPRB,& |
|---|
| 183 | & 0.100000000E+01_JPRB,& |
|---|
| 184 | & 0.000000000E-00_JPRB,& |
|---|
| 185 | & 0.000000000E-00_JPRB,& |
|---|
| 186 | & 0.000000000E-00_JPRB/) |
|---|
| 187 | ZBPAD2(1, 2, :) = (/& |
|---|
| 188 | & 0.5124541E+03_JPRB,& |
|---|
| 189 | & 0.3499516E+04_JPRB,& |
|---|
| 190 | & 0.1278186E+04_JPRB,& |
|---|
| 191 | & 0.100000000E+01_JPRB,& |
|---|
| 192 | & 0.000000000E-00_JPRB,& |
|---|
| 193 | & 0.000000000E-00_JPRB,& |
|---|
| 194 | & 0.000000000E-00_JPRB/) |
|---|
| 195 | ZBPAD2(1, 3, :) = (/& |
|---|
| 196 | & 0.4766857E-04_JPRB,& |
|---|
| 197 | & 0.4153343E-02_JPRB,& |
|---|
| 198 | & 0.1925810E+00_JPRB,& |
|---|
| 199 | & 0.1938938E+02_JPRB,& |
|---|
| 200 | & 0.6918191E+02_JPRB,& |
|---|
| 201 | & 0.2278795E+02_JPRB,& |
|---|
| 202 | & 0.100000000E+01_JPRB/) |
|---|
| 203 | |
|---|
| 204 | ZRAY2(1,:)= (/& |
|---|
| 205 | & .389107E-01_JPRB, .101608E+01_JPRB,-.354921E+01_JPRB,& |
|---|
| 206 | & .659510E+01_JPRB,-.592423E+01_JPRB, .202348E+01_JPRB/) |
|---|
| 207 | |
|---|
| 208 | !-- Original coefficients of ECMWF model (before 27 June 2000) |
|---|
| 209 | |
|---|
| 210 | !ZTDH2IS = 0.450_JPRB |
|---|
| 211 | !ZTDU2IS = 0.375_JPRB |
|---|
| 212 | !ZTH2IS = 296._JPRB |
|---|
| 213 | !ZTU2IS = 296._JPRB |
|---|
| 214 | !ZPDH2IS = 0.90_JPRB |
|---|
| 215 | !ZPDU2IS = 0.75_JPRB |
|---|
| 216 | !ZPRH2IS = 100000._JPRB |
|---|
| 217 | !ZPRU2IS = 100000._JPRB |
|---|
| 218 | |
|---|
| 219 | !ZD2(1,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /) |
|---|
| 220 | |
|---|
| 221 | !ZAPAD2(1, 1, :) = (/& |
|---|
| 222 | ! &0.184678379E+06_JPRB,& |
|---|
| 223 | ! &0.553080884E+05_JPRB,& |
|---|
| 224 | ! &0.248143712E+04_JPRB,& |
|---|
| 225 | ! &0.000000000E-00_JPRB,& |
|---|
| 226 | ! &0.000000000E-00_JPRB,& |
|---|
| 227 | ! &0.000000000E-00_JPRB,& |
|---|
| 228 | ! &0.000000000E-00_JPRB/) |
|---|
| 229 | !ZAPAD2(1, 2, :) = (/& |
|---|
| 230 | ! &0.715303869E+01_JPRB,& |
|---|
| 231 | ! &0.219386847E+03_JPRB,& |
|---|
| 232 | ! &0.830001089E+03_JPRB,& |
|---|
| 233 | ! &0.000000000E-00_JPRB,& |
|---|
| 234 | ! &0.000000000E-00_JPRB,& |
|---|
| 235 | ! &0.000000000E-00_JPRB,& |
|---|
| 236 | ! &0.000000000E-00_JPRB/) |
|---|
| 237 | !ZAPAD2(1, 3, :) = (/& |
|---|
| 238 | ! &0.925887084E-04_JPRB,& |
|---|
| 239 | ! &0.129353723E-01_JPRB,& |
|---|
| 240 | ! &0.800821928E+00_JPRB,& |
|---|
| 241 | ! &0.242715973E+02_JPRB,& |
|---|
| 242 | ! &0.878331486E+02_JPRB,& |
|---|
| 243 | ! &0.191559725E+02_JPRB,& |
|---|
| 244 | ! &0.000000000E+00_JPRB/) |
|---|
| 245 | |
|---|
| 246 | !ZBPAD2(1, 1, :) = (/& |
|---|
| 247 | ! &0.184678379E+06_JPRB,& |
|---|
| 248 | ! &0.555188347E+05_JPRB,& |
|---|
| 249 | ! &0.253257443E+04_JPRB,& |
|---|
| 250 | ! &0.100000000E+01_JPRB,& |
|---|
| 251 | ! &0.000000000E-00_JPRB,& |
|---|
| 252 | ! &0.000000000E-00_JPRB,& |
|---|
| 253 | ! &0.000000000E-00_JPRB/) |
|---|
| 254 | !ZBPAD2(1, 2, :) = (/& |
|---|
| 255 | ! &0.715303869E+01_JPRB,& |
|---|
| 256 | ! &0.219441875E+03_JPRB,& |
|---|
| 257 | ! &0.831119997E+03_JPRB,& |
|---|
| 258 | ! &0.100000000E+01_JPRB,& |
|---|
| 259 | ! &0.000000000E-00_JPRB,& |
|---|
| 260 | ! &0.000000000E-00_JPRB,& |
|---|
| 261 | ! &0.000000000E-00_JPRB/) |
|---|
| 262 | !ZBPAD2(1, 3, :) = (/& |
|---|
| 263 | ! &0.925887084E-04_JPRB,& |
|---|
| 264 | ! &0.131812683E-01_JPRB,& |
|---|
| 265 | ! &0.812706117E+00_JPRB,& |
|---|
| 266 | ! &0.249863591E+02_JPRB,& |
|---|
| 267 | ! &0.931071925E+02_JPRB,& |
|---|
| 268 | ! &0.252233437E+02_JPRB,& |
|---|
| 269 | ! &0.100000000E+01_JPRB/) |
|---|
| 270 | |
|---|
| 271 | ZRAY2(1,:)= (/& |
|---|
| 272 | & .428937E-01_JPRB, .890743E+00_JPRB,-.288555E+01_JPRB,& |
|---|
| 273 | & .522744E+01_JPRB,-.469173E+01_JPRB, .161645E+01_JPRB/) |
|---|
| 274 | |
|---|
| 275 | ZSWCE2(1) = 0.0_JPRB |
|---|
| 276 | ZSWCP2(1) = 0.0_JPRB |
|---|
| 277 | |
|---|
| 278 | ! ---------------------------------------------------------------- |
|---|
| 279 | |
|---|
| 280 | !* Near-Infrared (0.69 - 4.0 Microns) is not sub-divided |
|---|
| 281 | |
|---|
| 282 | ! ---------------------------------------------------------------- |
|---|
| 283 | |
|---|
| 284 | !* 0.68 - 4.00 Micron |
|---|
| 285 | |
|---|
| 286 | ZSUN2(2) = 0.540240_JPRB |
|---|
| 287 | |
|---|
| 288 | !-- New coefficients |
|---|
| 289 | |
|---|
| 290 | ZD2(2,:)= (/ 0.000000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB/) |
|---|
| 291 | |
|---|
| 292 | ZAPAD2(2, 1, :) = (/& |
|---|
| 293 | & 0.1345066E-03_JPRB,& |
|---|
| 294 | & 0.9276888E+00_JPRB,& |
|---|
| 295 | & 0.3449105E+03_JPRB,& |
|---|
| 296 | & 0.7143025E+04_JPRB,& |
|---|
| 297 | & 0.6531248E+04_JPRB,& |
|---|
| 298 | & 0.2210815E+03_JPRB,& |
|---|
| 299 | & 0.0000000E+00_JPRB/) |
|---|
| 300 | ZAPAD2(2, 2, :) = (/& |
|---|
| 301 | & 0.7776703E-05_JPRB,& |
|---|
| 302 | & 0.6851592E-02_JPRB,& |
|---|
| 303 | & 0.5533947E+00_JPRB,& |
|---|
| 304 | & 0.2538585E+01_JPRB,& |
|---|
| 305 | & 0.000000000E+00_JPRB,& |
|---|
| 306 | & 0.000000000E+00_JPRB,& |
|---|
| 307 | & 0.000000000E+00_JPRB/) |
|---|
| 308 | ZAPAD2(2, 3, :) = (/& |
|---|
| 309 | & 0.932012720E+02_JPRB,& |
|---|
| 310 | & 0.437866983E+03_JPRB,& |
|---|
| 311 | & 0.000000000E-00_JPRB,& |
|---|
| 312 | & 0.000000000E-00_JPRB,& |
|---|
| 313 | & 0.000000000E-00_JPRB,& |
|---|
| 314 | & 0.000000000E-00_JPRB,& |
|---|
| 315 | & 0.000000000E+00_JPRB/) |
|---|
| 316 | |
|---|
| 317 | ZBPAD2(2, 1, :) = (/& |
|---|
| 318 | & 0.1345066E-03_JPRB,& |
|---|
| 319 | & 0.9297728E+00_JPRB,& |
|---|
| 320 | & 0.3523821E+03_JPRB,& |
|---|
| 321 | & 0.7931590E+04_JPRB,& |
|---|
| 322 | & 0.8979086E+04_JPRB,& |
|---|
| 323 | & 0.5050864E+03_JPRB,& |
|---|
| 324 | & 0.1000000E+01_JPRB/) |
|---|
| 325 | ZBPAD2(2, 2, :) = (/& |
|---|
| 326 | & 0.7776703E-05_JPRB,& |
|---|
| 327 | & 0.6976301E-02_JPRB,& |
|---|
| 328 | & 0.6074521E+00_JPRB,& |
|---|
| 329 | & 0.3517521E+01_JPRB,& |
|---|
| 330 | & 0.100000000E+01_JPRB,& |
|---|
| 331 | & 0.000000000E+00_JPRB,& |
|---|
| 332 | & 0.000000000E+00_JPRB/) |
|---|
| 333 | ZBPAD2(2, 3, :) = (/& |
|---|
| 334 | & 0.932012720E+02_JPRB,& |
|---|
| 335 | & 0.439078599E+03_JPRB,& |
|---|
| 336 | & 0.100000000E+01_JPRB,& |
|---|
| 337 | & 0.000000000E+00_JPRB,& |
|---|
| 338 | & 0.000000000E+00_JPRB,& |
|---|
| 339 | & 0.000000000E+00_JPRB,& |
|---|
| 340 | & 0.000000000E+00_JPRB/) |
|---|
| 341 | |
|---|
| 342 | ZRAY2(2,:)= (/& |
|---|
| 343 | & .104635E-01_JPRB, .102535E+00_JPRB,-.536771E+00_JPRB,& |
|---|
| 344 | & .123640E+01_JPRB,-.128930E+01_JPRB, .502426E+00_JPRB/) |
|---|
| 345 | |
|---|
| 346 | !-- Original coefficients of ECMWF model (before 27 June 2000) |
|---|
| 347 | |
|---|
| 348 | !ZD2(2,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB/) |
|---|
| 349 | |
|---|
| 350 | !ZAPAD2(2, 1, :) = (/& |
|---|
| 351 | ! &0.577160808E-06_JPRB,& |
|---|
| 352 | ! &0.681760991E-02_JPRB,& |
|---|
| 353 | ! &0.502416370E+01_JPRB,& |
|---|
| 354 | ! &0.209771270E+03_JPRB,& |
|---|
| 355 | ! &0.644493139E+03_JPRB,& |
|---|
| 356 | ! &0.120654163E+03_JPRB,& |
|---|
| 357 | ! &0.000000000E+00_JPRB/) |
|---|
| 358 | !ZAPAD2(2, 2, :) = (/& |
|---|
| 359 | ! &0.422874202E-09_JPRB,& |
|---|
| 360 | ! &0.257499477E-05_JPRB,& |
|---|
| 361 | ! &0.349570458E-02_JPRB,& |
|---|
| 362 | ! &0.917758876E+00_JPRB,& |
|---|
| 363 | ! &0.327669989E+02_JPRB,& |
|---|
| 364 | ! &0.874155115E+02_JPRB,& |
|---|
| 365 | ! &0.000000000E+00_JPRB/) |
|---|
| 366 | !ZAPAD2(2, 3, :) = (/& |
|---|
| 367 | ! &0.932012720E+02_JPRB,& |
|---|
| 368 | ! &0.437866983E+03_JPRB,& |
|---|
| 369 | ! &0.000000000E-00_JPRB,& |
|---|
| 370 | ! &0.000000000E-00_JPRB,& |
|---|
| 371 | ! &0.000000000E-00_JPRB,& |
|---|
| 372 | ! &0.000000000E-00_JPRB,& |
|---|
| 373 | ! &0.000000000E+00_JPRB/) |
|---|
| 374 | |
|---|
| 375 | !ZBPAD2(2, 1, :) = (/& |
|---|
| 376 | ! &0.577160808E-06_JPRB,& |
|---|
| 377 | ! &0.682595589E-02_JPRB,& |
|---|
| 378 | ! &0.508797283E+01_JPRB,& |
|---|
| 379 | ! &0.223901000E+03_JPRB,& |
|---|
| 380 | ! &0.770640908E+03_JPRB,& |
|---|
| 381 | ! &0.186812480E+03_JPRB,& |
|---|
| 382 | ! &0.100000000E+01_JPRB/) |
|---|
| 383 | !ZBPAD2(2, 2, :) = (/& |
|---|
| 384 | ! &0.422874202E-09_JPRB,& |
|---|
| 385 | ! &0.257567242E-05_JPRB,& |
|---|
| 386 | ! &0.349912408E-02_JPRB,& |
|---|
| 387 | ! &0.920770481E+00_JPRB,& |
|---|
| 388 | ! &0.330990960E+02_JPRB,& |
|---|
| 389 | ! &0.903249595E+02_JPRB,& |
|---|
| 390 | ! &0.100000000E+01_JPRB/) |
|---|
| 391 | !ZBPAD2(2, 3, :) = (/& |
|---|
| 392 | ! &0.932012720E+02_JPRB,& |
|---|
| 393 | ! &0.439078599E+03_JPRB,& |
|---|
| 394 | ! &0.100000000E+01_JPRB,& |
|---|
| 395 | ! &0.000000000E+00_JPRB,& |
|---|
| 396 | ! &0.000000000E+00_JPRB,& |
|---|
| 397 | ! &0.000000000E+00_JPRB,& |
|---|
| 398 | ! &0.000000000E+00_JPRB/) |
|---|
| 399 | |
|---|
| 400 | !ZRAY2(2,:)= (/& |
|---|
| 401 | ! &.697200E-02_JPRB, .173297E-01_JPRB,-.850903E-01_JPRB,& |
|---|
| 402 | ! &.248261E+00_JPRB,-.302031E+00_JPRB, .129662E+00_JPRB/) |
|---|
| 403 | |
|---|
| 404 | ZSWCE2(2) = 0.0_JPRB |
|---|
| 405 | ZSWCP2(2) = 0.0_JPRB |
|---|
| 406 | |
|---|
| 407 | !===================================================================== |
|---|
| 408 | |
|---|
| 409 | !* 1.2 COEFFICIENTS FOR FOUR SPECTRAL INTERVALS |
|---|
| 410 | ! ---------------------------------------- |
|---|
| 411 | |
|---|
| 412 | !* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING |
|---|
| 413 | ! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9 |
|---|
| 414 | ! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2) |
|---|
| 415 | ! O3 : unchanged in interval 1, from HITRAN 92 in interval 4 |
|---|
| 416 | |
|---|
| 417 | ZADJUS4 = 1.0_JPRB |
|---|
| 418 | |
|---|
| 419 | !* 1st spectral interval: U.V. and Visible (0.25 - 0.69 Micron) |
|---|
| 420 | |
|---|
| 421 | ZSUN4(1) = 0.459760_JPRB |
|---|
| 422 | |
|---|
| 423 | !-- New coefficients |
|---|
| 424 | |
|---|
| 425 | !ZD4(1,:)= (/ 0.800000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB /) |
|---|
| 426 | |
|---|
| 427 | !ZAPAD4(1, 1, :) = (/& |
|---|
| 428 | ! &0.1095121E+03_JPRB,& |
|---|
| 429 | ! &0.8906647E+03_JPRB,& |
|---|
| 430 | ! &0.1030601E+03_JPRB,& |
|---|
| 431 | ! &0.0000000E-00_JPRB,& |
|---|
| 432 | ! &0.0000000E-00_JPRB,& |
|---|
| 433 | ! &0.0000000E-00_JPRB,& |
|---|
| 434 | ! &0.0000000E-00_JPRB/) |
|---|
| 435 | !ZAPAD4(1, 2, :) = (/& |
|---|
| 436 | ! &0.5124541E+03_JPRB,& |
|---|
| 437 | ! &0.3478800E+04_JPRB,& |
|---|
| 438 | ! &0.1244272E+04_JPRB,& |
|---|
| 439 | ! &0.0000000E-00_JPRB,& |
|---|
| 440 | ! &0.0000000E-00_JPRB,& |
|---|
| 441 | ! &0.0000000E-00_JPRB,& |
|---|
| 442 | ! &0.0000000E-00_JPRB/) |
|---|
| 443 | !ZAPAD4(1, 3, :) = (/& |
|---|
| 444 | ! &0.4766857E-04_JPRB,& |
|---|
| 445 | ! &0.4019806E-02_JPRB,& |
|---|
| 446 | ! &0.1936247E+00_JPRB,& |
|---|
| 447 | ! &0.1886522E+02_JPRB,& |
|---|
| 448 | ! &0.6498215E+02_JPRB,& |
|---|
| 449 | ! &0.1799387E+02_JPRB,& |
|---|
| 450 | ! &0.0000000E+00_JPRB/) |
|---|
| 451 | |
|---|
| 452 | !ZBPAD4(1, 1, :) = (/& |
|---|
| 453 | ! &0.1095121E+03_JPRB,& |
|---|
| 454 | ! &0.8919103E+03_JPRB,& |
|---|
| 455 | ! &0.1142218E+03_JPRB,& |
|---|
| 456 | ! &0.1000000E+01_JPRB,& |
|---|
| 457 | ! &0.0000000E-00_JPRB,& |
|---|
| 458 | ! &0.0000000E-00_JPRB,& |
|---|
| 459 | ! &0.0000000E-00_JPRB/) |
|---|
| 460 | !ZBPAD4(1, 2, :) = (/& |
|---|
| 461 | ! &0.5124541E+03_JPRB,& |
|---|
| 462 | ! &0.3499516E+04_JPRB,& |
|---|
| 463 | ! &0.1278186E+04_JPRB,& |
|---|
| 464 | ! &0.1000000E+01_JPRB,& |
|---|
| 465 | ! &0.0000000E-00_JPRB,& |
|---|
| 466 | ! &0.0000000E-00_JPRB,& |
|---|
| 467 | ! &0.0000000E-00_JPRB/) |
|---|
| 468 | !ZBPAD4(1, 3, :) = (/& |
|---|
| 469 | ! &0.4766857E-04_JPRB,& |
|---|
| 470 | ! &0.4153343E-02_JPRB,& |
|---|
| 471 | ! &0.1925810E+00_JPRB,& |
|---|
| 472 | ! &0.1938938E+02_JPRB,& |
|---|
| 473 | ! &0.6918191E+02_JPRB,& |
|---|
| 474 | ! &0.2278795E+02_JPRB,& |
|---|
| 475 | ! &0.1000000E+01_JPRB/) |
|---|
| 476 | |
|---|
| 477 | !ZRAY4(1,:)= (/& |
|---|
| 478 | ! &.389107E-01_JPRB, .101608E+01_JPRB,-.354921E+01_JPRB,& |
|---|
| 479 | ! &.659510E+01_JPRB,-.592423E+01_JPRB, .202348E+01_JPRB/) |
|---|
| 480 | |
|---|
| 481 | !-- Original coefficients in ECMWF scheme |
|---|
| 482 | |
|---|
| 483 | ZTDH4IS = 0.450_JPRB |
|---|
| 484 | ZTDU4IS = 0.375_JPRB |
|---|
| 485 | ZTH4IS = 296._JPRB |
|---|
| 486 | ZTU4IS = 296._JPRB |
|---|
| 487 | ZPDH4IS = 0.90_JPRB |
|---|
| 488 | ZPDU4IS = 0.75_JPRB |
|---|
| 489 | ZPRH4IS = 100000._JPRB |
|---|
| 490 | ZPRU4IS = 100000._JPRB |
|---|
| 491 | |
|---|
| 492 | ZD4(1,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /) |
|---|
| 493 | |
|---|
| 494 | ZAPAD4(1, 1, :) = (/& |
|---|
| 495 | & 0.184678379E+06_JPRB,& |
|---|
| 496 | & 0.553080884E+05_JPRB,& |
|---|
| 497 | & 0.248143712E+04_JPRB,& |
|---|
| 498 | & 0.000000000E-00_JPRB,& |
|---|
| 499 | & 0.000000000E-00_JPRB,& |
|---|
| 500 | & 0.000000000E-00_JPRB,& |
|---|
| 501 | & 0.000000000E-00_JPRB/) |
|---|
| 502 | ZAPAD4(1, 2, :) = (/& |
|---|
| 503 | & 0.715303869E+01_JPRB,& |
|---|
| 504 | & 0.219386847E+03_JPRB,& |
|---|
| 505 | & 0.830001089E+03_JPRB,& |
|---|
| 506 | & 0.000000000E-00_JPRB,& |
|---|
| 507 | & 0.000000000E-00_JPRB,& |
|---|
| 508 | & 0.000000000E-00_JPRB,& |
|---|
| 509 | & 0.000000000E-00_JPRB/) |
|---|
| 510 | ZAPAD4(1, 3, :) = (/& |
|---|
| 511 | & 0.925887084E-04_JPRB,& |
|---|
| 512 | & 0.129353723E-01_JPRB,& |
|---|
| 513 | & 0.800821928E+00_JPRB,& |
|---|
| 514 | & 0.242715973E+02_JPRB,& |
|---|
| 515 | & 0.878331486E+02_JPRB,& |
|---|
| 516 | & 0.191559725E+02_JPRB,& |
|---|
| 517 | & 0.000000000E+00_JPRB/) |
|---|
| 518 | |
|---|
| 519 | ZBPAD4(1, 1, :) = (/& |
|---|
| 520 | & 0.184678379E+06_JPRB,& |
|---|
| 521 | & 0.555188347E+05_JPRB,& |
|---|
| 522 | & 0.253257443E+04_JPRB,& |
|---|
| 523 | & 0.100000000E+01_JPRB,& |
|---|
| 524 | & 0.000000000E-00_JPRB,& |
|---|
| 525 | & 0.000000000E-00_JPRB,& |
|---|
| 526 | & 0.000000000E-00_JPRB/) |
|---|
| 527 | ZBPAD4(1, 2, :) = (/& |
|---|
| 528 | & 0.715303869E+01_JPRB,& |
|---|
| 529 | & 0.219441875E+03_JPRB,& |
|---|
| 530 | & 0.831119997E+03_JPRB,& |
|---|
| 531 | & 0.100000000E+01_JPRB,& |
|---|
| 532 | & 0.000000000E-00_JPRB,& |
|---|
| 533 | & 0.000000000E-00_JPRB,& |
|---|
| 534 | & 0.000000000E-00_JPRB/) |
|---|
| 535 | ZBPAD4(1, 3, :) = (/& |
|---|
| 536 | & 0.925887084E-04_JPRB,& |
|---|
| 537 | & 0.131812683E-01_JPRB,& |
|---|
| 538 | & 0.812706117E+00_JPRB,& |
|---|
| 539 | & 0.249863591E+02_JPRB,& |
|---|
| 540 | & 0.931071925E+02_JPRB,& |
|---|
| 541 | & 0.252233437E+02_JPRB,& |
|---|
| 542 | & 0.100000000E+01_JPRB/) |
|---|
| 543 | |
|---|
| 544 | ZRAY4(1,:)= (/& |
|---|
| 545 | & .428937E-01_JPRB, .890743E+00_JPRB,-.288555E+01_JPRB,& |
|---|
| 546 | & .522744E+01_JPRB,-.469173E+01_JPRB, .161645E+01_JPRB/) |
|---|
| 547 | |
|---|
| 548 | ZSWCE4(1) = 0.0_JPRB |
|---|
| 549 | ZSWCP4(1) = 0.0_JPRB |
|---|
| 550 | |
|---|
| 551 | ! ---------------------------------------------------------------- |
|---|
| 552 | |
|---|
| 553 | !* Near-Infrared (0.69 - 4.0 Microns) is sub-divided into: |
|---|
| 554 | |
|---|
| 555 | ! ---------------------------------------------------------------- |
|---|
| 556 | |
|---|
| 557 | !* 0.69 - 1.19 Micron |
|---|
| 558 | |
|---|
| 559 | ZSUN4(2) = 0.326158_JPRB |
|---|
| 560 | |
|---|
| 561 | !-- New coefficients |
|---|
| 562 | |
|---|
| 563 | !ZD4(2,:)= (/ 0.000000000_JPRB, 0.900000000_JPRB, 1.000000000_JPRB /) |
|---|
| 564 | |
|---|
| 565 | !ZAPAD4(2, 1, :) = (/& |
|---|
| 566 | ! &0.1335726E+02_JPRB,& |
|---|
| 567 | ! &0.2939136E+04_JPRB,& |
|---|
| 568 | ! &0.4010585E+05_JPRB,& |
|---|
| 569 | ! &0.7195030E+05_JPRB,& |
|---|
| 570 | ! &0.1648338E+05_JPRB,& |
|---|
| 571 | ! &0.3373738E+03_JPRB,& |
|---|
| 572 | ! &0.0000000E+00_JPRB/) |
|---|
| 573 | !ZAPAD4(2, 2, :) = (/& |
|---|
| 574 | ! &0.2001271E-01_JPRB,& |
|---|
| 575 | ! &0.2480831E+01_JPRB,& |
|---|
| 576 | ! &0.3444162E+02_JPRB,& |
|---|
| 577 | ! &0.4788946E+02_JPRB,& |
|---|
| 578 | ! &0.0000000E+00_JPRB,& |
|---|
| 579 | ! &0.0000000E+00_JPRB,& |
|---|
| 580 | ! &0.0000000E+00_JPRB/) |
|---|
| 581 | !ZAPAD4(2, 3, :) = (/& |
|---|
| 582 | ! &0.000000000E+00_JPRB,& |
|---|
| 583 | ! &0.000000000E+00_JPRB,& |
|---|
| 584 | ! &0.000000000E+00_JPRB,& |
|---|
| 585 | ! &0.000000000E+00_JPRB,& |
|---|
| 586 | ! &0.000000000E+00_JPRB,& |
|---|
| 587 | ! &0.000000000E+00_JPRB,& |
|---|
| 588 | ! &0.000000000E+00_JPRB/) |
|---|
| 589 | |
|---|
| 590 | !ZBPAD4(2, 1, :) = (/& |
|---|
| 591 | ! &0.1335726E+02_JPRB,& |
|---|
| 592 | ! &0.2942327E+04_JPRB,& |
|---|
| 593 | ! &0.4077237E+05_JPRB,& |
|---|
| 594 | ! &0.7749017E+05_JPRB,& |
|---|
| 595 | ! &0.2123132E+05_JPRB,& |
|---|
| 596 | ! &0.6659687E+03_JPRB,& |
|---|
| 597 | ! &0.1000000E+01_JPRB/) |
|---|
| 598 | !ZBPAD4(2, 2, :) = (/& |
|---|
| 599 | ! &0.2001271E-01_JPRB,& |
|---|
| 600 | ! &0.2549067E+01_JPRB,& |
|---|
| 601 | ! &0.3752433E+02_JPRB,& |
|---|
| 602 | ! &0.6276637E+02_JPRB,& |
|---|
| 603 | ! &0.1000000E+01_JPRB,& |
|---|
| 604 | ! &0.0000000E+00_JPRB,& |
|---|
| 605 | ! &0.0000000E+00_JPRB/) |
|---|
| 606 | !ZBPAD4(2, 3, :) = (/& |
|---|
| 607 | ! &1.000000000E+00_JPRB,& |
|---|
| 608 | ! &0.000000000E+00_JPRB,& |
|---|
| 609 | ! &0.000000000E+00_JPRB,& |
|---|
| 610 | ! &0.000000000E+00_JPRB,& |
|---|
| 611 | ! &0.000000000E+00_JPRB,& |
|---|
| 612 | ! &0.000000000E+00_JPRB,& |
|---|
| 613 | ! &0.000000000E+00_JPRB/) |
|---|
| 614 | |
|---|
| 615 | !DATA (ZRAY4(2,K),K=1,6) /& |
|---|
| 616 | ! & .519710E-01_JPRB,-.150667E+00_JPRB, .154170E+01_JPRB,& |
|---|
| 617 | ! &-.518901E+01_JPRB, .668011E+01_JPRB,-.294843E+01_JPRB/ |
|---|
| 618 | |
|---|
| 619 | !-- Original coefficients in the ECMWF model |
|---|
| 620 | |
|---|
| 621 | ZD4(2,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 1.000000000_JPRB /) |
|---|
| 622 | |
|---|
| 623 | ZAPAD4(2, 1, :) = (/& |
|---|
| 624 | & 0.690730834E-02_JPRB,& |
|---|
| 625 | & 0.151704275E+01_JPRB,& |
|---|
| 626 | & 0.751477543E+02_JPRB,& |
|---|
| 627 | & 0.759770236E+03_JPRB,& |
|---|
| 628 | & 0.109800326E+04_JPRB,& |
|---|
| 629 | & 0.148407574E+03_JPRB,& |
|---|
| 630 | & 0.000000000E+00_JPRB/) |
|---|
| 631 | ZAPAD4(2, 2, :) = (/& |
|---|
| 632 | & 0.863790752E-03_JPRB,& |
|---|
| 633 | & 0.448762291E+00_JPRB,& |
|---|
| 634 | & 0.332530367E+02_JPRB,& |
|---|
| 635 | & 0.190914146E+03_JPRB,& |
|---|
| 636 | & 0.000000000E+00_JPRB,& |
|---|
| 637 | & 0.000000000E+00_JPRB,& |
|---|
| 638 | & 0.000000000E+00_JPRB/) |
|---|
| 639 | ZAPAD4(2, 3, :) = (/& |
|---|
| 640 | & 0.000000000E+00_JPRB,& |
|---|
| 641 | & 0.000000000E+00_JPRB,& |
|---|
| 642 | & 0.000000000E+00_JPRB,& |
|---|
| 643 | & 0.000000000E+00_JPRB,& |
|---|
| 644 | & 0.000000000E+00_JPRB,& |
|---|
| 645 | & 0.000000000E+00_JPRB,& |
|---|
| 646 | & 0.000000000E+00_JPRB/) |
|---|
| 647 | |
|---|
| 648 | ZBPAD4(2, 1, :) = (/& |
|---|
| 649 | & 0.690730834E-02_JPRB,& |
|---|
| 650 | & 0.151954406E+01_JPRB,& |
|---|
| 651 | & 0.756512527E+02_JPRB,& |
|---|
| 652 | & 0.779384997E+03_JPRB,& |
|---|
| 653 | & 0.121113108E+04_JPRB,& |
|---|
| 654 | & 0.207678436E+03_JPRB,& |
|---|
| 655 | & 0.100000000E+01_JPRB/) |
|---|
| 656 | ZBPAD4(2, 2, :) = (/& |
|---|
| 657 | & 0.863790752E-03_JPRB,& |
|---|
| 658 | & 0.448948107E+00_JPRB,& |
|---|
| 659 | & 0.333186750E+02_JPRB,& |
|---|
| 660 | & 0.192727216E+03_JPRB,& |
|---|
| 661 | & 0.100000000E+01_JPRB,& |
|---|
| 662 | & 0.000000000E+00_JPRB,& |
|---|
| 663 | & 0.000000000E+00_JPRB/) |
|---|
| 664 | ZBPAD4(2, 3, :) = (/& |
|---|
| 665 | & 1.000000000E+00_JPRB,& |
|---|
| 666 | & 0.000000000E+00_JPRB,& |
|---|
| 667 | & 0.000000000E+00_JPRB,& |
|---|
| 668 | & 0.000000000E+00_JPRB,& |
|---|
| 669 | & 0.000000000E+00_JPRB,& |
|---|
| 670 | & 0.000000000E+00_JPRB,& |
|---|
| 671 | & 0.000000000E+00_JPRB/) |
|---|
| 672 | |
|---|
| 673 | DATA (ZRAY4(2,I_K),I_K=1,6) /& |
|---|
| 674 | & .164261E-01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,& |
|---|
| 675 | & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/ |
|---|
| 676 | |
|---|
| 677 | ZSWCE4(2) = 0.0_JPRB |
|---|
| 678 | ZSWCP4(2) = 0.0_JPRB |
|---|
| 679 | |
|---|
| 680 | ! ---------------------------------------------------------------- |
|---|
| 681 | |
|---|
| 682 | !* 1.19 - 2.38 Microns |
|---|
| 683 | |
|---|
| 684 | ZSUN4(3) = 0.180608_JPRB |
|---|
| 685 | |
|---|
| 686 | !-- New coefficients |
|---|
| 687 | |
|---|
| 688 | !ZD4(3,:)= (/ 0.000000000_JPRB, 0.800000000_JPRB, 1.000000000_JPRB /) |
|---|
| 689 | |
|---|
| 690 | !ZAPAD4(3, 1, :) = (/& |
|---|
| 691 | ! &0.3325841E-03_JPRB,& |
|---|
| 692 | ! &0.6194496E+00_JPRB,& |
|---|
| 693 | ! &0.1497138E+03_JPRB,& |
|---|
| 694 | ! &0.2314864E+04_JPRB,& |
|---|
| 695 | ! &0.2380109E+04_JPRB,& |
|---|
| 696 | ! &0.9553823E+02_JPRB,& |
|---|
| 697 | ! &0.0000000E+00_JPRB/) |
|---|
| 698 | !ZAPAD4(3, 2, :) = (/& |
|---|
| 699 | ! &0.4552471E-03_JPRB,& |
|---|
| 700 | ! &0.4084154E+00_JPRB,& |
|---|
| 701 | ! &0.6114905E+01_JPRB,& |
|---|
| 702 | ! &0.7102540E+01_JPRB,& |
|---|
| 703 | ! &0.0000000E+00_JPRB,& |
|---|
| 704 | ! &0.0000000E+00_JPRB,& |
|---|
| 705 | ! &0.0000000E+00_JPRB/) |
|---|
| 706 | !ZAPAD4(3, 3, :) = (/& |
|---|
| 707 | ! &0.000000000E+00_JPRB,& |
|---|
| 708 | ! &0.000000000E+00_JPRB,& |
|---|
| 709 | ! &0.000000000E+00_JPRB,& |
|---|
| 710 | ! &0.000000000E+00_JPRB,& |
|---|
| 711 | ! &0.000000000E+00_JPRB,& |
|---|
| 712 | ! &0.000000000E+00_JPRB,& |
|---|
| 713 | ! &0.000000000E+00_JPRB/) |
|---|
| 714 | |
|---|
| 715 | !ZBPAD4(3, 1, :) = (/& |
|---|
| 716 | ! &0.3325841E-03_JPRB,& |
|---|
| 717 | ! &0.6231947E+00_JPRB,& |
|---|
| 718 | ! &0.1553098E+03_JPRB,& |
|---|
| 719 | ! &0.2822458E+04_JPRB,& |
|---|
| 720 | ! &0.3885194E+04_JPRB,& |
|---|
| 721 | ! &0.2700235E+03_JPRB,& |
|---|
| 722 | ! &0.1000000E+01_JPRB/) |
|---|
| 723 | !ZBPAD4(3, 2, :) = (/& |
|---|
| 724 | ! &0.4552471E-03_JPRB,& |
|---|
| 725 | ! &0.4088242E+00_JPRB,& |
|---|
| 726 | ! &0.6411905E+01_JPRB,& |
|---|
| 727 | ! &0.9444439E+01_JPRB,& |
|---|
| 728 | ! &0.1000000E+01_JPRB,& |
|---|
| 729 | ! &0.0000000E+00_JPRB,& |
|---|
| 730 | ! &0.0000000E+00_JPRB/) |
|---|
| 731 | !ZBPAD4(3, 3, :) = (/& |
|---|
| 732 | ! &1.000000000E+00_JPRB,& |
|---|
| 733 | ! &0.000000000E+00_JPRB,& |
|---|
| 734 | ! &0.000000000E+00_JPRB,& |
|---|
| 735 | ! &0.000000000E+00_JPRB,& |
|---|
| 736 | ! &0.000000000E+00_JPRB,& |
|---|
| 737 | ! &0.000000000E+00_JPRB,& |
|---|
| 738 | ! &0.000000000E+00_JPRB/) |
|---|
| 739 | |
|---|
| 740 | !ZRAY4(3,:)= (/& |
|---|
| 741 | ! &.288474E-02_JPRB, .133984E-02_JPRB,-.152812E-01_JPRB,& |
|---|
| 742 | ! &.540845E-01_JPRB,-.747444E-01_JPRB, .364285E-01_JPRB/) |
|---|
| 743 | |
|---|
| 744 | ZD4(3,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 1.000000000_JPRB /) |
|---|
| 745 | |
|---|
| 746 | ZAPAD4(3, 1, :) = (/& |
|---|
| 747 | & 0.837531303E-05_JPRB,& |
|---|
| 748 | & 0.173886341E-01_JPRB,& |
|---|
| 749 | & 0.518852799E+01_JPRB,& |
|---|
| 750 | & 0.159078416E+03_JPRB,& |
|---|
| 751 | & 0.493273523E+03_JPRB,& |
|---|
| 752 | & 0.102567293E+03_JPRB,& |
|---|
| 753 | & 0.000000000E+00_JPRB/) |
|---|
| 754 | ZAPAD4(3, 2, :) = (/& |
|---|
| 755 | & 0.657978575E-02_JPRB,& |
|---|
| 756 | & 0.752617872E+00_JPRB,& |
|---|
| 757 | & 0.158209734E+02_JPRB,& |
|---|
| 758 | & 0.410274915E+02_JPRB,& |
|---|
| 759 | & 0.000000000E+00_JPRB,& |
|---|
| 760 | & 0.000000000E+00_JPRB,& |
|---|
| 761 | & 0.000000000E+00_JPRB/) |
|---|
| 762 | ZAPAD4(3, 3, :) = (/& |
|---|
| 763 | & 0.000000000E+00_JPRB,& |
|---|
| 764 | & 0.000000000E+00_JPRB,& |
|---|
| 765 | & 0.000000000E+00_JPRB,& |
|---|
| 766 | & 0.000000000E+00_JPRB,& |
|---|
| 767 | & 0.000000000E+00_JPRB,& |
|---|
| 768 | & 0.000000000E+00_JPRB,& |
|---|
| 769 | & 0.000000000E+00_JPRB/) |
|---|
| 770 | |
|---|
| 771 | ZBPAD4(3, 1, :) = (/& |
|---|
| 772 | & 0.837531303E-05_JPRB,& |
|---|
| 773 | & 0.174882536E-01_JPRB,& |
|---|
| 774 | & 0.534536580E+01_JPRB,& |
|---|
| 775 | & 0.180351767E+03_JPRB,& |
|---|
| 776 | & 0.673126838E+03_JPRB,& |
|---|
| 777 | & 0.182718543E+03_JPRB,& |
|---|
| 778 | & 0.100000000E+01_JPRB/) |
|---|
| 779 | ZBPAD4(3, 2, :) = (/& |
|---|
| 780 | & 0.657978575E-02_JPRB,& |
|---|
| 781 | & 0.753752065E+00_JPRB,& |
|---|
| 782 | & 0.159286262E+02_JPRB,& |
|---|
| 783 | & 0.424278450E+02_JPRB,& |
|---|
| 784 | & 0.100000000E+01_JPRB,& |
|---|
| 785 | & 0.000000000E+00_JPRB,& |
|---|
| 786 | & 0.000000000E+00_JPRB/) |
|---|
| 787 | ZBPAD4(3, 3, :) = (/& |
|---|
| 788 | & 1.000000000E+00_JPRB,& |
|---|
| 789 | & 0.000000000E+00_JPRB,& |
|---|
| 790 | & 0.000000000E+00_JPRB,& |
|---|
| 791 | & 0.000000000E+00_JPRB,& |
|---|
| 792 | & 0.000000000E+00_JPRB,& |
|---|
| 793 | & 0.000000000E+00_JPRB,& |
|---|
| 794 | & 0.000000000E+00_JPRB/) |
|---|
| 795 | |
|---|
| 796 | ZRAY4(3,:)= (/& |
|---|
| 797 | & .180438E-02_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,& |
|---|
| 798 | & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/) |
|---|
| 799 | |
|---|
| 800 | ZSWCE4(3) = 0.0_JPRB |
|---|
| 801 | ZSWCP4(3) = 0.0_JPRB |
|---|
| 802 | |
|---|
| 803 | ! ---------------------------------------------------------------- |
|---|
| 804 | |
|---|
| 805 | !* 2.38 - 4.00 Microns |
|---|
| 806 | |
|---|
| 807 | ZSUN4(4) = 0.033474_JPRB |
|---|
| 808 | |
|---|
| 809 | !-- New coefficients |
|---|
| 810 | |
|---|
| 811 | !ZD4(4,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /) |
|---|
| 812 | |
|---|
| 813 | !ZAPAD4(4, 1, :) = (/& |
|---|
| 814 | ! &0.2122889E-06_JPRB,& |
|---|
| 815 | ! &0.9030576E-03_JPRB,& |
|---|
| 816 | ! &0.2431282E+00_JPRB,& |
|---|
| 817 | ! &0.4901345E+01_JPRB,& |
|---|
| 818 | ! &0.3996347E+01_JPRB,& |
|---|
| 819 | ! &0.3910227E+01_JPRB,& |
|---|
| 820 | ! &0.0000000E+00_JPRB/) |
|---|
| 821 | !ZAPAD4(4, 2, :) = (/& |
|---|
| 822 | ! &0.1215163E-03_JPRB,& |
|---|
| 823 | ! &0.1222574E+00_JPRB,& |
|---|
| 824 | ! &0.9382420E+01_JPRB,& |
|---|
| 825 | ! &0.6875727E+02_JPRB,& |
|---|
| 826 | ! &0.2746421E+02_JPRB,& |
|---|
| 827 | ! &0.0000000E+00_JPRB,& |
|---|
| 828 | ! &0.0000000E+00_JPRB/) |
|---|
| 829 | !ZAPAD4(4, 3, :) = (/& |
|---|
| 830 | ! &0.263068898E+02_JPRB,& |
|---|
| 831 | ! &0.146425875E+03_JPRB,& |
|---|
| 832 | ! &0.860137809E+02_JPRB,& |
|---|
| 833 | ! &0.000000000E+00_JPRB,& |
|---|
| 834 | ! &0.000000000E+00_JPRB,& |
|---|
| 835 | ! &0.000000000E+00_JPRB,& |
|---|
| 836 | ! &0.000000000E+00_JPRB/) |
|---|
| 837 | |
|---|
| 838 | !ZBPAD4(4, 1, :) = (/& |
|---|
| 839 | ! &0.2122889E-06_JPRB,& |
|---|
| 840 | ! &0.9379083E-03_JPRB,& |
|---|
| 841 | ! &0.2957335E+00_JPRB,& |
|---|
| 842 | ! &0.8747190E+01_JPRB,& |
|---|
| 843 | ! &0.1015794E+02_JPRB,& |
|---|
| 844 | ! &0.1361277E+02_JPRB,& |
|---|
| 845 | ! &0.1000000E+01_JPRB/) |
|---|
| 846 | !ZBPAD4(4, 2, :) = (/& |
|---|
| 847 | ! &0.1215163E-03_JPRB,& |
|---|
| 848 | ! &0.1255648E+00_JPRB,& |
|---|
| 849 | ! &0.1060119E+02_JPRB,& |
|---|
| 850 | ! &0.8414439E+02_JPRB,& |
|---|
| 851 | ! &0.4299438E+02_JPRB,& |
|---|
| 852 | ! &0.1000000E+01_JPRB,& |
|---|
| 853 | ! &0.0000000E+00_JPRB/) |
|---|
| 854 | !ZBPAD4(4, 3, :) = (/& |
|---|
| 855 | ! &0.263068898E+02_JPRB,& |
|---|
| 856 | ! &0.152569217E+03_JPRB,& |
|---|
| 857 | ! &0.976791971E+02_JPRB,& |
|---|
| 858 | ! &0.100000000E+01_JPRB,& |
|---|
| 859 | ! &0.000000000E+00_JPRB,& |
|---|
| 860 | ! &0.000000000E+00_JPRB,& |
|---|
| 861 | ! &0.000000000E+00_JPRB/) |
|---|
| 862 | |
|---|
| 863 | !ZRAY4(4,:)= (/& |
|---|
| 864 | ! &.493545E-03_JPRB, .355696E-03_JPRB,-.284543E-02_JPRB,& |
|---|
| 865 | ! &.741248E-02_JPRB,-.773691E-02_JPRB, .370679E-02_JPRB/) |
|---|
| 866 | |
|---|
| 867 | ZD4(4,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /) |
|---|
| 868 | |
|---|
| 869 | ZAPAD4(4, 1, :) = (/& |
|---|
| 870 | & 0.122118185E-06_JPRB,& |
|---|
| 871 | & 0.154042531E-02_JPRB,& |
|---|
| 872 | & 0.141152193E+01_JPRB,& |
|---|
| 873 | & 0.685368761E+02_JPRB,& |
|---|
| 874 | & 0.216522281E+03_JPRB,& |
|---|
| 875 | & 0.421228746E+02_JPRB,& |
|---|
| 876 | & 0.000000000E+00_JPRB/) |
|---|
| 877 | ZAPAD4(4, 2, :) = (/& |
|---|
| 878 | & 0.364233560E-10_JPRB,& |
|---|
| 879 | & 0.217340835E-06_JPRB,& |
|---|
| 880 | & 0.292623386E-03_JPRB,& |
|---|
| 881 | & 0.797100631E-01_JPRB,& |
|---|
| 882 | & 0.319103672E+01_JPRB,& |
|---|
| 883 | & 0.110530283E+02_JPRB,& |
|---|
| 884 | & 0.000000000E+00_JPRB/) |
|---|
| 885 | ZAPAD4(4, 3, :) = (/& |
|---|
| 886 | & 0.263068898E+02_JPRB,& |
|---|
| 887 | & 0.146425875E+03_JPRB,& |
|---|
| 888 | & 0.860137809E+02_JPRB,& |
|---|
| 889 | & 0.000000000E+00_JPRB,& |
|---|
| 890 | & 0.000000000E+00_JPRB,& |
|---|
| 891 | & 0.000000000E+00_JPRB,& |
|---|
| 892 | & 0.000000000E+00_JPRB/) |
|---|
| 893 | |
|---|
| 894 | ZBPAD4(4, 1, :) = (/& |
|---|
| 895 | & 0.122118185E-06_JPRB,& |
|---|
| 896 | & 0.156287582E-02_JPRB,& |
|---|
| 897 | & 0.156994562E+01_JPRB,& |
|---|
| 898 | & 0.102304103E+03_JPRB,& |
|---|
| 899 | & 0.475479878E+03_JPRB,& |
|---|
| 900 | & 0.188714799E+03_JPRB,& |
|---|
| 901 | & 0.100000000E+01_JPRB/) |
|---|
| 902 | ZBPAD4(4, 2, :) = (/& |
|---|
| 903 | & 0.364233560E-10_JPRB,& |
|---|
| 904 | & 0.218265806E-06_JPRB,& |
|---|
| 905 | & 0.297085917E-03_JPRB,& |
|---|
| 906 | & 0.834253403E-01_JPRB,& |
|---|
| 907 | & 0.358290002E+01_JPRB,& |
|---|
| 908 | & 0.139206164E+02_JPRB,& |
|---|
| 909 | & 0.100000000E+01_JPRB/) |
|---|
| 910 | ZBPAD4(4, 3, :) = (/& |
|---|
| 911 | & 0.263068898E+02_JPRB,& |
|---|
| 912 | & 0.152569217E+03_JPRB,& |
|---|
| 913 | & 0.976791971E+02_JPRB,& |
|---|
| 914 | & 0.100000000E+01_JPRB,& |
|---|
| 915 | & 0.000000000E+00_JPRB,& |
|---|
| 916 | & 0.000000000E+00_JPRB,& |
|---|
| 917 | & 0.000000000E+00_JPRB/) |
|---|
| 918 | |
|---|
| 919 | ZRAY4(4,:)= (/& |
|---|
| 920 | & .136020E-03_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,& |
|---|
| 921 | & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/) |
|---|
| 922 | |
|---|
| 923 | ZSWCE4(4) = 0.0_JPRB |
|---|
| 924 | ZSWCP4(4) = 0.0_JPRB |
|---|
| 925 | |
|---|
| 926 | !===================================================================== |
|---|
| 927 | |
|---|
| 928 | !* 1.3 COEFFICIENTS FOR SIX SPECTRAL INTERVALS |
|---|
| 929 | ! --------------------------------------- |
|---|
| 930 | |
|---|
| 931 | !* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING |
|---|
| 932 | ! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9 |
|---|
| 933 | ! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2) |
|---|
| 934 | ! O3 : unchanged in interval 1, from HITRAN 92 in interval 4 |
|---|
| 935 | |
|---|
| 936 | ZADJUS6 = 0.9905709_JPRB |
|---|
| 937 | ZADJUS6 = 1.0_JPRB |
|---|
| 938 | |
|---|
| 939 | ZTDH6IS = 0.450_JPRB |
|---|
| 940 | ZTDU6IS = 0.375_JPRB |
|---|
| 941 | ZTH6IS = 273._JPRB |
|---|
| 942 | ZTU6IS = 273._JPRB |
|---|
| 943 | ZPDH6IS = 0.90_JPRB |
|---|
| 944 | ZPDU6IS = 0.75_JPRB |
|---|
| 945 | ZPRH6IS = 101300._JPRB |
|---|
| 946 | ZPRU6IS = 101300._JPRB |
|---|
| 947 | |
|---|
| 948 | !* 1st spectral interval: U.V. (0.18 - 0.25 Micron) |
|---|
| 949 | |
|---|
| 950 | ZSUN6(1) = 0.001917_JPRB |
|---|
| 951 | |
|---|
| 952 | ZD6(1,:)= (/ 1.000000000_JPRB, 1.000000000_JPRB, 0.000000000_JPRB /) |
|---|
| 953 | |
|---|
| 954 | ZAPAD6(1, 1, :) = (/& |
|---|
| 955 | & 0.000000000E-00_JPRB,& |
|---|
| 956 | & 0.000000000E-00_JPRB,& |
|---|
| 957 | & 0.000000000E-00_JPRB,& |
|---|
| 958 | & 0.000000000E-00_JPRB,& |
|---|
| 959 | & 0.000000000E-00_JPRB,& |
|---|
| 960 | & 0.000000000E-00_JPRB,& |
|---|
| 961 | & 0.000000000E-00_JPRB/) |
|---|
| 962 | ZAPAD6(1, 2, :) = (/& |
|---|
| 963 | & 0.000000000E-00_JPRB,& |
|---|
| 964 | & 0.000000000E-00_JPRB,& |
|---|
| 965 | & 0.000000000E-00_JPRB,& |
|---|
| 966 | & 0.000000000E-00_JPRB,& |
|---|
| 967 | & 0.000000000E-00_JPRB,& |
|---|
| 968 | & 0.000000000E-00_JPRB,& |
|---|
| 969 | & 0.000000000E-00_JPRB/) |
|---|
| 970 | ZAPAD6(1, 3, :) = (/& |
|---|
| 971 | & 0.000000000E-00_JPRB,& |
|---|
| 972 | & 0.000000000E-00_JPRB,& |
|---|
| 973 | & 0.000000000E-00_JPRB,& |
|---|
| 974 | & 0.000000000E-00_JPRB,& |
|---|
| 975 | & 0.000000000E-00_JPRB,& |
|---|
| 976 | & 0.000000000E-00_JPRB,& |
|---|
| 977 | & 0.000000000E-00_JPRB/) |
|---|
| 978 | |
|---|
| 979 | ZBPAD6(1, 1, :) = (/& |
|---|
| 980 | & 0.100000000E+01_JPRB,& |
|---|
| 981 | & 0.000000000E-00_JPRB,& |
|---|
| 982 | & 0.000000000E-00_JPRB,& |
|---|
| 983 | & 0.000000000E-00_JPRB,& |
|---|
| 984 | & 0.000000000E-00_JPRB,& |
|---|
| 985 | & 0.000000000E-00_JPRB,& |
|---|
| 986 | & 0.000000000E-00_JPRB/) |
|---|
| 987 | ZBPAD6(1, 2, :) = (/& |
|---|
| 988 | & 0.100000000E+01_JPRB,& |
|---|
| 989 | & 0.000000000E-00_JPRB,& |
|---|
| 990 | & 0.000000000E-00_JPRB,& |
|---|
| 991 | & 0.000000000E-00_JPRB,& |
|---|
| 992 | & 0.000000000E-00_JPRB,& |
|---|
| 993 | & 0.000000000E-00_JPRB,& |
|---|
| 994 | & 0.000000000E-00_JPRB/) |
|---|
| 995 | ZBPAD6(1, 3, :) = (/& |
|---|
| 996 | & 0.100000000E+01_JPRB,& |
|---|
| 997 | & 0.000000000E-00_JPRB,& |
|---|
| 998 | & 0.000000000E-00_JPRB,& |
|---|
| 999 | & 0.000000000E-00_JPRB,& |
|---|
| 1000 | & 0.000000000E-00_JPRB,& |
|---|
| 1001 | & 0.000000000E-00_JPRB,& |
|---|
| 1002 | & 0.000000000E-00_JPRB/) |
|---|
| 1003 | |
|---|
| 1004 | ZRAY6(1,:)= (/& |
|---|
| 1005 | & .334018E+00_JPRB,-.298776E+01_JPRB, .613911E+02_JPRB,& |
|---|
| 1006 | & -.164331E+03_JPRB, .197946E+03_JPRB,-.652206E+02_JPRB/) |
|---|
| 1007 | ZRBY6(1,:)= (/& |
|---|
| 1008 | & .43959002E+01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,& |
|---|
| 1009 | & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/) |
|---|
| 1010 | |
|---|
| 1011 | ZSWCE6(1) = 0.0_JPRB |
|---|
| 1012 | ZSWCP6(1) = 0.0_JPRB |
|---|
| 1013 | |
|---|
| 1014 | NEXPO3(1) = 7 |
|---|
| 1015 | REXPO3(1, 1, :) =(/& |
|---|
| 1016 | & 0.051395E+00_JPRB, 0.048250E+00_JPRB, 0.112339E+00_JPRB,& |
|---|
| 1017 | & 0.101426E+00_JPRB, 0.007700E+00_JPRB, 0.441320E+00_JPRB,& |
|---|
| 1018 | & 0.237571E+00_JPRB /) |
|---|
| 1019 | REXPO3(1, 2, :) =(/& |
|---|
| 1020 | & 0.100022E+02_JPRB, 0.851159E+02_JPRB, 0.346737E+03_JPRB,& |
|---|
| 1021 | & 0.158501E+02_JPRB, 0.724223E+01_JPRB, 0.177828E+03_JPRB,& |
|---|
| 1022 | & 0.467708E+02_JPRB /) |
|---|
| 1023 | |
|---|
| 1024 | !* 2nd spectral interval: U.V. (0.25 - 0.44 Micron) |
|---|
| 1025 | |
|---|
| 1026 | ZSUN6(2) = 0.135708_JPRB |
|---|
| 1027 | |
|---|
| 1028 | ZD6(2,:)= (/ 1.000000000_JPRB, 1.000000000_JPRB, 0.000000000_JPRB /) |
|---|
| 1029 | |
|---|
| 1030 | ZAPAD6(2, 1, :) = (/& |
|---|
| 1031 | & 0.000000000E-00_JPRB,& |
|---|
| 1032 | & 0.000000000E-00_JPRB,& |
|---|
| 1033 | & 0.000000000E-00_JPRB,& |
|---|
| 1034 | & 0.000000000E-00_JPRB,& |
|---|
| 1035 | & 0.000000000E-00_JPRB,& |
|---|
| 1036 | & 0.000000000E-00_JPRB,& |
|---|
| 1037 | & 0.000000000E-00_JPRB/) |
|---|
| 1038 | ZAPAD6(2, 2, :) = (/& |
|---|
| 1039 | & 0.000000000E-00_JPRB,& |
|---|
| 1040 | & 0.000000000E-00_JPRB,& |
|---|
| 1041 | & 0.000000000E-00_JPRB,& |
|---|
| 1042 | & 0.000000000E-00_JPRB,& |
|---|
| 1043 | & 0.000000000E-00_JPRB,& |
|---|
| 1044 | & 0.000000000E-00_JPRB,& |
|---|
| 1045 | & 0.000000000E-00_JPRB/) |
|---|
| 1046 | ZAPAD6(2, 3, :) = (/& |
|---|
| 1047 | & 0.000000000E-00_JPRB,& |
|---|
| 1048 | & 0.000000000E-00_JPRB,& |
|---|
| 1049 | & 0.000000000E-00_JPRB,& |
|---|
| 1050 | & 0.000000000E-00_JPRB,& |
|---|
| 1051 | & 0.000000000E-00_JPRB,& |
|---|
| 1052 | & 0.000000000E-00_JPRB,& |
|---|
| 1053 | & 0.000000000E-00_JPRB/) |
|---|
| 1054 | |
|---|
| 1055 | ZBPAD6(2, 1, :) = (/& |
|---|
| 1056 | & 0.100000000E+01_JPRB,& |
|---|
| 1057 | & 0.000000000E-00_JPRB,& |
|---|
| 1058 | & 0.000000000E-00_JPRB,& |
|---|
| 1059 | & 0.000000000E-00_JPRB,& |
|---|
| 1060 | & 0.000000000E-00_JPRB,& |
|---|
| 1061 | & 0.000000000E-00_JPRB,& |
|---|
| 1062 | & 0.000000000E-00_JPRB/) |
|---|
| 1063 | ZBPAD6(2, 2, :) = (/& |
|---|
| 1064 | & 0.100000000E+01_JPRB,& |
|---|
| 1065 | & 0.000000000E-00_JPRB,& |
|---|
| 1066 | & 0.000000000E-00_JPRB,& |
|---|
| 1067 | & 0.000000000E-00_JPRB,& |
|---|
| 1068 | & 0.000000000E-00_JPRB,& |
|---|
| 1069 | & 0.000000000E-00_JPRB,& |
|---|
| 1070 | & 0.000000000E-00_JPRB/) |
|---|
| 1071 | ZBPAD6(2, 3, :) = (/& |
|---|
| 1072 | & 0.100000000E+01_JPRB,& |
|---|
| 1073 | & 0.000000000E-00_JPRB,& |
|---|
| 1074 | & 0.000000000E-00_JPRB,& |
|---|
| 1075 | & 0.000000000E-00_JPRB,& |
|---|
| 1076 | & 0.000000000E-00_JPRB,& |
|---|
| 1077 | & 0.000000000E-00_JPRB,& |
|---|
| 1078 | & 0.000000000E-00_JPRB/) |
|---|
| 1079 | |
|---|
| 1080 | ZRAY6(2,:)= (/& |
|---|
| 1081 | & -.895408E-03_JPRB, .283130E+01_JPRB,-.867917E+01_JPRB,& |
|---|
| 1082 | & .148672E+02_JPRB,-.127282E+02_JPRB, .422212E+01_JPRB/) |
|---|
| 1083 | ZRBY6(2,:)= (/& |
|---|
| 1084 | & .55503070E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,& |
|---|
| 1085 | & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/) |
|---|
| 1086 | |
|---|
| 1087 | ZSWCE6(2) = 0.0_JPRB |
|---|
| 1088 | ZSWCP6(2) = 0.0_JPRB |
|---|
| 1089 | |
|---|
| 1090 | NEXPO3(2) = 7 |
|---|
| 1091 | REXPO3(2, 1, :) =(/& |
|---|
| 1092 | & 0.043801E+00_JPRB, 0.078893E+00_JPRB, 0.036839E+00_JPRB,& |
|---|
| 1093 | & 0.022503E+00_JPRB, 0.042333E+00_JPRB, 0.037870E+00_JPRB,& |
|---|
| 1094 | & 0.737762E+00_JPRB /) |
|---|
| 1095 | REXPO3(2, 2, :) =(/& |
|---|
| 1096 | & 0.234249E+01_JPRB, 0.125170E+00_JPRB, 0.549527E+02_JPRB,& |
|---|
| 1097 | & 0.257041E+03_JPRB, 0.476838E+00_JPRB, 0.911993E+01_JPRB,& |
|---|
| 1098 | & 0.000000E+00_JPRB /) |
|---|
| 1099 | |
|---|
| 1100 | !* 3rd spectral interval: Visible (0.44 - 0.69 Micron) |
|---|
| 1101 | |
|---|
| 1102 | ZSUN6(3) = 0.322135_JPRB |
|---|
| 1103 | |
|---|
| 1104 | ZD6(3,:)= (/ 0.800000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB /) |
|---|
| 1105 | |
|---|
| 1106 | ZAPAD6(3, 1, :) = (/& |
|---|
| 1107 | & 0.1762097E+03_JPRB,& |
|---|
| 1108 | & 0.1641762E+03_JPRB,& |
|---|
| 1109 | & 0.8687919E+02_JPRB,& |
|---|
| 1110 | & 0.0000000E-00_JPRB,& |
|---|
| 1111 | & 0.0000000E-00_JPRB,& |
|---|
| 1112 | & 0.0000000E-00_JPRB,& |
|---|
| 1113 | & 0.0000000E-00_JPRB/) |
|---|
| 1114 | ZAPAD6(3, 2, :) = (/& |
|---|
| 1115 | & 0.5581224E+00_JPRB,& |
|---|
| 1116 | & 0.1748430E+03_JPRB,& |
|---|
| 1117 | & 0.1134123E+04_JPRB,& |
|---|
| 1118 | & 0.3490429E+03_JPRB,& |
|---|
| 1119 | & 0.0000000E-00_JPRB,& |
|---|
| 1120 | & 0.0000000E-00_JPRB,& |
|---|
| 1121 | & 0.0000000E-00_JPRB/) |
|---|
| 1122 | ZAPAD6(3, 3, :) = (/& |
|---|
| 1123 | & 0.000000000E-00_JPRB,& |
|---|
| 1124 | & 0.000000000E-00_JPRB,& |
|---|
| 1125 | & 0.000000000E-00_JPRB,& |
|---|
| 1126 | & 0.000000000E-00_JPRB,& |
|---|
| 1127 | & 0.000000000E-00_JPRB,& |
|---|
| 1128 | & 0.000000000E-00_JPRB,& |
|---|
| 1129 | & 0.000000000E-00_JPRB/) |
|---|
| 1130 | |
|---|
| 1131 | ZBPAD6(3, 1, :) = (/& |
|---|
| 1132 | & 0.1762097E+03_JPRB,& |
|---|
| 1133 | & 0.1663950E+03_JPRB,& |
|---|
| 1134 | & 0.8939724E+02_JPRB,& |
|---|
| 1135 | & 0.1000000E+01_JPRB,& |
|---|
| 1136 | & 0.0000000E-00_JPRB,& |
|---|
| 1137 | & 0.0000000E-00_JPRB,& |
|---|
| 1138 | & 0.0000000E-00_JPRB/) |
|---|
| 1139 | ZBPAD6(3, 2, :) = (/& |
|---|
| 1140 | & 0.5581224E+00_JPRB,& |
|---|
| 1141 | & 0.1749251E+03_JPRB,& |
|---|
| 1142 | & 0.1159910E+04_JPRB,& |
|---|
| 1143 | & 0.3893268E+03_JPRB,& |
|---|
| 1144 | & 0.1000000E+01_JPRB,& |
|---|
| 1145 | & 0.0000000E-00_JPRB,& |
|---|
| 1146 | & 0.0000000E-00_JPRB/) |
|---|
| 1147 | ZBPAD6(3, 3, :) = (/& |
|---|
| 1148 | & 0.100000000E+01_JPRB,& |
|---|
| 1149 | & 0.000000000E-00_JPRB,& |
|---|
| 1150 | & 0.000000000E-00_JPRB,& |
|---|
| 1151 | & 0.000000000E-00_JPRB,& |
|---|
| 1152 | & 0.000000000E-00_JPRB,& |
|---|
| 1153 | & 0.000000000E-00_JPRB,& |
|---|
| 1154 | & 0.000000000E-00_JPRB/) |
|---|
| 1155 | |
|---|
| 1156 | ZRAY6(3,:)= (/& |
|---|
| 1157 | & .474852E-01_JPRB, .605110E+00_JPRB,-.248070E+01_JPRB,& |
|---|
| 1158 | & .485580E+01_JPRB,-.444043E+01_JPRB, .152369E+01_JPRB/) |
|---|
| 1159 | ZRBY6(3,:)= (/& |
|---|
| 1160 | & .10528199E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,& |
|---|
| 1161 | & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/) |
|---|
| 1162 | |
|---|
| 1163 | ZSWCE6(3) = 0.0_JPRB |
|---|
| 1164 | ZSWCP6(3) = 0.0_JPRB |
|---|
| 1165 | |
|---|
| 1166 | NEXPO3(3) = 6 |
|---|
| 1167 | REXPO3(3, 1, :) =(/& |
|---|
| 1168 | & 0.063442E+00_JPRB, 0.058550E+00_JPRB, 0.237534E+00_JPRB,& |
|---|
| 1169 | & 0.412292E+00_JPRB, 0.126141E+00_JPRB, 0.102041E+00_JPRB,& |
|---|
| 1170 | & 0.000000E+00_JPRB /) |
|---|
| 1171 | REXPO3(3, 2, :) =(/& |
|---|
| 1172 | & 0.125170E+00_JPRB, 0.119209E-01_JPRB, 0.119209E+00_JPRB,& |
|---|
| 1173 | & 0.417233E-01_JPRB, 0.894070E-01_JPRB, 0.000000E+00_JPRB,& |
|---|
| 1174 | & 0.000000E+00_JPRB /) |
|---|
| 1175 | |
|---|
| 1176 | ! ---------------------------------------------------------------- |
|---|
| 1177 | |
|---|
| 1178 | !* Near-Infrared (0.69 - 4.0 Microns) is sub-divided into: |
|---|
| 1179 | |
|---|
| 1180 | ! ---------------------------------------------------------------- |
|---|
| 1181 | |
|---|
| 1182 | !* 0.69 - 1.19 Micron |
|---|
| 1183 | |
|---|
| 1184 | !* UMG is O2 only |
|---|
| 1185 | |
|---|
| 1186 | ZSUN6(4) = 0.326158_JPRB |
|---|
| 1187 | |
|---|
| 1188 | ZD6(4,:)= (/ 0.000000000_JPRB, 0.900000000_JPRB, 1.000000000_JPRB /) |
|---|
| 1189 | |
|---|
| 1190 | ZAPAD6(4, 1, :) = (/& |
|---|
| 1191 | & 0.1335726E+02_JPRB,& |
|---|
| 1192 | & 0.2939136E+04_JPRB,& |
|---|
| 1193 | & 0.4010585E+05_JPRB,& |
|---|
| 1194 | & 0.7195030E+05_JPRB,& |
|---|
| 1195 | & 0.1648338E+05_JPRB,& |
|---|
| 1196 | & 0.3373738E+03_JPRB,& |
|---|
| 1197 | & 0.0000000E+00_JPRB/) |
|---|
| 1198 | ZAPAD6(4, 2, :) = (/& |
|---|
| 1199 | & 0.2001271E-01_JPRB,& |
|---|
| 1200 | & 0.2480831E+01_JPRB,& |
|---|
| 1201 | & 0.3444162E+02_JPRB,& |
|---|
| 1202 | & 0.4788946E+02_JPRB,& |
|---|
| 1203 | & 0.0000000E+00_JPRB,& |
|---|
| 1204 | & 0.0000000E+00_JPRB,& |
|---|
| 1205 | & 0.0000000E+00_JPRB/) |
|---|
| 1206 | ZAPAD6(4, 3, :) = (/& |
|---|
| 1207 | & 0.000000000E+00_JPRB,& |
|---|
| 1208 | & 0.000000000E+00_JPRB,& |
|---|
| 1209 | & 0.000000000E+00_JPRB,& |
|---|
| 1210 | & 0.000000000E+00_JPRB,& |
|---|
| 1211 | & 0.000000000E+00_JPRB,& |
|---|
| 1212 | & 0.000000000E+00_JPRB,& |
|---|
| 1213 | & 0.000000000E+00_JPRB/) |
|---|
| 1214 | |
|---|
| 1215 | ZBPAD6(4, 1, :) = (/& |
|---|
| 1216 | & 0.1335726E+02_JPRB,& |
|---|
| 1217 | & 0.2942327E+04_JPRB,& |
|---|
| 1218 | & 0.4077237E+05_JPRB,& |
|---|
| 1219 | & 0.7749017E+05_JPRB,& |
|---|
| 1220 | & 0.2123132E+05_JPRB,& |
|---|
| 1221 | & 0.6659687E+03_JPRB,& |
|---|
| 1222 | & 0.1000000E+01_JPRB/) |
|---|
| 1223 | ZBPAD6(4, 2, :) = (/& |
|---|
| 1224 | & 0.2001271E-01_JPRB,& |
|---|
| 1225 | & 0.2549067E+01_JPRB,& |
|---|
| 1226 | & 0.3752433E+02_JPRB,& |
|---|
| 1227 | & 0.6276637E+02_JPRB,& |
|---|
| 1228 | & 0.1000000E+01_JPRB,& |
|---|
| 1229 | & 0.0000000E+00_JPRB,& |
|---|
| 1230 | & 0.0000000E+00_JPRB/) |
|---|
| 1231 | ZBPAD6(4, 3, :) = (/& |
|---|
| 1232 | & 1.000000000E+00_JPRB,& |
|---|
| 1233 | & 0.000000000E+00_JPRB,& |
|---|
| 1234 | & 0.000000000E+00_JPRB,& |
|---|
| 1235 | & 0.000000000E+00_JPRB,& |
|---|
| 1236 | & 0.000000000E+00_JPRB,& |
|---|
| 1237 | & 0.000000000E+00_JPRB,& |
|---|
| 1238 | & 0.000000000E+00_JPRB/) |
|---|
| 1239 | |
|---|
| 1240 | ZRAY6(4,:)= (/& |
|---|
| 1241 | & .519710E-01_JPRB,-.150667E+00_JPRB, .154170E+01_JPRB,& |
|---|
| 1242 | & -.518901E+01_JPRB, .668011E+01_JPRB,-.294843E+01_JPRB/) |
|---|
| 1243 | ZRBY6(4,:)= (/& |
|---|
| 1244 | & .16436996E-01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,& |
|---|
| 1245 | & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/) |
|---|
| 1246 | |
|---|
| 1247 | ZSWCE6(4) = 0.0_JPRB |
|---|
| 1248 | ZSWCP6(4) = 0.0_JPRB |
|---|
| 1249 | |
|---|
| 1250 | NEXPO3(4) = 4 |
|---|
| 1251 | REXPO3(4, 1, :) =(/& |
|---|
| 1252 | & 0.000074E+00_JPRB, 0.320194E+00_JPRB, 0.082915E+00_JPRB,& |
|---|
| 1253 | & 0.596816E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,& |
|---|
| 1254 | & 0.000000E+00_JPRB /) |
|---|
| 1255 | REXPO3(4, 2, :) =(/& |
|---|
| 1256 | & 0.232458E+00_JPRB, 0.119209E-01_JPRB, 0.178814E-01_JPRB,& |
|---|
| 1257 | & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,& |
|---|
| 1258 | & 0.000000E+00_JPRB /) |
|---|
| 1259 | |
|---|
| 1260 | ! ---------------------------------------------------------------- |
|---|
| 1261 | |
|---|
| 1262 | !* 1.19 - 2.38 Microns |
|---|
| 1263 | |
|---|
| 1264 | !* UMG is CO2 only |
|---|
| 1265 | |
|---|
| 1266 | ZSUN6(5) = 0.180608_JPRB |
|---|
| 1267 | |
|---|
| 1268 | ZD6(5,:)= (/ 0.000000000_JPRB, 0.800000000_JPRB, 1.000000000_JPRB /) |
|---|
| 1269 | |
|---|
| 1270 | ZAPAD6(5, 1, :) = (/& |
|---|
| 1271 | & 0.3325841E-03_JPRB,& |
|---|
| 1272 | & 0.6194496E+00_JPRB,& |
|---|
| 1273 | & 0.1497138E+03_JPRB,& |
|---|
| 1274 | & 0.2314864E+04_JPRB,& |
|---|
| 1275 | & 0.2380109E+04_JPRB,& |
|---|
| 1276 | & 0.9553823E+02_JPRB,& |
|---|
| 1277 | & 0.0000000E+00_JPRB/) |
|---|
| 1278 | ZAPAD6(5, 2, :) = (/& |
|---|
| 1279 | & 0.4552471E-03_JPRB,& |
|---|
| 1280 | & 0.4084154E+00_JPRB,& |
|---|
| 1281 | & 0.6114905E+01_JPRB,& |
|---|
| 1282 | & 0.7102540E+01_JPRB,& |
|---|
| 1283 | & 0.0000000E+00_JPRB,& |
|---|
| 1284 | & 0.0000000E+00_JPRB,& |
|---|
| 1285 | & 0.0000000E+00_JPRB/) |
|---|
| 1286 | ZAPAD6(5, 3, :) = (/& |
|---|
| 1287 | & 0.000000000E+00_JPRB,& |
|---|
| 1288 | & 0.000000000E+00_JPRB,& |
|---|
| 1289 | & 0.000000000E+00_JPRB,& |
|---|
| 1290 | & 0.000000000E+00_JPRB,& |
|---|
| 1291 | & 0.000000000E+00_JPRB,& |
|---|
| 1292 | & 0.000000000E+00_JPRB,& |
|---|
| 1293 | & 0.000000000E+00_JPRB/) |
|---|
| 1294 | |
|---|
| 1295 | ZBPAD6(5, 1, :) = (/& |
|---|
| 1296 | & 0.3325841E-03_JPRB,& |
|---|
| 1297 | & 0.6231947E+00_JPRB,& |
|---|
| 1298 | & 0.1553098E+03_JPRB,& |
|---|
| 1299 | & 0.2822458E+04_JPRB,& |
|---|
| 1300 | & 0.3885194E+04_JPRB,& |
|---|
| 1301 | & 0.2700235E+03_JPRB,& |
|---|
| 1302 | & 0.1000000E+01_JPRB/) |
|---|
| 1303 | ZBPAD6(5, 2, :) = (/& |
|---|
| 1304 | & 0.4552471E-03_JPRB,& |
|---|
| 1305 | & 0.4088242E+00_JPRB,& |
|---|
| 1306 | & 0.6411905E+01_JPRB,& |
|---|
| 1307 | & 0.9444439E+01_JPRB,& |
|---|
| 1308 | & 0.1000000E+01_JPRB,& |
|---|
| 1309 | & 0.0000000E+00_JPRB,& |
|---|
| 1310 | & 0.0000000E+00_JPRB/) |
|---|
| 1311 | ZBPAD6(5, 3, :) = (/& |
|---|
| 1312 | & 1.000000000E+00_JPRB,& |
|---|
| 1313 | & 0.000000000E+00_JPRB,& |
|---|
| 1314 | & 0.000000000E+00_JPRB,& |
|---|
| 1315 | & 0.000000000E+00_JPRB,& |
|---|
| 1316 | & 0.000000000E+00_JPRB,& |
|---|
| 1317 | & 0.000000000E+00_JPRB,& |
|---|
| 1318 | & 0.000000000E+00_JPRB/) |
|---|
| 1319 | |
|---|
| 1320 | ZRAY6(5,:)= (/& |
|---|
| 1321 | & .288474E-02_JPRB, .133984E-02_JPRB,-.152812E-01_JPRB,& |
|---|
| 1322 | & .540845E-01_JPRB,-.747444E-01_JPRB, .364285E-01_JPRB/) |
|---|
| 1323 | ZRBY6(5,:)= (/& |
|---|
| 1324 | & .18073079E-02_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,& |
|---|
| 1325 | & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/) |
|---|
| 1326 | |
|---|
| 1327 | ZSWCE6(5) = 0.0_JPRB |
|---|
| 1328 | ZSWCP6(5) = 0.0_JPRB |
|---|
| 1329 | |
|---|
| 1330 | NEXPO3(5) = 0 |
|---|
| 1331 | REXPO3(5, 1, :) =(/& |
|---|
| 1332 | & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,& |
|---|
| 1333 | & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,& |
|---|
| 1334 | & 0.000000E+00_JPRB /) |
|---|
| 1335 | REXPO3(5, 2, :) =(/& |
|---|
| 1336 | & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,& |
|---|
| 1337 | & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,& |
|---|
| 1338 | & 0.000000E+00_JPRB /) |
|---|
| 1339 | |
|---|
| 1340 | ! ---------------------------------------------------------------- |
|---|
| 1341 | |
|---|
| 1342 | !* 2.38 - 4.00 Microns |
|---|
| 1343 | |
|---|
| 1344 | ZSUN6(6) = 0.033474_JPRB |
|---|
| 1345 | |
|---|
| 1346 | ZD6(6,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /) |
|---|
| 1347 | |
|---|
| 1348 | ZAPAD6(6, 1, :) = (/& |
|---|
| 1349 | & 0.2122889E-06_JPRB,& |
|---|
| 1350 | & 0.9030576E-03_JPRB,& |
|---|
| 1351 | & 0.2431282E+00_JPRB,& |
|---|
| 1352 | & 0.4901345E+01_JPRB,& |
|---|
| 1353 | & 0.3996347E+01_JPRB,& |
|---|
| 1354 | & 0.3910227E+01_JPRB,& |
|---|
| 1355 | & 0.0000000E+00_JPRB/) |
|---|
| 1356 | ZAPAD6(6, 2, :) = (/& |
|---|
| 1357 | & 0.1215163E-03_JPRB,& |
|---|
| 1358 | & 0.1222574E+00_JPRB,& |
|---|
| 1359 | & 0.9382420E+01_JPRB,& |
|---|
| 1360 | & 0.6875727E+02_JPRB,& |
|---|
| 1361 | & 0.2746421E+02_JPRB,& |
|---|
| 1362 | & 0.0000000E+00_JPRB,& |
|---|
| 1363 | & 0.0000000E+00_JPRB/) |
|---|
| 1364 | ZAPAD6(6, 3, :) = (/& |
|---|
| 1365 | & 0.263068898E+02_JPRB,& |
|---|
| 1366 | & 0.146425875E+03_JPRB,& |
|---|
| 1367 | & 0.860137809E+02_JPRB,& |
|---|
| 1368 | & 0.000000000E+00_JPRB,& |
|---|
| 1369 | & 0.000000000E+00_JPRB,& |
|---|
| 1370 | & 0.000000000E+00_JPRB,& |
|---|
| 1371 | & 0.000000000E+00_JPRB/) |
|---|
| 1372 | |
|---|
| 1373 | ZBPAD6(6, 1, :) = (/& |
|---|
| 1374 | & 0.2122889E-06_JPRB,& |
|---|
| 1375 | & 0.9379083E-03_JPRB,& |
|---|
| 1376 | & 0.2957335E+00_JPRB,& |
|---|
| 1377 | & 0.8747190E+01_JPRB,& |
|---|
| 1378 | & 0.1015794E+02_JPRB,& |
|---|
| 1379 | & 0.1361277E+02_JPRB,& |
|---|
| 1380 | & 0.1000000E+01_JPRB/) |
|---|
| 1381 | ZBPAD6(6, 2, :) = (/& |
|---|
| 1382 | & 0.1215163E-03_JPRB,& |
|---|
| 1383 | & 0.1255648E+00_JPRB,& |
|---|
| 1384 | & 0.1060119E+02_JPRB,& |
|---|
| 1385 | & 0.8414439E+02_JPRB,& |
|---|
| 1386 | & 0.4299438E+02_JPRB,& |
|---|
| 1387 | & 0.1000000E+01_JPRB,& |
|---|
| 1388 | & 0.0000000E+00_JPRB/) |
|---|
| 1389 | ZBPAD6(6, 3, :) = (/& |
|---|
| 1390 | & 0.263068898E+02_JPRB,& |
|---|
| 1391 | & 0.152569217E+03_JPRB,& |
|---|
| 1392 | & 0.976791971E+02_JPRB,& |
|---|
| 1393 | & 0.100000000E+01_JPRB,& |
|---|
| 1394 | & 0.000000000E+00_JPRB,& |
|---|
| 1395 | & 0.000000000E+00_JPRB,& |
|---|
| 1396 | & 0.000000000E+00_JPRB/) |
|---|
| 1397 | |
|---|
| 1398 | ZRAY6(6,:)= (/& |
|---|
| 1399 | & .493545E-03_JPRB, .355696E-03_JPRB,-.284543E-02_JPRB,& |
|---|
| 1400 | & .741248E-02_JPRB,-.773691E-02_JPRB, .370679E-02_JPRB/) |
|---|
| 1401 | ZRBY6(6,:)= (/& |
|---|
| 1402 | & .13618247E-03_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,& |
|---|
| 1403 | & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/) |
|---|
| 1404 | |
|---|
| 1405 | ZSWCE6(6) = 0.0_JPRB |
|---|
| 1406 | ZSWCP6(6) = 0.0_JPRB |
|---|
| 1407 | |
|---|
| 1408 | NEXPO3(6) = 0 |
|---|
| 1409 | REXPO3(6, 1, :) =(/& |
|---|
| 1410 | & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,& |
|---|
| 1411 | & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,& |
|---|
| 1412 | & 0.000000E+00_JPRB /) |
|---|
| 1413 | REXPO3(6, 2, :) =(/& |
|---|
| 1414 | & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,& |
|---|
| 1415 | & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,& |
|---|
| 1416 | & 0.000000E+00_JPRB /) |
|---|
| 1417 | |
|---|
| 1418 | !===================================================================== |
|---|
| 1419 | !* 2.4 OPTICAL PARAMETERS FOR RAIN DROPS |
|---|
| 1420 | ! Savijarvi et al. (1996) |
|---|
| 1421 | |
|---|
| 1422 | ZRTO1 = 0.003_JPRB |
|---|
| 1423 | ZRTO2 = -0.22_JPRB |
|---|
| 1424 | ! CAUTION JUST TEMPORARY PARAMETERS |
|---|
| 1425 | ZROMA2(1:2)= (/ 0.00008_JPRB , 0.264_JPRB /) |
|---|
| 1426 | ZROMB2(1:2)= (/ 0.23_JPRB , 0.09_JPRB /) |
|---|
| 1427 | ZRASY2(1:2)= (/ 0.88_JPRB , 0.94_JPRB /) |
|---|
| 1428 | |
|---|
| 1429 | ZRA2(1:2)= (/ 1.5_JPRB , 1.5_JPRB /) |
|---|
| 1430 | ZRB2(1:2)= (/ 0.50_JPRB , 1.13_JPRB /) |
|---|
| 1431 | ZRC2(1:2)= (/ 5.58E-7_JPRB , 8.55E-4_JPRB/) |
|---|
| 1432 | ZRD2(1:2)= (/ 1.25E-7_JPRB , 1.28E-3_JPRB/) |
|---|
| 1433 | ZRE2(1:2)= (/ 0.841_JPRB , 0.786_JPRB /) |
|---|
| 1434 | ZRF2(1:2)= (/ 2.08E-3_JPRB , 5.32E-3_JPRB/) |
|---|
| 1435 | |
|---|
| 1436 | ZROMA4(1:4)= (/ 0.00008_JPRB , 0.0105_JPRB , 0.264_JPRB , 0.465_JPRB /) |
|---|
| 1437 | ZROMB4(1:4)= (/ 0.23_JPRB , 0.22_JPRB , 0.09_JPRB , 0.001_JPRB /) |
|---|
| 1438 | ZRASY4(1:4)= (/ 0.88_JPRB , 0.89_JPRB , 0.94_JPRB , 0.97_JPRB /) |
|---|
| 1439 | |
|---|
| 1440 | ZRA4(1:4)= (/ 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB /) |
|---|
| 1441 | ZRB4(1:4)= (/ 0.50_JPRB , 0.78_JPRB , 1.13_JPRB , 2.00_JPRB /) |
|---|
| 1442 | ZRC4(1:4)= (/ 5.58E-7_JPRB , 2.18E-5_JPRB, 8.55E-4_JPRB, 1.94E-1_JPRB /) |
|---|
| 1443 | ZRD4(1:4)= (/ 1.25E-7_JPRB , 2.25E-5_JPRB, 1.28E-3_JPRB, 8.04E-3_JPRB /) |
|---|
| 1444 | ZRE4(1:4)= (/ 0.841_JPRB , 0.821_JPRB , 0.786_JPRB , 0.820_JPRB /) |
|---|
| 1445 | ZRF4(1:4)= (/ 2.08E-3_JPRB , 3.06E-3_JPRB, 5.32E-3_JPRB, 5.59E-3_JPRB /) |
|---|
| 1446 | |
|---|
| 1447 | ZROMA6(1:6)= (/ 0.00008_JPRB , 0.00008_JPRB , 0.00008_JPRB , 0.0105_JPRB , 0.264_JPRB , 0.465_JPRB /) |
|---|
| 1448 | ZROMB6(1:6)= (/ 0.23_JPRB , 0.23_JPRB , 0.23_JPRB , 0.22_JPRB , 0.09_JPRB , 0.001_JPRB /) |
|---|
| 1449 | ZRASY6(1:6)= (/ 0.88_JPRB , 0.88_JPRB , 0.88_JPRB , 0.89_JPRB , 0.94_JPRB , 0.97_JPRB /) |
|---|
| 1450 | |
|---|
| 1451 | ZRA6(1:6)= (/ 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB /) |
|---|
| 1452 | ZRB6(1:6)= (/ 0.50_JPRB , 0.50_JPRB , 0.50_JPRB , 0.78_JPRB , 1.13_JPRB , 2.00_JPRB /) |
|---|
| 1453 | ZRC6(1:6)= (/ 5.58E-7_JPRB , 5.58E-7_JPRB , 5.58E-7_JPRB , 2.18E-5_JPRB, 8.55E-4_JPRB, 1.94E-1_JPRB /) |
|---|
| 1454 | ZRD6(1:6)= (/ 1.25E-7_JPRB , 1.25E-7_JPRB , 1.25E-7_JPRB , 2.25E-5_JPRB, 1.28E-3_JPRB, 8.04E-3_JPRB /) |
|---|
| 1455 | ZRE6(1:6)= (/ 0.841_JPRB , 0.841_JPRB , 0.841_JPRB , 0.821_JPRB , 0.786_JPRB , 0.820_JPRB /) |
|---|
| 1456 | ZRF6(1:6)= (/ 2.08E-3_JPRB , 2.08E-3_JPRB , 2.08E-3_JPRB , 3.06E-3_JPRB, 5.32E-3_JPRB, 5.59E-3_JPRB /) |
|---|
| 1457 | |
|---|
| 1458 | !===================================================================== |
|---|
| 1459 | ! 3. FOURTEEN SPECTRAL INTERVALS FOR RRTM_SW |
|---|
| 1460 | |
|---|
| 1461 | ! SW : 14 spectral intervals : mapping to SW6 spectral intervals |
|---|
| 1462 | ! 0.185-0.25-0.44-0.69-1.19-2.38-4.00 |
|---|
| 1463 | ! 3.846 - 3.077 6 |
|---|
| 1464 | ! 3.077 - 2.500 6 |
|---|
| 1465 | ! 2.500 - 2.150 5 |
|---|
| 1466 | ! 2.150 - 1.942 5 |
|---|
| 1467 | ! 1.942 - 1.626 5 |
|---|
| 1468 | ! 1.626 - 1.299 5 |
|---|
| 1469 | ! 1.299 - 1.242 5 |
|---|
| 1470 | ! 1.242 - 0.7782 4 |
|---|
| 1471 | ! 0.7782- 0.6250 4 |
|---|
| 1472 | ! 0.6250- 0.4415 3 |
|---|
| 1473 | ! 0.4415- 0.3448 2 |
|---|
| 1474 | ! 0.3448- 0.2632 2 |
|---|
| 1475 | ! 0.2632- 0.2000 1 |
|---|
| 1476 | ! 12.195 - 3.846 6 |
|---|
| 1477 | |
|---|
| 1478 | IMPSRTM(1:14)= (/ 6,6,5,5,5,5,5,4,4,3,2,2,1,6 /) |
|---|
| 1479 | !===================================================================== |
|---|
| 1480 | |
|---|
| 1481 | !* 4. SET VALUES. |
|---|
| 1482 | ! ----------- |
|---|
| 1483 | |
|---|
| 1484 | IF (LONEWSW) THEN |
|---|
| 1485 | ZPDH2O = ZPDHNIS |
|---|
| 1486 | ZPDUMG = ZPDUNIS |
|---|
| 1487 | ZPRH2O = ZPRHNIS |
|---|
| 1488 | ZPRUMG = ZPRUNIS |
|---|
| 1489 | RTDH2O = ZTDHNIS |
|---|
| 1490 | RTDUMG = ZTDUNIS |
|---|
| 1491 | RTH2O = ZTHNIS |
|---|
| 1492 | RTUMG = ZTUNIS |
|---|
| 1493 | ENDIF |
|---|
| 1494 | |
|---|
| 1495 | RPDH1=ZPDH2O+1.0_JPRB |
|---|
| 1496 | RPDU1=ZPDUMG+1.0_JPRB |
|---|
| 1497 | ZH2O=1.0_JPRB/( 10._JPRB* RG * RPDH1 ) |
|---|
| 1498 | ZUMG=1.0_JPRB/( 10._JPRB* RG * RPDU1 ) |
|---|
| 1499 | RPNU = ZUMG/(ZPRUMG**ZPDUMG) |
|---|
| 1500 | RPNH = ZH2O/(ZPRH2O**ZPDH2O) |
|---|
| 1501 | |
|---|
| 1502 | IF(.NOT.ALLOCATED(RSUN)) ALLOCATE (RSUN(NTSW)) |
|---|
| 1503 | |
|---|
| 1504 | DO JW=1,ISW6 |
|---|
| 1505 | RSUN(JW)=0.0_JPRB |
|---|
| 1506 | RSWCE(JW)=0.0_JPRB |
|---|
| 1507 | RSWCP(JW)=0.0_JPRB |
|---|
| 1508 | |
|---|
| 1509 | DO JC3=1,3 |
|---|
| 1510 | D(JW,JC3)=0.0_JPRB |
|---|
| 1511 | ENDDO |
|---|
| 1512 | DO JC6=1,6 |
|---|
| 1513 | RRAY(JW,JC6)=0.0_JPRB |
|---|
| 1514 | ENDDO |
|---|
| 1515 | DO JI=1,3 |
|---|
| 1516 | DO JJ=1,7 |
|---|
| 1517 | APAD(JW,JI,JJ)=0.0_JPRB |
|---|
| 1518 | BPAD(JW,JI,JJ)=0.0_JPRB |
|---|
| 1519 | ENDDO |
|---|
| 1520 | ENDDO |
|---|
| 1521 | ENDDO |
|---|
| 1522 | |
|---|
| 1523 | IF (KSW == 2) THEN |
|---|
| 1524 | ZPDH2O = ZPDH2IS |
|---|
| 1525 | ZPDUMG = ZPDU2IS |
|---|
| 1526 | ZPRH2O = ZPRH2IS |
|---|
| 1527 | ZPRUMG = ZPRU2IS |
|---|
| 1528 | RTDH2O = ZTDH2IS |
|---|
| 1529 | RTDUMG = ZTDU2IS |
|---|
| 1530 | RTH2O = ZTH2IS |
|---|
| 1531 | RTUMG = ZTU2IS |
|---|
| 1532 | |
|---|
| 1533 | RADJUST=ZADJUS2 |
|---|
| 1534 | RHSRTA=ZRTO1 |
|---|
| 1535 | RHSRTB=ZRTO2 |
|---|
| 1536 | DO JW=1,KSW |
|---|
| 1537 | RSUN (JW)=ZSUN2(JW) |
|---|
| 1538 | |
|---|
| 1539 | RSWCE(JW)=ZSWCE2(JW) |
|---|
| 1540 | RSWCP(JW)=ZSWCP2(JW) |
|---|
| 1541 | |
|---|
| 1542 | RROMA(JW)=ZROMA2(JW) |
|---|
| 1543 | RROMB(JW)=ZROMB2(JW) |
|---|
| 1544 | RRASY(JW)=ZRASY2(JW) |
|---|
| 1545 | RHSRA(JW)=ZRA2(JW) |
|---|
| 1546 | RHSRB(JW)=ZRB2(JW) |
|---|
| 1547 | RHSRC(JW)=ZRC2(JW) |
|---|
| 1548 | RHSRD(JW)=ZRD2(JW) |
|---|
| 1549 | RHSRE(JW)=ZRE2(JW) |
|---|
| 1550 | RHSRF(JW)=ZRF2(JW) |
|---|
| 1551 | |
|---|
| 1552 | DO JC3=1,3 |
|---|
| 1553 | D(JW,JC3)=ZD2(JW,JC3) |
|---|
| 1554 | ENDDO |
|---|
| 1555 | DO JC6=1,6 |
|---|
| 1556 | RRAY(JW,JC6)=ZRAY2(JW,JC6) |
|---|
| 1557 | ENDDO |
|---|
| 1558 | DO JI=1,3 |
|---|
| 1559 | DO JJ=1,7 |
|---|
| 1560 | APAD(JW,JI,JJ)=ZAPAD2(JW,JI,JJ) |
|---|
| 1561 | BPAD(JW,JI,JJ)=ZBPAD2(JW,JI,JJ) |
|---|
| 1562 | ENDDO |
|---|
| 1563 | ENDDO |
|---|
| 1564 | ENDDO |
|---|
| 1565 | |
|---|
| 1566 | ELSEIF (KSW == 4) THEN |
|---|
| 1567 | ZPDH2O = ZPDH4IS |
|---|
| 1568 | ZPDUMG = ZPDU4IS |
|---|
| 1569 | ZPRH2O = ZPRH4IS |
|---|
| 1570 | ZPRUMG = ZPRU4IS |
|---|
| 1571 | RTDH2O = ZTDH4IS |
|---|
| 1572 | RTDUMG = ZTDU4IS |
|---|
| 1573 | RTH2O = ZTH4IS |
|---|
| 1574 | RTUMG = ZTU4IS |
|---|
| 1575 | |
|---|
| 1576 | RADJUST=ZADJUS4 |
|---|
| 1577 | RHSRTA=ZRTO1 |
|---|
| 1578 | RHSRTB=ZRTO2 |
|---|
| 1579 | DO JW=1,KSW |
|---|
| 1580 | RSUN (JW)=ZSUN4(JW) |
|---|
| 1581 | |
|---|
| 1582 | RSWCE(JW)=ZSWCE4(JW) |
|---|
| 1583 | RSWCP(JW)=ZSWCP4(JW) |
|---|
| 1584 | |
|---|
| 1585 | RROMA(JW)=ZROMA4(JW) |
|---|
| 1586 | RROMB(JW)=ZROMB4(JW) |
|---|
| 1587 | RRASY(JW)=ZRASY4(JW) |
|---|
| 1588 | RHSRA(JW)=ZRA4(JW) |
|---|
| 1589 | RHSRB(JW)=ZRB4(JW) |
|---|
| 1590 | RHSRC(JW)=ZRC4(JW) |
|---|
| 1591 | RHSRD(JW)=ZRD4(JW) |
|---|
| 1592 | RHSRE(JW)=ZRE4(JW) |
|---|
| 1593 | RHSRF(JW)=ZRF4(JW) |
|---|
| 1594 | |
|---|
| 1595 | DO JC3=1,3 |
|---|
| 1596 | D(JW,JC3)=ZD4(JW,JC3) |
|---|
| 1597 | ENDDO |
|---|
| 1598 | DO JC6=1,6 |
|---|
| 1599 | RRAY(JW,JC6)=ZRAY4(JW,JC6) |
|---|
| 1600 | ENDDO |
|---|
| 1601 | DO JI=1,3 |
|---|
| 1602 | DO JJ=1,7 |
|---|
| 1603 | APAD(JW,JI,JJ)=ZAPAD4(JW,JI,JJ) |
|---|
| 1604 | BPAD(JW,JI,JJ)=ZBPAD4(JW,JI,JJ) |
|---|
| 1605 | ENDDO |
|---|
| 1606 | ENDDO |
|---|
| 1607 | ENDDO |
|---|
| 1608 | |
|---|
| 1609 | ELSEIF (KSW == 6) THEN |
|---|
| 1610 | ZPDH2O = ZPDH6IS |
|---|
| 1611 | ZPDUMG = ZPDU6IS |
|---|
| 1612 | ZPRH2O = ZPRH6IS |
|---|
| 1613 | ZPRUMG = ZPRU6IS |
|---|
| 1614 | RTDH2O = ZTDH6IS |
|---|
| 1615 | RTDUMG = ZTDU6IS |
|---|
| 1616 | RTH2O = ZTH6IS |
|---|
| 1617 | RTUMG = ZTU6IS |
|---|
| 1618 | |
|---|
| 1619 | RADJUST=ZADJUS6 |
|---|
| 1620 | RHSRTA=ZRTO1 |
|---|
| 1621 | RHSRTB=ZRTO2 |
|---|
| 1622 | DO JW=1,KSW |
|---|
| 1623 | RSUN (JW)=ZSUN6(JW) |
|---|
| 1624 | |
|---|
| 1625 | RSWCE(JW)=ZSWCE6(JW) |
|---|
| 1626 | RSWCP(JW)=ZSWCP6(JW) |
|---|
| 1627 | |
|---|
| 1628 | RROMA(JW)=ZROMA6(JW) |
|---|
| 1629 | RROMB(JW)=ZROMB6(JW) |
|---|
| 1630 | RRASY(JW)=ZRASY6(JW) |
|---|
| 1631 | RHSRA(JW)=ZRA6(JW) |
|---|
| 1632 | RHSRB(JW)=ZRB6(JW) |
|---|
| 1633 | RHSRC(JW)=ZRC6(JW) |
|---|
| 1634 | RHSRD(JW)=ZRD6(JW) |
|---|
| 1635 | RHSRE(JW)=ZRE6(JW) |
|---|
| 1636 | RHSRF(JW)=ZRF6(JW) |
|---|
| 1637 | |
|---|
| 1638 | DO JC3=1,3 |
|---|
| 1639 | D(JW,JC3)=ZD6(JW,JC3) |
|---|
| 1640 | ENDDO |
|---|
| 1641 | IF (LRAYL) THEN |
|---|
| 1642 | DO JC6=1,6 |
|---|
| 1643 | RRAY(JW,JC6)=ZRAY6(JW,JC6) |
|---|
| 1644 | ENDDO |
|---|
| 1645 | ELSE |
|---|
| 1646 | DO JC6=1,6 |
|---|
| 1647 | RRAY(JW,JC6)=ZRBY6(JW,JC6) |
|---|
| 1648 | ENDDO |
|---|
| 1649 | ENDIF |
|---|
| 1650 | DO JI=1,3 |
|---|
| 1651 | DO JJ=1,7 |
|---|
| 1652 | APAD(JW,JI,JJ)=ZAPAD6(JW,JI,JJ) |
|---|
| 1653 | BPAD(JW,JI,JJ)=ZBPAD6(JW,JI,JJ) |
|---|
| 1654 | ENDDO |
|---|
| 1655 | ENDDO |
|---|
| 1656 | ENDDO |
|---|
| 1657 | |
|---|
| 1658 | ELSEIF (KTSW == 14 .AND. KSW == 14) THEN |
|---|
| 1659 | PRINT *,'SUSWN: 14-SPECTRAL INTERVALS --> RRTM_SW' |
|---|
| 1660 | |
|---|
| 1661 | ELSE |
|---|
| 1662 | PRINT *,'SUSWN: NTSW, NSW: ',KTSW,KSW |
|---|
| 1663 | CALL ABOR1('SUSWN: WRONG NUMBER OF SW INTERVALS') |
|---|
| 1664 | ENDIF |
|---|
| 1665 | |
|---|
| 1666 | ! ---------------------------------------------------------------- |
|---|
| 1667 | |
|---|
| 1668 | !* 5. PRINT VALUES. |
|---|
| 1669 | ! ------------ |
|---|
| 1670 | |
|---|
| 1671 | ! PRINT 9700,RPDH1,RPDU1,RPNU,RPNH |
|---|
| 1672 | !9700 FORMAT(1X,'RPDH1,RPDU1,RPNU,RPNH ',4F10.5) |
|---|
| 1673 | ! DO JW=1,KSW |
|---|
| 1674 | ! PRINT *,' SW INTERVAL ',JW |
|---|
| 1675 | ! PRINT 9701,JW,RSUN(JW),RSWCE(JW),RSWCP(JW) |
|---|
| 1676 | !9701 FORMAT(1X,I3,F10.7,2F10.5) |
|---|
| 1677 | ! PRINT 9702,JW,(D(JW,JC3),JC3=1,3) |
|---|
| 1678 | !9702 FORMAT(1X,I3,3F10.5) |
|---|
| 1679 | ! PRINT 9703,JW,(RRAY(JW,JC6),JC6=1,6) |
|---|
| 1680 | !9703 FORMAT(1X,I3,6E15.8) |
|---|
| 1681 | ! PRINT 9704,JW,(JABS,(APAD(JW,JABS,JJ),JJ=1,7),JABS=1,3) |
|---|
| 1682 | !9704 FORMAT(1X,I3,(I3,7E15.8)) |
|---|
| 1683 | ! PRINT 9706,JW,RROMA(JW),RROMB(JW),RRASY(JW),RHSRTA,RHSRTB |
|---|
| 1684 | !9706 FORMAT(1X,I3,' Rain Optics from Rain Rate: ',F8.5,F6.3,3F6.3) |
|---|
| 1685 | ! PRINT 9707,JW,RHSRA(JW),RHSRB(JW),RHSRC(JW),RHSRD(JW) |
|---|
| 1686 | ! & ,RHSRE(JW),RHSRF(JW) |
|---|
| 1687 | !9707 FORMAT(1X,I3,' Rain Optics as function of Reff: ',6E11.4) |
|---|
| 1688 | ! ENDDO |
|---|
| 1689 | |
|---|
| 1690 | ! ---------------------------------------------------------------- |
|---|
| 1691 | |
|---|
| 1692 | IF (LHOOK) CALL DR_HOOK('SUSWN',1,ZHOOK_HANDLE) |
|---|
| 1693 | END SUBROUTINE SUSWN |
|---|
