| [2089] | 1 | SUBROUTINE SUSW |
|---|
| 2 | |
|---|
| 3 | !**** *SUSW* - INITIALIZE COMMON YOESW |
|---|
| 4 | |
|---|
| 5 | ! PURPOSE. |
|---|
| 6 | ! -------- |
|---|
| 7 | ! INITIALIZE YOESW, THE COMMON THAT CONTAINS COEFFICIENTS |
|---|
| 8 | ! NEEDED TO RUN THE SHORTWAVE RADIATION SUBROUTINES |
|---|
| 9 | |
|---|
| 10 | !** INTERFACE. |
|---|
| 11 | ! ---------- |
|---|
| 12 | ! *CALL* *SUSW |
|---|
| 13 | |
|---|
| 14 | ! EXPLICIT ARGUMENTS : |
|---|
| 15 | ! -------------------- |
|---|
| 16 | ! NONE |
|---|
| 17 | |
|---|
| 18 | ! IMPLICIT ARGUMENTS : |
|---|
| 19 | ! -------------------- |
|---|
| 20 | ! COMMON YOESW |
|---|
| 21 | |
|---|
| 22 | ! METHOD. |
|---|
| 23 | ! ------- |
|---|
| 24 | ! SEE DOCUMENTATION |
|---|
| 25 | |
|---|
| 26 | ! EXTERNALS. |
|---|
| 27 | ! ---------- |
|---|
| 28 | |
|---|
| 29 | ! REFERENCE. |
|---|
| 30 | ! ---------- |
|---|
| 31 | ! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS |
|---|
| 32 | |
|---|
| 33 | ! AUTHOR. |
|---|
| 34 | ! ------- |
|---|
| 35 | ! JEAN-JACQUES MORCRETTE *ECMWF* |
|---|
| 36 | |
|---|
| 37 | ! MODIFICATIONS. |
|---|
| 38 | ! -------------- |
|---|
| 39 | ! ORIGINAL : 88-12-15 |
|---|
| 40 | |
|---|
| 41 | ! ------------------------------------------------------------------ |
|---|
| 42 | |
|---|
| 43 | #include "tsmbkind.h" |
|---|
| 44 | |
|---|
| 45 | USE YOMCST , ONLY : RG |
|---|
| 46 | USE YOESW , ONLY : APAD ,BPAD ,RRAY ,RSUN ,& |
|---|
| 47 | &RPDH1 ,RPDU1 ,RPNH ,RPNU ,RSWCE ,& |
|---|
| 48 | &RSWCP ,RTDH2O ,RTDUMG ,RTH2O ,RTUMG ,D |
|---|
| 49 | |
|---|
| 50 | ! ---------------------------------------------------------------- |
|---|
| 51 | |
|---|
| 52 | !* 1. SET VALUES. |
|---|
| 53 | ! ----------- |
|---|
| 54 | |
|---|
| 55 | |
|---|
| 56 | |
|---|
| 57 | IMPLICIT NONE |
|---|
| 58 | REAL_B :: ZPDH2O=0.8_JPRB |
|---|
| 59 | REAL_B :: ZPDUMG=0.75_JPRB |
|---|
| 60 | REAL_B :: ZPRH2O=30000._JPRB |
|---|
| 61 | REAL_B :: ZPRUMG=30000._JPRB |
|---|
| 62 | |
|---|
| 63 | ! LOCAL REAL SCALARS |
|---|
| 64 | REAL_B :: ZH2O, ZUMG |
|---|
| 65 | |
|---|
| 66 | RPDH1=ZPDH2O+_ONE_ |
|---|
| 67 | RPDU1=ZPDUMG+_ONE_ |
|---|
| 68 | ZH2O=_ONE_/( 10._JPRB* RG * RPDH1 ) |
|---|
| 69 | ZUMG=_ONE_/( 10._JPRB* RG * RPDU1 ) |
|---|
| 70 | RPNU = ZUMG/(ZPRUMG**ZPDUMG) |
|---|
| 71 | RPNH = ZH2O/(ZPRH2O**ZPDH2O) |
|---|
| 72 | |
|---|
| 73 | ! RSWCP=0.002*RSWCE |
|---|
| 74 | |
|---|
| 75 | |
|---|
| 76 | RTDH2O = 0.40_JPRB |
|---|
| 77 | RTDUMG = 0.375_JPRB |
|---|
| 78 | RTH2O = 240._JPRB |
|---|
| 79 | RTUMG = 240._JPRB |
|---|
| 80 | RSWCE = 0._JPRB |
|---|
| 81 | RSWCP = 0._JPRB |
|---|
| 82 | |
|---|
| 83 | RSUN(1) = 0.441676_JPRB |
|---|
| 84 | D(1,:) = (/ 0.00_JPRB, 0.00_JPRB, 0.00_JPRB /) |
|---|
| 85 | !* DERIVED FROM HITRAN APRIL 1991 |
|---|
| 86 | ! H2O: Pref=300 hPa, Tref=240K, Pdep=0.8 |
|---|
| 87 | ! O3 : unchanged |
|---|
| 88 | |
|---|
| 89 | APAD(1, 1, :) = (/& |
|---|
| 90 | &0.912418292E+05_JPRB,& |
|---|
| 91 | &0.723613782E+05_JPRB,& |
|---|
| 92 | &0.596037057E+04_JPRB,& |
|---|
| 93 | &0.000000000E-00_JPRB,& |
|---|
| 94 | &0.000000000E-00_JPRB,& |
|---|
| 95 | &0.000000000E-00_JPRB,& |
|---|
| 96 | &0.000000000E-00_JPRB/) |
|---|
| 97 | APAD(1, 2, :) = (/& |
|---|
| 98 | &0.000000000E-00_JPRB,& |
|---|
| 99 | &0.000000000E-00_JPRB,& |
|---|
| 100 | &0.000000000E-00_JPRB,& |
|---|
| 101 | &0.000000000E-00_JPRB,& |
|---|
| 102 | &0.000000000E-00_JPRB,& |
|---|
| 103 | &0.000000000E-00_JPRB,& |
|---|
| 104 | &0.000000000E-00_JPRB/) |
|---|
| 105 | APAD(1, 3, :) = (/& |
|---|
| 106 | &0.925887084E-04_JPRB,& |
|---|
| 107 | &0.129353723E-01_JPRB,& |
|---|
| 108 | &0.800821928E+00_JPRB,& |
|---|
| 109 | &0.242715973E+02_JPRB,& |
|---|
| 110 | &0.878331486E+02_JPRB,& |
|---|
| 111 | &0.191559725E+02_JPRB,& |
|---|
| 112 | &0.000000000E+00_JPRB/) |
|---|
| 113 | |
|---|
| 114 | BPAD(1, 1, :) = (/& |
|---|
| 115 | &0.912418292E+05_JPRB,& |
|---|
| 116 | &0.724555318E+05_JPRB,& |
|---|
| 117 | &0.602593328E+04_JPRB,& |
|---|
| 118 | &0.100000000E+01_JPRB,& |
|---|
| 119 | &0.000000000E-00_JPRB,& |
|---|
| 120 | &0.000000000E-00_JPRB,& |
|---|
| 121 | &0.000000000E-00_JPRB/) |
|---|
| 122 | BPAD(1, 2, :) = (/& |
|---|
| 123 | &0.000000000E-00_JPRB,& |
|---|
| 124 | &0.000000000E-00_JPRB,& |
|---|
| 125 | &0.000000000E-00_JPRB,& |
|---|
| 126 | &0.000000000E-00_JPRB,& |
|---|
| 127 | &0.000000000E-00_JPRB,& |
|---|
| 128 | &0.000000000E-00_JPRB,& |
|---|
| 129 | &0.000000000E-00_JPRB/) |
|---|
| 130 | BPAD(1, 3, :) = (/& |
|---|
| 131 | &0.925887084E-04_JPRB,& |
|---|
| 132 | &0.131812683E-01_JPRB,& |
|---|
| 133 | &0.812706117E+00_JPRB,& |
|---|
| 134 | &0.249863591E+02_JPRB,& |
|---|
| 135 | &0.931071925E+02_JPRB,& |
|---|
| 136 | &0.252233437E+02_JPRB,& |
|---|
| 137 | &0.100000000E+01_JPRB/) |
|---|
| 138 | |
|---|
| 139 | RRAY(1,1:6) = (/& |
|---|
| 140 | &.428937E-01_JPRB, .890743E+00_JPRB,-.288555E+01_JPRB,& |
|---|
| 141 | &.522744E+01_JPRB,-.469173E+01_JPRB, .161645E+01_JPRB/) |
|---|
| 142 | |
|---|
| 143 | RSUN(2) = 0.558324_JPRB |
|---|
| 144 | |
|---|
| 145 | D(2,1:3) = (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.800000000_JPRB /) |
|---|
| 146 | |
|---|
| 147 | !* INTERVAL 2: 0.68 - 4.00 MICRONS |
|---|
| 148 | !* DERIVED FROM HITRAN APRIL 1991 |
|---|
| 149 | ! H2O: Pref=300 hPa, Tref=240K, Pdep=0.80 |
|---|
| 150 | ! UMG: Pref=300 hPa, Tref=240K, Pdep=0.75 (CO2+O2+CH4+N2O+CO) |
|---|
| 151 | ! O3 : unchanged |
|---|
| 152 | |
|---|
| 153 | APAD(2, 1, :) = (/& |
|---|
| 154 | &0.376655383E-08_JPRB,& |
|---|
| 155 | &0.978576773E-04_JPRB,& |
|---|
| 156 | &0.387714006E+00_JPRB,& |
|---|
| 157 | &0.118461660E+03_JPRB,& |
|---|
| 158 | &0.119079797E+04_JPRB,& |
|---|
| 159 | &0.293353397E+03_JPRB,& |
|---|
| 160 | &0.000000000E+00_JPRB/) |
|---|
| 161 | APAD(2, 2, :) = (/& |
|---|
| 162 | &0.739646016E-08_JPRB,& |
|---|
| 163 | &0.131849595E-03_JPRB,& |
|---|
| 164 | &0.437772681E+00_JPRB,& |
|---|
| 165 | &0.151345118E+03_JPRB,& |
|---|
| 166 | &0.233628890E+04_JPRB,& |
|---|
| 167 | &0.797219934E+03_JPRB,& |
|---|
| 168 | &0.000000000E+00_JPRB/) |
|---|
| 169 | APAD(2, 3, :) = (/& |
|---|
| 170 | &0.410177786E+03_JPRB,& |
|---|
| 171 | &0.672595424E+02_JPRB,& |
|---|
| 172 | &0.000000000E-00_JPRB,& |
|---|
| 173 | &0.000000000E-00_JPRB,& |
|---|
| 174 | &0.000000000E-00_JPRB,& |
|---|
| 175 | &0.000000000E-00_JPRB,& |
|---|
| 176 | &0.000000000E+00_JPRB/) |
|---|
| 177 | |
|---|
| 178 | BPAD(2, 1, :) = (/& |
|---|
| 179 | &0.376655383E-08_JPRB,& |
|---|
| 180 | &0.979023421E-04_JPRB,& |
|---|
| 181 | &0.388611139E+00_JPRB,& |
|---|
| 182 | &0.120291383E+03_JPRB,& |
|---|
| 183 | &0.130531005E+04_JPRB,& |
|---|
| 184 | &0.415049409E+03_JPRB,& |
|---|
| 185 | &0.100000000E+01_JPRB/) |
|---|
| 186 | BPAD(2, 2, :) = (/& |
|---|
| 187 | &0.739646016E-08_JPRB,& |
|---|
| 188 | &0.131861712E-03_JPRB,& |
|---|
| 189 | &0.437949001E+00_JPRB,& |
|---|
| 190 | &0.151692730E+03_JPRB,& |
|---|
| 191 | &0.237071130E+04_JPRB,& |
|---|
| 192 | &0.867914360E+03_JPRB,& |
|---|
| 193 | &0.100000000E+01_JPRB/) |
|---|
| 194 | BPAD(2, 3, :) = (/& |
|---|
| 195 | &0.410177786E+03_JPRB,& |
|---|
| 196 | &0.731185438E+02_JPRB,& |
|---|
| 197 | &0.100000000E+01_JPRB,& |
|---|
| 198 | &0.000000000E+00_JPRB,& |
|---|
| 199 | &0.000000000E+00_JPRB,& |
|---|
| 200 | &0.000000000E+00_JPRB,& |
|---|
| 201 | &0.000000000E+00_JPRB/) |
|---|
| 202 | |
|---|
| 203 | RRAY(2,1:6) = (/& |
|---|
| 204 | &.697200E-02_JPRB, .173297E-01_JPRB,-.850903E-01_JPRB,& |
|---|
| 205 | &.248261E+00_JPRB,-.302031E+00_JPRB, .129662E+00_JPRB/) |
|---|
| 206 | RETURN |
|---|
| 207 | END SUBROUTINE SUSW |
|---|
