[3880] | 1 | SUBROUTINE SUSRTM |
---|
| 2 | |
---|
| 3 | ! Adapted from E.J. Mlawer, J. Delamere, Atmospheric & Environmental Research. |
---|
| 4 | ! by JJMorcrette, ECMWF |
---|
| 5 | ! Modified to add arrays relevant to mapping for g-point reduction, |
---|
| 6 | ! M.J. Iacono, Atmospheric & Environmental Research, Inc. |
---|
| 7 | ! JJMorcrette 20010610 Flexible configuration for number of g-points |
---|
| 8 | ! ------------------------------------------------------------------ |
---|
| 9 | |
---|
| 10 | USE PARKIND1 ,ONLY : JPRB , JPIM |
---|
| 11 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
---|
| 12 | |
---|
| 13 | USE YOESRTM , ONLY : JPGPT, NGBSW, NGN |
---|
| 14 | USE YOESRTWN , ONLY : NG , NSPA, NSPB , NMPSRTM, & |
---|
| 15 | & PREF , PREFLOG , TREF , & |
---|
| 16 | & NGM , WT , NGC , NGS |
---|
| 17 | ! & NGM , WT , NGC , NGS , NGN , NGBSW |
---|
| 18 | ! & WAVENUM1, WAVENUM2, DELWAVE, PREF, PREFLOG, TREF , & |
---|
| 19 | |
---|
| 20 | ! ------------------------------------------------------------------ |
---|
| 21 | |
---|
| 22 | IMPLICIT NONE |
---|
| 23 | |
---|
| 24 | INTEGER(KIND=JPIM) :: IGC56(14), IGC112(14) , IGC224(14) |
---|
| 25 | INTEGER(KIND=JPIM) :: IGS56(14), IGS112(14) , IGS224(14) |
---|
| 26 | |
---|
| 27 | INTEGER(KIND=JPIM) :: IGM56(224),IGM112(224), IGM224(224) |
---|
| 28 | |
---|
| 29 | INTEGER(KIND=JPIM) :: IGN56(56), IGN112(112), IGN224(224) |
---|
| 30 | INTEGER(KIND=JPIM) :: IGB56(56), IGB112(112), IGB224(224) |
---|
| 31 | |
---|
| 32 | REAL(KIND=JPRB) :: ZHOOK_HANDLE |
---|
| 33 | !----------------------------------------------------------------------- |
---|
| 34 | IF (LHOOK) CALL DR_HOOK('SUSRTM',0,ZHOOK_HANDLE) |
---|
| 35 | |
---|
| 36 | NG(:) =(/ 16,16,16,16,16,16,16,16,16,16,16,16,16,16 /) |
---|
| 37 | NSPA(:) =(/ 9, 9, 9, 9, 1, 9, 9, 1, 9, 1, 0, 1, 9, 1 /) |
---|
| 38 | NSPB(:) =(/ 1, 5, 1, 1, 1, 5, 1, 0, 1, 0, 0, 1, 5, 1 /) |
---|
| 39 | NMPSRTM(:)=(/ 6, 6, 5, 5, 5, 5, 5, 4, 4, 3, 2, 2, 1, 6 /) |
---|
| 40 | |
---|
| 41 | !WAVENUM1( :) = (/& |
---|
| 42 | ! & 2600._JPRB, 3250._JPRB, 4000._JPRB, 4650._JPRB, 5150._JPRB, 6150._JPRB, 7700._JPRB & |
---|
| 43 | ! & , 8050._JPRB,12850._JPRB,16000._JPRB,22650._JPRB,29000._JPRB,38000._JPRB, 820._JPRB /) |
---|
| 44 | !WAVENUM2( :) = (/& |
---|
| 45 | ! & 3250._JPRB, 4000._JPRB, 4650._JPRB, 5150._JPRB, 6150._JPRB, 7700._JPRB, 8050._JPRB & |
---|
| 46 | ! & ,12850._JPRB,16000._JPRB,22650._JPRB,29000._JPRB,38000._JPRB,50000._JPRB, 2600._JPRB /) |
---|
| 47 | !DELWAVE( :) = (/& |
---|
| 48 | ! & 650._JPRB, 750._JPRB, 650._JPRB, 500._JPRB, 1000._JPRB, 1550._JPRB, 350._JPRB & |
---|
| 49 | ! & , 4800._JPRB, 3150._JPRB, 6650._JPRB, 6350._JPRB, 9000._JPRB,12000._JPRB, 1780._JPRB /) |
---|
| 50 | |
---|
| 51 | !===================================================================== |
---|
| 52 | ! Set arrays needed for the g-point reduction from 224 to |
---|
| 53 | ! - either 112 for the high-resolution forecast model configuration |
---|
| 54 | ! - or 56 for the EPS-type configuration |
---|
| 55 | ! in the 14 SW bands: |
---|
| 56 | |
---|
| 57 | ! NB: This mapping from 224 to 112 points has been carefully selected to |
---|
| 58 | ! minimize the effect on the resulting fluxes and cooling rates, and |
---|
| 59 | ! caution should be used if the mapping is modified. |
---|
| 60 | ! The further reduction to 56 for EPS configuration is considered |
---|
| 61 | ! acceptable, only because of the random perturbations introduced on |
---|
| 62 | ! the total heating rates produced by the physical parametrization package. |
---|
| 63 | ! While a reduction to 56 obviously speeds up the model, it as obviously |
---|
| 64 | ! reduces the accuracy that could be expected from the radiation scheme. |
---|
| 65 | |
---|
| 66 | ! JPGPT The total number of new g-points (NGPT) |
---|
| 67 | ! NGC The number of new g-points in each band (14) |
---|
| 68 | ! NGS The cumulative sum of new g-points for each band (14) |
---|
| 69 | ! NGM The index of each new g-point relative to the original |
---|
| 70 | ! 16 g-points for each band. |
---|
| 71 | ! NGN The number of original g-points that are combined to make |
---|
| 72 | ! each new g-point in each band. |
---|
| 73 | ! NGB The band index for each new g-point. |
---|
| 74 | ! WT RRTM weights for 16 g-points. (16) |
---|
| 75 | |
---|
| 76 | !-- ECMWF EPS model RRTM_SW configuration with 56 g-points |
---|
| 77 | IGC56(:) = (/ 3, 6, 4, 4, 5, 5, 1, 5, 4, 3, 3, 4, 3, 6 /) |
---|
| 78 | IGS56(:) = (/ 3, 9,13,17,22,27,28,33,37,40,43,47, 50, 56 /) |
---|
| 79 | |
---|
| 80 | IGM56(:) = (/ 1,1,1,1,2,2,2,2,3,3,3,3,3,3,3,3, & ! Band 16 |
---|
| 81 | & 1,1,2,2,3,3,3,4,4,4,5,5,5,6,6,6, & ! Band 17 |
---|
| 82 | & 1,1,2,2,3,3,3,3,4,4,4,4,4,4,4,4, & ! Band 18 |
---|
| 83 | & 1,1,2,2,3,3,3,3,4,4,4,4,4,4,4,4, & ! Band 19 |
---|
| 84 | & 1,1,2,2,3,3,4,4,5,5,5,5,5,5,5,5, & ! Band 20 |
---|
| 85 | & 1,1,2,2,3,3,4,4,5,5,5,5,5,5,5,5, & ! Band 21 |
---|
| 86 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 22 |
---|
| 87 | & 1,1,1,1,2,2,3,3,4,4,5,5,5,5,5,5, & ! Band 23 |
---|
| 88 | & 1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4, & ! Band 24 |
---|
| 89 | & 1,1,2,2,2,2,3,3,3,3,3,3,3,3,3,3, & ! Band 25 |
---|
| 90 | & 1,1,2,2,2,2,3,3,3,3,3,3,3,3,3,3, & ! Band 26 |
---|
| 91 | & 1,1,2,2,3,3,4,4,4,4,4,4,4,4,4,4, & ! Band 27 |
---|
| 92 | & 1,1,2,2,2,2,3,3,3,3,3,3,3,3,3,3, & ! Band 28 |
---|
| 93 | & 1,1,2,2,3,3,3,3,4,4,4,4,5,5,6,6 /) ! Band 29 |
---|
| 94 | |
---|
| 95 | IGN56(:) = (/ 4,4,8, & ! Band 16 |
---|
| 96 | & 2,2,3,3,3,3, & ! Band 17 |
---|
| 97 | & 2,2,4,8, & ! Band 18 |
---|
| 98 | & 2,2,4,8, & ! Band 19 |
---|
| 99 | & 2,2,2,2,8, & ! Band 20 |
---|
| 100 | & 2,2,2,2,8, & ! Band 21 |
---|
| 101 | & 16, & ! Band 22 |
---|
| 102 | & 4,2,2,2,6, & ! Band 23 |
---|
| 103 | & 4,4,4,4, & ! Band 24 |
---|
| 104 | & 2,4,10, & ! Band 25 |
---|
| 105 | & 2,4,10, & ! Band 26 |
---|
| 106 | & 2,2,2,10, & ! Band 27 |
---|
| 107 | & 2,4,10, & ! Band 28 |
---|
| 108 | & 2,2,4,4,2,2 /) ! Band 29 |
---|
| 109 | |
---|
| 110 | IGB56(:) = (/ 16,16,16, & ! Band 16 |
---|
| 111 | & 17,17,17,17,17,17, & ! Band 17 |
---|
| 112 | & 18,18,18,18, & ! Band 18 |
---|
| 113 | & 19,19,19,19, & ! Band 19 |
---|
| 114 | & 20,20,20,20,20, & ! Band 20 |
---|
| 115 | & 21,21,21,21,21, & ! Band 21 |
---|
| 116 | & 22, & ! Band 22 |
---|
| 117 | & 23,23,23,23,23, & ! Band 23 |
---|
| 118 | & 24,24,24,24, & ! Band 24 |
---|
| 119 | & 25,25,25, & ! Band 25 |
---|
| 120 | & 26,26,26, & ! Band 26 |
---|
| 121 | & 27,27,27,27, & ! Band 27 |
---|
| 122 | & 28,28,28, & ! Band 28 |
---|
| 123 | & 29,29,29,29,29,29 /) ! Band 29 |
---|
| 124 | |
---|
| 125 | !------------------------------------------------------------------------------- |
---|
| 126 | !-- ECMWF high-resolution model RRTM_SW configuration with 112 g-points |
---|
| 127 | ! Use this NGC, NGS, NGM, and NGN for reduced (112) g-point set |
---|
| 128 | ! (A related code change is required in modules parsrtm.F90 and yoesrtwn.F90) |
---|
| 129 | |
---|
| 130 | IGC112(:) = (/ 6,12, 8, 8,10,10, 2,10, 8, 6, 6, 8, 6,12 /) |
---|
| 131 | IGS112(:) = (/ 6,18,26,34,44,54,56,66,74,80,86,94,100,112 /) |
---|
| 132 | |
---|
| 133 | !NGM(:) |
---|
| 134 | IGM112(:) = (/ 1,1,2,2,3,3,4,4,5,5,5,5,6,6,6,6, & ! Band 16 |
---|
| 135 | & 1,2,3,4,5,6,6,7,8,8,9,10,10,11,12,12, & ! Band 17 |
---|
| 136 | & 1,2,3,4,5,5,6,6,7,7,7,7,8,8,8,8, & ! Band 18 |
---|
| 137 | & 1,2,3,4,5,5,6,6,7,7,7,7,8,8,8,8, & ! Band 19 |
---|
| 138 | & 1,2,3,4,5,6,7,8,9,9,10,10,10,10,10,10, & ! Band 20 |
---|
| 139 | & 1,2,3,4,5,6,7,8,9,9,10,10,10,10,10,10, & ! Band 21 |
---|
| 140 | & 1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2, & ! Band 22 |
---|
| 141 | & 1,1,2,2,3,4,5,6,7,8,9,9,10,10,10,10, & ! Band 23 |
---|
| 142 | & 1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8, & ! Band 24 |
---|
| 143 | & 1,2,3,3,4,4,5,5,5,5,6,6,6,6,6,6, & ! Band 25 |
---|
| 144 | & 1,2,3,3,4,4,5,5,5,5,6,6,6,6,6,6, & ! Band 26 |
---|
| 145 | & 1,2,3,4,5,6,7,7,7,7,8,8,8,8,8,8, & ! Band 27 |
---|
| 146 | & 1,2,3,3,4,4,5,5,5,5,6,6,6,6,6,6, & ! Band 28 |
---|
| 147 | & 1,2,3,4,5,5,6,6,7,7,8,8,9,10,11,12 /) ! Band 29 |
---|
| 148 | !NGN(:) |
---|
| 149 | IGN112(:) = (/ 2,2,2,2,4,4, & ! Band 16 |
---|
| 150 | & 1,1,1,1,1,2,1,2,1,2,1,2, & ! Band 17 |
---|
| 151 | & 1,1,1,1,2,2,4,4, & ! Band 18 |
---|
| 152 | & 1,1,1,1,2,2,4,4, & ! Band 19 |
---|
| 153 | & 1,1,1,1,1,1,1,1,2,6, & ! Band 20 |
---|
| 154 | & 1,1,1,1,1,1,1,1,2,6, & ! Band 21 |
---|
| 155 | & 8,8, & ! Band 22 |
---|
| 156 | & 2,2,1,1,1,1,1,1,2,4, & ! Band 23 |
---|
| 157 | & 2,2,2,2,2,2,2,2, & ! Band 24 |
---|
| 158 | & 1,1,2,2,4,6, & ! Band 25 |
---|
| 159 | & 1,1,2,2,4,6, & ! Band 26 |
---|
| 160 | & 1,1,1,1,1,1,4,6, & ! Band 27 |
---|
| 161 | & 1,1,2,2,4,6, & ! Band 28 |
---|
| 162 | & 1,1,1,1,2,2,2,2,1,1,1,1 /) ! Band 29 |
---|
| 163 | !NGBSW(:) |
---|
| 164 | IGB112(:) = (/ 16,16,16,16,16,16, & ! Band 16 |
---|
| 165 | & 17,17,17,17,17,17,17,17,17,17,17,17, & ! Band 17 |
---|
| 166 | & 18,18,18,18,18,18,18,18, & ! Band 18 |
---|
| 167 | & 19,19,19,19,19,19,19,19, & ! Band 19 |
---|
| 168 | & 20,20,20,20,20,20,20,20,20,20, & ! Band 20 |
---|
| 169 | & 21,21,21,21,21,21,21,21,21,21, & ! Band 21 |
---|
| 170 | & 22,22, & ! Band 22 |
---|
| 171 | & 23,23,23,23,23,23,23,23,23,23, & ! Band 23 |
---|
| 172 | & 24,24,24,24,24,24,24,24, & ! Band 24 |
---|
| 173 | & 25,25,25,25,25,25, & ! Band 25 |
---|
| 174 | & 26,26,26,26,26,26, & ! Band 26 |
---|
| 175 | & 27,27,27,27,27,27,27,27, & ! Band 27 |
---|
| 176 | & 28,28,28,28,28,28, & ! Band 28 |
---|
| 177 | & 29,29,29,29,29,29,29,29,29,29,29,29 /) ! Band 29 |
---|
| 178 | |
---|
| 179 | !------------------------------------------------------------------------------- |
---|
| 180 | !-- original RRTM_SW configuration with 224 (14*16 g-points) |
---|
| 181 | ! Use this NGC, NGS, NGM, and NGN for full (224) g-point set |
---|
| 182 | ! (A related code change is required in modules parsrtm.F90 and yoesrtwn.F90) |
---|
| 183 | IGC224(:) = (/ 16,16,16,16,16,16,16,16,16,16,16,16,16,16 /) |
---|
| 184 | IGS224(:) = (/ 16,32,48,64,80,96,112,128,144,160,176,192,208,224 /) |
---|
| 185 | |
---|
| 186 | IGM224(:) = (/ 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 16 |
---|
| 187 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 17 |
---|
| 188 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 18 |
---|
| 189 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 19 |
---|
| 190 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 20 |
---|
| 191 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 21 |
---|
| 192 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 22 |
---|
| 193 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 23 |
---|
| 194 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 24 |
---|
| 195 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 25 |
---|
| 196 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 26 |
---|
| 197 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 27 |
---|
| 198 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! Band 28 |
---|
| 199 | & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16 /) ! Band 29 |
---|
| 200 | |
---|
| 201 | IGN224(:) = (/ 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 16 |
---|
| 202 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 17 |
---|
| 203 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 18 |
---|
| 204 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 19 |
---|
| 205 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 20 |
---|
| 206 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 21 |
---|
| 207 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 22 |
---|
| 208 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 23 |
---|
| 209 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 24 |
---|
| 210 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 25 |
---|
| 211 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 26 |
---|
| 212 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 27 |
---|
| 213 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! Band 28 |
---|
| 214 | & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 /) ! Band 29 |
---|
| 215 | |
---|
| 216 | IGB224(:) = (/ 16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16, & ! Band 16 |
---|
| 217 | & 17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17, & ! Band 17 |
---|
| 218 | & 18,18,18,18,18,18,18,18,18,18,18,18,18,18,18,18, & ! Band 18 |
---|
| 219 | & 19,19,19,19,19,19,19,19,19,19,19,19,19,19,19,19, & ! Band 19 |
---|
| 220 | & 20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20, & ! Band 20 |
---|
| 221 | & 21,21,21,21,21,21,21,21,21,21,21,21,21,21,21,21, & ! Band 21 |
---|
| 222 | & 22,22,22,22,22,22,22,22,22,22,22,22,22,22,22,22, & ! Band 22 |
---|
| 223 | & 23,23,23,23,23,23,23,23,23,23,23,23,23,23,23,23, & ! Band 23 |
---|
| 224 | & 24,24,24,24,24,24,24,24,24,24,24,24,24,24,24,24, & ! Band 24 |
---|
| 225 | & 25,25,25,25,25,25,25,25,25,25,25,25,25,25,25,25, & ! Band 25 |
---|
| 226 | & 26,26,26,26,26,26,26,26,26,26,26,26,26,26,26,26, & ! Band 26 |
---|
| 227 | & 27,27,27,27,27,27,27,27,27,27,27,27,27,27,27,27, & ! Band 27 |
---|
| 228 | & 28,28,28,28,28,28,28,28,28,28,28,28,28,28,28,28, & ! Band 28 |
---|
| 229 | & 29,29,29,29,29,29,29,29,29,29,29,29,29,29,29,28 /) ! Band 29 |
---|
| 230 | |
---|
| 231 | !============================================================================= |
---|
| 232 | |
---|
| 233 | WT(:) = (/ 0.1527534276_JPRB, 0.1491729617_JPRB, 0.1420961469_JPRB, & |
---|
| 234 | & 0.1316886544_JPRB, 0.1181945205_JPRB, 0.1019300893_JPRB, & |
---|
| 235 | & 0.0832767040_JPRB, 0.0626720116_JPRB, 0.0424925000_JPRB, & |
---|
| 236 | & 0.0046269894_JPRB, 0.0038279891_JPRB, 0.0030260086_JPRB, & |
---|
| 237 | & 0.0022199750_JPRB, 0.0014140010_JPRB, 0.0005330000_JPRB, & |
---|
| 238 | & 0.0000750000_JPRB /) |
---|
| 239 | |
---|
| 240 | !============================================================================= |
---|
| 241 | |
---|
| 242 | ! These pressures are chosen such that the ln of the first pressure |
---|
| 243 | ! has only a few non-zero digits (i.e. ln(PREF(1)) = 6.96000) and |
---|
| 244 | ! each subsequent ln(pressure) differs from the previous one by 0.2. |
---|
| 245 | PREF = (/ & |
---|
| 246 | & 1.05363E+03_JPRB,8.62642E+02_JPRB,7.06272E+02_JPRB,5.78246E+02_JPRB,4.73428E+02_JPRB, & |
---|
| 247 | & 3.87610E+02_JPRB,3.17348E+02_JPRB,2.59823E+02_JPRB,2.12725E+02_JPRB,1.74164E+02_JPRB, & |
---|
| 248 | & 1.42594E+02_JPRB,1.16746E+02_JPRB,9.55835E+01_JPRB,7.82571E+01_JPRB,6.40715E+01_JPRB, & |
---|
| 249 | & 5.24573E+01_JPRB,4.29484E+01_JPRB,3.51632E+01_JPRB,2.87892E+01_JPRB,2.35706E+01_JPRB, & |
---|
| 250 | & 1.92980E+01_JPRB,1.57998E+01_JPRB,1.29358E+01_JPRB,1.05910E+01_JPRB,8.67114E+00_JPRB, & |
---|
| 251 | & 7.09933E+00_JPRB,5.81244E+00_JPRB,4.75882E+00_JPRB,3.89619E+00_JPRB,3.18993E+00_JPRB, & |
---|
| 252 | & 2.61170E+00_JPRB,2.13828E+00_JPRB,1.75067E+00_JPRB,1.43333E+00_JPRB,1.17351E+00_JPRB, & |
---|
| 253 | & 9.60789E-01_JPRB,7.86628E-01_JPRB,6.44036E-01_JPRB,5.27292E-01_JPRB,4.31710E-01_JPRB, & |
---|
| 254 | & 3.53455E-01_JPRB,2.89384E-01_JPRB,2.36928E-01_JPRB,1.93980E-01_JPRB,1.58817E-01_JPRB, & |
---|
| 255 | & 1.30029E-01_JPRB,1.06458E-01_JPRB,8.71608E-02_JPRB,7.13612E-02_JPRB,5.84256E-02_JPRB, & |
---|
| 256 | & 4.78349E-02_JPRB,3.91639E-02_JPRB,3.20647E-02_JPRB,2.62523E-02_JPRB,2.14936E-02_JPRB, & |
---|
| 257 | & 1.75975E-02_JPRB,1.44076E-02_JPRB,1.17959E-02_JPRB,9.65769E-03_JPRB /) |
---|
| 258 | PREFLOG = (/ & |
---|
| 259 | & 6.9600E+00_JPRB, 6.7600E+00_JPRB, 6.5600E+00_JPRB, 6.3600E+00_JPRB, 6.1600E+00_JPRB, & |
---|
| 260 | & 5.9600E+00_JPRB, 5.7600E+00_JPRB, 5.5600E+00_JPRB, 5.3600E+00_JPRB, 5.1600E+00_JPRB, & |
---|
| 261 | & 4.9600E+00_JPRB, 4.7600E+00_JPRB, 4.5600E+00_JPRB, 4.3600E+00_JPRB, 4.1600E+00_JPRB, & |
---|
| 262 | & 3.9600E+00_JPRB, 3.7600E+00_JPRB, 3.5600E+00_JPRB, 3.3600E+00_JPRB, 3.1600E+00_JPRB, & |
---|
| 263 | & 2.9600E+00_JPRB, 2.7600E+00_JPRB, 2.5600E+00_JPRB, 2.3600E+00_JPRB, 2.1600E+00_JPRB, & |
---|
| 264 | & 1.9600E+00_JPRB, 1.7600E+00_JPRB, 1.5600E+00_JPRB, 1.3600E+00_JPRB, 1.1600E+00_JPRB, & |
---|
| 265 | & 9.6000E-01_JPRB, 7.6000E-01_JPRB, 5.6000E-01_JPRB, 3.6000E-01_JPRB, 1.6000E-01_JPRB, & |
---|
| 266 | & -4.0000E-02_JPRB,-2.4000E-01_JPRB,-4.4000E-01_JPRB,-6.4000E-01_JPRB,-8.4000E-01_JPRB, & |
---|
| 267 | & -1.0400E+00_JPRB,-1.2400E+00_JPRB,-1.4400E+00_JPRB,-1.6400E+00_JPRB,-1.8400E+00_JPRB, & |
---|
| 268 | & -2.0400E+00_JPRB,-2.2400E+00_JPRB,-2.4400E+00_JPRB,-2.6400E+00_JPRB,-2.8400E+00_JPRB, & |
---|
| 269 | & -3.0400E+00_JPRB,-3.2400E+00_JPRB,-3.4400E+00_JPRB,-3.6400E+00_JPRB,-3.8400E+00_JPRB, & |
---|
| 270 | & -4.0400E+00_JPRB,-4.2400E+00_JPRB,-4.4400E+00_JPRB,-4.6400E+00_JPRB /) |
---|
| 271 | ! These are the temperatures associated with the respective |
---|
| 272 | ! pressures for the MLS standard atmosphere. |
---|
| 273 | TREF = (/ & |
---|
| 274 | & 2.9420E+02_JPRB, 2.8799E+02_JPRB, 2.7894E+02_JPRB, 2.6925E+02_JPRB, 2.5983E+02_JPRB, & |
---|
| 275 | & 2.5017E+02_JPRB, 2.4077E+02_JPRB, 2.3179E+02_JPRB, 2.2306E+02_JPRB, 2.1578E+02_JPRB, & |
---|
| 276 | & 2.1570E+02_JPRB, 2.1570E+02_JPRB, 2.1570E+02_JPRB, 2.1706E+02_JPRB, 2.1858E+02_JPRB, & |
---|
| 277 | & 2.2018E+02_JPRB, 2.2174E+02_JPRB, 2.2328E+02_JPRB, 2.2479E+02_JPRB, 2.2655E+02_JPRB, & |
---|
| 278 | & 2.2834E+02_JPRB, 2.3113E+02_JPRB, 2.3401E+02_JPRB, 2.3703E+02_JPRB, 2.4022E+02_JPRB, & |
---|
| 279 | & 2.4371E+02_JPRB, 2.4726E+02_JPRB, 2.5085E+02_JPRB, 2.5457E+02_JPRB, 2.5832E+02_JPRB, & |
---|
| 280 | & 2.6216E+02_JPRB, 2.6606E+02_JPRB, 2.6999E+02_JPRB, 2.7340E+02_JPRB, 2.7536E+02_JPRB, & |
---|
| 281 | & 2.7568E+02_JPRB, 2.7372E+02_JPRB, 2.7163E+02_JPRB, 2.6955E+02_JPRB, 2.6593E+02_JPRB, & |
---|
| 282 | & 2.6211E+02_JPRB, 2.5828E+02_JPRB, 2.5360E+02_JPRB, 2.4854E+02_JPRB, 2.4348E+02_JPRB, & |
---|
| 283 | & 2.3809E+02_JPRB, 2.3206E+02_JPRB, 2.2603E+02_JPRB, 2.2000E+02_JPRB, 2.1435E+02_JPRB, & |
---|
| 284 | & 2.0887E+02_JPRB, 2.0340E+02_JPRB, 1.9792E+02_JPRB, 1.9290E+02_JPRB, 1.8809E+02_JPRB, & |
---|
| 285 | & 1.8329E+02_JPRB, 1.7849E+02_JPRB, 1.7394E+02_JPRB, 1.7212E+02_JPRB /) |
---|
| 286 | ! ----------------------------------------------------------------- |
---|
| 287 | |
---|
| 288 | IF (JPGPT == 56) THEN |
---|
| 289 | |
---|
| 290 | !- 14 |
---|
| 291 | NGC(:)=IGC56(:) |
---|
| 292 | NGS(:)=IGS56(:) |
---|
| 293 | !- 14*16=224 |
---|
| 294 | NGM(:)=IGM56(:) |
---|
| 295 | |
---|
| 296 | NGN(1:56)=IGN56(1:56) |
---|
| 297 | NGBSW(1:56)=IGB56(1:56) |
---|
| 298 | |
---|
| 299 | ELSEIF (JPGPT == 112) THEN |
---|
| 300 | !- 14 |
---|
| 301 | NGC(:)=IGC112(:) |
---|
| 302 | NGS(:)=IGS112(:) |
---|
| 303 | !- 14*16=224 |
---|
| 304 | NGM(:)=IGM112(:) |
---|
| 305 | |
---|
| 306 | NGN(1:112)=IGN112(1:112) |
---|
| 307 | NGBSW(1:112)=IGB112(1:112) |
---|
| 308 | |
---|
| 309 | ELSEIF (JPGPT == 224) THEN |
---|
| 310 | !- 14 |
---|
| 311 | NGC(:)=IGC224(:) |
---|
| 312 | NGS(:)=IGS224(:) |
---|
| 313 | !- 14*16=224 |
---|
| 314 | NGM(:)=IGM224(:) |
---|
| 315 | |
---|
| 316 | NGN(1:224)=IGN224(1:224) |
---|
| 317 | NGBSW(1:224)=IGB224(1:224) |
---|
| 318 | |
---|
| 319 | ENDIF |
---|
| 320 | |
---|
| 321 | ! ----------------------------------------------------------------- |
---|
| 322 | IF (LHOOK) CALL DR_HOOK('SUSRTM',1,ZHOOK_HANDLE) |
---|
| 323 | END SUBROUTINE SUSRTM |
---|
| 324 | |
---|