| Line |
Branch |
Exec |
Source |
| 1 |
|
✗ |
SUBROUTINE SUSRTAER |
| 2 |
|
|
|
| 3 |
|
|
!**** *SUSRTAER* - INITIALIZE COMMON YOESRTAER |
| 4 |
|
|
|
| 5 |
|
|
! PURPOSE. |
| 6 |
|
|
! -------- |
| 7 |
|
|
! INITIALIZE YOESRTAER, THE COMMON THAT CONTAINS THE |
| 8 |
|
|
! RADIATIVE CHARACTERISTICS OF THE AEROSOLS |
| 9 |
|
|
|
| 10 |
|
|
!** INTERFACE. |
| 11 |
|
|
! ---------- |
| 12 |
|
|
! ----- ----- |
| 13 |
|
|
|
| 14 |
|
|
! EXPLICIT ARGUMENTS : |
| 15 |
|
|
! -------------------- |
| 16 |
|
|
! NONE |
| 17 |
|
|
|
| 18 |
|
|
! IMPLICIT ARGUMENTS : |
| 19 |
|
|
! -------------------- |
| 20 |
|
|
! COMMON YOESRTAER |
| 21 |
|
|
|
| 22 |
|
|
! METHOD. |
| 23 |
|
|
! ------- |
| 24 |
|
|
! SEE DOCUMENTATION |
| 25 |
|
|
|
| 26 |
|
|
! EXTERNALS. |
| 27 |
|
|
! ---------- |
| 28 |
|
|
|
| 29 |
|
|
! REFERENCE. |
| 30 |
|
|
! ---------- |
| 31 |
|
|
! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE "IFS MODEL" |
| 32 |
|
|
|
| 33 |
|
|
! AUTHOR. |
| 34 |
|
|
! ------- |
| 35 |
|
|
! JEAN-JACQUES MORCRETTE *ECMWF* |
| 36 |
|
|
|
| 37 |
|
|
! MODIFICATIONS. |
| 38 |
|
|
! -------------- |
| 39 |
|
|
! ORIGINAL : 03-03-06 |
| 40 |
|
|
! 03-03-06 JJMorcrette 14 spectral intervals for RRTM_SW |
| 41 |
|
|
! M.Hamrud 01-Oct-2003 CY28 Cleaning |
| 42 |
|
|
|
| 43 |
|
|
! ------------------------------------------------------------------ |
| 44 |
|
|
|
| 45 |
|
|
USE PARKIND1 ,ONLY : JPIM ,JPRB |
| 46 |
|
|
USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
| 47 |
|
|
|
| 48 |
|
|
USE YOESRTAER , ONLY : RSRTAUA ,RSRPIZA ,RSRASYA |
| 49 |
|
|
|
| 50 |
|
|
! ---------------------------------------------------------------- |
| 51 |
|
|
|
| 52 |
|
|
IMPLICIT NONE |
| 53 |
|
|
|
| 54 |
|
|
! DUMMY INTEGER REAL |
| 55 |
|
|
REAL(KIND=JPRB) :: ZTAUA14(14,6),ZPIZA14(14,6),ZCGA14(14,6) |
| 56 |
|
|
|
| 57 |
|
|
INTEGER(KIND=JPIM) :: JAER, JNU |
| 58 |
|
|
REAL(KIND=JPRB) :: ZHOOK_HANDLE |
| 59 |
|
|
|
| 60 |
|
|
! ---------------------------------------------------------------- |
| 61 |
|
|
|
| 62 |
|
|
!* 1. SHORTWAVE COEFFICIENTS |
| 63 |
|
|
! ---------------------- |
| 64 |
|
|
!======================================================================= |
| 65 |
|
|
!-- The (old) five aerosol types were respectively: |
| 66 |
|
|
|
| 67 |
|
|
! 1/ continental average (+desert) 2/ maritime |
| 68 |
|
|
! 3/ urban 4/ volcanic active |
| 69 |
|
|
! 5/ stratospheric background |
| 70 |
|
|
|
| 71 |
|
|
!-- old values were not spectrally defined: |
| 72 |
|
|
! ZTAU2 = .730719, .912819, .725059, .745405, .682188 |
| 73 |
|
|
! ZPIZA2 = .872212, .982545, .623143, .944887, .997975 |
| 74 |
|
|
! ZCGA2 = .647596, .739002, .580845, .662657, .624246 |
| 75 |
|
|
!======================================================================= |
| 76 |
|
|
|
| 77 |
|
|
!-- The six aerosol types are respectively: |
| 78 |
|
|
|
| 79 |
|
|
! 1/ continental average 2/ maritime |
| 80 |
|
|
! 3/ desert 4/ urban |
| 81 |
|
|
! 5/ volcanic active 6/ stratospheric background |
| 82 |
|
|
|
| 83 |
|
|
! The quantities given are: |
| 84 |
|
|
! TAU : ratio of average optical thickness in interval to that at 0.55 |
| 85 |
|
|
! micron |
| 86 |
|
|
! PIZA: average single scattering albedo |
| 87 |
|
|
! CGA : average asymmetry factor |
| 88 |
|
|
|
| 89 |
|
|
! computed from Hess and Koepke (con, mar, des, urb) |
| 90 |
|
|
! from Bonnel et al. (vol, str) |
| 91 |
|
|
|
| 92 |
|
|
! 1.1 FOURTEEN SPECTRAL INTERVALS FOR RRTM_SW |
| 93 |
|
|
|
| 94 |
|
|
! SW : 14 spectral intervals |
| 95 |
|
|
! 3.846 - 3.077 |
| 96 |
|
|
! 3.077 - 2.500 |
| 97 |
|
|
! 2.500 - 2.150 |
| 98 |
|
|
! 2.150 - 1.942 |
| 99 |
|
|
! 1.942 - 1.626 |
| 100 |
|
|
! 1.626 - 1.299 |
| 101 |
|
|
! 1.299 - 1.242 |
| 102 |
|
|
! 1.242 - 0.7782 |
| 103 |
|
|
! 0.7782- 0.6250 |
| 104 |
|
|
! 0.6250- 0.4415 |
| 105 |
|
|
! 0.4415- 0.3448 |
| 106 |
|
|
! 0.3448- 0.2632 |
| 107 |
|
|
! 0.2632- 0.2000 |
| 108 |
|
|
! 12.195 - 3.846 |
| 109 |
|
|
|
| 110 |
|
✗ |
IF (LHOOK) CALL DR_HOOK('SUSRTAER',0,ZHOOK_HANDLE) |
| 111 |
|
|
ZTAUA14( 1, :)= (/& |
| 112 |
|
✗ |
& 0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /) |
| 113 |
|
|
ZTAUA14( 2, :)= (/& |
| 114 |
|
✗ |
& 0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /) |
| 115 |
|
|
ZTAUA14( 3, :)= (/& |
| 116 |
|
✗ |
& 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) |
| 117 |
|
|
ZTAUA14( 4, :)= (/& |
| 118 |
|
✗ |
& 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) |
| 119 |
|
|
ZTAUA14( 5, :)= (/& |
| 120 |
|
✗ |
& 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) |
| 121 |
|
|
ZTAUA14( 6, :)= (/& |
| 122 |
|
✗ |
& 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) |
| 123 |
|
|
ZTAUA14( 7, :)= (/& |
| 124 |
|
✗ |
& 0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /) |
| 125 |
|
|
ZTAUA14( 8, :)= (/& |
| 126 |
|
✗ |
& 0.52838_JPRB , 0.93285_JPRB , 0.93449_JPRB , 0.53078_JPRB , 0.67148_JPRB , 0.46608_JPRB /) |
| 127 |
|
|
ZTAUA14( 9, :)= (/& |
| 128 |
|
✗ |
& 0.52838_JPRB , 0.93285_JPRB , 0.93449_JPRB , 0.53078_JPRB , 0.67148_JPRB , 0.46608_JPRB /) |
| 129 |
|
|
ZTAUA14(10, :)= (/& |
| 130 |
|
✗ |
& 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) |
| 131 |
|
|
ZTAUA14(11, :)= (/& |
| 132 |
|
✗ |
& 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) |
| 133 |
|
|
ZTAUA14(12, :)= (/& |
| 134 |
|
✗ |
& 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) |
| 135 |
|
|
ZTAUA14(13, :)= (/& |
| 136 |
|
✗ |
& 1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /) |
| 137 |
|
|
ZTAUA14(14, :)= (/& |
| 138 |
|
✗ |
& 0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /) |
| 139 |
|
|
|
| 140 |
|
|
ZPIZA14( 1, :)= (/& |
| 141 |
|
✗ |
& .5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/) |
| 142 |
|
|
ZPIZA14( 2, :)= (/& |
| 143 |
|
✗ |
& .5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/) |
| 144 |
|
|
ZPIZA14( 3, :)= (/& |
| 145 |
|
✗ |
& .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) |
| 146 |
|
|
ZPIZA14( 4, :)= (/& |
| 147 |
|
✗ |
& .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) |
| 148 |
|
|
ZPIZA14( 5, :)= (/& |
| 149 |
|
✗ |
& .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) |
| 150 |
|
|
ZPIZA14( 6, :)= (/& |
| 151 |
|
✗ |
& .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) |
| 152 |
|
|
ZPIZA14( 7, :)= (/& |
| 153 |
|
✗ |
& .8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/) |
| 154 |
|
|
ZPIZA14( 8, :)= (/& |
| 155 |
|
✗ |
& .8970131_JPRB, .9984940_JPRB, .9245594_JPRB, .7768385_JPRB, .9532763_JPRB, .9999999_JPRB/) |
| 156 |
|
|
ZPIZA14( 9, :)= (/& |
| 157 |
|
✗ |
& .8970131_JPRB, .9984940_JPRB, .9245594_JPRB, .7768385_JPRB, .9532763_JPRB, .9999999_JPRB/) |
| 158 |
|
|
ZPIZA14(10, :)= (/& |
| 159 |
|
✗ |
& .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) |
| 160 |
|
|
ZPIZA14(11, :)= (/& |
| 161 |
|
✗ |
& .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) |
| 162 |
|
|
ZPIZA14(12, :)= (/& |
| 163 |
|
✗ |
& .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) |
| 164 |
|
|
ZPIZA14(13, :)= (/& |
| 165 |
|
✗ |
& .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/) |
| 166 |
|
|
ZPIZA14(14, :)= (/& |
| 167 |
|
✗ |
& .5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/) |
| 168 |
|
|
|
| 169 |
|
|
ZCGA14( 1, :)= (/& |
| 170 |
|
✗ |
& 0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/) |
| 171 |
|
|
ZCGA14( 2, :)= (/& |
| 172 |
|
✗ |
& 0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/) |
| 173 |
|
|
ZCGA14( 3, :)= (/& |
| 174 |
|
✗ |
& 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) |
| 175 |
|
|
ZCGA14( 4, :)= (/& |
| 176 |
|
✗ |
& 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) |
| 177 |
|
|
ZCGA14( 5, :)= (/& |
| 178 |
|
✗ |
& 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) |
| 179 |
|
|
ZCGA14( 6, :)= (/& |
| 180 |
|
✗ |
& 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) |
| 181 |
|
|
ZCGA14( 7, :)= (/& |
| 182 |
|
✗ |
& 0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/) |
| 183 |
|
|
ZCGA14( 8, :)= (/& |
| 184 |
|
✗ |
& 0.668431_JPRB, 0.788530_JPRB, 0.698682_JPRB, 0.657422_JPRB, .6735182_JPRB, .6519706_JPRB/) |
| 185 |
|
|
ZCGA14( 9, :)= (/& |
| 186 |
|
✗ |
& 0.668431_JPRB, 0.788530_JPRB, 0.698682_JPRB, 0.657422_JPRB, .6735182_JPRB, .6519706_JPRB/) |
| 187 |
|
|
ZCGA14(10, :)= (/& |
| 188 |
|
✗ |
& 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) |
| 189 |
|
|
ZCGA14(11, :)= (/& |
| 190 |
|
✗ |
& 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) |
| 191 |
|
|
ZCGA14(12, :)= (/& |
| 192 |
|
✗ |
& 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) |
| 193 |
|
|
ZCGA14(13, :)= (/& |
| 194 |
|
✗ |
& 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/) |
| 195 |
|
|
ZCGA14(14, :)= (/& |
| 196 |
|
✗ |
& 0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/) |
| 197 |
|
|
|
| 198 |
|
|
! ---------------------------------------------------------------- |
| 199 |
|
|
|
| 200 |
|
✗ |
DO JNU=1,14 |
| 201 |
|
✗ |
DO JAER=1,6 |
| 202 |
|
✗ |
RSRTAUA(JNU,JAER)=ZTAUA14(JNU,JAER) |
| 203 |
|
✗ |
RSRPIZA(JNU,JAER)=ZPIZA14(JNU,JAER) |
| 204 |
|
✗ |
RSRASYA(JNU,JAER)=ZCGA14 (JNU,JAER) |
| 205 |
|
|
ENDDO |
| 206 |
|
|
ENDDO |
| 207 |
|
|
|
| 208 |
|
|
! ---------------------------------------------------------------- |
| 209 |
|
|
|
| 210 |
|
✗ |
IF (LHOOK) CALL DR_HOOK('SUSRTAER',1,ZHOOK_HANDLE) |
| 211 |
|
✗ |
END SUBROUTINE SUSRTAER |
| 212 |
|
|
|
| 213 |
|
|
|