| 1 | SUBROUTINE SRTM_KGB18 |
|---|
| 2 | |
|---|
| 3 | ! Originally by J.Delamere, Atmospheric & Environmental Research. |
|---|
| 4 | ! Revision: 2.4 |
|---|
| 5 | ! BAND 18: 4000-4650 cm-1 (low - H2O,CH4; high - CH4) |
|---|
| 6 | ! Reformatted for F90 by JJMorcrette, ECMWF |
|---|
| 7 | ! R. Elkhatib 12-10-2005 Split for faster and more robust compilation. |
|---|
| 8 | ! G.Mozdzynski March 2011 read constants from files |
|---|
| 9 | ! T. Wilhelmsson and K. Yessad (Oct 2013) Geometry and setup refactoring. |
|---|
| 10 | ! F. Vana 05-Mar-2015 Support for single precision |
|---|
| 11 | ! ------------------------------------------------------------------ |
|---|
| 12 | |
|---|
| 13 | USE PARKIND1 , ONLY : JPRB |
|---|
| 14 | USE YOMHOOK , ONLY : LHOOK, DR_HOOK, JPHOOK |
|---|
| 15 | USE YOMLUN , ONLY : NULRAD |
|---|
| 16 | USE YOMMP0 , ONLY : NPROC, MYPROC |
|---|
| 17 | USE MPL_MODULE, ONLY : MPL_BROADCAST |
|---|
| 18 | USE YOMTAG , ONLY : MTAGRAD |
|---|
| 19 | USE YOESRTA18 , ONLY : KA, KB, SELFREF, FORREF, SFLUXREF, RAYL, STRRAT, LAYREFFR, & |
|---|
| 20 | & KA_D, KB_D |
|---|
| 21 | |
|---|
| 22 | ! ------------------------------------------------------------------ |
|---|
| 23 | |
|---|
| 24 | IMPLICIT NONE |
|---|
| 25 | |
|---|
| 26 | ! KURUCZ |
|---|
| 27 | REAL(KIND=JPHOOK) :: ZHOOK_HANDLE |
|---|
| 28 | |
|---|
| 29 | #include "abor1.intfb.h" |
|---|
| 30 | |
|---|
| 31 | IF (LHOOK) CALL DR_HOOK('SRTM_KGB18',0,ZHOOK_HANDLE) |
|---|
| 32 | |
|---|
| 33 | IF( MYPROC==1 )THEN |
|---|
| 34 | READ(NULRAD,ERR=1001) KA_D,KB_D |
|---|
| 35 | KA = REAL(KA_D,JPRB) |
|---|
| 36 | KB = REAL(KB_D,JPRB) |
|---|
| 37 | ENDIF |
|---|
| 38 | IF( NPROC>1 )THEN |
|---|
| 39 | CALL MPL_BROADCAST (KA,MTAGRAD,1,CDSTRING='SRTM_KGB18:') |
|---|
| 40 | CALL MPL_BROADCAST (KB,MTAGRAD,1,CDSTRING='SRTM_KGB18:') |
|---|
| 41 | ENDIF |
|---|
| 42 | |
|---|
| 43 | SFLUXREF(:,1) = (/ & |
|---|
| 44 | & 3.65840_JPRB , 3.54375_JPRB , 3.34481_JPRB , 3.10534_JPRB , & |
|---|
| 45 | & 2.79879_JPRB , 2.42841_JPRB , 1.98748_JPRB , 1.49377_JPRB , & |
|---|
| 46 | & 1.00196_JPRB , 0.108342_JPRB , 8.95099E-02_JPRB, 7.05199E-02_JPRB, & |
|---|
| 47 | & 5.16432E-02_JPRB, 3.27635E-02_JPRB, 1.25133E-02_JPRB, 1.73001E-03_JPRB /) |
|---|
| 48 | SFLUXREF(:,2) = (/ & |
|---|
| 49 | & 3.86372_JPRB , 3.48521_JPRB , 3.30790_JPRB , 3.08103_JPRB , & |
|---|
| 50 | & 2.77552_JPRB , 2.40722_JPRB , 1.97307_JPRB , 1.48023_JPRB , & |
|---|
| 51 | & 0.993055_JPRB , 0.107691_JPRB , 8.84430E-02_JPRB, 6.99354E-02_JPRB, & |
|---|
| 52 | & 5.07881E-02_JPRB, 3.24121E-02_JPRB, 1.19442E-02_JPRB, 1.57612E-03_JPRB /) |
|---|
| 53 | SFLUXREF(:,3) = (/ & |
|---|
| 54 | & 3.90370_JPRB , 3.50657_JPRB , 3.30629_JPRB , 3.06046_JPRB , & |
|---|
| 55 | & 2.76982_JPRB , 2.39907_JPRB , 1.96358_JPRB , 1.47458_JPRB , & |
|---|
| 56 | & 0.988475_JPRB , 0.106698_JPRB , 8.75242E-02_JPRB, 6.85898E-02_JPRB, & |
|---|
| 57 | & 5.04798E-02_JPRB, 3.13718E-02_JPRB, 1.09533E-02_JPRB, 1.57612E-03_JPRB /) |
|---|
| 58 | SFLUXREF(:,4) = (/ & |
|---|
| 59 | & 3.93165_JPRB , 3.52058_JPRB , 3.31346_JPRB , 3.04944_JPRB , & |
|---|
| 60 | & 2.76074_JPRB , 2.39433_JPRB , 1.95556_JPRB , 1.46712_JPRB , & |
|---|
| 61 | & 0.984056_JPRB , 0.105885_JPRB , 8.73062E-02_JPRB, 6.84054E-02_JPRB, & |
|---|
| 62 | & 4.87443E-02_JPRB, 2.99295E-02_JPRB, 1.09533E-02_JPRB, 1.57612E-03_JPRB /) |
|---|
| 63 | SFLUXREF(:,5) = (/ & |
|---|
| 64 | & 3.94082_JPRB , 3.55221_JPRB , 3.31863_JPRB , 3.04730_JPRB , & |
|---|
| 65 | & 2.74918_JPRB , 2.38328_JPRB , 1.95212_JPRB , 1.45889_JPRB , & |
|---|
| 66 | & 0.978888_JPRB , 0.105102_JPRB , 8.65732E-02_JPRB, 6.74563E-02_JPRB, & |
|---|
| 67 | & 4.76592E-02_JPRB, 2.91017E-02_JPRB, 1.09533E-02_JPRB, 1.57612E-03_JPRB /) |
|---|
| 68 | SFLUXREF(:,6) = (/ & |
|---|
| 69 | & 3.94198_JPRB , 3.58743_JPRB , 3.32106_JPRB , 3.05866_JPRB , & |
|---|
| 70 | & 2.74115_JPRB , 2.36939_JPRB , 1.94305_JPRB , 1.45180_JPRB , & |
|---|
| 71 | & 0.971784_JPRB , 1.04045E-01_JPRB, 8.53731E-02_JPRB, 6.60654E-02_JPRB, & |
|---|
| 72 | & 4.63228E-02_JPRB, 2.91016E-02_JPRB, 1.09552E-02_JPRB, 1.57612E-03_JPRB /) |
|---|
| 73 | SFLUXREF(:,7) = (/ & |
|---|
| 74 | & 3.93596_JPRB , 3.63366_JPRB , 3.33144_JPRB , 3.06252_JPRB , & |
|---|
| 75 | & 2.74054_JPRB , 2.35492_JPRB , 1.92769_JPRB , 1.44300_JPRB , & |
|---|
| 76 | & 0.961809_JPRB , 1.02867E-01_JPRB, 8.34164E-02_JPRB, 6.41005E-02_JPRB, & |
|---|
| 77 | & 4.61826E-02_JPRB, 2.91006E-02_JPRB, 1.09553E-02_JPRB, 1.57612E-03_JPRB /) |
|---|
| 78 | SFLUXREF(:,8) = (/ & |
|---|
| 79 | & 3.92520_JPRB , 3.69078_JPRB , 3.35656_JPRB , 3.07055_JPRB , & |
|---|
| 80 | & 2.73862_JPRB , 2.34430_JPRB , 1.90187_JPRB , 1.42242_JPRB , & |
|---|
| 81 | & 0.946676_JPRB , 9.96302E-02_JPRB, 8.14421E-02_JPRB, 6.38622E-02_JPRB, & |
|---|
| 82 | & 4.61794E-02_JPRB, 2.91017E-02_JPRB, 1.09553E-02_JPRB, 1.57612E-03_JPRB /) |
|---|
| 83 | SFLUXREF(:,9) = (/ & |
|---|
| 84 | & 3.80721_JPRB , 3.74437_JPRB , 3.50205_JPRB , 3.18009_JPRB , & |
|---|
| 85 | & 2.75757_JPRB , 2.29188_JPRB , 1.84382_JPRB , 1.35694_JPRB , & |
|---|
| 86 | & 0.914040_JPRB , 9.86811E-02_JPRB, 8.14321E-02_JPRB, 6.38541E-02_JPRB, & |
|---|
| 87 | & 4.61795E-02_JPRB, 2.90960E-02_JPRB, 1.09613E-02_JPRB, 1.57612E-03_JPRB /) |
|---|
| 88 | |
|---|
| 89 | ! Rayleigh extinction coefficient at v = 4325 cm-1. |
|---|
| 90 | RAYL = 1.39E-09_JPRB |
|---|
| 91 | |
|---|
| 92 | STRRAT = 38.9589_JPRB |
|---|
| 93 | |
|---|
| 94 | LAYREFFR = 6 |
|---|
| 95 | |
|---|
| 96 | ! ------------------------------------------------------------------ |
|---|
| 97 | |
|---|
| 98 | ! The array KA contains absorption coefs at the 16 chosen g-values |
|---|
| 99 | ! for a range of pressure levels> ~100mb, temperatures, and binary |
|---|
| 100 | ! species parameters (see taumol.f for definition). The first |
|---|
| 101 | ! index in the array, JS, runs from 1 to 9, and corresponds to |
|---|
| 102 | ! different values of the binary species parameter. For instance, |
|---|
| 103 | ! JS=1 refers to dry air, JS = 2 corresponds to the paramter value 1/8, |
|---|
| 104 | ! JS = 3 corresponds to the parameter value 2/8, etc. The second index |
|---|
| 105 | ! in the array, JT, which runs from 1 to 5, corresponds to different |
|---|
| 106 | ! temperatures. More specifically, JT = 3 means that the data are for |
|---|
| 107 | ! the reference temperature TREF for this pressure level, JT = 2 refers |
|---|
| 108 | ! to TREF-15, JT = 1 is for TREF-30, JT = 4 is for TREF+15, and JT = 5 |
|---|
| 109 | ! is for TREF+30. The third index, JP, runs from 1 to 13 and refers |
|---|
| 110 | ! to the JPth reference pressure level (see taumol.f for these levels |
|---|
| 111 | ! in mb). The fourth index, IG, goes from 1 to 16, and indicates |
|---|
| 112 | ! which g-interval the absorption coefficients are for. |
|---|
| 113 | ! ----------------------------------------------------------------- |
|---|
| 114 | |
|---|
| 115 | ! ----------------------------------------------------------------- |
|---|
| 116 | ! The array KB contains absorption coefs at the 16 chosen g-values |
|---|
| 117 | ! for a range of pressure levels < ~100mb and temperatures. The first |
|---|
| 118 | ! index in the array, JT, which runs from 1 to 5, corresponds to |
|---|
| 119 | ! different temperatures. More specifically, JT = 3 means that the |
|---|
| 120 | ! data are for the reference temperature TREF for this pressure |
|---|
| 121 | ! level, JT = 2 refers to the temperature TREF-15, JT = 1 is for |
|---|
| 122 | ! TREF-30, JT = 4 is for TREF+15, and JT = 5 is for TREF+30. |
|---|
| 123 | ! The second index, JP, runs from 13 to 59 and refers to the JPth |
|---|
| 124 | ! reference pressure level (see taumol.f for the value of these |
|---|
| 125 | ! pressure levels in mb). The third index, IG, goes from 1 to 16, |
|---|
| 126 | ! and tells us which g-interval the absorption coefficients are for. |
|---|
| 127 | ! ----------------------------------------------------------------- |
|---|
| 128 | |
|---|
| 129 | |
|---|
| 130 | FORREF(:, 1) = (/ 0.860560E-06_JPRB, 0.130439E-05_JPRB, 0.382378E-05_JPRB /) |
|---|
| 131 | FORREF(:, 2) = (/ 0.817926E-06_JPRB, 0.158599E-05_JPRB, 0.658771E-04_JPRB /) |
|---|
| 132 | FORREF(:, 3) = (/ 0.129369E-05_JPRB, 0.824406E-05_JPRB, 0.952778E-04_JPRB /) |
|---|
| 133 | FORREF(:, 4) = (/ 0.438918E-05_JPRB, 0.375356E-04_JPRB, 0.119111E-03_JPRB /) |
|---|
| 134 | FORREF(:, 5) = (/ 0.306057E-04_JPRB, 0.622798E-04_JPRB, 0.100740E-03_JPRB /) |
|---|
| 135 | FORREF(:, 6) = (/ 0.891934E-04_JPRB, 0.856393E-04_JPRB, 0.635583E-04_JPRB /) |
|---|
| 136 | FORREF(:, 7) = (/ 0.171959E-03_JPRB, 0.173431E-03_JPRB, 0.611721E-04_JPRB /) |
|---|
| 137 | FORREF(:, 8) = (/ 0.357795E-03_JPRB, 0.247261E-03_JPRB, 0.488864E-04_JPRB /) |
|---|
| 138 | FORREF(:, 9) = (/ 0.326623E-03_JPRB, 0.289471E-03_JPRB, 0.548834E-04_JPRB /) |
|---|
| 139 | FORREF(:,10) = (/ 0.345103E-03_JPRB, 0.320898E-03_JPRB, 0.633214E-04_JPRB /) |
|---|
| 140 | FORREF(:,11) = (/ 0.392567E-03_JPRB, 0.325153E-03_JPRB, 0.744479E-04_JPRB /) |
|---|
| 141 | FORREF(:,12) = (/ 0.349277E-03_JPRB, 0.345610E-03_JPRB, 0.916479E-04_JPRB /) |
|---|
| 142 | FORREF(:,13) = (/ 0.425161E-03_JPRB, 0.348452E-03_JPRB, 0.125788E-03_JPRB /) |
|---|
| 143 | FORREF(:,14) = (/ 0.407594E-03_JPRB, 0.435836E-03_JPRB, 0.287583E-03_JPRB /) |
|---|
| 144 | FORREF(:,15) = (/ 0.521605E-03_JPRB, 0.486596E-03_JPRB, 0.483511E-03_JPRB /) |
|---|
| 145 | FORREF(:,16) = (/ 0.773790E-03_JPRB, 0.737247E-03_JPRB, 0.665939E-03_JPRB /) |
|---|
| 146 | |
|---|
| 147 | ! ----------------------------------------------------------------- |
|---|
| 148 | ! The array SELFREF contains the coefficient of the water vapor |
|---|
| 149 | ! self-continuum (including the energy term). The first index |
|---|
| 150 | ! refers to temperature in 7.2 degree increments. For instance, |
|---|
| 151 | ! JT = 1 refers to a temperature of 245.6, JT = 2 refers to 252.8, |
|---|
| 152 | ! etc. The second index runs over the g-channel (1 to 16). |
|---|
| 153 | |
|---|
| 154 | SELFREF(:, 1) = (/ & |
|---|
| 155 | & 0.750370E-03_JPRB, 0.644938E-03_JPRB, 0.554321E-03_JPRB, 0.476436E-03_JPRB, 0.409494E-03_JPRB, & |
|---|
| 156 | & 0.351957E-03_JPRB, 0.302505E-03_JPRB, 0.260002E-03_JPRB, 0.223470E-03_JPRB, 0.192071E-03_JPRB /) |
|---|
| 157 | SELFREF(:, 2) = (/ & |
|---|
| 158 | & 0.136135E-02_JPRB, 0.113187E-02_JPRB, 0.941076E-03_JPRB, 0.782440E-03_JPRB, 0.650546E-03_JPRB, & |
|---|
| 159 | & 0.540885E-03_JPRB, 0.449709E-03_JPRB, 0.373902E-03_JPRB, 0.310874E-03_JPRB, 0.258471E-03_JPRB /) |
|---|
| 160 | SELFREF(:, 3) = (/ & |
|---|
| 161 | & 0.333950E-02_JPRB, 0.256391E-02_JPRB, 0.196845E-02_JPRB, 0.151129E-02_JPRB, 0.116030E-02_JPRB, & |
|---|
| 162 | & 0.890824E-03_JPRB, 0.683934E-03_JPRB, 0.525093E-03_JPRB, 0.403143E-03_JPRB, 0.309515E-03_JPRB /) |
|---|
| 163 | SELFREF(:, 4) = (/ & |
|---|
| 164 | & 0.793392E-02_JPRB, 0.589865E-02_JPRB, 0.438548E-02_JPRB, 0.326048E-02_JPRB, 0.242408E-02_JPRB, & |
|---|
| 165 | & 0.180223E-02_JPRB, 0.133991E-02_JPRB, 0.996186E-03_JPRB, 0.740636E-03_JPRB, 0.550642E-03_JPRB /) |
|---|
| 166 | SELFREF(:, 5) = (/ & |
|---|
| 167 | & 0.828169E-02_JPRB, 0.703139E-02_JPRB, 0.596984E-02_JPRB, 0.506856E-02_JPRB, 0.430335E-02_JPRB, & |
|---|
| 168 | & 0.365366E-02_JPRB, 0.310206E-02_JPRB, 0.263374E-02_JPRB, 0.223612E-02_JPRB, 0.189852E-02_JPRB /) |
|---|
| 169 | SELFREF(:, 6) = (/ & |
|---|
| 170 | & 0.834190E-02_JPRB, 0.780225E-02_JPRB, 0.729750E-02_JPRB, 0.682541E-02_JPRB, 0.638386E-02_JPRB, & |
|---|
| 171 | & 0.597087E-02_JPRB, 0.558460E-02_JPRB, 0.522332E-02_JPRB, 0.488541E-02_JPRB, 0.456936E-02_JPRB /) |
|---|
| 172 | SELFREF(:, 7) = (/ & |
|---|
| 173 | & 0.119082E-01_JPRB, 0.112566E-01_JPRB, 0.106406E-01_JPRB, 0.100583E-01_JPRB, 0.950785E-02_JPRB, & |
|---|
| 174 | & 0.898755E-02_JPRB, 0.849571E-02_JPRB, 0.803080E-02_JPRB, 0.759132E-02_JPRB, 0.717590E-02_JPRB /) |
|---|
| 175 | SELFREF(:, 8) = (/ & |
|---|
| 176 | & 0.144004E-01_JPRB, 0.141762E-01_JPRB, 0.139554E-01_JPRB, 0.137381E-01_JPRB, 0.135241E-01_JPRB, & |
|---|
| 177 | & 0.133135E-01_JPRB, 0.131062E-01_JPRB, 0.129021E-01_JPRB, 0.127011E-01_JPRB, 0.125033E-01_JPRB /) |
|---|
| 178 | SELFREF(:, 9) = (/ & |
|---|
| 179 | & 0.186171E-01_JPRB, 0.175281E-01_JPRB, 0.165027E-01_JPRB, 0.155373E-01_JPRB, 0.146284E-01_JPRB, & |
|---|
| 180 | & 0.137726E-01_JPRB, 0.129670E-01_JPRB, 0.122084E-01_JPRB, 0.114942E-01_JPRB, 0.108218E-01_JPRB /) |
|---|
| 181 | SELFREF(:,10) = (/ & |
|---|
| 182 | & 0.209396E-01_JPRB, 0.195077E-01_JPRB, 0.181737E-01_JPRB, 0.169309E-01_JPRB, 0.157731E-01_JPRB, & |
|---|
| 183 | & 0.146945E-01_JPRB, 0.136897E-01_JPRB, 0.127535E-01_JPRB, 0.118814E-01_JPRB, 0.110689E-01_JPRB /) |
|---|
| 184 | SELFREF(:,11) = (/ & |
|---|
| 185 | & 0.203661E-01_JPRB, 0.193311E-01_JPRB, 0.183487E-01_JPRB, 0.174163E-01_JPRB, 0.165312E-01_JPRB, & |
|---|
| 186 | & 0.156911E-01_JPRB, 0.148937E-01_JPRB, 0.141368E-01_JPRB, 0.134184E-01_JPRB, 0.127365E-01_JPRB /) |
|---|
| 187 | SELFREF(:,12) = (/ & |
|---|
| 188 | & 0.226784E-01_JPRB, 0.210210E-01_JPRB, 0.194848E-01_JPRB, 0.180608E-01_JPRB, 0.167409E-01_JPRB, & |
|---|
| 189 | & 0.155174E-01_JPRB, 0.143834E-01_JPRB, 0.133322E-01_JPRB, 0.123579E-01_JPRB, 0.114547E-01_JPRB /) |
|---|
| 190 | SELFREF(:,13) = (/ & |
|---|
| 191 | & 0.221773E-01_JPRB, 0.210306E-01_JPRB, 0.199433E-01_JPRB, 0.189122E-01_JPRB, 0.179344E-01_JPRB, & |
|---|
| 192 | & 0.170071E-01_JPRB, 0.161278E-01_JPRB, 0.152939E-01_JPRB, 0.145032E-01_JPRB, 0.137533E-01_JPRB /) |
|---|
| 193 | SELFREF(:,14) = (/ & |
|---|
| 194 | & 0.275920E-01_JPRB, 0.252595E-01_JPRB, 0.231241E-01_JPRB, 0.211693E-01_JPRB, 0.193797E-01_JPRB, & |
|---|
| 195 | & 0.177415E-01_JPRB, 0.162417E-01_JPRB, 0.148687E-01_JPRB, 0.136117E-01_JPRB, 0.124610E-01_JPRB /) |
|---|
| 196 | SELFREF(:,15) = (/ & |
|---|
| 197 | & 0.288687E-01_JPRB, 0.269968E-01_JPRB, 0.252462E-01_JPRB, 0.236092E-01_JPRB, 0.220783E-01_JPRB, & |
|---|
| 198 | & 0.206466E-01_JPRB, 0.193078E-01_JPRB, 0.180559E-01_JPRB, 0.168851E-01_JPRB, 0.157902E-01_JPRB /) |
|---|
| 199 | SELFREF(:,16) = (/ & |
|---|
| 200 | & 0.371842E-01_JPRB, 0.347595E-01_JPRB, 0.324929E-01_JPRB, 0.303741E-01_JPRB, 0.283934E-01_JPRB, & |
|---|
| 201 | & 0.265419E-01_JPRB, 0.248112E-01_JPRB, 0.231933E-01_JPRB, 0.216809E-01_JPRB, 0.202671E-01_JPRB /) |
|---|
| 202 | |
|---|
| 203 | IF (LHOOK) CALL DR_HOOK('SRTM_KGB18',1,ZHOOK_HANDLE) |
|---|
| 204 | RETURN |
|---|
| 205 | |
|---|
| 206 | 1001 CONTINUE |
|---|
| 207 | CALL ABOR1("SRTM_KGB18:ERROR READING FILE RADSRTM") |
|---|
| 208 | |
|---|
| 209 | END SUBROUTINE SRTM_KGB18 |
|---|
