| [2089] | 1 | subroutine Atm_RT_INI |
|---|
| 2 | |
|---|
| 3 | !--------------------------------------------------------------------------+ |
|---|
| 4 | ! Mon 24-Jun-2013 MAR | |
|---|
| 5 | ! subroutine Atm_RT_INI initializes MAR PHYsics interface | |
|---|
| 6 | ! with ECMWF Solar/Infrared Radiative Transfer Scheme | |
|---|
| 7 | ! | |
|---|
| 8 | ! version 3.p.4.1 created by H. Gallee, Tue 2-Apr-2013 | |
|---|
| 9 | ! Last Modification by H. Gallee, Mon 24-Jun-2013 | |
|---|
| 10 | ! | |
|---|
| 11 | !--------------------------------------------------------------------------+ |
|---|
| 12 | ! | |
|---|
| 13 | ! Content: CALL of - ECMWF Code initializing the Radiation Transfert | |
|---|
| 14 | ! - ECMWF Radiation Transfert | |
|---|
| 15 | ! | |
|---|
| 16 | ! ECMWF Code Source: J.-J. Morcrette, 28 nov 2002 | |
|---|
| 17 | ! | |
|---|
| 18 | !--------------------------------------------------------------------------+ |
|---|
| 19 | |
|---|
| 20 | use Mod_Real |
|---|
| 21 | use Mod_PHY_RT_dat |
|---|
| 22 | use Mod_PHY____grd |
|---|
| 23 | use Mod_PHY_S0_grd |
|---|
| 24 | use Mod_PHY_RT_grd |
|---|
| 25 | use Mod_PHY_DY_kkl |
|---|
| 26 | use Mod_SISVAT_gpt |
|---|
| 27 | |
|---|
| 28 | |
|---|
| 29 | |
|---|
| 30 | #include "tsmbkind.h" |
|---|
| 31 | |
|---|
| 32 | |
|---|
| 33 | |
|---|
| 34 | |
|---|
| 35 | ! Global Variables (ECMWF) |
|---|
| 36 | ! ======================== |
|---|
| 37 | |
|---|
| 38 | USE PARRTM1D , ONLY : JP_LON ,JP_IDIA ,JP_FDIA ,JP_TDIA ,& |
|---|
| 39 | & JP_LEV ,JP_LW ,JP_SW ,JP_NUA ,JP_MODE ,& |
|---|
| 40 | & JP_AER ,JP_LEVP1 |
|---|
| 41 | USE YOMCST , ONLY : RD ,RG ,RTT ,RSIGMA ,& |
|---|
| 42 | & RCPD ,RPI ,RDAY ,REA ,RI0 ,& |
|---|
| 43 | & REPSM ,RMD ,RKBOL ,RNAVO ,R ,& |
|---|
| 44 | & RLVTT ,RLSTT |
|---|
| 45 | USE YOERAD , ONLY : NSW ,NTSW ,NRADFR ,& |
|---|
| 46 | & LRRTM ,LINHOM ,LOIFUEC ,LTEMPDS ,LOWASYF ,& |
|---|
| 47 | & LOWHSSS,LONEWSW ,LNEWAER ,LHVOLCA ,& |
|---|
| 48 | & NRADIP ,NRADLP ,NOZOCL ,& |
|---|
| 49 | & NICEOPT,NLIQOPT ,NOVLP ,NHOWINH ,RMINICE |
|---|
| 50 | USE YOEAERD , ONLY : CVDAES ,CVDAEL ,CVDAEU ,CVDAED ,& |
|---|
| 51 | & RCTRBGA,RCVOBGA ,RCSTBGA ,RCAEOPS ,RCAEOPL ,RCAEOPU ,& |
|---|
| 52 | & RCAEOPD,RCTRPT ,RCAEADK ,RCAEADM ,RCAEROS |
|---|
| 53 | USE YOERDI , ONLY : RCH4 ,RN2O ,RO3 ,RCFC11 ,& |
|---|
| 54 | & RCFC12 |
|---|
| 55 | USE YOERDU , ONLY : RCDAY ,R10E ,DIFF ,REPLOG ,& |
|---|
| 56 | & REPSC ,REPSCO ,REPSCQ ,REPSCT ,REPSCW ,& |
|---|
| 57 | & NTRAER |
|---|
| 58 | USE YOETHF , ONLY : R2ES ,R3LES ,R3IES ,R4LES ,& |
|---|
| 59 | & R4IES ,R5LES ,R5IES ,R5ALVCP ,R5ALSCP ,& |
|---|
| 60 | & RALVDCP,RALSDCP ,RTWAT ,RTICE ,RTICECU |
|---|
| 61 | |
|---|
| 62 | |
|---|
| 63 | |
|---|
| 64 | |
|---|
| 65 | IMPLICIT NONE |
|---|
| 66 | |
|---|
| 67 | |
|---|
| 68 | |
|---|
| 69 | |
|---|
| 70 | ! INTERNAL VARIABLES |
|---|
| 71 | ! ================== |
|---|
| 72 | |
|---|
| 73 | integer :: i ,j ,ikl |
|---|
| 74 | |
|---|
| 75 | ! For Use in radlsw |
|---|
| 76 | ! ^^^^^^^^^^^^^^^^^ |
|---|
| 77 | INTEGER_M :: KIDIA ,KFDIA ,KTDIA ,KLON ,KLEV |
|---|
| 78 | INTEGER_M :: KMODE ,KAER ,KSW |
|---|
| 79 | |
|---|
| 80 | INTEGER_M :: KBOX ,NBOX |
|---|
| 81 | INTEGER_M :: NDUMP ,ILWRAD |
|---|
| 82 | ! |
|---|
| 83 | REAL_B :: PRII05 |
|---|
| 84 | |
|---|
| 85 | REAL_B :: PAER5(JP_LON,JP_AER,JP_LEV) ! Aerosol Optical Depth |
|---|
| 86 | REAL_B :: PALBD5(JP_LON,JP_SW) |
|---|
| 87 | REAL_B :: PALBP5(JP_LON,JP_SW) |
|---|
| 88 | ! |
|---|
| 89 | REAL_B :: PAPH5(JP_LON,JP_LEVP1) |
|---|
| 90 | REAL_B :: PAP5(JP_LON,JP_LEV) |
|---|
| 91 | ! |
|---|
| 92 | REAL_B :: PCCO25 |
|---|
| 93 | REAL_B :: PCLFR5(JP_LON,JP_LEV) |
|---|
| 94 | REAL_B :: PDP5(JP_LON,JP_LEV) |
|---|
| 95 | REAL_B :: PEMIS5(JP_LON) |
|---|
| 96 | REAL_B :: PEMIW5(JP_LON) |
|---|
| 97 | REAL_B :: PLSM5(JP_LON) |
|---|
| 98 | REAL_B :: PMU05(JP_LON) |
|---|
| 99 | REAL_B :: POZON5(JP_LON,JP_LEV) |
|---|
| 100 | REAL_B :: PQ5(JP_LON,JP_LEV) |
|---|
| 101 | ! |
|---|
| 102 | REAL_B :: PQIWP5(JP_LON,JP_LEV) |
|---|
| 103 | REAL_B :: PQLWP5(JP_LON,JP_LEV) ! Dropplets Concentration |
|---|
| 104 | REAL_B :: PSQIW5(JP_LON,JP_LEV) ! Ice Crystals Concentration |
|---|
| 105 | REAL_B :: PSQLW5(JP_LON,JP_LEV) ! |
|---|
| 106 | REAL_B :: PQS5(JP_LON,JP_LEV) |
|---|
| 107 | REAL_B :: PQRAIN5(JP_LON,JP_LEV) |
|---|
| 108 | REAL_B :: PRAINT5(JP_LON,JP_LEV) |
|---|
| 109 | REAL_B :: PRLVRI5(JP_LON,JP_LEV) |
|---|
| 110 | REAL_B :: PRLVRL5(JP_LON,JP_LEV) |
|---|
| 111 | REAL_B :: PTH5(JP_LON,JP_LEVP1) |
|---|
| 112 | REAL_B :: PT5(JP_LON,JP_LEV) |
|---|
| 113 | REAL_B :: PTS5(JP_LON) |
|---|
| 114 | REAL_B :: PNBAS5(JP_LON) |
|---|
| 115 | REAL_B :: PNTOP5(JP_LON) |
|---|
| 116 | |
|---|
| 117 | ! For Use in SUCLD |
|---|
| 118 | ! ^^^^^^^^^^^^^^^^ |
|---|
| 119 | REAL_B :: ZETA(JP_LEV) |
|---|
| 120 | REAL_B :: ZETAH(JP_LEVP1) |
|---|
| 121 | |
|---|
| 122 | ! For Use in SUOVLP |
|---|
| 123 | ! ^^^^^^^^^^^^^^^^^ |
|---|
| 124 | REAL_B :: ZTVIR |
|---|
| 125 | REAL_B :: ZFACT |
|---|
| 126 | REAL_B :: ZAZ(JP_LEV) |
|---|
| 127 | REAL_B :: ZAZH(JP_LEVP1) |
|---|
| 128 | |
|---|
| 129 | |
|---|
| 130 | |
|---|
| 131 | ! Local Variables |
|---|
| 132 | ! ---------------- |
|---|
| 133 | |
|---|
| 134 | REAL_B :: RTIMTR |
|---|
| 135 | |
|---|
| 136 | INTEGER_M :: JL , JK , JAER ,JNU |
|---|
| 137 | INTEGER_M :: KULOUT , NINDAT , NSSSSS ,KPRTLEV |
|---|
| 138 | INTEGER_M :: IYR , MONTH , IDAY ,IMINUT |
|---|
| 139 | INTEGER_M :: KPRINT , SKSHIFT |
|---|
| 140 | |
|---|
| 141 | |
|---|
| 142 | |
|---|
| 143 | |
|---|
| 144 | ! Load External Functions |
|---|
| 145 | ! ======================= |
|---|
| 146 | |
|---|
| 147 | #include "fctast.h" |
|---|
| 148 | #include "fcttim.h" |
|---|
| 149 | #include "fcttre.h" |
|---|
| 150 | |
|---|
| 151 | |
|---|
| 152 | KPRTLEV= 1 |
|---|
| 153 | KPRINT = 1 |
|---|
| 154 | SKSHIFT= 0 |
|---|
| 155 | |
|---|
| 156 | |
|---|
| 157 | |
|---|
| 158 | |
|---|
| 159 | ! Time Base |
|---|
| 160 | ! ========= |
|---|
| 161 | |
|---|
| 162 | NINDAT = min(yearTU,2004)*10000+mon_TU*100+Day_TU ! Date in the form yyyyMMdd |
|---|
| 163 | NSSSSS = HourTU * 3600+minuTU* 60+sec_TU ! Nb of second since day Begin |
|---|
| 164 | IYR = NAA(NINDAT) |
|---|
| 165 | MONTH = NMM(NINDAT) |
|---|
| 166 | IDAY = NDD(NINDAT) |
|---|
| 167 | RTIMTR = RTIME(IYR,MONTH,IDAY,NSSSSS) |
|---|
| 168 | IMINUT = INT(FLOAT(NSSSSS)/60.) |
|---|
| 169 | |
|---|
| 170 | |
|---|
| 171 | |
|---|
| 172 | |
|---|
| 173 | ! Basic Initialization |
|---|
| 174 | ! ==================== |
|---|
| 175 | |
|---|
| 176 | ! Dimensions (auxiliary variables) |
|---|
| 177 | ! -------------------------------- |
|---|
| 178 | |
|---|
| 179 | KIDIA = 1 ! DO'NT CHANGE |
|---|
| 180 | KFDIA = JP_LON ! Nb Columns |
|---|
| 181 | KTDIA = 1 ! |
|---|
| 182 | KLON = JP_LON ! Nb Columns |
|---|
| 183 | KLEV = JP_LEV ! Nb Levels |
|---|
| 184 | KMODE = JP_MODE ! Used in Planck Fcts Specification |
|---|
| 185 | KAER = JP_AER ! |
|---|
| 186 | |
|---|
| 187 | ! Nb of Solar Spectral Intervals |
|---|
| 188 | ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ |
|---|
| 189 | KSW = JP_SW ! SW Nb of Spectral Intervals (max is JP_SW=6) |
|---|
| 190 | NSW = JP_SW ! SW Nb of Spectral Intervals (max is JP_SW=6) |
|---|
| 191 | NTSW = JP_SW ! SW Nb of Spectral Intervals (max is JP_SW=6) |
|---|
| 192 | |
|---|
| 193 | KBOX = 0 ! \VER |
|---|
| 194 | NBOX = 1 ! \VER |
|---|
| 195 | |
|---|
| 196 | ILWRAD = 1 ! 0: Morcrette, 1991 operational before 20000627 |
|---|
| 197 | ! 1: Mlawer et al., 1997 now ECMWF-operational |
|---|
| 198 | ! 2: Morcrette, 1991 original as in ERA-15' |
|---|
| 199 | |
|---|
| 200 | NDUMP = 3 ! No Print |
|---|
| 201 | ! NDUMP = 2 ! 1D Results |
|---|
| 202 | ! NDUMP = 1 ! Debug |
|---|
| 203 | ! NDUMP = 0 ! ALL |
|---|
| 204 | KULOUT = 6 ! Output Device for SUCST |
|---|
| 205 | |
|---|
| 206 | |
|---|
| 207 | ! Verification of the Dimensions |
|---|
| 208 | ! ------------------------------ |
|---|
| 209 | |
|---|
| 210 | write(6,*) 'INITIALISATION OF ECMWF RADIATIVE TRANSFERT: BEGIN' |
|---|
| 211 | |
|---|
| 212 | write(6,*) 'CONTROL' |
|---|
| 213 | write(6,*) 'JP_LON=',JP_LON |
|---|
| 214 | write(6,*) 'kcolp',kcolp |
|---|
| 215 | IF (kcolp .GT.JP_LON) THEN |
|---|
| 216 | write(6,6001) kcolp ,JP_LON |
|---|
| 217 | 6001 format(' @!#& BAD dimensions SET-UP (kcolp > JP_LON); = ('& |
|---|
| 218 | & ,i3,',',i3,')') |
|---|
| 219 | STOP |
|---|
| 220 | END IF |
|---|
| 221 | IF (mzp .NE.JP_LEV) THEN |
|---|
| 222 | write(6,6002) mzp ,JP_LEV |
|---|
| 223 | 6002 format(' @!#& BAD dimensions SET-UP (mzp /= JP_LEV); = ('& |
|---|
| 224 | & ,i3,',',i3,')') |
|---|
| 225 | STOP |
|---|
| 226 | END IF |
|---|
| 227 | IF (naero .NE.JP_AER) THEN |
|---|
| 228 | write(6,6003) naero ,JP_AER |
|---|
| 229 | 6003 format(' @!#& BAD dimensions SET-UP (naero/= JP_AER); = ('& |
|---|
| 230 | & ,i3,',',i3,')') |
|---|
| 231 | STOP |
|---|
| 232 | END IF |
|---|
| 233 | |
|---|
| 234 | |
|---|
| 235 | ! Mathematical Constants |
|---|
| 236 | ! ---------------------- |
|---|
| 237 | |
|---|
| 238 | REPLOG = 1.E-12 ! Minimum Logarithm Argument |
|---|
| 239 | REPSC = 1.E-12 |
|---|
| 240 | REPSCO = 1.E-12 |
|---|
| 241 | REPSCQ = 1.E-12 |
|---|
| 242 | REPSCT = 1.E-12 |
|---|
| 243 | REPSCW = 1.E-12 |
|---|
| 244 | |
|---|
| 245 | |
|---|
| 246 | ! Switches (general) |
|---|
| 247 | ! ------------------ |
|---|
| 248 | |
|---|
| 249 | LONEWSW= .TRUE. ! .TRUE. SWSN radiative routine is used |
|---|
| 250 | IF (ILWRAD.EQ.1) THEN |
|---|
| 251 | LRRTM = .TRUE. ! .TRUE. RRTM radiative routine is used |
|---|
| 252 | ELSE |
|---|
| 253 | LRRTM = .FALSE. ! .FALSE.RRTM radiative routine is NOT used |
|---|
| 254 | END IF |
|---|
| 255 | LTEMPDS= .FALSE. ! .TRUE. ALLOWS FOR SURF. T DISCONTIN. IN RAD.COMPUT. |
|---|
| 256 | |
|---|
| 257 | NTRAER = 19 ! NUMBER OF TRANSMISSION FUNCTIONS W OR W/O AEROSOLS |
|---|
| 258 | |
|---|
| 259 | |
|---|
| 260 | ! Switches (Clouds Optical Properties) |
|---|
| 261 | ! ------------------------------------ |
|---|
| 262 | |
|---|
| 263 | LINHOM = .FALSE. ! Tiedke (1995) correct. factor (0.7) of tau not used |
|---|
| 264 | NHOWINH= 2 ! Tau correction factor: exp(-(sig/tau)^2) |
|---|
| 265 | ! (used if LINHOM = .TRUE.) |
|---|
| 266 | LOIFUEC= .FALSE. ! .FALSE. IF ICE CLOUDS AS EBERT-CURRY (LW & SW) |
|---|
| 267 | LOWASYF= .FALSE. ! .FALSE. IF WATER CLOUDS AS FOUQUART (SW) |
|---|
| 268 | LOWHSSS= .FALSE. ! .FALSE. IF WATER CLOUDS AS SMITH-SHI (LW) |
|---|
| 269 | NRADIP = 3 ! Ice effective Radius: |
|---|
| 270 | ! 0 fixed 40 microns |
|---|
| 271 | ! 1 f(T) 40 - 130 microns |
|---|
| 272 | ! 2 f(T) 30 - 60 microns Jakob-Klein |
|---|
| 273 | ! 3 f(T,IWC) Sun-Rikus, 1999 |
|---|
| 274 | RMINICE= 15 ! Minimum Diameter for Ice Particles (micronm) |
|---|
| 275 | ! Needed only if NRADIP = 3 |
|---|
| 276 | ! (see von Walden et al., 2003 (Oct) Tab. 2 p.1393) |
|---|
| 277 | NRADLP = 0 ! Liquid effective Radius: f(Pressure) |
|---|
| 278 | ! 0 effective radius - liquid as f(Pressure) |
|---|
| 279 | ! 1 fixed 10 microns over land, 13 over ocean |
|---|
| 280 | ! 2 computed from LWC Martin et al, 1994 |
|---|
| 281 | NLIQOPT= 1 ! Cloud Optical Properties (Water): 1=ECMWF Operat. |
|---|
| 282 | ! 0 LW: Smith-Shi, 1992; SW: Fouquart, 1987 |
|---|
| 283 | ! 1 LW: Savijarvi, 1997; SW: Slingo , 1989 |
|---|
| 284 | ! 2 LW: Lindner,Li, 2000; SW: Slingo , 1989 |
|---|
| 285 | NICEOPT= 1 ! Cloud Optical Properties (Ice) : 1=ECMWF Operat. |
|---|
| 286 | ! 0 LW: Smith,Shi , 1992; SW: Ebert-Curry, 1992 |
|---|
| 287 | ! 1 LW: Ebert,Curry, 1992; SW: Ebert-Curry, 1992 |
|---|
| 288 | ! 2 LW: Fu,Liou , 1993; SW: Fu,Liou , 1993 |
|---|
| 289 | ! 3 LW: Fu et al. , 1998; SW: Fu , 1996 |
|---|
| 290 | NOVLP = 2 ! CLOUD OVERLAP CONFIGURATION: |
|---|
| 291 | ! 1=MRN, 2=MAX, 3=RAN, 4=Hogan |
|---|
| 292 | |
|---|
| 293 | |
|---|
| 294 | ! Switches (Aerosols/O3) |
|---|
| 295 | ! ---------------------- |
|---|
| 296 | |
|---|
| 297 | LNEWAER= .TRUE. ! Climatology of Aerosols: TEGEN ET AL. 1997 / GADS |
|---|
| 298 | LHVOLCA= .TRUE. ! .TRUE. IF GISS HISTORY OF VOLCANIC AEROSOLS IS ON |
|---|
| 299 | ! NOZOCL = -1 ! TESTS the vertical quadrature (NO absorber) |
|---|
| 300 | ! NOZOCL = 0 ! whatever is read for O3 as input profile |
|---|
| 301 | ! NOZOCL = 1 ! OLD ECMWF O3 climatology and aerosols |
|---|
| 302 | NOZOCL = 2 ! Fortuin-Langematz O3 climatology and aerosols |
|---|
| 303 | ! NOZOCL = 3 ! OLD ECMWF O3 climatology and NO aerosols |
|---|
| 304 | ! NOZOCL = 4 ! Fortuin-Langematz O3 climatology and NO aerosols |
|---|
| 305 | |
|---|
| 306 | |
|---|
| 307 | ! BASIC CONSTANTS |
|---|
| 308 | ! --------------- |
|---|
| 309 | |
|---|
| 310 | ! ***** |
|---|
| 311 | CALL SUCST (KULOUT, NINDAT, NSSSSS, KPRTLEV) ! Initialize common YOMCST |
|---|
| 312 | ! ***** ! (Basic Constants) |
|---|
| 313 | ! Initialize common YOMRIP |
|---|
| 314 | ! (only date and time) |
|---|
| 315 | |
|---|
| 316 | ! YOENCST - THERMODYNAMIC TRANSITION OF PHASE |
|---|
| 317 | ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ |
|---|
| 318 | RTT=273.16 ! |
|---|
| 319 | |
|---|
| 320 | ! YOETHF - DERIVED CONSTANTS SPECIFIC TO ECMWF THERMODYNAMICS |
|---|
| 321 | ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ |
|---|
| 322 | RTWAT=RTT ! |
|---|
| 323 | RTICE=RTT-23. ! |
|---|
| 324 | |
|---|
| 325 | ! YOERDU - CONTROL, PARAMETERS AND SECURITY IN RADIATION |
|---|
| 326 | ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ |
|---|
| 327 | RCDAY = RDAY * RG / RCPD ! |
|---|
| 328 | R10E = 0.4342945 ! DECIMAL / NATURAL |
|---|
| 329 | ! LOG. FACTOR |
|---|
| 330 | DIFF = 1.66 ! DIFFUSIVITY FACTOR |
|---|
| 331 | |
|---|
| 332 | |
|---|
| 333 | ! Space/Time Independant Coefficients |
|---|
| 334 | ! ----------------------------------- |
|---|
| 335 | |
|---|
| 336 | CALL SURDI ! ECMWF Surface Albedo, Emissivity |
|---|
| 337 | CALL SULWN ! Initialize common YOELW (new LW Coeff.) |
|---|
| 338 | CALL SUOLW ! Initialize common YOELW (old LW Coeff.) |
|---|
| 339 | |
|---|
| 340 | ! Initialization routine for RRTM |
|---|
| 341 | ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ |
|---|
| 342 | CALL SURRTAB ! AER'S RRTM LW RADIATION |
|---|
| 343 | CALL SURRTPK ! Initialize common YOERRTWN |
|---|
| 344 | ! (k-coefficients in spectral intervals) |
|---|
| 345 | CALL SURRTRF ! Initialize common YOERRTRF |
|---|
| 346 | ! (RRTM Reference Atmosphere) |
|---|
| 347 | CALL SURRTFTR ! Initialize common YOERRTRF |
|---|
| 348 | |
|---|
| 349 | ! RRTM routine ! BAND [cm-1] ! low ! high |
|---|
| 350 | ! ----------------------+---------------+---------------+--------- |
|---|
| 351 | CALL RRTM_KGB1 ! 1: 10- 250 ! H2O ! H2O |
|---|
| 352 | CALL RRTM_KGB2 ! 2: 250- 500 ! H2O ! H2O |
|---|
| 353 | CALL RRTM_KGB3 ! 3: 500- 630 ! H2O,CO2 ! H2O,CO2 |
|---|
| 354 | CALL RRTM_KGB4 ! 4: 630- 700 ! H2O,CO2 ! O3,CO2 |
|---|
| 355 | CALL RRTM_KGB5 ! 5: 700- 820 ! H2O,CO2 ! O3,CO2 |
|---|
| 356 | CALL RRTM_KGB6 ! 6: 820- 980 ! H2O ! nothing |
|---|
| 357 | CALL RRTM_KGB7 ! 7: 980-1080 ! H2O,O3 ! O3 |
|---|
| 358 | CALL RRTM_KGB8 ! 8: 1080-1180 ! (i.e.>~300mb) ! O3 |
|---|
| 359 | ! ! H2O ! |
|---|
| 360 | CALL RRTM_KGB9 ! 9: 1180-1390 ! H2O,CH4 ! CH4 |
|---|
| 361 | CALL RRTM_KGB10 ! 10: 1390-1480 ! H2O ! H2O |
|---|
| 362 | CALL RRTM_KGB11 ! 11: 1480-1800 ! H2O ! H2O |
|---|
| 363 | CALL RRTM_KGB12 ! 12: 1800-2080 ! H2O,CO2 ! nothing |
|---|
| 364 | CALL RRTM_KGB13 ! 13: 2080-2250 ! H2O,N2O ! nothing |
|---|
| 365 | CALL RRTM_KGB14 ! 14: 2250-2380 ! CO2 ! CO2 |
|---|
| 366 | CALL RRTM_KGB15 ! 15: 2380-2600 ! N2O,CO2 ! nothing |
|---|
| 367 | CALL RRTM_KGB16 ! 16: 2600-3000 ! H2O,CH4 ! nothing |
|---|
| 368 | |
|---|
| 369 | ! Reduce absorption coefficient data from 256 to 140 g-points |
|---|
| 370 | ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ |
|---|
| 371 | CALL RRTM_INIT_140GP |
|---|
| 372 | |
|---|
| 373 | ! Initialization routine for SW (6 spectral interval resolution) |
|---|
| 374 | ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ |
|---|
| 375 | CALL SUSWN ( NTSW , KSW ) ! Initialize common YOESW |
|---|
| 376 | |
|---|
| 377 | |
|---|
| 378 | write(6,*) 'INITIALISATION OF ECMWF RADIATIVE TRANSFERT: END ' |
|---|
| 379 | |
|---|
| 380 | |
|---|
| 381 | |
|---|
| 382 | |
|---|
| 383 | ! Radiation: Global (Time dependant) Parameters |
|---|
| 384 | ! ============================================= |
|---|
| 385 | |
|---|
| 386 | PRII05 = RI0 ! INSOLATION |
|---|
| 387 | ! PRII05 = RI0/(dST_UA*dST_UA) ! INSOLATION (dST_UA: Distance Soleil-Terre [UA]) |
|---|
| 388 | PCCO25 = 360.E-06*44./29. ! CONCENTRATION IN CO2 [Pa/Pa] |
|---|
| 389 | |
|---|
| 390 | |
|---|
| 391 | |
|---|
| 392 | |
|---|
| 393 | ! Surface Properties |
|---|
| 394 | ! ================== |
|---|
| 395 | |
|---|
| 396 | ! Land/sea Mask |
|---|
| 397 | ! ------------- |
|---|
| 398 | |
|---|
| 399 | ! Martin CONTROL |
|---|
| 400 | !PRINT*, 'CONTROL dans Atm_RT_INI' |
|---|
| 401 | !PRINT*,'SIZE(psa_DY)=)',SIZE(psa_DY) |
|---|
| 402 | !PRINT*, 'psa_DY=',psa_DY |
|---|
| 403 | ! Martin CONTROL |
|---|
| 404 | |
|---|
| 405 | DO jl=1,KLON |
|---|
| 406 | ikl = min(jl,kcolp) |
|---|
| 407 | i = ii__AP(ikl) |
|---|
| 408 | j = jj__AP(ikl) |
|---|
| 409 | PLSM5 (JL) = 1 - MaskSV_gpt(ikl) |
|---|
| 410 | |
|---|
| 411 | |
|---|
| 412 | |
|---|
| 413 | |
|---|
| 414 | ! Atmospheric Thermodynamics (Time and Space dependant) |
|---|
| 415 | ! ===================================================== |
|---|
| 416 | |
|---|
| 417 | ! Pressure Distribution |
|---|
| 418 | ! -------------------------- |
|---|
| 419 | |
|---|
| 420 | JK=1+KLEV |
|---|
| 421 | PAPH5 (JL,JK) = (psa_DY(ikl) + pt__DY) * 1000. ! Pressure (Layer Interface)[Pa] |
|---|
| 422 | DO JK=1,KLEV |
|---|
| 423 | PAPH5 (JL,JK) = (psa_DY(ikl) * sigma(JK) + pt__DY) * 1000. ! Pressure (Layer) [Pa] |
|---|
| 424 | PAP5 (JL,JK) = (psa_DY(ikl) * sigmi(JK) + pt__DY) * 1000. ! Pressure (Layer Interface)[Pa] |
|---|
| 425 | PDP5 (JL,JK) = psa_DY(ikl) *dsigmi(JK) * 1000. ! Pressure (Layer Thickness)[Pa] |
|---|
| 426 | |
|---|
| 427 | |
|---|
| 428 | ! Temperature Distribution |
|---|
| 429 | ! -------------------------- |
|---|
| 430 | |
|---|
| 431 | PT5 (JL,JK) = Ta__DY(ikl,JK) |
|---|
| 432 | |
|---|
| 433 | |
|---|
| 434 | ! Water Species Distribution |
|---|
| 435 | ! -------------------------- |
|---|
| 436 | |
|---|
| 437 | PQ5 (JL,JK) = qv__DY(ikl,JK) ! Water Vapor Concentr. [kg/kg] |
|---|
| 438 | |
|---|
| 439 | ENDDO |
|---|
| 440 | ENDDO |
|---|
| 441 | |
|---|
| 442 | |
|---|
| 443 | |
|---|
| 444 | |
|---|
| 445 | ! Initialization (Climatologies, Time independant) |
|---|
| 446 | ! ================================================ |
|---|
| 447 | |
|---|
| 448 | ! Aerosols Radiative Characteristics (YOEAER) |
|---|
| 449 | ! ---------------------------------- |
|---|
| 450 | |
|---|
| 451 | ! ******* |
|---|
| 452 | CALL SUAERL ! Aerosols LW Radiative Charact. |
|---|
| 453 | CALL SUAERSN ( NTSW , KSW ) ! Aerosols SW Radiative Charact. |
|---|
| 454 | ! ******* |
|---|
| 455 | |
|---|
| 456 | |
|---|
| 457 | ! Aerosols Optical Thickness Horizontal Distribution (model grid |
|---|
| 458 | ! -------------------------------------------------- independant) |
|---|
| 459 | |
|---|
| 460 | ! ******** |
|---|
| 461 | CALL SUAERH ! |
|---|
| 462 | CALL SUECAEBC ! BLACK CARBON (URBAN/FOREST FIRE ORIGIN) |
|---|
| 463 | CALL SUECAEOR ! ORGANIC-TYPE |
|---|
| 464 | CALL SUECAESD ! SOIL-DUST ORIGIN |
|---|
| 465 | CALL SUECAESS ! SEA -SALT ORIGIN |
|---|
| 466 | CALL SUECAESU ! SULFATE-TYPE |
|---|
| 467 | ! ******** |
|---|
| 468 | |
|---|
| 469 | |
|---|
| 470 | ! Clouds (YOECLD) |
|---|
| 471 | ! --------------- |
|---|
| 472 | |
|---|
| 473 | DO jk=1,klev |
|---|
| 474 | ZETA(jk) =PAP5(1,jk) /PAPH5(1,klev+1) |
|---|
| 475 | ENDDO |
|---|
| 476 | DO jk=1,klev+1 |
|---|
| 477 | ZETAH(jk)=PAPH5(1,jk)/PAPH5(1,klev+1) |
|---|
| 478 | ENDDO |
|---|
| 479 | ! Martin Control |
|---|
| 480 | !PRINT*, 'PAPH5(1,:)=',PAPH5(1,:) |
|---|
| 481 | !PRINT*, 'ZETAH=',ZETAH |
|---|
| 482 | ! Martin Control |
|---|
| 483 | |
|---|
| 484 | ! ***** |
|---|
| 485 | CALL SUCLD ( KLEV , ZETA ) |
|---|
| 486 | ! ***** |
|---|
| 487 | |
|---|
| 488 | |
|---|
| 489 | ! Cloud Optical Parameters SW/LW (all parameterizations) |
|---|
| 490 | ! ------------------------------------------------------ |
|---|
| 491 | |
|---|
| 492 | ! ******* |
|---|
| 493 | CALL SUCLOPN ( NTSW , KSW , KLEV ) ! Initialize YOECLOP |
|---|
| 494 | ! ******* |
|---|
| 495 | |
|---|
| 496 | |
|---|
| 497 | ! Radar Reflectivity |
|---|
| 498 | ! ------------------ |
|---|
| 499 | |
|---|
| 500 | ZAZH(KLEV+1)= 0. |
|---|
| 501 | ZAZ(1) =100000. |
|---|
| 502 | JL=KIDIA |
|---|
| 503 | DO jk=KLEV,2,-1 |
|---|
| 504 | ZTVIR = PT5(JL,jk) /(1.-0.608*PQ5(jl,jk)) |
|---|
| 505 | ZFACT = LOG(PAPH5(jl,jk+1))- LOG(PAPH5(jl,jk)) |
|---|
| 506 | ZAZH(jk) = ZAZH(jk+1) + R * ZTVIR/(RMD*RG)*ZFACT |
|---|
| 507 | ZAZ(jk) = 0.5*(ZAZH(jk+1)+ZAZH(jk))*1000. |
|---|
| 508 | END DO |
|---|
| 509 | |
|---|
| 510 | ! ****** |
|---|
| 511 | CALL SUOVLP ( KLEV , ZAZ ) ! Initialize ALPHA1 (%radar refl.) |
|---|
| 512 | ! ****** ! (Hogan & Illingsworth, 1999) |
|---|
| 513 | |
|---|
| 514 | |
|---|
| 515 | ! NO Absorber |
|---|
| 516 | ! ----------- |
|---|
| 517 | |
|---|
| 518 | IF (NOZOCL.EQ.-1) THEN |
|---|
| 519 | RCH4 = 1.E-18 |
|---|
| 520 | RN2O = 1.E-18 |
|---|
| 521 | RO3 = 1.E-18 |
|---|
| 522 | RCFC11 = 1.E-18 |
|---|
| 523 | RCFC12 = 1.E-18 |
|---|
| 524 | PCCO25 = 1.E-18 |
|---|
| 525 | DO jk=1,klev |
|---|
| 526 | DO jl=KIDIA,KFDIA |
|---|
| 527 | POZON5(JL,JK) = 0. |
|---|
| 528 | ENDDO |
|---|
| 529 | ENDDO |
|---|
| 530 | DO JK=1,KLEV |
|---|
| 531 | DO JAER=1,KAER |
|---|
| 532 | DO JL=KIDIA,KFDIA |
|---|
| 533 | PAER5(JL,JAER,JK)= ZEPAER |
|---|
| 534 | ENDDO |
|---|
| 535 | ENDDO |
|---|
| 536 | ENDDO |
|---|
| 537 | END IF |
|---|
| 538 | |
|---|
| 539 | |
|---|
| 540 | |
|---|
| 541 | return |
|---|
| 542 | end subroutine Atm_RT_INI |
|---|
