[2089] | 1 | SUBROUTINE SUCLOPN (KTSW,KSW , KLEV) |
---|
| 2 | |
---|
| 3 | !**** *SUCLOP* - INITIALIZE COMMON YOECLOP |
---|
| 4 | |
---|
| 5 | ! PURPOSE. |
---|
| 6 | ! -------- |
---|
| 7 | ! INITIALIZE YOMCLOP, WITH CLOUD OPTICAL PARAMETERS |
---|
| 8 | |
---|
| 9 | !** INTERFACE. |
---|
| 10 | ! ---------- |
---|
| 11 | ! *CALL* SUCLOPN |
---|
| 12 | ! FROM *SUECRAD* |
---|
| 13 | |
---|
| 14 | ! EXPLICIT ARGUMENTS : |
---|
| 15 | ! -------------------- |
---|
| 16 | ! NONE |
---|
| 17 | |
---|
| 18 | ! IMPLICIT ARGUMENTS : |
---|
| 19 | ! -------------------- |
---|
| 20 | ! COMMON YOECLOP |
---|
| 21 | |
---|
| 22 | ! METHOD. |
---|
| 23 | ! ------- |
---|
| 24 | ! SEE DOCUMENTATION |
---|
| 25 | |
---|
| 26 | ! EXTERNALS. |
---|
| 27 | ! ---------- |
---|
| 28 | ! NONE |
---|
| 29 | |
---|
| 30 | ! REFERENCE. |
---|
| 31 | ! ---------- |
---|
| 32 | ! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE |
---|
| 33 | ! "INTEGRATED FORECASTING SYSTEM" |
---|
| 34 | |
---|
| 35 | ! Fouquart, 1987: NATO ASI, 223-284 |
---|
| 36 | ! A. Slingo, 1989: J. Atmos. Sci., 46, 1419-1427 |
---|
| 37 | ! Ebert and Curry, 1992: J. Geophys. Res., 97D, 3831-3836 |
---|
| 38 | ! Sun and Shine, 1994: Quart. J. Roy. Meteor. Soc., 120, 111-138 |
---|
| 39 | ! Fu and Liou, 1993: J. Atmos. Sci., 50, 2008-2025 |
---|
| 40 | ! Fu, 1996: J. Climate, 9, 2058-2082 |
---|
| 41 | ! Fu et al., 1998: J. Climate, 11, 2223-2237 |
---|
| 42 | |
---|
| 43 | ! AUTHOR. |
---|
| 44 | ! ------- |
---|
| 45 | ! JEAN-JACQUES MORCRETTE *ECMWF* |
---|
| 46 | |
---|
| 47 | ! MODIFICATIONS. |
---|
| 48 | ! -------------- |
---|
| 49 | ! ORIGINAL : 92-02-29 |
---|
| 50 | ! 97-04-16 JJ Morcrette 2 and 4 spectral intervals |
---|
| 51 | ! 00-10-25 JJMorcrette 6 spectral intervals |
---|
| 52 | ! 01-01-16 JJ Morcrette ice cloud from Fu 96, Fu et al 98 |
---|
| 53 | |
---|
| 54 | ! ------------------------------------------------------------------ |
---|
| 55 | |
---|
| 56 | #include "tsmbkind.h" |
---|
| 57 | |
---|
| 58 | USE YOESW , ONLY : RYFWCA ,RYFWCB ,RYFWCC ,RYFWCD ,& |
---|
| 59 | &RYFWCE ,RYFWCF ,REBCUA ,REBCUB ,REBCUC ,& |
---|
| 60 | &REBCUD ,REBCUE ,REBCUF ,REBCUG ,REBCUH ,& |
---|
| 61 | &REBCUI ,REBCUJ ,RASWCA ,RASWCB ,RASWCC ,& |
---|
| 62 | &RASWCD ,RASWCE ,RASWCF ,RSUSHE ,RSUSHF ,& |
---|
| 63 | &RSUSHH ,RSUSHK ,RSUSHA ,RSUSHG ,RSUSHFA ,& |
---|
| 64 | &RSUSHC ,RSUSHD ,REFFIA ,REFFIB ,RHSAVI ,& |
---|
| 65 | &RFUAA0 ,RFUAA1 ,RFUBB0 ,RFUBB1 ,RFUBB2 ,& |
---|
| 66 | &RFUBB3 ,RFUCC0 ,RFUCC1 ,RFUCC2 ,RFUCC3 ,& |
---|
| 67 | &RFUETA ,RFULIO ,RFLAA0 ,RFLAA1 ,RFLBB0 ,& |
---|
| 68 | &RFLBB1 ,RFLBB2 ,RFLBB3 ,RFLCC0 ,RFLCC1 ,& |
---|
| 69 | &RFLCC2 ,RFLCC3 ,RFLDD0 ,RFLDD1 ,RFLDD2 ,& |
---|
| 70 | &RFLDD3 ,RLINLI ,RTIW ,RRIW |
---|
| 71 | USE YOERAD , ONLY : RAOVLP ,RBOVLP |
---|
| 72 | |
---|
| 73 | IMPLICIT NONE |
---|
| 74 | |
---|
| 75 | |
---|
| 76 | ! DUMMY INTEGER SCALARS |
---|
| 77 | INTEGER_M :: KSW |
---|
| 78 | INTEGER_M :: KTSW |
---|
| 79 | INTEGER_M :: KLEV |
---|
| 80 | INTEGER_M :: JC, JNU |
---|
| 81 | |
---|
| 82 | |
---|
| 83 | |
---|
| 84 | |
---|
| 85 | ! ----------------------------------------------------------------- |
---|
| 86 | REAL_B :: ZEBCUA2(2) ,ZEBCUB2(2) ,ZEBCUC2(2) ,ZEBCUD2(2)& |
---|
| 87 | &, ZEBCUE2(2) ,ZEBCUF2(2) ,ZYFWCA2(2) ,ZYFWCB2(2)& |
---|
| 88 | &, ZYFWCC2(2) ,ZYFWCD2(2) ,ZYFWCE2(2) ,ZYFWCF2(2)& |
---|
| 89 | &, ZASWCA2(2) ,ZASWCB2(2) ,ZASWCC2(2) ,ZASWCD2(2)& |
---|
| 90 | &, ZASWCE2(2) ,ZASWCF2(2) ,ZSUSHE2(2) ,ZSUSHF2(2)& |
---|
| 91 | &, ZSUSHH2(2) ,ZSUSHK2(2) ,ZSUSHA2(2) ,ZSUSHG2(2)& |
---|
| 92 | &, ZFLAA02(2) ,ZFLAA12(2) ,ZFLBB02(2) ,ZFLBB12(2)& |
---|
| 93 | &, ZFLBB22(2) ,ZFLBB32(2) ,ZFLCC02(2) ,ZFLCC12(2)& |
---|
| 94 | &, ZFLCC22(2) ,ZFLCC32(2) ,ZFLDD02(2) ,ZFLDD12(2)& |
---|
| 95 | &, ZFLDD22(2) ,ZFLDD32(2)& |
---|
| 96 | &, ZFUAA02(2) ,ZFUAA12(2) ,ZFUBB02(2) ,ZFUBB12(2)& |
---|
| 97 | &, ZFUBB22(2) ,ZFUBB32(2) ,ZFUCC02(2) ,ZFUCC12(2)& |
---|
| 98 | &, ZFUCC22(2) ,ZFUCC32(2) ,ZFUDD02(2) ,ZFUDD12(2)& |
---|
| 99 | &, ZFUDD22(2) ,ZFUDD32(2) |
---|
| 100 | |
---|
| 101 | REAL_B :: ZEBCUA4(4) ,ZEBCUB4(4) ,ZEBCUC4(4) ,ZEBCUD4(4)& |
---|
| 102 | &, ZEBCUE4(4) ,ZEBCUF4(4) ,ZYFWCA4(4) ,ZYFWCB4(4)& |
---|
| 103 | &, ZYFWCC4(4) ,ZYFWCD4(4) ,ZYFWCE4(4) ,ZYFWCF4(4)& |
---|
| 104 | &, ZASWCA4(4) ,ZASWCB4(4) ,ZASWCC4(4) ,ZASWCD4(4)& |
---|
| 105 | &, ZASWCE4(4) ,ZASWCF4(4) ,ZSUSHE4(4) ,ZSUSHF4(4)& |
---|
| 106 | &, ZSUSHH4(4) ,ZSUSHK4(4) ,ZSUSHA4(4) ,ZSUSHG4(4)& |
---|
| 107 | &, ZFLAA04(4) ,ZFLAA14(4) ,ZFLBB04(4) ,ZFLBB14(4)& |
---|
| 108 | &, ZFLBB24(4) ,ZFLBB34(4) ,ZFLCC04(4) ,ZFLCC14(4)& |
---|
| 109 | &, ZFLCC24(4) ,ZFLCC34(4) ,ZFLDD04(4) ,ZFLDD14(4)& |
---|
| 110 | &, ZFLDD24(4) ,ZFLDD34(4)& |
---|
| 111 | &, ZFUAA04(4) ,ZFUAA14(4) ,ZFUBB04(4) ,ZFUBB14(4)& |
---|
| 112 | &, ZFUBB24(4) ,ZFUBB34(4) ,ZFUCC04(4) ,ZFUCC14(4)& |
---|
| 113 | &, ZFUCC24(4) ,ZFUCC34(4) ,ZFUDD04(4) ,ZFUDD14(4)& |
---|
| 114 | &, ZFUDD24(4) ,ZFUDD34(4) |
---|
| 115 | |
---|
| 116 | REAL_B :: ZEBCUA6(6) ,ZEBCUB6(6) ,ZEBCUC6(6) ,ZEBCUD6(6)& |
---|
| 117 | &, ZEBCUE6(6) ,ZEBCUF6(6) ,ZYFWCA6(6) ,ZYFWCB6(6)& |
---|
| 118 | &, ZYFWCC6(6) ,ZYFWCD6(6) ,ZYFWCE6(6) ,ZYFWCF6(6)& |
---|
| 119 | &, ZASWCA6(6) ,ZASWCB6(6) ,ZASWCC6(6) ,ZASWCD6(6)& |
---|
| 120 | &, ZASWCE6(6) ,ZASWCF6(6) ,ZSUSHE6(6) ,ZSUSHF6(6)& |
---|
| 121 | &, ZSUSHH6(6) ,ZSUSHK6(6) ,ZSUSHA6(6) ,ZSUSHG6(6)& |
---|
| 122 | &, ZFLAA06(6) ,ZFLAA16(6) ,ZFLBB06(6) ,ZFLBB16(6)& |
---|
| 123 | &, ZFLBB26(6) ,ZFLBB36(6) ,ZFLCC06(6) ,ZFLCC16(6)& |
---|
| 124 | &, ZFLCC26(6) ,ZFLCC36(6) ,ZFLDD06(6) ,ZFLDD16(6)& |
---|
| 125 | &, ZFLDD26(6) ,ZFLDD36(6)& |
---|
| 126 | &, ZFUAA06(6) ,ZFUAA16(6) ,ZFUBB06(6) ,ZFUBB16(6)& |
---|
| 127 | &, ZFUBB26(6) ,ZFUBB36(6) ,ZFUCC06(6) ,ZFUCC16(6)& |
---|
| 128 | &, ZFUCC26(6) ,ZFUCC36(6) ,ZFUDD06(6) ,ZFUDD16(6)& |
---|
| 129 | &, ZFUDD26(6) ,ZFUDD36(6) |
---|
| 130 | |
---|
| 131 | |
---|
| 132 | REAL_B :: ZAOVLP(3), ZBOVLP(3) |
---|
| 133 | |
---|
| 134 | ! ----------------------------------------------------------------- |
---|
| 135 | |
---|
| 136 | !* 1. SHORTWAVE CLOUD OPTICAL PROPERTIES |
---|
| 137 | ! ---------------------------------- |
---|
| 138 | |
---|
| 139 | ! ------------------------------------------------------------------ |
---|
| 140 | |
---|
| 141 | !* 1.1 TWO SPECTRAL INTERVALS |
---|
| 142 | ! ---------------------- |
---|
| 143 | |
---|
| 144 | ! SW : 0.25 - 0.68 - 4.00 microns |
---|
| 145 | |
---|
| 146 | !* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992 |
---|
| 147 | |
---|
| 148 | ! optical properties |
---|
| 149 | ZEBCUA2 = (/ 3.448E-04_JPRB , 3.448E-04_JPRB /) |
---|
| 150 | ZEBCUB2 = (/ 2.431_JPRB , 2.431_JPRB /) |
---|
| 151 | ZEBCUC2 = (/ 0.00001_JPRB , 0.024366_JPRB /) |
---|
| 152 | ZEBCUD2 = (/ _ZERO_ , 2.487E-04_JPRB /) |
---|
| 153 | ZEBCUE2 = (/ 0.7661_JPRB , 0.7866_JPRB /) |
---|
| 154 | ZEBCUF2 = (/ 5.851E-04_JPRB , 5.937E-04_JPRB /) |
---|
| 155 | |
---|
| 156 | ! optical properties |
---|
| 157 | ! ZEBCUA2 = (/ 3.448E-04 , 3.448E-04 /) |
---|
| 158 | ! ZEBCUB2 = (/ 2.431 , 2.431 /) |
---|
| 159 | ! ZEBCUC2 = (/ 0.00001 , 0.035589 /) |
---|
| 160 | ! ZEBCUD2 = (/ 0. , 2.757E-04 /) |
---|
| 161 | ! ZEBCUE2 = (/ 0.7661 , 0.7921 /) |
---|
| 162 | ! ZEBCUF2 = (/ 5.851E-04 , 5.893E-04 /) |
---|
| 163 | |
---|
| 164 | !* Water cloud properties - from Fouquart (1987) |
---|
| 165 | |
---|
| 166 | ZYFWCA2 = (/ _ZERO_ , _ZERO_ /) |
---|
| 167 | ZYFWCB2 = (/ 1.5_JPRB , 1.5_JPRB /) |
---|
| 168 | ZYFWCC2 = (/ 0.9999_JPRB , 0.9988_JPRB /) |
---|
| 169 | ZYFWCD2 = (/ 5.000E-04_JPRB , 2.500E-03_JPRB /) |
---|
| 170 | ZYFWCE2 = (/ _HALF_ , 0.05_JPRB /) |
---|
| 171 | ZYFWCF2 = (/ 0.865_JPRB , 0.910_JPRB /) |
---|
| 172 | |
---|
| 173 | !* Water cloud properties - from Slingo (1989) |
---|
| 174 | |
---|
| 175 | ZASWCA2 = (/ 2.817_JPRB , 2.4584_JPRB /) |
---|
| 176 | ZASWCB2 = (/ 1.305_JPRB , 1.3994_JPRB /) |
---|
| 177 | ZASWCC2 = (/-5.62E-08_JPRB , 1.111E-02_JPRB /) |
---|
| 178 | ZASWCD2 = (/ 1.63E-07_JPRB , 8.613E-04_JPRB /) |
---|
| 179 | ZASWCE2 = (/ 0.829_JPRB , 0.7819_JPRB /) |
---|
| 180 | ZASWCF2 = (/ 2.482_JPRB , 5.0150_JPRB /) |
---|
| 181 | |
---|
| 182 | !* Ice cloud properties - from Sun and Shine (1995) |
---|
| 183 | |
---|
| 184 | ZSUSHE2 = (/ _ZERO_ , 8.6822_JPRB /) |
---|
| 185 | ZSUSHF2 = (/ _ZERO_ , 9.6277_JPRB /) |
---|
| 186 | ZSUSHH2 = (/ 0.8522_JPRB , 0.8819_JPRB /) |
---|
| 187 | ZSUSHK2 = (/ 0.1620_JPRB , 0.1630_JPRB /) |
---|
| 188 | ZSUSHA2 = (/ _ZERO_ , 23.204_JPRB /) |
---|
| 189 | ZSUSHG2 = (/ 0.3270_JPRB , 0.4180_JPRB /) |
---|
| 190 | |
---|
| 191 | !* Ice cloud properties - from Fu and Liou (1993) |
---|
| 192 | |
---|
| 193 | ZFLAA02 = (/-6.656E-3_JPRB ,-6.656E-3_JPRB /) |
---|
| 194 | ZFLAA12 = (/ 3.686_JPRB , 3.686_JPRB /) |
---|
| 195 | ZFLBB02 = (/ .10998E-5_JPRB , .21136E-1_JPRB /) |
---|
| 196 | ZFLBB12 = (/-.26101E-7_JPRB , .39150E-3_JPRB /) |
---|
| 197 | ZFLBB22 = (/ .18096E-8_JPRB ,-.20740E-6_JPRB /) |
---|
| 198 | ZFLBB32 = (/-.47387E-11_JPRB,-.28829E-8_JPRB /) |
---|
| 199 | ZFLCC02 = (/ .22110E+1_JPRB , .22498E+1_JPRB /) |
---|
| 200 | ZFLCC12 = (/-.10398E-2_JPRB , .23656E-3_JPRB /) |
---|
| 201 | ZFLCC22 = (/ .65199E-4_JPRB , .51948E-4_JPRB /) |
---|
| 202 | ZFLCC32 = (/-.34498E-6_JPRB ,-.29768E-6_JPRB /) |
---|
| 203 | ZFLDD02 = (/ .12495_JPRB , .11716_JPRB /) |
---|
| 204 | ZFLDD12 = (/-.43582E-3_JPRB ,-.45208E-3_JPRB /) |
---|
| 205 | ZFLDD22 = (/ .14092E-4_JPRB , .12772E-4_JPRB /) |
---|
| 206 | ZFLDD32 = (/-.69565E-7_JPRB ,-.62779E-7_JPRB /) |
---|
| 207 | |
---|
| 208 | !* Ice cloud properties - from Fu (1996) |
---|
| 209 | |
---|
| 210 | !ZFUAA02 = (/-.291721E-04_JPRB , |
---|
| 211 | !ZFUAA12 = (/ .251925E+01_JPRB , |
---|
| 212 | !ZFUBB02 = (/ .135403E-06_JPRB , |
---|
| 213 | !ZFUBB12 = (/ .992822E-07_JPRB , |
---|
| 214 | !ZFUBB22 = (/-.738432E-10_JPRB , |
---|
| 215 | !ZFUBB32 = (/ .331119E-12_JPRB , |
---|
| 216 | !ZFUCC02 = (/ .748127E+00_JPRB , |
---|
| 217 | !ZFUCC12 = (/ .956845E-03_JPRB , |
---|
| 218 | !ZFUCC22 = (/-.111517E-05_JPRB , |
---|
| 219 | !ZFUCC32 = (/-.815573E-08_JPRB , |
---|
| 220 | !ZFUDD02 = (/ .115730E+00_JPRB , |
---|
| 221 | !ZFUDD12 = (/ .256481E-03_JPRB , |
---|
| 222 | !ZFUDD22 = (/ .191313E-05_JPRB , |
---|
| 223 | !ZFUDD32 = (/-.124603E-07_JPRB , |
---|
| 224 | ! ------------------------------------------------------------------ |
---|
| 225 | |
---|
| 226 | !* 1.2 FOUR SPECTRAL INTERVALS |
---|
| 227 | ! ----------------------- |
---|
| 228 | |
---|
| 229 | ! SW : 4 spectral intervals (0.25 - 0.69 - 1.19 - 2.38 - 4.00) |
---|
| 230 | |
---|
| 231 | !* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992 |
---|
| 232 | |
---|
| 233 | ! ZEBCUA4 = (/ 3.448E-03 , 3.448E-03 , 3.448E-03 , 3.448E-03 /) |
---|
| 234 | ! ZEBCUB4 = (/ 2.431 , 2.431 , 2.431 , 2.431 /) |
---|
| 235 | ! ZEBCUC4 = (/ 0.00001 , 0.00011 , 0.01861 , 0.46658 /) |
---|
| 236 | ! ZEBCUD4 = (/ 0. , 1.405E-05 , 8.328E-4 , 2.050E-05 /) |
---|
| 237 | ! ZEBCUE4 = (/ 0.7661 , 0.7730 , 0.7940 , 0.9595 /) |
---|
| 238 | ! ZEBCUF4 = (/ 5.851E-04 , 5.665E-04 , 7.267E-04 , 1.076E-04 /) |
---|
| 239 | |
---|
| 240 | ZEBCUA4 = (/ 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB /) |
---|
| 241 | ZEBCUB4 = (/ 2.431_JPRB , 2.431_JPRB , 2.431_JPRB , 2.431_JPRB /) |
---|
| 242 | ZEBCUC4 = (/ 0.00001_JPRB , 0.00011_JPRB , 0.0197796_JPRB , 0.46658_JPRB /) |
---|
| 243 | ZEBCUD4 = (/ _ZERO_ , 1.405E-05_JPRB , 7.95513E-4_JPRB, 2.050E-05_JPRB /) |
---|
| 244 | ZEBCUE4 = (/ 0.7661_JPRB , 0.7730_JPRB , 0.795653_JPRB , 0.9595_JPRB /) |
---|
| 245 | ZEBCUF4 = (/ 5.851E-04_JPRB , 5.665E-04_JPRB , 7.267E-04_JPRB , 1.076E-04_JPRB /) |
---|
| 246 | |
---|
| 247 | !* Water cloud properties - from Fouquart (1987) |
---|
| 248 | |
---|
| 249 | ZYFWCA4 = (/ _ZERO_ , _ZERO_ , _ZERO_ , _ZERO_ /) |
---|
| 250 | ZYFWCB4 = (/ 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB /) |
---|
| 251 | ZYFWCC4 = (/ 0.9999_JPRB , 0.9988_JPRB , 0.9988_JPRB , 0.9988_JPRB /) |
---|
| 252 | ZYFWCD4 = (/ 5.000E-04_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB, 2.500E-03_JPRB /) |
---|
| 253 | ZYFWCE4 = (/ _HALF_ , 0.05_JPRB , 0.05_JPRB , 0.05_JPRB /) |
---|
| 254 | ZYFWCF4 = (/ 0.865_JPRB , 0.910_JPRB , 0.910_JPRB , 0.910_JPRB /) |
---|
| 255 | |
---|
| 256 | !* Water cloud properties - from Slingo (1989) |
---|
| 257 | |
---|
| 258 | ZASWCA4 = (/ 2.817_JPRB , 2.682_JPRB , 2.264_JPRB , 1.281_JPRB /) |
---|
| 259 | ZASWCB4 = (/ 1.305_JPRB , 1.346_JPRB , 1.455_JPRB , 1.641_JPRB /) |
---|
| 260 | ZASWCC4 = (/-5.62E-08_JPRB ,-6.94E-06_JPRB , 4.64E-04_JPRB , 2.01E-01_JPRB /) |
---|
| 261 | ZASWCD4 = (/ 1.63E-07_JPRB , 2.35E-05_JPRB , 1.24E-03_JPRB , 7.56E-03_JPRB /) |
---|
| 262 | ZASWCE4 = (/ 0.829_JPRB , 0.794_JPRB , 0.754_JPRB , 0.826_JPRB /) |
---|
| 263 | ZASWCF4 = (/ 2.482_JPRB , 4.226_JPRB , 6.560_JPRB , 4.353_JPRB /) |
---|
| 264 | |
---|
| 265 | !* Ice cloud properties - from Sun and Shine (1995) |
---|
| 266 | |
---|
| 267 | ZSUSHE4 = (/ _ZERO_ , 7.2471E-02_JPRB , 17.5933_JPRB , 48.7166_JPRB /) |
---|
| 268 | ZSUSHF4 = (/ _ZERO_ , 4.01511_JPRB , 21.1249_JPRB , 1.20890_JPRB /) |
---|
| 269 | ZSUSHH4 = (/ 0.8522_JPRB , 0.85841_JPRB , 0.90778_JPRB ,0.982046_JPRB /) |
---|
| 270 | ZSUSHK4 = (/ 0.1620_JPRB , 0.160048_JPRB , 0.188521_JPRB ,0.0411446_JPRB /) |
---|
| 271 | ZSUSHA4 = (/ _ZERO_ , 0.273455_JPRB , 41.7675_JPRB , 161.104_JPRB /) |
---|
| 272 | ZSUSHG4 = (/ 0.3270_JPRB , 0.343668_JPRB , 0.526192_JPRB , 0.574040_JPRB /) |
---|
| 273 | |
---|
| 274 | !* Ice cloud properties - from Fu and Liou (1993) |
---|
| 275 | |
---|
| 276 | ZFLAA04 = (/-6.656E-03_JPRB ,-6.656E-03_JPRB,-6.656E-03_JPRB,-6.656E-03_JPRB /) |
---|
| 277 | ZFLAA14 = (/ 3.686_JPRB , 3.686_JPRB , 3.686_JPRB , 3.686_JPRB /) |
---|
| 278 | ZFLBB04 = (/ .10998E-5_JPRB , .20208E-4_JPRB, .51557E-3_JPRB, .39517E+0_JPRB /) |
---|
| 279 | ZFLBB14 = (/-.26101E-7_JPRB , .96483E-5_JPRB, .10731E-2_JPRB, .15787E-2_JPRB /) |
---|
| 280 | ZFLBB24 = (/ .18096E-8_JPRB , .83009E-7_JPRB, .17753E-5_JPRB,-.14337E-4_JPRB /) |
---|
| 281 | ZFLBB34 = (/-.47387E-11_JPRB,-.32217E-9_JPRB,-.18379E-7_JPRB, .46942E-7_JPRB /) |
---|
| 282 | ZFLCC04 = (/ .22110E+1_JPRB , .22151E+1_JPRB, .22534E+1_JPRB, .26653E+1_JPRB /) |
---|
| 283 | ZFLCC14 = (/-.10398E-2_JPRB ,-.77982E-3_JPRB, .16163E-2_JPRB, .56935E-2_JPRB /) |
---|
| 284 | ZFLCC24 = (/ .65199E-4_JPRB , .63750E-4_JPRB, .44037E-4_JPRB,-.54316E-4_JPRB /) |
---|
| 285 | ZFLCC34 = (/-.34498E-6_JPRB ,-.34466E-6_JPRB,-.27627E-6_JPRB, .17858E-6_JPRB /) |
---|
| 286 | ZFLDD04 = (/ .12495_JPRB , .12363_JPRB , .11983_JPRB , .21834E-1_JPRB /) |
---|
| 287 | ZFLDD14 = (/-.43582E-3_JPRB ,-.44419E-3_JPRB,-.50108E-3_JPRB,-.29204E-3_JPRB /) |
---|
| 288 | ZFLDD24 = (/ .14092E-4_JPRB , .14038E-4_JPRB, .11843E-4_JPRB, .18060E-5_JPRB /) |
---|
| 289 | ZFLDD34 = (/-.69565E-7_JPRB ,-.68851E-7_JPRB,-.59367E-7_JPRB,-.46257E-8_JPRB /) |
---|
| 290 | |
---|
| 291 | !* Ice cloud properties - from Fu (1996) as tabulated in Sun & Rikus (1999) |
---|
| 292 | |
---|
| 293 | ZFUAA04 = (/-1.30817E-04_JPRB,-6.39479E-05_JPRB,-6.74730E-06_JPRB, 1.62674E-04_JPRB /) |
---|
| 294 | ZFUAA14 = (/ 2.52883E+00_JPRB, 2.52393E+00_JPRB, 2.52056E+00_JPRB, 2.49823E+00_JPRB /) |
---|
| 295 | ZFUBB04 = (/-1.55357E-08_JPRB,-7.90657E-07_JPRB, 9.64842E-04_JPRB, 2.25112E-01_JPRB /) |
---|
| 296 | ZFUBB14 = (/ 1.95793E-07_JPRB, 7.79991E-06_JPRB, 9.09809E-04_JPRB, 3.05017E-03_JPRB /) |
---|
| 297 | ZFUBB24 = (/-2.31234E-10_JPRB, 2.90894E-10_JPRB,-3.57557E-06_JPRB,-2.54236E-05_JPRB /) |
---|
| 298 | ZFUBB34 = (/ 1.12247E-12_JPRB,-2.02818E-12_JPRB, 1.00197E-08_JPRB, 8.49116E-08_JPRB /) |
---|
| 299 | ZFUCC04 = (/ 7.39781E-01_JPRB, 7.52335E-01_JPRB, 7.56307E-01_JPRB, 8.30812E-01_JPRB /) |
---|
| 300 | ZFUCC14 = (/ 9.10564E-04_JPRB, 1.06211E-03_JPRB, 1.73364E-03_JPRB, 2.62788E-03_JPRB /) |
---|
| 301 | ZFUCC24 = (/-4.62479E-07_JPRB,-2.45770E-06_JPRB,-8.92191E-06_JPRB,-2.43196E-05_JPRB /) |
---|
| 302 | ZFUCC34 = (/-1.05910E-08_JPRB,-3.03712E-09_JPRB, 1.97757E-08_JPRB, 8.23543E-08_JPRB /) |
---|
| 303 | |
---|
| 304 | ! ------------------------------------------------------------------ |
---|
| 305 | |
---|
| 306 | !* 1.3 SIX SPECTRAL INTERVALS |
---|
| 307 | ! ---------------------- |
---|
| 308 | |
---|
| 309 | ! SW : 6 spectral intervals (0.185-0.25-0.44-0.69-1.19-2.38-4.00) |
---|
| 310 | |
---|
| 311 | !* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992 |
---|
| 312 | |
---|
| 313 | ZEBCUA6 = (/ 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB /) |
---|
| 314 | ZEBCUB6 = (/ 2.431_JPRB , 2.431_JPRB , 2.431_JPRB , 2.431_JPRB , 2.431_JPRB , 2.431_JPRB /) |
---|
| 315 | ZEBCUC6 = (/ 0.00001_JPRB , 0.00001_JPRB , 0.00001_JPRB , 0.00011_JPRB , 0.0197796_JPRB , 0.46658_JPRB /) |
---|
| 316 | ZEBCUD6 = (/ _ZERO_ , _ZERO_ , _ZERO_ , 1.405E-05_JPRB , 7.95513E-4_JPRB, 2.050E-05_JPRB /) |
---|
| 317 | ZEBCUE6 = (/ 0.7661_JPRB , 0.7661_JPRB , 0.7661_JPRB , 0.7730_JPRB , 0.795653_JPRB , 0.9595_JPRB /) |
---|
| 318 | ZEBCUF6 = (/ 5.851E-04_JPRB , 5.851E-04_JPRB , 5.851E-04_JPRB , 5.665E-04_JPRB , 7.267E-04_JPRB , 1.076E-04_JPRB /) |
---|
| 319 | |
---|
| 320 | !* Water cloud properties - from Fouquart (1987) |
---|
| 321 | |
---|
| 322 | ZYFWCA6 = (/ _ZERO_ , _ZERO_ , _ZERO_ , _ZERO_ , _ZERO_ , _ZERO_ /) |
---|
| 323 | ZYFWCB6 = (/ 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB /) |
---|
| 324 | ZYFWCC6 = (/ 0.9999_JPRB , 0.9999_JPRB , 0.9999_JPRB , 0.9988_JPRB , 0.9988_JPRB , 0.9988_JPRB /) |
---|
| 325 | ZYFWCD6 = (/ 5.000E-04_JPRB , 5.000E-04_JPRB , 5.000E-04_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB, 2.500E-03_JPRB /) |
---|
| 326 | ZYFWCE6 = (/ _HALF_ , _HALF_ , _HALF_ , 0.05_JPRB , 0.05_JPRB , 0.05_JPRB /) |
---|
| 327 | ZYFWCF6 = (/ 0.865_JPRB , 0.865_JPRB , 0.865_JPRB , 0.910_JPRB , 0.910_JPRB , 0.910_JPRB /) |
---|
| 328 | |
---|
| 329 | !* Water cloud properties - from Slingo (1989) |
---|
| 330 | |
---|
| 331 | ZASWCA6 = (/ 2.817_JPRB , 2.817_JPRB , 2.817_JPRB , 2.682_JPRB , 2.264_JPRB , 1.281_JPRB /) |
---|
| 332 | ZASWCB6 = (/ 1.305_JPRB , 1.305_JPRB , 1.305_JPRB , 1.346_JPRB , 1.455_JPRB , 1.641_JPRB /) |
---|
| 333 | ZASWCC6 = (/-5.62E-08_JPRB ,-5.62E-08_JPRB ,-5.62E-08_JPRB ,-6.94E-06_JPRB , 4.64E-04_JPRB , 2.01E-01_JPRB /) |
---|
| 334 | ZASWCD6 = (/ 1.63E-07_JPRB , 1.63E-07_JPRB , 1.63E-07_JPRB , 2.35E-05_JPRB , 1.24E-03_JPRB , 7.56E-03_JPRB /) |
---|
| 335 | ZASWCE6 = (/ 0.829_JPRB , 0.829_JPRB , 0.829_JPRB , 0.794_JPRB , 0.754_JPRB , 0.826_JPRB /) |
---|
| 336 | ZASWCF6 = (/ 2.482_JPRB , 2.482_JPRB , 2.482_JPRB , 4.226_JPRB , 6.560_JPRB , 4.353_JPRB /) |
---|
| 337 | |
---|
| 338 | !* Ice cloud properties - from Sun and Shine (1995) |
---|
| 339 | |
---|
| 340 | ZSUSHE6 = (/ _ZERO_ , _ZERO_ , _ZERO_ , 7.2471E-02_JPRB , 17.5933_JPRB , 48.7166_JPRB /) |
---|
| 341 | ZSUSHF6 = (/ _ZERO_ , _ZERO_ , _ZERO_ , 4.01511_JPRB , 21.1249_JPRB , 1.20890_JPRB /) |
---|
| 342 | ZSUSHH6 = (/ 0.8522_JPRB , 0.8522_JPRB , 0.8522_JPRB , 0.85841_JPRB , 0.90778_JPRB ,0.982046_JPRB /) |
---|
| 343 | ZSUSHK6 = (/ 0.1620_JPRB , 0.1620_JPRB , 0.1620_JPRB , 0.160048_JPRB , 0.188521_JPRB ,0.0411446_JPRB /) |
---|
| 344 | ZSUSHA6 = (/ _ZERO_ , _ZERO_ , _ZERO_ , 0.273455_JPRB , 41.7675_JPRB , 161.104_JPRB /) |
---|
| 345 | ZSUSHG6 = (/ 0.3270_JPRB , 0.3270_JPRB , 0.3270_JPRB , 0.343668_JPRB , 0.526192_JPRB , 0.574040_JPRB /) |
---|
| 346 | |
---|
| 347 | !* Ice cloud properties - from Fu and Liou (1993) |
---|
| 348 | |
---|
| 349 | ZFLAA06 = (/-6.656E-03_JPRB ,-6.656E-03_JPRB ,-6.656E-03_JPRB ,-6.656E-03_JPRB,-6.656E-03_JPRB,-6.656E-03_JPRB /) |
---|
| 350 | ZFLAA16 = (/ 3.686_JPRB , 3.686_JPRB , 3.686_JPRB , 3.686_JPRB , 3.686_JPRB , 3.686_JPRB /) |
---|
| 351 | ZFLBB06 = (/ .10998E-5_JPRB , .10998E-5_JPRB , .10998E-5_JPRB , .20208E-4_JPRB, .51557E-3_JPRB, .39517E+0_JPRB /) |
---|
| 352 | ZFLBB16 = (/-.26101E-7_JPRB ,-.26101E-7_JPRB ,-.26101E-7_JPRB , .96483E-5_JPRB, .10731E-2_JPRB, .15787E-2_JPRB /) |
---|
| 353 | ZFLBB26 = (/ .18096E-8_JPRB , .18096E-8_JPRB , .18096E-8_JPRB , .83009E-7_JPRB, .17753E-5_JPRB,-.14337E-4_JPRB /) |
---|
| 354 | ZFLBB36 = (/-.47387E-11_JPRB,-.47387E-11_JPRB,-.47387E-11_JPRB,-.32217E-9_JPRB,-.18379E-7_JPRB, .46942E-7_JPRB /) |
---|
| 355 | ZFLCC06 = (/ .22110E+1_JPRB , .22110E+1_JPRB , .22110E+1_JPRB , .22151E+1_JPRB, .22534E+1_JPRB, .26653E+1_JPRB /) |
---|
| 356 | ZFLCC16 = (/-.10398E-2_JPRB ,-.10398E-2_JPRB ,-.10398E-2_JPRB ,-.77982E-3_JPRB, .16163E-2_JPRB, .56935E-2_JPRB /) |
---|
| 357 | ZFLCC26 = (/ .65199E-4_JPRB , .65199E-4_JPRB , .65199E-4_JPRB , .63750E-4_JPRB, .44037E-4_JPRB,-.54316E-4_JPRB /) |
---|
| 358 | ZFLCC36 = (/-.34498E-6_JPRB ,-.34498E-6_JPRB ,-.34498E-6_JPRB ,-.34466E-6_JPRB,-.27627E-6_JPRB, .17858E-6_JPRB /) |
---|
| 359 | ZFLDD06 = (/ .12495_JPRB , .12495_JPRB , .12495_JPRB , .12363_JPRB , .11983_JPRB , .21834E-1_JPRB /) |
---|
| 360 | ZFLDD16 = (/-.43582E-3_JPRB ,-.43582E-3_JPRB ,-.43582E-3_JPRB ,-.44419E-3_JPRB,-.50108E-3_JPRB,-.29204E-3_JPRB /) |
---|
| 361 | ZFLDD26 = (/ .14092E-4_JPRB , .14092E-4_JPRB , .14092E-4_JPRB , .14038E-4_JPRB, .11843E-4_JPRB, .18060E-5_JPRB /) |
---|
| 362 | ZFLDD36 = (/-.69565E-7_JPRB ,-.69565E-7_JPRB ,-.69565E-7_JPRB ,-.68851E-7_JPRB,-.59367E-7_JPRB,-.46257E-8_JPRB /) |
---|
| 363 | |
---|
| 364 | !* Ice cloud properties - from Fu (1996) as tabulated in Sun & Rikus (1999) |
---|
| 365 | |
---|
| 366 | ZFUAA06 = (/-1.30817E-04_JPRB,-1.30817E-04_JPRB,-1.30817E-04_JPRB,-6.39479E-05_JPRB,-6.74730E-06_JPRB, 1.62674E-04_JPRB /) |
---|
| 367 | ZFUAA16 = (/ 2.52883E+00_JPRB, 2.52883E+00_JPRB, 2.52883E+00_JPRB, 2.52393E+00_JPRB, 2.52056E+00_JPRB, 2.49823E+00_JPRB /) |
---|
| 368 | ZFUBB06 = (/-1.55357E-08_JPRB,-1.55357E-08_JPRB,-1.55357E-08_JPRB,-7.90657E-07_JPRB, 9.64842E-04_JPRB, 2.25112E-01_JPRB /) |
---|
| 369 | ZFUBB16 = (/ 1.95793E-07_JPRB, 1.95793E-07_JPRB, 1.95793E-07_JPRB, 7.79991E-06_JPRB, 9.09809E-04_JPRB, 3.05017E-03_JPRB /) |
---|
| 370 | ZFUBB26 = (/-2.31234E-10_JPRB,-2.31234E-10_JPRB,-2.31234E-10_JPRB, 2.90894E-10_JPRB,-3.57557E-06_JPRB,-2.54236E-05_JPRB /) |
---|
| 371 | ZFUBB36 = (/ 1.12247E-12_JPRB, 1.12247E-12_JPRB, 1.12247E-12_JPRB,-2.02818E-12_JPRB, 1.00197E-08_JPRB, 8.49116E-08_JPRB /) |
---|
| 372 | ZFUCC06 = (/ 7.39781E-01_JPRB, 7.39781E-01_JPRB, 7.39781E-01_JPRB, 7.52335E-01_JPRB, 7.56307E-01_JPRB, 8.30812E-01_JPRB /) |
---|
| 373 | ZFUCC16 = (/ 9.10564E-04_JPRB, 9.10564E-04_JPRB, 9.10564E-04_JPRB, 1.06211E-03_JPRB, 1.73364E-03_JPRB, 2.62788E-03_JPRB /) |
---|
| 374 | ZFUCC26 = (/-4.62479E-07_JPRB,-4.62479E-07_JPRB,-4.62479E-07_JPRB,-2.45770E-06_JPRB,-8.92191E-06_JPRB,-2.43196E-05_JPRB /) |
---|
| 375 | ZFUCC36 = (/-1.05910E-08_JPRB,-1.05910E-08_JPRB,-1.05910E-08_JPRB,-3.03712E-09_JPRB, 1.97757E-08_JPRB, 8.23543E-08_JPRB /) |
---|
| 376 | |
---|
| 377 | ! ------------------------------------------------------------------ |
---|
| 378 | |
---|
| 379 | ! LW : absorption coefficient as a function of effective radius in RRTM |
---|
| 380 | |
---|
| 381 | ! water clouds from Savijarvi |
---|
| 382 | |
---|
| 383 | RHSAVI( 1, :) = (/& |
---|
| 384 | &0.1651082_JPRB , -0.003494839_JPRB , _ZERO_ /) |
---|
| 385 | RHSAVI( 2, :) = (/& |
---|
| 386 | &0.327820597_JPRB, -0.0219634383_JPRB , 0.000506783898_JPRB /) |
---|
| 387 | RHSAVI( 3, :) = (/& |
---|
| 388 | &0.504805453_JPRB, -0.0478602354_JPRB , 0.00141521102_JPRB /) |
---|
| 389 | RHSAVI( 4, :) = (/& |
---|
| 390 | &0.513169093_JPRB, -0.0508960145_JPRB , 0.00155498711_JPRB /) |
---|
| 391 | RHSAVI( 5, :) = (/& |
---|
| 392 | &0.448042082_JPRB, -0.0431857592_JPRB , 0.00130848978_JPRB /) |
---|
| 393 | RHSAVI( 6, :) = (/& |
---|
| 394 | &0.249547237_JPRB, -0.0185273835_JPRB , 0.00050361258_JPRB /) |
---|
| 395 | RHSAVI( 7, :) = (/& |
---|
| 396 | &0.135486796_JPRB, -0.00585852322_JPRB, 0.00011873119_JPRB /) |
---|
| 397 | RHSAVI( 8, :) = (/& |
---|
| 398 | &0.126668819_JPRB, -0.00497949082_JPRB, 0.0000927679172_JPRB/) |
---|
| 399 | RHSAVI( 9, :) = (/& |
---|
| 400 | &0.130938752_JPRB, -0.00531814674_JPRB, 0.0000984953029_JPRB/) |
---|
| 401 | RHSAVI(10, :) = (/& |
---|
| 402 | &0.147024519_JPRB, -0.00671655774_JPRB, 0.000130885091_JPRB /) |
---|
| 403 | RHSAVI(11, :) = (/& |
---|
| 404 | &0.241463914_JPRB, -0.0174486461_JPRB , 0.000457756556_JPRB /) |
---|
| 405 | RHSAVI(12, :) = (/& |
---|
| 406 | &0.07576579_JPRB , -0.001695588_JPRB , _ZERO_ /) |
---|
| 407 | RHSAVI(13, :) = (/& |
---|
| 408 | &0.1032178_JPRB , -0.00293412_JPRB , _ZERO_ /) |
---|
| 409 | RHSAVI(14, :) = (/& |
---|
| 410 | &0.07342832_JPRB , -0.001775135_JPRB , _ZERO_ /) |
---|
| 411 | RHSAVI(15, :) = (/& |
---|
| 412 | &0.04649514_JPRB , -0.0009165462_JPRB , _ZERO_ /) |
---|
| 413 | RHSAVI(16, :) = (/& |
---|
| 414 | &0.06893519_JPRB , -0.0001245402_JPRB , _ZERO_ /) |
---|
| 415 | |
---|
| 416 | ! water clouds from Lindner & Li (2000) |
---|
| 417 | |
---|
| 418 | RLINLI( 1, :) = (/& |
---|
| 419 | & 0.88116E-01_JPRB,-0.12857E-02_JPRB, 0.81658E+00_JPRB,-0.39428E+01_JPRB, 0.46652E+01_JPRB /) |
---|
| 420 | RLINLI( 2, :) = (/& |
---|
| 421 | & 0.41307E-03_JPRB,-0.59631E-04_JPRB, 0.24275E+01_JPRB,-0.90838E+01_JPRB, 0.96069E+01_JPRB /) |
---|
| 422 | RLINLI( 3, :) = (/& |
---|
| 423 | &-0.57709E-01_JPRB, 0.99071E-03_JPRB, 0.31118E+01_JPRB,-0.95540E+01_JPRB, 0.90189E+01_JPRB /) |
---|
| 424 | RLINLI( 4, :) = (/& |
---|
| 425 | &-0.53069E-01_JPRB, 0.99992E-03_JPRB, 0.28045E+01_JPRB,-0.72836E+01_JPRB, 0.62573E+01_JPRB /) |
---|
| 426 | RLINLI( 5, :) = (/& |
---|
| 427 | &-0.23627E-01_JPRB, 0.55291E-03_JPRB, 0.21785E+01_JPRB,-0.54664E+01_JPRB, 0.47379E+01_JPRB /) |
---|
| 428 | RLINLI( 6, :) = (/& |
---|
| 429 | & 0.29022E-01_JPRB,-0.39657E-03_JPRB, 0.14902E+01_JPRB,-0.50777E+01_JPRB, 0.52170E+01_JPRB /) |
---|
| 430 | RLINLI( 7, :) = (/& |
---|
| 431 | &-0.24901E-01_JPRB, 0.16195E-03_JPRB, 0.29375E+02_JPRB,-0.11437E+02_JPRB, 0.12273E+02_JPRB /) |
---|
| 432 | RLINLI( 8, :) = (/& |
---|
| 433 | &-0.14269E+00_JPRB, 0.22282E-02_JPRB, 0.46478E+01_JPRB,-0.16369E+02_JPRB, 0.16533E+02_JPRB /) |
---|
| 434 | RLINLI( 9, :) = (/& |
---|
| 435 | &-0.20398E+00_JPRB, 0.34708E-02_JPRB, 0.52858E+01_JPRB,-0.16603E+02_JPRB, 0.15392E+02_JPRB /) |
---|
| 436 | RLINLI(10, :) = (/& |
---|
| 437 | &-0.18318E+00_JPRB, 0.33080E-02_JPRB, 0.46120E+01_JPRB,-0.11550E+02_JPRB, 0.87086E+01_JPRB /) |
---|
| 438 | RLINLI(11, :) = (/& |
---|
| 439 | &-0.20420E+00_JPRB, 0.37167E-02_JPRB, 0.48566E+01_JPRB,-0.11972E+02_JPRB, 0.86344E+01_JPRB /) |
---|
| 440 | RLINLI(12, :) = (/& |
---|
| 441 | &-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /) |
---|
| 442 | RLINLI(13, :) = (/& |
---|
| 443 | &-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /) |
---|
| 444 | RLINLI(14, :) = (/& |
---|
| 445 | &-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /) |
---|
| 446 | RLINLI(15, :) = (/& |
---|
| 447 | &-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /) |
---|
| 448 | RLINLI(16, :) = (/& |
---|
| 449 | &-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /) |
---|
| 450 | |
---|
| 451 | ! ice clouds From Fu & Liou |
---|
| 452 | |
---|
| 453 | RFULIO( 1, :) = (/& |
---|
| 454 | & -7.752E-03_JPRB, 4.624_JPRB, -42.01_JPRB /) |
---|
| 455 | RFULIO( 2, :) = (/& |
---|
| 456 | & -1.741E-02_JPRB, 5.541_JPRB, -58.42_JPRB /) |
---|
| 457 | RFULIO( 3, :) = (/& |
---|
| 458 | & -1.704E-02_JPRB, 4.830_JPRB, 16.27_JPRB /) |
---|
| 459 | RFULIO( 4, :) = (/& |
---|
| 460 | & -1.151E-02_JPRB, 4.182_JPRB, 31.13_JPRB /) |
---|
| 461 | RFULIO( 5, :) = (/& |
---|
| 462 | & -1.026E-02_JPRB, 4.105_JPRB, 16.36_JPRB /) |
---|
| 463 | RFULIO( 6, :) = (/& |
---|
| 464 | & -8.294E-03_JPRB, 3.925_JPRB, 1.315_JPRB /) |
---|
| 465 | RFULIO( 7, :) = (/& |
---|
| 466 | & -1.153E-02_JPRB, 4.109_JPRB, 17.32_JPRB /) |
---|
| 467 | RFULIO( 8, :) = (/& |
---|
| 468 | & -9.609E-03_JPRB, 3.768_JPRB, 34.11_JPRB /) |
---|
| 469 | RFULIO( 9, :) = (/& |
---|
| 470 | & -9.061E-03_JPRB, 3.741_JPRB, 26.48_JPRB /) |
---|
| 471 | RFULIO(10, :) = (/& |
---|
| 472 | & -8.441E-03_JPRB, 3.715_JPRB, 19.48_JPRB /) |
---|
| 473 | RFULIO(11, :) = (/& |
---|
| 474 | & -8.088E-03_JPRB, 3.717_JPRB, 17.17_JPRB /) |
---|
| 475 | RFULIO(12, :) = (/& |
---|
| 476 | & -8.088E-03_JPRB, 3.717_JPRB, 17.17_JPRB /) |
---|
| 477 | RFULIO(13, :) = (/& |
---|
| 478 | & -7.770E-03_JPRB, 3.734_JPRB, 11.85_JPRB /) |
---|
| 479 | RFULIO(14, :) = (/& |
---|
| 480 | & -6.656E-03_JPRB, 3.686_JPRB, _ZERO_ /) |
---|
| 481 | RFULIO(15, :) = (/& |
---|
| 482 | & -6.656E-03_JPRB, 3.686_JPRB, _ZERO_ /) |
---|
| 483 | RFULIO(16, :) = (/& |
---|
| 484 | & -6.656E-03_JPRB, 3.686_JPRB, _ZERO_ /) |
---|
| 485 | |
---|
| 486 | ! ice clouds from Fu et al. (1998) |
---|
| 487 | |
---|
| 488 | RFUETA( 1, :) = (/& |
---|
| 489 | & 4.919685E-03_JPRB, 2.327741E+00_JPRB,-1.390858E+01_JPRB /) |
---|
| 490 | RFUETA( 2, :) = (/& |
---|
| 491 | & 3.325756E-03_JPRB, 2.601360E+00_JPRB,-1.909602E+01_JPRB /) |
---|
| 492 | RFUETA( 3, :) = (/& |
---|
| 493 | &-1.334860E-02_JPRB, 4.043808E+00_JPRB,-2.171029E+01_JPRB /) |
---|
| 494 | RFUETA( 4, :) = (/& |
---|
| 495 | &-9.524174E-03_JPRB, 3.587742E+00_JPRB,-1.068895E+01_JPRB /) |
---|
| 496 | RFUETA( 5, :) = (/& |
---|
| 497 | &-4.159424E-03_JPRB, 3.047325E+00_JPRB,-5.061568E+00_JPRB /) |
---|
| 498 | RFUETA( 6, :) = (/& |
---|
| 499 | &-1.691632E-03_JPRB, 2.765756E+00_JPRB,-8.331033E+00_JPRB /) |
---|
| 500 | RFUETA( 7, :) = (/& |
---|
| 501 | &-8.372696E-03_JPRB, 3.455018E+00_JPRB,-1.516692E+01_JPRB /) |
---|
| 502 | RFUETA( 8, :) = (/& |
---|
| 503 | &-8.178608E-03_JPRB, 3.401245E+00_JPRB,-8.812820E+00_JPRB /) |
---|
| 504 | RFUETA( 9, :) = (/& |
---|
| 505 | &-4.936610E-03_JPRB, 3.087764E+00_JPRB,-3.884262E+00_JPRB /) |
---|
| 506 | RFUETA(10, :) = (/& |
---|
| 507 | &-3.034573E-03_JPRB, 2.900043E+00_JPRB,-1.849911E+00_JPRB /) |
---|
| 508 | RFUETA(11, :) = (/& |
---|
| 509 | &-3.034573E-03_JPRB, 2.900043E+00_JPRB,-1.849911E+00_JPRB /) |
---|
| 510 | RFUETA(12, :) = (/& |
---|
| 511 | &-2.465236E-03_JPRB, 2.833187E+00_JPRB,-4.227573E-01_JPRB /) |
---|
| 512 | RFUETA(13, :) = (/& |
---|
| 513 | &-2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /) |
---|
| 514 | RFUETA(14, :) = (/& |
---|
| 515 | &-2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /) |
---|
| 516 | RFUETA(15, :) = (/& |
---|
| 517 | &-2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /) |
---|
| 518 | RFUETA(16, :) = (/& |
---|
| 519 | &-2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /) |
---|
| 520 | |
---|
| 521 | |
---|
| 522 | ! ---------------------------------------------------------------- |
---|
| 523 | ! Ebert-Curry |
---|
| 524 | |
---|
| 525 | ! LW : spectrally defined for EC-OPE |
---|
| 526 | |
---|
| 527 | REBCUI = (/ 1.136_JPRB, 1.338_JPRB, 1.166_JPRB, 1.166_JPRB, 1.118_JPRB, & |
---|
| 528 | &0.600_JPRB /) |
---|
| 529 | REBCUJ = (/ 0.0036_JPRB, 0.0003_JPRB, 0.0016_JPRB, 0.0016_JPRB, 0.0020_JPRB, & |
---|
| 530 | &0.0068_JPRB /) |
---|
| 531 | |
---|
| 532 | ! LW : spectrally defined for RRTM |
---|
| 533 | ! mass-absorption coefficients for vertical path: no diffusivity factor |
---|
| 534 | |
---|
| 535 | REBCUG = (/ 0.718_JPRB, 0.726_JPRB, 1.136_JPRB, 1.320_JPRB, 1.505_JPRB, & |
---|
| 536 | & 1.290_JPRB, 0.911_JPRB, 0.949_JPRB, 1.021_JPRB, 1.193_JPRB, & |
---|
| 537 | & 1.279_JPRB, 0.626_JPRB, 0.647_JPRB, 0.668_JPRB, 0.690_JPRB, & |
---|
| 538 | & 0.690_JPRB /) |
---|
| 539 | |
---|
| 540 | |
---|
| 541 | REBCUH = (/ 0.0069_JPRB, 0.0060_JPRB, 0.0024_JPRB, 0.0004_JPRB,-0.0016_JPRB, & |
---|
| 542 | & 0.0003_JPRB, 0.0043_JPRB, 0.0038_JPRB, 0.0030_JPRB, 0.0013_JPRB, & |
---|
| 543 | & 0.0005_JPRB, 0.0054_JPRB, 0.0052_JPRB, 0.0050_JPRB, 0.0048_JPRB, & |
---|
| 544 | & 0.0048_JPRB /) |
---|
| 545 | |
---|
| 546 | |
---|
| 547 | ! Sun-Shine |
---|
| 548 | |
---|
| 549 | RSUSHFA = (/ 1.047_JPRB, -0.913E-04_JPRB, 0.203E-03_JPRB, -0.106E-04_JPRB /) |
---|
| 550 | |
---|
| 551 | ! ------------------------------------------------------------------ |
---|
| 552 | |
---|
| 553 | !* 2. |
---|
| 554 | |
---|
| 555 | !* Liquid/Solid water transition |
---|
| 556 | |
---|
| 557 | RTIW= 263._JPRB |
---|
| 558 | RRIW= 20._JPRB |
---|
| 559 | |
---|
| 560 | ! Ice particle Effective Radius as a function of LWC |
---|
| 561 | |
---|
| 562 | REFFIA= 40._JPRB |
---|
| 563 | REFFIB= 0._JPRB |
---|
| 564 | |
---|
| 565 | ! Sun-Shine |
---|
| 566 | |
---|
| 567 | RSUSHC= 0.0306_JPRB |
---|
| 568 | RSUSHD= 0.2548_JPRB |
---|
| 569 | |
---|
| 570 | ! ------------------------------------------------------------------ |
---|
| 571 | |
---|
| 572 | ! SW : absorption coefficients |
---|
| 573 | |
---|
| 574 | IF (KSW == 2) THEN |
---|
| 575 | DO JNU=1,KSW |
---|
| 576 | RASWCA(JNU)=ZASWCA2(JNU)*1.E-02_JPRB |
---|
| 577 | RASWCB(JNU)=ZASWCB2(JNU) |
---|
| 578 | RASWCC(JNU)=ZASWCC2(JNU) |
---|
| 579 | RASWCD(JNU)=ZASWCD2(JNU) |
---|
| 580 | RASWCE(JNU)=ZASWCE2(JNU) |
---|
| 581 | RASWCF(JNU)=ZASWCF2(JNU)*1.E-03_JPRB |
---|
| 582 | |
---|
| 583 | REBCUA(JNU)=ZEBCUA2(JNU) |
---|
| 584 | REBCUB(JNU)=ZEBCUB2(JNU) |
---|
| 585 | REBCUC(JNU)=ZEBCUC2(JNU) |
---|
| 586 | REBCUD(JNU)=ZEBCUD2(JNU) |
---|
| 587 | REBCUE(JNU)=ZEBCUE2(JNU) |
---|
| 588 | REBCUF(JNU)=ZEBCUF2(JNU) |
---|
| 589 | |
---|
| 590 | RYFWCA(JNU)=ZYFWCA2(JNU) |
---|
| 591 | RYFWCB(JNU)=ZYFWCB2(JNU) |
---|
| 592 | RYFWCC(JNU)=ZYFWCC2(JNU) |
---|
| 593 | RYFWCD(JNU)=ZYFWCD2(JNU) |
---|
| 594 | RYFWCE(JNU)=ZYFWCE2(JNU) |
---|
| 595 | RYFWCF(JNU)=ZYFWCF2(JNU) |
---|
| 596 | |
---|
| 597 | RSUSHE(JNU)=ZSUSHE2(JNU)*1.E-02_JPRB |
---|
| 598 | RSUSHF(JNU)=ZSUSHF2(JNU)*1.E-02_JPRB |
---|
| 599 | RSUSHH(JNU)=ZSUSHH2(JNU) |
---|
| 600 | RSUSHK(JNU)=ZSUSHK2(JNU)*1.E-01_JPRB |
---|
| 601 | RSUSHA(JNU)=ZSUSHA2(JNU)*1.E-03_JPRB |
---|
| 602 | RSUSHG(JNU)=ZSUSHG2(JNU)*1.E-01_JPRB |
---|
| 603 | |
---|
| 604 | RFLAA0(JNU)=ZFLAA02(JNU) |
---|
| 605 | RFLAA1(JNU)=ZFLAA12(JNU) |
---|
| 606 | RFLBB0(JNU)=ZFLBB02(JNU) |
---|
| 607 | RFLBB1(JNU)=ZFLBB12(JNU) |
---|
| 608 | RFLBB2(JNU)=ZFLBB22(JNU) |
---|
| 609 | RFLBB3(JNU)=ZFLBB32(JNU) |
---|
| 610 | RFLCC0(JNU)=ZFLCC02(JNU) |
---|
| 611 | RFLCC1(JNU)=ZFLCC12(JNU) |
---|
| 612 | RFLCC2(JNU)=ZFLCC22(JNU) |
---|
| 613 | RFLCC3(JNU)=ZFLCC32(JNU) |
---|
| 614 | RFLDD0(JNU)=ZFLDD02(JNU) |
---|
| 615 | RFLDD1(JNU)=ZFLDD12(JNU) |
---|
| 616 | RFLDD2(JNU)=ZFLDD22(JNU) |
---|
| 617 | RFLDD3(JNU)=ZFLDD32(JNU) |
---|
| 618 | |
---|
| 619 | ! RFUAA0(JNU)=ZFUAA02(JNU) |
---|
| 620 | ! RFUAA1(JNU)=ZFUAA12(JNU) |
---|
| 621 | ! RFUBB0(JNU)=ZFUBB02(JNU) |
---|
| 622 | ! RFUBB1(JNU)=ZFUBB12(JNU) |
---|
| 623 | ! RFUBB2(JNU)=ZFUBB22(JNU) |
---|
| 624 | ! RFUBB3(JNU)=ZFUBB32(JNU) |
---|
| 625 | ! RFUCC0(JNU)=ZFUCC02(JNU) |
---|
| 626 | ! RFUCC1(JNU)=ZFUCC12(JNU) |
---|
| 627 | ! RFUCC2(JNU)=ZFUCC22(JNU) |
---|
| 628 | ! RFUCC3(JNU)=ZFUCC32(JNU) |
---|
| 629 | ! RFUDD0(JNU)=ZFUDD02(JNU) |
---|
| 630 | ! RFUDD1(JNU)=ZFUDD12(JNU) |
---|
| 631 | ! RFUDD2(JNU)=ZFUDD22(JNU) |
---|
| 632 | ! RFUDD3(JNU)=ZFUDD32(JNU) |
---|
| 633 | |
---|
| 634 | ENDDO |
---|
| 635 | ELSEIF (KSW == 4) THEN |
---|
| 636 | DO JNU=1,KSW |
---|
| 637 | RASWCA(JNU)=ZASWCA4(JNU)*1.E-02_JPRB |
---|
| 638 | RASWCB(JNU)=ZASWCB4(JNU) |
---|
| 639 | RASWCC(JNU)=ZASWCC4(JNU) |
---|
| 640 | RASWCD(JNU)=ZASWCD4(JNU) |
---|
| 641 | RASWCE(JNU)=ZASWCE4(JNU) |
---|
| 642 | RASWCF(JNU)=ZASWCF4(JNU)*1.E-03_JPRB |
---|
| 643 | |
---|
| 644 | REBCUA(JNU)=ZEBCUA4(JNU) |
---|
| 645 | REBCUB(JNU)=ZEBCUB4(JNU) |
---|
| 646 | REBCUC(JNU)=ZEBCUC4(JNU) |
---|
| 647 | REBCUD(JNU)=ZEBCUD4(JNU) |
---|
| 648 | REBCUE(JNU)=ZEBCUE4(JNU) |
---|
| 649 | REBCUF(JNU)=ZEBCUF4(JNU) |
---|
| 650 | |
---|
| 651 | RYFWCA(JNU)=ZYFWCA4(JNU) |
---|
| 652 | RYFWCB(JNU)=ZYFWCB4(JNU) |
---|
| 653 | RYFWCC(JNU)=ZYFWCC4(JNU) |
---|
| 654 | RYFWCD(JNU)=ZYFWCD4(JNU) |
---|
| 655 | RYFWCE(JNU)=ZYFWCE4(JNU) |
---|
| 656 | RYFWCF(JNU)=ZYFWCF4(JNU) |
---|
| 657 | |
---|
| 658 | RSUSHE(JNU)=ZSUSHE4(JNU)*1.E-02_JPRB |
---|
| 659 | RSUSHF(JNU)=ZSUSHF4(JNU)*1.E-02_JPRB |
---|
| 660 | RSUSHH(JNU)=ZSUSHH4(JNU) |
---|
| 661 | RSUSHK(JNU)=ZSUSHK4(JNU)*1.E-01_JPRB |
---|
| 662 | RSUSHA(JNU)=ZSUSHA4(JNU)*1.E-03_JPRB |
---|
| 663 | RSUSHG(JNU)=ZSUSHG4(JNU)*1.E-01_JPRB |
---|
| 664 | |
---|
| 665 | RFLAA0(JNU)=ZFLAA04(JNU) |
---|
| 666 | RFLAA1(JNU)=ZFLAA14(JNU) |
---|
| 667 | RFLBB0(JNU)=ZFLBB04(JNU) |
---|
| 668 | RFLBB1(JNU)=ZFLBB14(JNU) |
---|
| 669 | RFLBB2(JNU)=ZFLBB24(JNU) |
---|
| 670 | RFLBB3(JNU)=ZFLBB34(JNU) |
---|
| 671 | RFLCC0(JNU)=ZFLCC04(JNU) |
---|
| 672 | RFLCC1(JNU)=ZFLCC14(JNU) |
---|
| 673 | RFLCC2(JNU)=ZFLCC24(JNU) |
---|
| 674 | RFLCC3(JNU)=ZFLCC34(JNU) |
---|
| 675 | RFLDD0(JNU)=ZFLDD04(JNU) |
---|
| 676 | RFLDD1(JNU)=ZFLDD14(JNU) |
---|
| 677 | RFLDD2(JNU)=ZFLDD24(JNU) |
---|
| 678 | RFLDD3(JNU)=ZFLDD34(JNU) |
---|
| 679 | |
---|
| 680 | RFUAA0(JNU)=ZFUAA04(JNU) |
---|
| 681 | RFUAA1(JNU)=ZFUAA14(JNU) |
---|
| 682 | RFUBB0(JNU)=ZFUBB04(JNU) |
---|
| 683 | RFUBB1(JNU)=ZFUBB14(JNU) |
---|
| 684 | RFUBB2(JNU)=ZFUBB24(JNU) |
---|
| 685 | RFUBB3(JNU)=ZFUBB34(JNU) |
---|
| 686 | RFUCC0(JNU)=ZFUCC04(JNU) |
---|
| 687 | RFUCC1(JNU)=ZFUCC14(JNU) |
---|
| 688 | RFUCC2(JNU)=ZFUCC24(JNU) |
---|
| 689 | RFUCC3(JNU)=ZFUCC34(JNU) |
---|
| 690 | ! RFUDD0(JNU)=ZFUDD04(JNU) |
---|
| 691 | ! RFUDD1(JNU)=ZFUDD14(JNU) |
---|
| 692 | ! RFUDD2(JNU)=ZFUDD24(JNU) |
---|
| 693 | ! RFUDD3(JNU)=ZFUDD34(JNU) |
---|
| 694 | |
---|
| 695 | ENDDO |
---|
| 696 | ELSEIF (KSW == 6) THEN |
---|
| 697 | DO JNU=1,KSW |
---|
| 698 | RASWCA(JNU)=ZASWCA6(JNU)*1.E-02_JPRB |
---|
| 699 | RASWCB(JNU)=ZASWCB6(JNU) |
---|
| 700 | RASWCC(JNU)=ZASWCC6(JNU) |
---|
| 701 | RASWCD(JNU)=ZASWCD6(JNU) |
---|
| 702 | RASWCE(JNU)=ZASWCE6(JNU) |
---|
| 703 | RASWCF(JNU)=ZASWCF6(JNU)*1.E-03_JPRB |
---|
| 704 | |
---|
| 705 | REBCUA(JNU)=ZEBCUA6(JNU) |
---|
| 706 | REBCUB(JNU)=ZEBCUB6(JNU) |
---|
| 707 | REBCUC(JNU)=ZEBCUC6(JNU) |
---|
| 708 | REBCUD(JNU)=ZEBCUD6(JNU) |
---|
| 709 | REBCUE(JNU)=ZEBCUE6(JNU) |
---|
| 710 | REBCUF(JNU)=ZEBCUF6(JNU) |
---|
| 711 | |
---|
| 712 | RYFWCA(JNU)=ZYFWCA6(JNU) |
---|
| 713 | RYFWCB(JNU)=ZYFWCB6(JNU) |
---|
| 714 | RYFWCC(JNU)=ZYFWCC6(JNU) |
---|
| 715 | RYFWCD(JNU)=ZYFWCD6(JNU) |
---|
| 716 | RYFWCE(JNU)=ZYFWCE6(JNU) |
---|
| 717 | RYFWCF(JNU)=ZYFWCF6(JNU) |
---|
| 718 | |
---|
| 719 | RSUSHE(JNU)=ZSUSHE6(JNU)*1.E-02_JPRB |
---|
| 720 | RSUSHF(JNU)=ZSUSHF6(JNU)*1.E-02_JPRB |
---|
| 721 | RSUSHH(JNU)=ZSUSHH6(JNU) |
---|
| 722 | RSUSHK(JNU)=ZSUSHK6(JNU)*1.E-01_JPRB |
---|
| 723 | RSUSHA(JNU)=ZSUSHA6(JNU)*1.E-03_JPRB |
---|
| 724 | RSUSHG(JNU)=ZSUSHG6(JNU)*1.E-01_JPRB |
---|
| 725 | |
---|
| 726 | RFLAA0(JNU)=ZFLAA06(JNU) |
---|
| 727 | RFLAA1(JNU)=ZFLAA16(JNU) |
---|
| 728 | RFLBB0(JNU)=ZFLBB06(JNU) |
---|
| 729 | RFLBB1(JNU)=ZFLBB16(JNU) |
---|
| 730 | RFLBB2(JNU)=ZFLBB26(JNU) |
---|
| 731 | RFLBB3(JNU)=ZFLBB36(JNU) |
---|
| 732 | RFLCC0(JNU)=ZFLCC06(JNU) |
---|
| 733 | RFLCC1(JNU)=ZFLCC16(JNU) |
---|
| 734 | RFLCC2(JNU)=ZFLCC26(JNU) |
---|
| 735 | RFLCC3(JNU)=ZFLCC36(JNU) |
---|
| 736 | RFLDD0(JNU)=ZFLDD06(JNU) |
---|
| 737 | RFLDD1(JNU)=ZFLDD16(JNU) |
---|
| 738 | RFLDD2(JNU)=ZFLDD26(JNU) |
---|
| 739 | RFLDD3(JNU)=ZFLDD36(JNU) |
---|
| 740 | |
---|
| 741 | RFUAA0(JNU)=ZFUAA06(JNU) |
---|
| 742 | RFUAA1(JNU)=ZFUAA16(JNU) |
---|
| 743 | RFUBB0(JNU)=ZFUBB06(JNU) |
---|
| 744 | RFUBB1(JNU)=ZFUBB16(JNU) |
---|
| 745 | RFUBB2(JNU)=ZFUBB26(JNU) |
---|
| 746 | RFUBB3(JNU)=ZFUBB36(JNU) |
---|
| 747 | RFUCC0(JNU)=ZFUCC06(JNU) |
---|
| 748 | RFUCC1(JNU)=ZFUCC16(JNU) |
---|
| 749 | RFUCC2(JNU)=ZFUCC26(JNU) |
---|
| 750 | RFUCC3(JNU)=ZFUCC36(JNU) |
---|
| 751 | ! RFUDD0(JNU)=ZFUDD06(JNU) |
---|
| 752 | ! RFUDD1(JNU)=ZFUDD16(JNU) |
---|
| 753 | ! RFUDD2(JNU)=ZFUDD26(JNU) |
---|
| 754 | ! RFUDD3(JNU)=ZFUDD36(JNU) |
---|
| 755 | |
---|
| 756 | ENDDO |
---|
| 757 | ELSE |
---|
| 758 | ! CALL ABOR1('SUCLOPN: WRONG SW SPECTRAL RESOLUTION') |
---|
| 759 | STOP 'SUCLOPN: WRONG SW SPECTRAL RESOLUTION' |
---|
| 760 | ENDIF |
---|
| 761 | |
---|
| 762 | ! ------------------------------------------------------------------ |
---|
| 763 | |
---|
| 764 | !* 2. CLOUD OVERLAP PARAMETERS |
---|
| 765 | ! ------------------------ |
---|
| 766 | |
---|
| 767 | ZAOVLP = (/ -2.250E-05_JPRB,-7.316E-06_JPRB,-1.966E-05_JPRB /) |
---|
| 768 | ZBOVLP = (/ 0.7865_JPRB , 0.8186_JPRB , 0.8900_JPRB /) |
---|
| 769 | |
---|
| 770 | IF (KLEV == 19) THEN |
---|
| 771 | RAOVLP=ZAOVLP(1) |
---|
| 772 | RBOVLP=ZBOVLP(1) |
---|
| 773 | ELSE IF (KLEV == 31) THEN |
---|
| 774 | RAOVLP=ZAOVLP(2) |
---|
| 775 | RBOVLP=ZBOVLP(2) |
---|
| 776 | ELSE IF (KLEV == 60) THEN |
---|
| 777 | RAOVLP=ZAOVLP(3) |
---|
| 778 | RBOVLP=ZBOVLP(3) |
---|
| 779 | ELSE |
---|
| 780 | RAOVLP=ZAOVLP(3) |
---|
| 781 | RBOVLP=ZBOVLP(3) |
---|
| 782 | END IF |
---|
| 783 | |
---|
| 784 | ! ------------------------------------------------------------------ |
---|
| 785 | |
---|
| 786 | RETURN |
---|
| 787 | END SUBROUTINE SUCLOPN |
---|