| 1 | MODULE chemparam_mod |
|---|
| 2 | |
|---|
| 3 | !MODULE qui définit les indices des traceurs et leurs masses molaires. |
|---|
| 4 | ! utilise aussi pour variables communes nuages/photochimie |
|---|
| 5 | |
|---|
| 6 | IMPLICIT NONE |
|---|
| 7 | |
|---|
| 8 | !---------------------------------------------------------------------------- |
|---|
| 9 | INTEGER, SAVE :: i_co2, i_co, i_h2, i_h2o, i_o1d, & |
|---|
| 10 | i_o, i_o2, i_o2dg, i_o3, i_h, & |
|---|
| 11 | i_oh, i_ho2, i_h2o2, i_cl, i_clo, & |
|---|
| 12 | i_cl2, i_hcl, i_hocl, i_clco, i_clco3, & |
|---|
| 13 | i_cocl2, i_s, i_so, i_so2, i_so3, & |
|---|
| 14 | i_s2o2, i_ocs, i_hso3, i_h2so4, i_s2, & |
|---|
| 15 | i_clso2, i_oscl, i_h2oliq, i_h2so4liq |
|---|
| 16 | |
|---|
| 17 | REAL, DIMENSION(:), SAVE, ALLOCATABLE :: M_tr |
|---|
| 18 | |
|---|
| 19 | |
|---|
| 20 | !---------------------------------------------------------------------------- |
|---|
| 21 | ! number of clouds layer modelized |
|---|
| 22 | ! INTEGER, PARAMETER :: nbr_cloud = 1 |
|---|
| 23 | INTEGER :: i_cloud |
|---|
| 24 | INTEGER, SAVE :: cloudmax |
|---|
| 25 | INTEGER, SAVE :: cloudmin |
|---|
| 26 | REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: R_MEDIAN |
|---|
| 27 | REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: STDDEV |
|---|
| 28 | REAL, SAVE :: RMI |
|---|
| 29 | REAL, SAVE :: RMA |
|---|
| 30 | REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: NBRTOT |
|---|
| 31 | REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: WH2SO4 |
|---|
| 32 | REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: rho_droplet |
|---|
| 33 | !---------------------------------------------------------------------------- |
|---|
| 34 | |
|---|
| 35 | !---------------------------------------------------------------------------- |
|---|
| 36 | CONTAINS |
|---|
| 37 | !---------------------------------------------------------------------------- |
|---|
| 38 | |
|---|
| 39 | SUBROUTINE cloud_ini(nbr_lon,nbr_lev) |
|---|
| 40 | |
|---|
| 41 | !============================================================= |
|---|
| 42 | ! cloud_ini definit le champ 3D des caracteristiques du nuage |
|---|
| 43 | |
|---|
| 44 | INTEGER :: nbr_lon,nbr_lev,i_lev |
|---|
| 45 | |
|---|
| 46 | ALLOCATE(NBRTOT(nbr_lon,nbr_lev)) |
|---|
| 47 | ALLOCATE(R_MEDIAN(nbr_lon,nbr_lev)) |
|---|
| 48 | ALLOCATE(STDDEV(nbr_lon,nbr_lev)) |
|---|
| 49 | ALLOCATE(WH2SO4(nbr_lon,nbr_lev)) |
|---|
| 50 | ALLOCATE(rho_droplet(nbr_lon,nbr_lev)) |
|---|
| 51 | |
|---|
| 52 | PRINT*,'==========================' |
|---|
| 53 | PRINT*,'Initialisation cloud layer' |
|---|
| 54 | PRINT*,'==========================' |
|---|
| 55 | PRINT*,'nbr_lon',nbr_lon |
|---|
| 56 | PRINT*,'nbr_lev',nbr_lev |
|---|
| 57 | |
|---|
| 58 | !============================================================= |
|---|
| 59 | ! Initialisation cloud layer 1 |
|---|
| 60 | !============================================================= |
|---|
| 61 | ! cloudmin et cloudmax niveaux du GCM |
|---|
| 62 | cloudmin= 15 |
|---|
| 63 | cloudmax= 50 |
|---|
| 64 | ! radius min et max en microns, *e-6 dans cloud_venus -> SETBIN |
|---|
| 65 | ! RMI=0.001 |
|---|
| 66 | ! RMA=100.0 |
|---|
| 67 | ! NBRTOT= 0. |
|---|
| 68 | |
|---|
| 69 | ! radius R_MEDIAN en m (donc *e-6 pour microns) |
|---|
| 70 | |
|---|
| 71 | R_MEDIAN(:,:)=0.0 ! Geometric Average Radius |
|---|
| 72 | STDDEV(:,:)=0.0 ! Geometric Std Deviation |
|---|
| 73 | |
|---|
| 74 | ! =============================================== |
|---|
| 75 | ! Knollenberg & Hunten, 1980 and James et al 1997 |
|---|
| 76 | ! =============================================== |
|---|
| 77 | |
|---|
| 78 | ! Lower Haze: mode 1 |
|---|
| 79 | DO i_lev=cloudmin,20 |
|---|
| 80 | R_MEDIAN(:,i_lev)=0.2e-6 |
|---|
| 81 | PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev) |
|---|
| 82 | STDDEV(:,i_lev)=1.56 |
|---|
| 83 | PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev) |
|---|
| 84 | END DO |
|---|
| 85 | |
|---|
| 86 | ! Lower Cloud: mode 2 |
|---|
| 87 | ! DO i_lev=21,23 |
|---|
| 88 | ! R_MEDIAN(:,i_lev)=1.4e-6 |
|---|
| 89 | ! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev) |
|---|
| 90 | ! STDDEV(:,i_lev)=1.23 |
|---|
| 91 | ! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev) |
|---|
| 92 | ! END DO |
|---|
| 93 | |
|---|
| 94 | ! Lower Cloud: mode 3 |
|---|
| 95 | DO i_lev=21,23 |
|---|
| 96 | R_MEDIAN(:,i_lev)=3.65e-6 |
|---|
| 97 | PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev) |
|---|
| 98 | STDDEV(:,i_lev)=1.28 |
|---|
| 99 | PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev) |
|---|
| 100 | END DO |
|---|
| 101 | |
|---|
| 102 | ! Middle Cloud: mode 2 prime |
|---|
| 103 | DO i_lev=24,28 |
|---|
| 104 | R_MEDIAN(:,i_lev)=1.4e-6 |
|---|
| 105 | PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev) |
|---|
| 106 | STDDEV(:,i_lev)=1.23 |
|---|
| 107 | PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev) |
|---|
| 108 | END DO |
|---|
| 109 | |
|---|
| 110 | ! Middle Cloud: mode 3 |
|---|
| 111 | ! DO i_lev=24,28 |
|---|
| 112 | ! R_MEDIAN(:,i_lev)=3.65e-6 |
|---|
| 113 | ! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev) |
|---|
| 114 | ! STDDEV(:,i_lev)=1.28 |
|---|
| 115 | ! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev) |
|---|
| 116 | ! END DO |
|---|
| 117 | |
|---|
| 118 | ! Middle Cloud: mode 4 |
|---|
| 119 | ! DO i_lev=24,28 |
|---|
| 120 | ! R_MEDIAN(:,i_lev)=7.0e-6 |
|---|
| 121 | ! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev) |
|---|
| 122 | ! STDDEV(:,i_lev)=1.3 |
|---|
| 123 | ! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev) |
|---|
| 124 | ! END DO |
|---|
| 125 | |
|---|
| 126 | ! Upper Cloud: mode 4 |
|---|
| 127 | ! DO i_lev=29,35 |
|---|
| 128 | ! R_MEDIAN(:,i_lev)=7.0e-6 |
|---|
| 129 | ! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev) |
|---|
| 130 | ! STDDEV(:,i_lev)=1.3 |
|---|
| 131 | ! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev) |
|---|
| 132 | ! END DO |
|---|
| 133 | |
|---|
| 134 | ! Upper Cloud: mode 3 |
|---|
| 135 | ! DO i_lev=29,35 |
|---|
| 136 | ! R_MEDIAN(:,i_lev)=3.65e-6 |
|---|
| 137 | ! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev) |
|---|
| 138 | ! STDDEV(:,i_lev)=1.28 |
|---|
| 139 | ! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev) |
|---|
| 140 | ! END DO |
|---|
| 141 | |
|---|
| 142 | ! Upper Cloud: mode 2 |
|---|
| 143 | DO i_lev=29,35 |
|---|
| 144 | R_MEDIAN(:,i_lev)=1.0e-6 |
|---|
| 145 | PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev) |
|---|
| 146 | STDDEV(:,i_lev)=1.29 |
|---|
| 147 | PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev) |
|---|
| 148 | END DO |
|---|
| 149 | |
|---|
| 150 | ! Upper Haze: mode 1 |
|---|
| 151 | DO i_lev=36, cloudmax |
|---|
| 152 | R_MEDIAN(:,i_lev)=0.2e-6 |
|---|
| 153 | PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev) |
|---|
| 154 | STDDEV(:,i_lev)=2.16 |
|---|
| 155 | PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev) |
|---|
| 156 | END DO |
|---|
| 157 | |
|---|
| 158 | !============================================================= |
|---|
| 159 | |
|---|
| 160 | PRINT*,'===============================' |
|---|
| 161 | PRINT*,'FIN Initialisation cloud layer' |
|---|
| 162 | PRINT*,'===============================' |
|---|
| 163 | |
|---|
| 164 | END SUBROUTINE |
|---|
| 165 | |
|---|
| 166 | SUBROUTINE chemparam_ini |
|---|
| 167 | USE infotrac |
|---|
| 168 | IMPLICIT NONE |
|---|
| 169 | INTEGER :: i |
|---|
| 170 | |
|---|
| 171 | ALLOCATE(M_tr(nqtot)) |
|---|
| 172 | |
|---|
| 173 | DO i=1, nqtot |
|---|
| 174 | |
|---|
| 175 | PRINT*,'i',i |
|---|
| 176 | PRINT*,'tname(i)',tname(i) |
|---|
| 177 | |
|---|
| 178 | SELECT CASE(tname(i)) |
|---|
| 179 | CASE('co2') |
|---|
| 180 | i_co2=i |
|---|
| 181 | PRINT*,'co2',i_co2 |
|---|
| 182 | M_tr(i_co2) = 44.0095 |
|---|
| 183 | CASE('co') |
|---|
| 184 | i_co=i |
|---|
| 185 | PRINT*,'co',i_co |
|---|
| 186 | M_tr(i_co)=28.0101 |
|---|
| 187 | CASE('h2') |
|---|
| 188 | i_h2=i |
|---|
| 189 | PRINT*,'h2',i_h2 |
|---|
| 190 | M_tr(i_h2)= 18.0153 |
|---|
| 191 | CASE('h2o') |
|---|
| 192 | i_h2o=i |
|---|
| 193 | PRINT*,'h2o',i_h2o |
|---|
| 194 | M_tr(i_h2o)=18.0153 |
|---|
| 195 | CASE('o1d') |
|---|
| 196 | i_o1d=i |
|---|
| 197 | PRINT*,'o1d',i_o1d |
|---|
| 198 | M_tr(i_o1d)=15.994 |
|---|
| 199 | CASE('o') |
|---|
| 200 | i_o=i |
|---|
| 201 | PRINT*,'o',i_o |
|---|
| 202 | M_tr(i_o)=15.994 |
|---|
| 203 | CASE('o2') |
|---|
| 204 | i_o2=i |
|---|
| 205 | PRINT*,'o2',i_o2 |
|---|
| 206 | M_tr(i_o2)=31.9988 |
|---|
| 207 | CASE('o2dg') |
|---|
| 208 | i_o2dg=i |
|---|
| 209 | PRINT*,'o2dg',i_o2dg |
|---|
| 210 | M_tr(i_o2dg)=31.9988 |
|---|
| 211 | CASE('o3') |
|---|
| 212 | i_o3=i |
|---|
| 213 | PRINT*,'o3',i_o3 |
|---|
| 214 | M_tr(i_o3)= 47.9982 |
|---|
| 215 | CASE('h') |
|---|
| 216 | i_h=i |
|---|
| 217 | PRINT*,'h',i_h |
|---|
| 218 | M_tr(i_h)= 1.00794 |
|---|
| 219 | CASE('oh') |
|---|
| 220 | i_oh=i |
|---|
| 221 | PRINT*,'oh',i_oh |
|---|
| 222 | M_tr(i_oh)=17.0073 |
|---|
| 223 | CASE('ho2') |
|---|
| 224 | i_ho2=i |
|---|
| 225 | PRINT*,'ho2',i_ho2 |
|---|
| 226 | M_tr(i_ho2)=33.0067 |
|---|
| 227 | CASE('h2o2') |
|---|
| 228 | i_h2o2=i |
|---|
| 229 | PRINT*,'h2o2',i_h2o2 |
|---|
| 230 | M_tr(i_h2o2)=34.0147 |
|---|
| 231 | CASE('cl') |
|---|
| 232 | i_cl=i |
|---|
| 233 | PRINT*,'cl',i_cl |
|---|
| 234 | M_tr(i_cl)=35.453 |
|---|
| 235 | CASE('clo') |
|---|
| 236 | i_clo=i |
|---|
| 237 | PRINT*,'clo',i_clo |
|---|
| 238 | M_tr(i_clo)=51.452 |
|---|
| 239 | CASE('cl2') |
|---|
| 240 | i_cl2=i |
|---|
| 241 | PRINT*,'cl2',i_cl2 |
|---|
| 242 | M_tr(i_cl2)=70.906 |
|---|
| 243 | CASE('hcl') |
|---|
| 244 | i_hcl=i |
|---|
| 245 | PRINT*,'hcl',i_hcl |
|---|
| 246 | M_tr(i_hcl)=36.461 |
|---|
| 247 | CASE('hocl') |
|---|
| 248 | i_hocl=i |
|---|
| 249 | PRINT*,'hocl',i_hocl |
|---|
| 250 | M_tr(i_hocl)=52.46 |
|---|
| 251 | CASE('clco') |
|---|
| 252 | i_clco=i |
|---|
| 253 | PRINT*,'clco',i_clco |
|---|
| 254 | M_tr(i_clco)=63.463 |
|---|
| 255 | CASE('clco3') |
|---|
| 256 | i_clco3=i |
|---|
| 257 | PRINT*,'clco3',i_clco3 |
|---|
| 258 | M_tr(i_clco3)=95.462 |
|---|
| 259 | CASE('cocl2') |
|---|
| 260 | i_cocl2=i |
|---|
| 261 | PRINT*,'cocl2',i_cocl2 |
|---|
| 262 | M_tr(i_cocl2)=98.916 |
|---|
| 263 | CASE('s') |
|---|
| 264 | i_s=i |
|---|
| 265 | PRINT*,'s',i_s |
|---|
| 266 | M_tr(i_s)=32.065 |
|---|
| 267 | CASE('so') |
|---|
| 268 | i_so=i |
|---|
| 269 | PRINT*,'so',i_so |
|---|
| 270 | M_tr(i_so)=48.0644 |
|---|
| 271 | CASE('so2') |
|---|
| 272 | i_so2=i |
|---|
| 273 | PRINT*,'so2',i_so2 |
|---|
| 274 | M_tr(i_so2)=64.064 |
|---|
| 275 | CASE('so3') |
|---|
| 276 | i_so3=i |
|---|
| 277 | PRINT*,'so3',i_so3 |
|---|
| 278 | M_tr(i_so3)=80.063 |
|---|
| 279 | CASE('s2o2') |
|---|
| 280 | i_s2o2=i |
|---|
| 281 | PRINT*,'s2o2',i_s2o2 |
|---|
| 282 | M_tr(i_s2o2)= 96.1288 |
|---|
| 283 | CASE('ocs') |
|---|
| 284 | i_ocs=i |
|---|
| 285 | PRINT*,'ocs',i_ocs |
|---|
| 286 | M_tr(i_ocs)=60.0751 |
|---|
| 287 | CASE('hso3') |
|---|
| 288 | i_hso3=i |
|---|
| 289 | PRINT*,'hso3',i_hso3 |
|---|
| 290 | M_tr(i_hso3)=81.071 |
|---|
| 291 | CASE('h2so4') |
|---|
| 292 | i_h2so4=i |
|---|
| 293 | PRINT*,'h2so4',i_h2so4 |
|---|
| 294 | M_tr(i_h2so4)=98.078 |
|---|
| 295 | CASE('s2') |
|---|
| 296 | i_s2=i |
|---|
| 297 | PRINT*,'s2',i_s2 |
|---|
| 298 | M_tr(i_s2)=64.13 |
|---|
| 299 | CASE('clso2') |
|---|
| 300 | i_clso2=i |
|---|
| 301 | PRINT*,'clso2',i_clso2 |
|---|
| 302 | M_tr(i_clso2)=99.517 |
|---|
| 303 | CASE('oscl') |
|---|
| 304 | i_oscl=i |
|---|
| 305 | PRINT*,'oscl',i_oscl |
|---|
| 306 | M_tr(i_oscl)=83.517 |
|---|
| 307 | CASE('h2oliq') |
|---|
| 308 | i_h2oliq=i |
|---|
| 309 | PRINT*,'h2oliq',i_h2oliq |
|---|
| 310 | M_tr(i_h2oliq)=18.0153 |
|---|
| 311 | CASE('h2so4liq') |
|---|
| 312 | i_h2so4liq=i |
|---|
| 313 | PRINT*,'h2so4liq',i_h2so4liq |
|---|
| 314 | M_tr(i_h2so4liq)=98.078 |
|---|
| 315 | END SELECT |
|---|
| 316 | |
|---|
| 317 | ! PRINT*,'M_tr(i)',M_tr(i) |
|---|
| 318 | END DO |
|---|
| 319 | |
|---|
| 320 | END SUBROUTINE chemparam_ini |
|---|
| 321 | |
|---|
| 322 | END MODULE chemparam_mod |
|---|
| 323 | |
|---|
