| 1 | SUBROUTINE SRTM_KGB24 |
|---|
| 2 | |
|---|
| 3 | ! Originally by J.Delamere, Atmospheric & Environmental Research. |
|---|
| 4 | ! Revision: 2.4 |
|---|
| 5 | ! BAND 24: 12850-16000 cm-1 (low - H2O,O2; high - O2) |
|---|
| 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 YOESRTA24 , ONLY : KA, KB, SELFREF, FORREF, SFLUXREF, RAYLA, RAYLB, & |
|---|
| 20 | & ABSO3A, ABSO3B, STRRAT, LAYREFFR , 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_KGB24',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_KGB24:') |
|---|
| 40 | CALL MPL_BROADCAST (KB,MTAGRAD,1,CDSTRING='SRTM_KGB24:') |
|---|
| 41 | ENDIF |
|---|
| 42 | |
|---|
| 43 | SFLUXREF(:,1) = (/ & |
|---|
| 44 | & 34.3610_JPRB , 33.1240_JPRB, 31.3948_JPRB, 28.7248_JPRB, & |
|---|
| 45 | & 24.7884_JPRB , 21.4892_JPRB, 17.3972_JPRB, 13.7928_JPRB, & |
|---|
| 46 | & 9.54462_JPRB , 1.05002_JPRB,0.867332_JPRB,0.685753_JPRB, & |
|---|
| 47 | & 0.504718_JPRB,0.323112_JPRB,0.122183_JPRB, 1.70288E-02_JPRB /) |
|---|
| 48 | SFLUXREF(:,2) = (/ & |
|---|
| 49 | & 34.2367_JPRB , 32.4327_JPRB, 30.0863_JPRB, 28.2085_JPRB, & |
|---|
| 50 | & 25.6533_JPRB , 22.3412_JPRB, 18.3112_JPRB, 13.8521_JPRB, & |
|---|
| 51 | & 9.51035_JPRB , 1.04138_JPRB,0.863493_JPRB,0.682790_JPRB, & |
|---|
| 52 | & 0.504721_JPRB,0.323102_JPRB,0.122193_JPRB, 1.70288E-02_JPRB /) |
|---|
| 53 | SFLUXREF(:,3) = (/ & |
|---|
| 54 | & 34.1883_JPRB , 32.2479_JPRB, 30.2650_JPRB, 28.2914_JPRB, & |
|---|
| 55 | & 25.6626_JPRB , 22.3163_JPRB, 18.3327_JPRB, 13.8508_JPRB, & |
|---|
| 56 | & 9.49190_JPRB , 1.03672_JPRB,0.858272_JPRB,0.681485_JPRB, & |
|---|
| 57 | & 0.501363_JPRB,0.323110_JPRB,0.122183_JPRB, 1.70288E-02_JPRB /) |
|---|
| 58 | SFLUXREF(:,4) = (/ & |
|---|
| 59 | & 34.1365_JPRB , 32.2316_JPRB, 30.3325_JPRB, 28.3305_JPRB, & |
|---|
| 60 | & 25.6420_JPRB , 22.3223_JPRB, 18.3411_JPRB, 13.8471_JPRB, & |
|---|
| 61 | & 9.47492_JPRB , 1.03376_JPRB,0.855380_JPRB,0.679085_JPRB, & |
|---|
| 62 | & 0.497998_JPRB,0.323053_JPRB,0.122183_JPRB, 1.70288E-02_JPRB /) |
|---|
| 63 | SFLUXREF(:,5) = (/ & |
|---|
| 64 | & 34.0460_JPRB , 32.2795_JPRB, 30.4147_JPRB, 28.3123_JPRB, & |
|---|
| 65 | & 25.6438_JPRB , 22.3238_JPRB, 18.3441_JPRB, 13.8528_JPRB, & |
|---|
| 66 | & 9.45222_JPRB , 1.03058_JPRB,0.854037_JPRB,0.675554_JPRB, & |
|---|
| 67 | & 0.498344_JPRB,0.320072_JPRB,0.122193_JPRB, 1.70288E-02_JPRB /) |
|---|
| 68 | SFLUXREF(:,6) = (/ & |
|---|
| 69 | & 33.9909_JPRB , 32.3127_JPRB, 30.4854_JPRB, 28.3005_JPRB, & |
|---|
| 70 | & 25.6310_JPRB , 22.3294_JPRB, 18.3459_JPRB, 13.8488_JPRB, & |
|---|
| 71 | & 9.43336_JPRB , 1.02901_JPRB,0.852728_JPRB,0.672322_JPRB, & |
|---|
| 72 | & 0.498056_JPRB,0.317753_JPRB,0.122183_JPRB, 1.70288E-02_JPRB /) |
|---|
| 73 | SFLUXREF(:,7) = (/ & |
|---|
| 74 | & 33.9225_JPRB , 32.4097_JPRB, 30.5125_JPRB, 28.2810_JPRB, & |
|---|
| 75 | & 25.6387_JPRB , 22.3080_JPRB, 18.3715_JPRB, 13.8248_JPRB, & |
|---|
| 76 | & 9.41834_JPRB , 1.02735_JPRB,0.850807_JPRB,0.671379_JPRB, & |
|---|
| 77 | & 0.496975_JPRB,0.317158_JPRB,0.119297_JPRB, 1.70207E-02_JPRB /) |
|---|
| 78 | SFLUXREF(:,8) = (/ & |
|---|
| 79 | & 33.8940_JPRB , 32.4951_JPRB, 30.5494_JPRB, 28.2788_JPRB, & |
|---|
| 80 | & 25.5975_JPRB , 22.3225_JPRB, 18.3358_JPRB, 13.8199_JPRB, & |
|---|
| 81 | & 9.40283_JPRB , 1.02751_JPRB,0.850729_JPRB,0.670152_JPRB, & |
|---|
| 82 | & 0.494294_JPRB,0.315829_JPRB,0.116195_JPRB, 1.64138E-02_JPRB /) |
|---|
| 83 | SFLUXREF(:,9) = (/ & |
|---|
| 84 | & 34.6501_JPRB , 32.6690_JPRB, 30.2872_JPRB, 28.0955_JPRB, & |
|---|
| 85 | & 25.4662_JPRB , 22.1446_JPRB, 18.2754_JPRB, 13.7573_JPRB, & |
|---|
| 86 | & 9.36645_JPRB , 1.02356_JPRB,0.847154_JPRB,0.668519_JPRB, & |
|---|
| 87 | & 0.489186_JPRB,0.313790_JPRB,0.117074_JPRB, 1.60943E-02_JPRB /) |
|---|
| 88 | |
|---|
| 89 | ! Rayleigh extinction coefficient at all v |
|---|
| 90 | RAYLA(:,1) = (/ & |
|---|
| 91 | & 1.28405E-07_JPRB,1.45501E-07_JPRB,1.67272E-07_JPRB,1.94856E-07_JPRB, & |
|---|
| 92 | & 2.15248E-07_JPRB,2.34920E-07_JPRB,2.48558E-07_JPRB,1.80004E-07_JPRB, & |
|---|
| 93 | & 1.46504E-07_JPRB,1.31355E-07_JPRB,1.33562E-07_JPRB,1.35618E-07_JPRB, & |
|---|
| 94 | & 1.22412E-07_JPRB,1.19842E-07_JPRB,1.19924E-07_JPRB,1.20264E-07_JPRB /) |
|---|
| 95 | RAYLA(:,2) = (/ & |
|---|
| 96 | & 1.41622E-07_JPRB,1.93436E-07_JPRB,2.25057E-07_JPRB,2.01025E-07_JPRB, & |
|---|
| 97 | & 1.85138E-07_JPRB,1.72672E-07_JPRB,1.64771E-07_JPRB,1.59312E-07_JPRB, & |
|---|
| 98 | & 1.44961E-07_JPRB,1.37448E-07_JPRB,1.37506E-07_JPRB,1.38081E-07_JPRB, & |
|---|
| 99 | & 1.22432E-07_JPRB,1.19844E-07_JPRB,1.19921E-07_JPRB,1.20287E-07_JPRB /) |
|---|
| 100 | RAYLA(:,3) = (/ & |
|---|
| 101 | & 1.45382E-07_JPRB,1.97020E-07_JPRB,2.22781E-07_JPRB,1.96062E-07_JPRB, & |
|---|
| 102 | & 1.83495E-07_JPRB,1.72495E-07_JPRB,1.64910E-07_JPRB,1.58797E-07_JPRB, & |
|---|
| 103 | & 1.46208E-07_JPRB,1.42274E-07_JPRB,1.40445E-07_JPRB,1.39496E-07_JPRB, & |
|---|
| 104 | & 1.26940E-07_JPRB,1.19844E-07_JPRB,1.19921E-07_JPRB,1.20287E-07_JPRB /) |
|---|
| 105 | RAYLA(:,4) = (/ & |
|---|
| 106 | & 1.48247E-07_JPRB,1.99958E-07_JPRB,2.18048E-07_JPRB,1.93896E-07_JPRB, & |
|---|
| 107 | & 1.83125E-07_JPRB,1.73244E-07_JPRB,1.64320E-07_JPRB,1.58298E-07_JPRB, & |
|---|
| 108 | & 1.48428E-07_JPRB,1.44769E-07_JPRB,1.43704E-07_JPRB,1.38498E-07_JPRB, & |
|---|
| 109 | & 1.31732E-07_JPRB,1.22299E-07_JPRB,1.19921E-07_JPRB,1.20287E-07_JPRB /) |
|---|
| 110 | RAYLA(:,5) = (/ & |
|---|
| 111 | & 1.51343E-07_JPRB,1.99621E-07_JPRB,2.14563E-07_JPRB,1.93824E-07_JPRB, & |
|---|
| 112 | & 1.82992E-07_JPRB,1.73143E-07_JPRB,1.64587E-07_JPRB,1.57355E-07_JPRB, & |
|---|
| 113 | & 1.51198E-07_JPRB,1.46373E-07_JPRB,1.45438E-07_JPRB,1.38095E-07_JPRB, & |
|---|
| 114 | & 1.35026E-07_JPRB,1.27504E-07_JPRB,1.19921E-07_JPRB,1.20287E-07_JPRB /) |
|---|
| 115 | RAYLA(:,6) = (/ & |
|---|
| 116 | & 1.54462E-07_JPRB,1.97610E-07_JPRB,2.11992E-07_JPRB,1.93831E-07_JPRB, & |
|---|
| 117 | & 1.83900E-07_JPRB,1.73125E-07_JPRB,1.64093E-07_JPRB,1.57651E-07_JPRB, & |
|---|
| 118 | & 1.53158E-07_JPRB,1.46843E-07_JPRB,1.44733E-07_JPRB,1.40611E-07_JPRB, & |
|---|
| 119 | & 1.37320E-07_JPRB,1.33932E-07_JPRB,1.20423E-07_JPRB,1.20287E-07_JPRB /) |
|---|
| 120 | RAYLA(:,7) = (/ & |
|---|
| 121 | & 1.59068E-07_JPRB,1.92757E-07_JPRB,2.09865E-07_JPRB,1.95132E-07_JPRB, & |
|---|
| 122 | & 1.83641E-07_JPRB,1.73778E-07_JPRB,1.63215E-07_JPRB,1.59462E-07_JPRB, & |
|---|
| 123 | & 1.54331E-07_JPRB,1.46177E-07_JPRB,1.45819E-07_JPRB,1.43177E-07_JPRB, & |
|---|
| 124 | & 1.39797E-07_JPRB,1.36780E-07_JPRB,1.33385E-07_JPRB,1.20287E-07_JPRB /) |
|---|
| 125 | RAYLA(:,8) = (/ & |
|---|
| 126 | & 1.62066E-07_JPRB,1.87529E-07_JPRB,2.07191E-07_JPRB,1.97788E-07_JPRB, & |
|---|
| 127 | & 1.84920E-07_JPRB,1.72951E-07_JPRB,1.65450E-07_JPRB,1.60344E-07_JPRB, & |
|---|
| 128 | & 1.54403E-07_JPRB,1.47679E-07_JPRB,1.47287E-07_JPRB,1.44951E-07_JPRB, & |
|---|
| 129 | & 1.42517E-07_JPRB,1.41107E-07_JPRB,1.48688E-07_JPRB,1.51127E-07_JPRB /) |
|---|
| 130 | RAYLA(:,9) = (/ & |
|---|
| 131 | & 1.19177E-07_JPRB,1.86522E-07_JPRB,2.20324E-07_JPRB,2.13543E-07_JPRB, & |
|---|
| 132 | & 1.92198E-07_JPRB,1.81641E-07_JPRB,1.70092E-07_JPRB,1.65072E-07_JPRB, & |
|---|
| 133 | & 1.59804E-07_JPRB,1.56745E-07_JPRB,1.51235E-07_JPRB,1.51400E-07_JPRB, & |
|---|
| 134 | & 1.49635E-07_JPRB,1.48056E-07_JPRB,1.49046E-07_JPRB,1.51010E-07_JPRB /) |
|---|
| 135 | |
|---|
| 136 | RAYLB = (/ & |
|---|
| 137 | & 1.23766E-07_JPRB,1.40524E-07_JPRB,1.61610E-07_JPRB,1.83232E-07_JPRB, & |
|---|
| 138 | & 2.02951E-07_JPRB,2.21367E-07_JPRB,2.38367E-07_JPRB,2.53019E-07_JPRB, & |
|---|
| 139 | & 2.12202E-07_JPRB,1.36977E-07_JPRB,1.39118E-07_JPRB,1.37097E-07_JPRB, & |
|---|
| 140 | & 1.33223E-07_JPRB,1.38695E-07_JPRB,1.19868E-07_JPRB,1.20062E-07_JPRB /) |
|---|
| 141 | |
|---|
| 142 | ABSO3A = (/ & |
|---|
| 143 | & 8.03067E-02_JPRB,0.180926_JPRB ,0.227484_JPRB ,0.168015_JPRB , & |
|---|
| 144 | & 0.138284_JPRB ,0.114537_JPRB ,9.50114E-02_JPRB,8.06816E-02_JPRB, & |
|---|
| 145 | & 6.76406E-02_JPRB,5.69802E-02_JPRB,5.63283E-02_JPRB,4.57592E-02_JPRB, & |
|---|
| 146 | & 4.21862E-02_JPRB,3.47949E-02_JPRB,2.65731E-02_JPRB,2.67628E-02_JPRB /) |
|---|
| 147 | |
|---|
| 148 | ABSO3B = (/ & |
|---|
| 149 | & 2.94848E-02_JPRB,4.33642E-02_JPRB,6.70197E-02_JPRB,0.104990_JPRB , & |
|---|
| 150 | & 0.156180_JPRB ,0.214638_JPRB ,0.266281_JPRB ,0.317941_JPRB , & |
|---|
| 151 | & 0.355327_JPRB ,0.371241_JPRB ,0.374396_JPRB ,0.326847_JPRB , & |
|---|
| 152 | & 0.126497_JPRB ,6.95264E-02_JPRB,2.58175E-02_JPRB,2.52862E-02_JPRB /) |
|---|
| 153 | |
|---|
| 154 | STRRAT = 0.124692_JPRB |
|---|
| 155 | |
|---|
| 156 | LAYREFFR = 1 |
|---|
| 157 | |
|---|
| 158 | ! ------------------------------------------------------------------ |
|---|
| 159 | |
|---|
| 160 | ! The array KA contains absorption coefs at the 16 chosen g-values |
|---|
| 161 | ! for a range of pressure levels> ~100mb, temperatures, and binary |
|---|
| 162 | ! species parameters (see taumol.f for definition). The first |
|---|
| 163 | ! index in the array, JS, runs from 1 to 9, and corresponds to |
|---|
| 164 | ! different values of the binary species parameter. For instance, |
|---|
| 165 | ! JS=1 refers to dry air, JS = 2 corresponds to the paramter value 1/8, |
|---|
| 166 | ! JS = 3 corresponds to the parameter value 2/8, etc. The second index |
|---|
| 167 | ! in the array, JT, which runs from 1 to 5, corresponds to different |
|---|
| 168 | ! temperatures. More specifically, JT = 3 means that the data are for |
|---|
| 169 | ! the reference temperature TREF for this pressure level, JT = 2 refers |
|---|
| 170 | ! to TREF-15, JT = 1 is for TREF-30, JT = 4 is for TREF+15, and JT = 5 |
|---|
| 171 | ! is for TREF+30. The third index, JP, runs from 1 to 13 and refers |
|---|
| 172 | ! to the JPth reference pressure level (see taumol.f for these levels |
|---|
| 173 | ! in mb). The fourth index, IG, goes from 1 to 16, and indicates |
|---|
| 174 | ! which g-interval the absorption coefficients are for. |
|---|
| 175 | ! ----------------------------------------------------------------- |
|---|
| 176 | |
|---|
| 177 | ! ----------------------------------------------------------------- |
|---|
| 178 | ! The array KB contains absorption coefs at the 16 chosen g-values |
|---|
| 179 | ! for a range of pressure levels < ~100mb and temperatures. The first |
|---|
| 180 | ! index in the array, JT, which runs from 1 to 5, corresponds to |
|---|
| 181 | ! different temperatures. More specifically, JT = 3 means that the |
|---|
| 182 | ! data are for the reference temperature TREF for this pressure |
|---|
| 183 | ! level, JT = 2 refers to the temperature TREF-15, JT = 1 is for |
|---|
| 184 | ! TREF-30, JT = 4 is for TREF+15, and JT = 5 is for TREF+30. |
|---|
| 185 | ! The second index, JP, runs from 13 to 59 and refers to the JPth |
|---|
| 186 | ! reference pressure level (see taumol.f for the value of these |
|---|
| 187 | ! pressure levels in mb). The third index, IG, goes from 1 to 16, |
|---|
| 188 | ! and tells us which g-interval the absorption coefficients are for. |
|---|
| 189 | ! ----------------------------------------------------------------- |
|---|
| 190 | |
|---|
| 191 | |
|---|
| 192 | FORREF(:, 1) = (/ 0.515619E-08_JPRB, 0.131078E-06_JPRB, 0.349038E-06_JPRB /) |
|---|
| 193 | FORREF(:, 2) = (/ 0.329605E-07_JPRB, 0.430497E-06_JPRB, 0.458569E-05_JPRB /) |
|---|
| 194 | FORREF(:, 3) = (/ 0.188244E-06_JPRB, 0.792931E-06_JPRB, 0.267176E-05_JPRB /) |
|---|
| 195 | FORREF(:, 4) = (/ 0.611237E-06_JPRB, 0.798868E-06_JPRB, 0.411583E-06_JPRB /) |
|---|
| 196 | FORREF(:, 5) = (/ 0.111903E-05_JPRB, 0.914895E-06_JPRB, 0.444828E-06_JPRB /) |
|---|
| 197 | FORREF(:, 6) = (/ 0.235399E-05_JPRB, 0.269099E-05_JPRB, 0.739855E-06_JPRB /) |
|---|
| 198 | FORREF(:, 7) = (/ 0.400131E-05_JPRB, 0.378135E-05_JPRB, 0.231265E-06_JPRB /) |
|---|
| 199 | FORREF(:, 8) = (/ 0.464257E-05_JPRB, 0.371927E-05_JPRB, 0.460611E-06_JPRB /) |
|---|
| 200 | FORREF(:, 9) = (/ 0.476792E-05_JPRB, 0.311841E-05_JPRB, 0.934811E-06_JPRB /) |
|---|
| 201 | FORREF(:,10) = (/ 0.555683E-05_JPRB, 0.238129E-05_JPRB, 0.400334E-07_JPRB /) |
|---|
| 202 | FORREF(:,11) = (/ 0.569068E-05_JPRB, 0.196039E-05_JPRB, 0.374476E-07_JPRB /) |
|---|
| 203 | FORREF(:,12) = (/ 0.554154E-05_JPRB, 0.131724E-05_JPRB, 0.399720E-07_JPRB /) |
|---|
| 204 | FORREF(:,13) = (/ 0.462684E-05_JPRB, 0.238826E-07_JPRB, 0.325793E-07_JPRB /) |
|---|
| 205 | FORREF(:,14) = (/ 0.808644E-06_JPRB, 0.105126E-11_JPRB, 0.148691E-07_JPRB /) |
|---|
| 206 | FORREF(:,15) = (/ 0.865024E-12_JPRB, 0.822434E-12_JPRB, 0.825756E-12_JPRB /) |
|---|
| 207 | FORREF(:,16) = (/ 0.945747E-12_JPRB, 0.802065E-12_JPRB, 0.724732E-12_JPRB /) |
|---|
| 208 | |
|---|
| 209 | ! ----------------------------------------------------------------- |
|---|
| 210 | ! The array SELFREF contains the coefficient of the water vapor |
|---|
| 211 | ! self-continuum (including the energy term). The first index |
|---|
| 212 | ! refers to temperature in 7.2 degree increments. For instance, |
|---|
| 213 | ! JT = 1 refers to a temperature of 245.6, JT = 2 refers to 252.8, |
|---|
| 214 | ! etc. The second index runs over the g-channel (1 to 16). |
|---|
| 215 | |
|---|
| 216 | SELFREF(:, 1) = (/ & |
|---|
| 217 | & 0.637755E-05_JPRB, 0.403921E-05_JPRB, 0.255823E-05_JPRB, 0.162025E-05_JPRB, 0.102618E-05_JPRB, & |
|---|
| 218 | & 0.649930E-06_JPRB, 0.411632E-06_JPRB, 0.260707E-06_JPRB, 0.165118E-06_JPRB, 0.104577E-06_JPRB /) |
|---|
| 219 | SELFREF(:, 2) = (/ & |
|---|
| 220 | & 0.180887E-04_JPRB, 0.108890E-04_JPRB, 0.655493E-05_JPRB, 0.394592E-05_JPRB, 0.237536E-05_JPRB, & |
|---|
| 221 | & 0.142991E-05_JPRB, 0.860774E-06_JPRB, 0.518167E-06_JPRB, 0.311925E-06_JPRB, 0.187772E-06_JPRB /) |
|---|
| 222 | SELFREF(:, 3) = (/ & |
|---|
| 223 | & 0.212261E-04_JPRB, 0.150697E-04_JPRB, 0.106989E-04_JPRB, 0.759581E-05_JPRB, 0.539274E-05_JPRB, & |
|---|
| 224 | & 0.382864E-05_JPRB, 0.271819E-05_JPRB, 0.192981E-05_JPRB, 0.137009E-05_JPRB, 0.972711E-06_JPRB /) |
|---|
| 225 | SELFREF(:, 4) = (/ & |
|---|
| 226 | & 0.132497E-04_JPRB, 0.118071E-04_JPRB, 0.105216E-04_JPRB, 0.937599E-05_JPRB, 0.835516E-05_JPRB, & |
|---|
| 227 | & 0.744547E-05_JPRB, 0.663482E-05_JPRB, 0.591243E-05_JPRB, 0.526870E-05_JPRB, 0.469506E-05_JPRB /) |
|---|
| 228 | SELFREF(:, 5) = (/ & |
|---|
| 229 | & 0.124069E-04_JPRB, 0.120785E-04_JPRB, 0.117589E-04_JPRB, 0.114477E-04_JPRB, 0.111447E-04_JPRB, & |
|---|
| 230 | & 0.108498E-04_JPRB, 0.105626E-04_JPRB, 0.102831E-04_JPRB, 0.100109E-04_JPRB, 0.974601E-05_JPRB /) |
|---|
| 231 | SELFREF(:, 6) = (/ & |
|---|
| 232 | & 0.411994E-04_JPRB, 0.372560E-04_JPRB, 0.336901E-04_JPRB, 0.304654E-04_JPRB, 0.275494E-04_JPRB, & |
|---|
| 233 | & 0.249126E-04_JPRB, 0.225281E-04_JPRB, 0.203718E-04_JPRB, 0.184219E-04_JPRB, 0.166587E-04_JPRB /) |
|---|
| 234 | SELFREF(:, 7) = (/ & |
|---|
| 235 | & 0.537376E-04_JPRB, 0.501002E-04_JPRB, 0.467090E-04_JPRB, 0.435473E-04_JPRB, 0.405996E-04_JPRB, & |
|---|
| 236 | & 0.378515E-04_JPRB, 0.352893E-04_JPRB, 0.329006E-04_JPRB, 0.306736E-04_JPRB, 0.285974E-04_JPRB /) |
|---|
| 237 | SELFREF(:, 8) = (/ & |
|---|
| 238 | & 0.494279E-04_JPRB, 0.475365E-04_JPRB, 0.457175E-04_JPRB, 0.439681E-04_JPRB, 0.422857E-04_JPRB, & |
|---|
| 239 | & 0.406676E-04_JPRB, 0.391114E-04_JPRB, 0.376148E-04_JPRB, 0.361755E-04_JPRB, 0.347912E-04_JPRB /) |
|---|
| 240 | SELFREF(:, 9) = (/ & |
|---|
| 241 | & 0.377444E-04_JPRB, 0.378199E-04_JPRB, 0.378956E-04_JPRB, 0.379715E-04_JPRB, 0.380475E-04_JPRB, & |
|---|
| 242 | & 0.381236E-04_JPRB, 0.381999E-04_JPRB, 0.382763E-04_JPRB, 0.383529E-04_JPRB, 0.384297E-04_JPRB /) |
|---|
| 243 | SELFREF(:,10) = (/ & |
|---|
| 244 | & 0.245916E-04_JPRB, 0.267183E-04_JPRB, 0.290289E-04_JPRB, 0.315394E-04_JPRB, 0.342669E-04_JPRB, & |
|---|
| 245 | & 0.372304E-04_JPRB, 0.404501E-04_JPRB, 0.439483E-04_JPRB, 0.477490E-04_JPRB, 0.518784E-04_JPRB /) |
|---|
| 246 | SELFREF(:,11) = (/ & |
|---|
| 247 | & 0.186528E-04_JPRB, 0.211417E-04_JPRB, 0.239628E-04_JPRB, 0.271603E-04_JPRB, 0.307845E-04_JPRB, & |
|---|
| 248 | & 0.348923E-04_JPRB, 0.395482E-04_JPRB, 0.448254E-04_JPRB, 0.508068E-04_JPRB, 0.575863E-04_JPRB /) |
|---|
| 249 | SELFREF(:,12) = (/ & |
|---|
| 250 | & 0.109896E-04_JPRB, 0.133794E-04_JPRB, 0.162890E-04_JPRB, 0.198312E-04_JPRB, 0.241438E-04_JPRB, & |
|---|
| 251 | & 0.293942E-04_JPRB, 0.357864E-04_JPRB, 0.435686E-04_JPRB, 0.530432E-04_JPRB, 0.645781E-04_JPRB /) |
|---|
| 252 | SELFREF(:,13) = (/ & |
|---|
| 253 | & 0.183885E-06_JPRB, 0.391019E-06_JPRB, 0.831472E-06_JPRB, 0.176806E-05_JPRB, 0.375966E-05_JPRB, & |
|---|
| 254 | & 0.799463E-05_JPRB, 0.170000E-04_JPRB, 0.361492E-04_JPRB, 0.768686E-04_JPRB, 0.163455E-03_JPRB /) |
|---|
| 255 | SELFREF(:,14) = (/ & |
|---|
| 256 | & 0.466057E-07_JPRB, 0.937419E-07_JPRB, 0.188551E-06_JPRB, 0.379248E-06_JPRB, 0.762813E-06_JPRB, & |
|---|
| 257 | & 0.153431E-05_JPRB, 0.308608E-05_JPRB, 0.620729E-05_JPRB, 0.124852E-04_JPRB, 0.251126E-04_JPRB /) |
|---|
| 258 | SELFREF(:,15) = (/ & |
|---|
| 259 | & 0.248961E-06_JPRB, 0.216780E-06_JPRB, 0.188758E-06_JPRB, 0.164358E-06_JPRB, 0.143113E-06_JPRB, & |
|---|
| 260 | & 0.124613E-06_JPRB, 0.108505E-06_JPRB, 0.944795E-07_JPRB, 0.822667E-07_JPRB, 0.716326E-07_JPRB /) |
|---|
| 261 | SELFREF(:,16) = (/ & |
|---|
| 262 | & 0.252246E-06_JPRB, 0.220335E-06_JPRB, 0.192462E-06_JPRB, 0.168114E-06_JPRB, 0.146847E-06_JPRB, & |
|---|
| 263 | & 0.128270E-06_JPRB, 0.112043E-06_JPRB, 0.978688E-07_JPRB, 0.854878E-07_JPRB, 0.746731E-07_JPRB /) |
|---|
| 264 | |
|---|
| 265 | IF (LHOOK) CALL DR_HOOK('SRTM_KGB24',1,ZHOOK_HANDLE) |
|---|
| 266 | RETURN |
|---|
| 267 | |
|---|
| 268 | 1001 CONTINUE |
|---|
| 269 | CALL ABOR1("SRTM_KGB24:ERROR READING FILE RADSRTM") |
|---|
| 270 | |
|---|
| 271 | END SUBROUTINE SRTM_KGB24 |
|---|
