[3880] | 1 | SUBROUTINE RRTM_KGB2 |
---|
| 2 | |
---|
| 3 | ! Originally by Eli J. Mlawer, Atmospheric & Environmental Research. |
---|
| 4 | ! BAND 2: 250-500 cm-1 (low - H2O; high - H2O) |
---|
| 5 | ! Reformatted for F90 by JJMorcrette, ECMWF |
---|
| 6 | ! R. Elkhatib 12-10-2005 Split for faster and more robust compilation. |
---|
| 7 | ! G.Mozdzynski March 2011 read constants from files |
---|
| 8 | ! ABozzo May 2013 update to RRTMG v4.85 |
---|
| 9 | ! band 2: 350-500 cm-1 |
---|
| 10 | ! T. Wilhelmsson and K. Yessad (Oct 2013) Geometry and setup refactoring. |
---|
| 11 | ! ------------------------------------------------------------------ |
---|
| 12 | |
---|
| 13 | USE PARKIND1 ,ONLY : JPRB |
---|
| 14 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
---|
| 15 | USE YOMLUN ,ONLY : NULRAD |
---|
| 16 | USE MPL_MODULE,ONLY : MPL_BROADCAST |
---|
| 17 | USE YOMTAG ,ONLY : MTAGRAD |
---|
| 18 | |
---|
| 19 | USE YOERRTO2 , ONLY : KAO ,KBO ,KAO_D,KBO_D,SELFREFO ,FRACREFAO ,& |
---|
| 20 | & FRACREFBO ,FORREFO |
---|
| 21 | USE YOMMP0 , ONLY : NPROC, MYPROC |
---|
| 22 | |
---|
| 23 | ! ------------------------------------------------------------------ |
---|
| 24 | |
---|
| 25 | IMPLICIT NONE |
---|
| 26 | REAL(KIND=JPRB) :: ZHOOK_HANDLE |
---|
| 27 | |
---|
| 28 | #include "abor1.intfb.h" |
---|
| 29 | |
---|
| 30 | IF (LHOOK) CALL DR_HOOK('RRTM_KGB2',0,ZHOOK_HANDLE) |
---|
| 31 | |
---|
| 32 | IF( MYPROC==1 )THEN |
---|
| 33 | READ(NULRAD,ERR=1001) KAO_D,KBO_D |
---|
| 34 | ! Convert the data into model actual precision. |
---|
| 35 | KAO = REAL(KAO_D,JPRB) |
---|
| 36 | KBO = REAL(KBO_D,JPRB) |
---|
| 37 | ENDIF |
---|
| 38 | IF( NPROC>1 )THEN |
---|
| 39 | CALL MPL_BROADCAST (KAO,MTAGRAD,1,CDSTRING='RRTM_KGB2:') |
---|
| 40 | CALL MPL_BROADCAST (KBO,MTAGRAD,1,CDSTRING='RRTM_KGB2:') |
---|
| 41 | ENDIF |
---|
| 42 | |
---|
| 43 | |
---|
| 44 | ! Planck fraction mapping level: P = 1053.630 mbar, T = 294.2 K |
---|
| 45 | FRACREFAO(:) = (/ & |
---|
| 46 | & 1.6388e-01_JPRB, 1.5241e-01_JPRB, 1.4290e-01_JPRB, 1.2864e-01_JPRB, & |
---|
| 47 | & 1.1615e-01_JPRB, 1.0047e-01_JPRB, 8.0013e-02_JPRB, 6.0445e-02_JPRB, & |
---|
| 48 | & 4.0530e-02_JPRB, 4.3879e-03_JPRB, 3.5726e-03_JPRB, 2.7669e-03_JPRB, & |
---|
| 49 | & 2.0078e-03_JPRB, 1.2864e-03_JPRB, 4.7630e-04_JPRB, 6.9109e-05_JPRB/) |
---|
| 50 | |
---|
| 51 | ! Planck fraction mapping level: P = 3.206e-2 mb, T = 197.92 K |
---|
| 52 | FRACREFBO(:) = (/ & |
---|
| 53 | & 1.4697e-01_JPRB, 1.4826e-01_JPRB, 1.4278e-01_JPRB, 1.3320e-01_JPRB, & |
---|
| 54 | & 1.1965e-01_JPRB, 1.0297e-01_JPRB, 8.4170e-02_JPRB, 6.3282e-02_JPRB, & |
---|
| 55 | & 4.2868e-02_JPRB, 4.6644e-03_JPRB, 3.8619e-03_JPRB, 3.0533e-03_JPRB, & |
---|
| 56 | & 2.2359e-03_JPRB, 1.4226e-03_JPRB, 5.3642e-04_JPRB, 7.6316e-05_JPRB/) |
---|
| 57 | |
---|
| 58 | ! The array FORREFO contains the coefficient of the water vapor |
---|
| 59 | ! foreign-continuum (including the energy term). The first |
---|
| 60 | ! index refers to reference temperature (296,260,224,260) and |
---|
| 61 | ! pressure (970,475,219,3 mbar) levels. The second index |
---|
| 62 | ! runs over the g-channel (1 to 16). |
---|
| 63 | |
---|
| 64 | FORREFO(1,:) = (/ & |
---|
| 65 | & 2.8549e-03_JPRB,4.8281e-03_JPRB,6.2570e-03_JPRB,8.2731e-03_JPRB,7.9056e-03_JPRB,7.7840e-03_JPRB, & |
---|
| 66 | & 1.0115e-02_JPRB,9.6599e-03_JPRB,1.0153e-02_JPRB,1.0921e-02_JPRB,1.2408e-02_JPRB,1.3496e-02_JPRB, & |
---|
| 67 | & 1.5059e-02_JPRB,1.4636e-02_JPRB,1.6483e-02_JPRB,1.2394e-02_JPRB/) |
---|
| 68 | FORREFO(2,:) = (/ & |
---|
| 69 | & 3.0036e-03_JPRB,5.1093e-03_JPRB,5.7317e-03_JPRB,9.2246e-03_JPRB,8.9829e-03_JPRB,8.6477e-03_JPRB, & |
---|
| 70 | & 1.1448e-02_JPRB,1.0391e-02_JPRB,1.0211e-02_JPRB,1.2921e-02_JPRB,1.2726e-02_JPRB,1.2426e-02_JPRB, & |
---|
| 71 | & 1.4609e-02_JPRB,1.5783e-02_JPRB,1.6617e-02_JPRB,1.6858e-02_JPRB/) |
---|
| 72 | FORREFO(3,:) = (/ & |
---|
| 73 | & 3.0771e-03_JPRB,5.1206e-03_JPRB,5.8426e-03_JPRB,9.5727e-03_JPRB,1.0338e-02_JPRB,9.3737e-03_JPRB, & |
---|
| 74 | & 1.2805e-02_JPRB,1.1272e-02_JPRB,1.1353e-02_JPRB,1.1837e-02_JPRB,1.1550e-02_JPRB,1.3020e-02_JPRB, & |
---|
| 75 | & 1.3536e-02_JPRB,1.6226e-02_JPRB,1.6039e-02_JPRB,2.2578e-02_JPRB/) |
---|
| 76 | FORREFO(4,:) = (/ & |
---|
| 77 | & 3.3072e-03_JPRB,5.0240e-03_JPRB,6.8474e-03_JPRB,8.2736e-03_JPRB,8.6151e-03_JPRB,8.6762e-03_JPRB, & |
---|
| 78 | & 1.1476e-02_JPRB,1.0246e-02_JPRB,1.0819e-02_JPRB,1.0640e-02_JPRB,1.0545e-02_JPRB,1.0533e-02_JPRB, & |
---|
| 79 | & 1.0496e-02_JPRB,1.0142e-02_JPRB,9.7979e-03_JPRB,1.5255e-02_JPRB/) |
---|
| 80 | |
---|
| 81 | |
---|
| 82 | ! The following are parameters related to the reference water vapor |
---|
| 83 | ! mixing ratios by REFPARAM(I) = REFH2O(I) / (.002+REFH2O(I)). |
---|
| 84 | ! These parameters are used for the Planck function interpolation. |
---|
| 85 | !REFPARAM( :) = (/& |
---|
| 86 | ! & 0.903661_JPRB , 0.859386_JPRB , 0.746542_JPRB , 0.580496_JPRB , 0.412889_JPRB ,& |
---|
| 87 | ! & 0.275283_JPRB , 0.162745_JPRB , 7.63929E-02_JPRB, 1.82553E-02_JPRB, 3.72432E-03_JPRB, & |
---|
| 88 | ! & 2.14946E-03_JPRB, 1.66320E-03_JPRB, 1.59940E-03_JPRB/) |
---|
| 89 | |
---|
| 90 | ! ------------------------------------------------------------------ |
---|
| 91 | |
---|
| 92 | ! The array KAO contains absorption coefs at the 16 chosen g-values |
---|
| 93 | ! for a range of pressure levels > ~100mb and temperatures. The first |
---|
| 94 | ! index in the array, JT, which runs from 1 to 5, corresponds to |
---|
| 95 | ! different temperatures. More specifically, JT = 3 means that the |
---|
| 96 | ! data are for the corresponding TREF for this pressure level, |
---|
| 97 | ! JT = 2 refers to the temperatureTREF-15, JT = 1 is for TREF-30, |
---|
| 98 | ! JT = 4 is for TREF+15, and JT = 5 is for TREF+30. The second |
---|
| 99 | ! index, JP, runs from 1 to 13 and refers to the corresponding |
---|
| 100 | ! pressure level in PREF (e.g. JP = 1 is for a pressure of 1053.63 mb). |
---|
| 101 | ! The third index, IG, goes from 1 to 16, and tells us which |
---|
| 102 | ! g-interval the absorption coefficients are for. |
---|
| 103 | |
---|
| 104 | |
---|
| 105 | |
---|
| 106 | ! The array KBO contains absorption coefs at the 16 chosen g-values |
---|
| 107 | ! for a range of pressure levels < ~100mb and temperatures. The first |
---|
| 108 | ! index in the array, JT, which runs from 1 to 5, corresponds to |
---|
| 109 | ! different temperatures. More specifically, JT = 3 means that the |
---|
| 110 | ! data are for the reference temperature TREF for this pressure |
---|
| 111 | ! level, JT = 2 refers to the temperature TREF-15, JT = 1 is for |
---|
| 112 | ! TREF-30, JT = 4 is for TREF+15, and JT = 5 is for TREF+30. |
---|
| 113 | ! The second index, JP, runs from 13 to 59 and refers to the JPth |
---|
| 114 | ! reference pressure level (see taumol.f for the value of these |
---|
| 115 | ! pressure levels in mb). The third index, IG, goes from 1 to 16, |
---|
| 116 | ! and tells us which g-interval the absorption coefficients are for. |
---|
| 117 | |
---|
| 118 | |
---|
| 119 | |
---|
| 120 | ! The array SELFREFO contains the coefficient of the water vapor |
---|
| 121 | ! self-continuum (including the energy term). The first index |
---|
| 122 | ! refers to temperature in 7.2 degree increments. For instance, |
---|
| 123 | ! JT = 1 refers to a temperature of 245.6, JT = 2 refers to 252.8, |
---|
| 124 | ! etc. The second index runs over the g-channel (1 to 16). |
---|
| 125 | |
---|
| 126 | SELFREFO(:, 1) = (/ & |
---|
| 127 | & 7.25695e-01_JPRB, 6.53591e-01_JPRB, 5.88650e-01_JPRB, 5.30162e-01_JPRB, 4.77485e-01_JPRB, & |
---|
| 128 | & 4.30042e-01_JPRB, 3.87313e-01_JPRB, 3.48830e-01_JPRB, 3.14170e-01_JPRB, 2.82954e-01_JPRB/) |
---|
| 129 | SELFREFO(:, 2) = (/ & |
---|
| 130 | & 9.61996e-01_JPRB, 8.77853e-01_JPRB, 8.01070e-01_JPRB, 7.31003e-01_JPRB, 6.67064e-01_JPRB, & |
---|
| 131 | & 6.08718e-01_JPRB, 5.55476e-01_JPRB, 5.06890e-01_JPRB, 4.62554e-01_JPRB, 4.22096e-01_JPRB/) |
---|
| 132 | SELFREFO(:, 3) = (/ & |
---|
| 133 | & 9.72584e-01_JPRB, 9.02658e-01_JPRB, 8.37760e-01_JPRB, 7.77527e-01_JPRB, 7.21626e-01_JPRB, & |
---|
| 134 | & 6.69743e-01_JPRB, 6.21591e-01_JPRB, 5.76900e-01_JPRB, 5.35423e-01_JPRB, 4.96927e-01_JPRB/) |
---|
| 135 | SELFREFO(:, 4) = (/ & |
---|
| 136 | & 1.24790e+00_JPRB, 1.14353e+00_JPRB, 1.04790e+00_JPRB, 9.60263e-01_JPRB, 8.79956e-01_JPRB, & |
---|
| 137 | & 8.06364e-01_JPRB, 7.38927e-01_JPRB, 6.77130e-01_JPRB, 6.20501e-01_JPRB, 5.68608e-01_JPRB/) |
---|
| 138 | SELFREFO(:, 5) = (/ & |
---|
| 139 | & 1.23574e+00_JPRB, 1.12928e+00_JPRB, 1.03200e+00_JPRB, 9.43096e-01_JPRB, 8.61851e-01_JPRB, & |
---|
| 140 | & 7.87605e-01_JPRB, 7.19755e-01_JPRB, 6.57750e-01_JPRB, 6.01087e-01_JPRB, 5.49305e-01_JPRB/) |
---|
| 141 | SELFREFO(:, 6) = (/ & |
---|
| 142 | & 1.20921e+00_JPRB, 1.10660e+00_JPRB, 1.01270e+00_JPRB, 9.26766e-01_JPRB, 8.48124e-01_JPRB, & |
---|
| 143 | & 7.76155e-01_JPRB, 7.10293e-01_JPRB, 6.50020e-01_JPRB, 5.94861e-01_JPRB, 5.44384e-01_JPRB/) |
---|
| 144 | SELFREFO(:, 7) = (/ & |
---|
| 145 | & 1.38112e+00_JPRB, 1.26727e+00_JPRB, 1.16280e+00_JPRB, 1.06694e+00_JPRB, 9.78990e-01_JPRB, & |
---|
| 146 | & 8.98287e-01_JPRB, 8.24236e-01_JPRB, 7.56290e-01_JPRB, 6.93945e-01_JPRB, 6.36739e-01_JPRB/) |
---|
| 147 | SELFREFO(:, 8) = (/ & |
---|
| 148 | & 1.30321e+00_JPRB, 1.20127e+00_JPRB, 1.10730e+00_JPRB, 1.02068e+00_JPRB, 9.40840e-01_JPRB, & |
---|
| 149 | & 8.67243e-01_JPRB, 7.99403e-01_JPRB, 7.36870e-01_JPRB, 6.79229e-01_JPRB, 6.26096e-01_JPRB/) |
---|
| 150 | SELFREFO(:, 9) = (/ & |
---|
| 151 | & 1.26713e+00_JPRB, 1.17927e+00_JPRB, 1.09750e+00_JPRB, 1.02140e+00_JPRB, 9.50575e-01_JPRB, & |
---|
| 152 | & 8.84662e-01_JPRB, 8.23319e-01_JPRB, 7.66230e-01_JPRB, 7.13099e-01_JPRB, 6.63653e-01_JPRB/) |
---|
| 153 | SELFREFO(:,10) = (/ & |
---|
| 154 | & 1.49824e+00_JPRB, 1.37053e+00_JPRB, 1.25370e+00_JPRB, 1.14683e+00_JPRB, 1.04908e+00_JPRB, & |
---|
| 155 | & 9.59651e-01_JPRB, 8.77849e-01_JPRB, 8.03020e-01_JPRB, 7.34569e-01_JPRB, 6.71954e-01_JPRB/) |
---|
| 156 | SELFREFO(:,11) = (/ & |
---|
| 157 | & 1.44786e+00_JPRB, 1.34594e+00_JPRB, 1.25120e+00_JPRB, 1.16313e+00_JPRB, 1.08125e+00_JPRB, & |
---|
| 158 | & 1.00514e+00_JPRB, 9.34392e-01_JPRB, 8.68620e-01_JPRB, 8.07477e-01_JPRB, 7.50639e-01_JPRB/) |
---|
| 159 | SELFREFO(:,12) = (/ & |
---|
| 160 | & 1.38460e+00_JPRB, 1.30437e+00_JPRB, 1.22880e+00_JPRB, 1.15760e+00_JPRB, 1.09053e+00_JPRB, & |
---|
| 161 | & 1.02735e+00_JPRB, 9.67825e-01_JPRB, 9.11750e-01_JPRB, 8.58924e-01_JPRB, 8.09159e-01_JPRB/) |
---|
| 162 | SELFREFO(:,13) = (/ & |
---|
| 163 | & 1.51953e+00_JPRB, 1.42822e+00_JPRB, 1.34240e+00_JPRB, 1.26173e+00_JPRB, 1.18592e+00_JPRB, & |
---|
| 164 | & 1.11465e+00_JPRB, 1.04768e+00_JPRB, 9.84720e-01_JPRB, 9.25548e-01_JPRB, 8.69932e-01_JPRB/) |
---|
| 165 | SELFREFO(:,14) = (/ & |
---|
| 166 | & 1.62608e+00_JPRB, 1.51021e+00_JPRB, 1.40260e+00_JPRB, 1.30266e+00_JPRB, 1.20983e+00_JPRB, & |
---|
| 167 | & 1.12363e+00_JPRB, 1.04356e+00_JPRB, 9.69200e-01_JPRB, 9.00138e-01_JPRB, 8.35998e-01_JPRB/) |
---|
| 168 | SELFREFO(:,15) = (/ & |
---|
| 169 | & 1.65383e+00_JPRB, 1.54808e+00_JPRB, 1.44910e+00_JPRB, 1.35644e+00_JPRB, 1.26971e+00_JPRB, & |
---|
| 170 | & 1.18853e+00_JPRB, 1.11254e+00_JPRB, 1.04140e+00_JPRB, 9.74813e-01_JPRB, 9.12484e-01_JPRB/) |
---|
| 171 | SELFREFO(:,16) = (/ & |
---|
| 172 | & 1.78105e+00_JPRB, 1.61421e+00_JPRB, 1.46300e+00_JPRB, 1.32595e+00_JPRB, 1.20174e+00_JPRB, & |
---|
| 173 | & 1.08917e+00_JPRB, 9.87141e-01_JPRB, 8.94670e-01_JPRB, 8.10861e-01_JPRB, 7.34904e-01_JPRB/) |
---|
| 174 | |
---|
| 175 | |
---|
| 176 | IF (LHOOK) CALL DR_HOOK('RRTM_KGB2',1,ZHOOK_HANDLE) |
---|
| 177 | RETURN |
---|
| 178 | |
---|
| 179 | 1001 CONTINUE |
---|
| 180 | CALL ABOR1("RRTM_KGB2:ERROR READING FILE RADRRTM") |
---|
| 181 | |
---|
| 182 | END SUBROUTINE RRTM_KGB2 |
---|