| [1969] | 1 | MODULE updatereffrad_mod |
|---|
| 2 | |
|---|
| 3 | IMPLICIT NONE |
|---|
| 4 | |
|---|
| 5 | CONTAINS |
|---|
| 6 | |
|---|
| [38] | 7 | SUBROUTINE updatereffrad(ngrid,nlayer, |
|---|
| [1974] | 8 | & rdust,rstormdust,rice,nuice, |
|---|
| [38] | 9 | & reffrad,nueffrad, |
|---|
| [744] | 10 | & pq,tauscaling,tau,pplay) |
|---|
| [1969] | 11 | USE updaterad, ONLY: updaterdust, updaterice_micro, |
|---|
| 12 | & updaterice_typ |
|---|
| [1036] | 13 | use tracer_mod, only: nqmx, igcm_dust_mass, igcm_dust_number, |
|---|
| 14 | & igcm_h2o_ice, igcm_ccn_mass, radius, |
|---|
| 15 | & igcm_ccn_number, nuice_ref, varian, |
|---|
| [1974] | 16 | & ref_r0, igcm_dust_submicron, |
|---|
| 17 | & igcm_stormdust_mass,igcm_stormdust_number |
|---|
| [1246] | 18 | USE dimradmars_mod, only: nueffdust,naerkind, |
|---|
| 19 | & name_iaer, |
|---|
| 20 | & iaer_dust_conrath,iaer_dust_doubleq, |
|---|
| [1974] | 21 | & iaer_dust_submicron,iaer_h2o_ice, |
|---|
| 22 | & iaer_stormdust_doubleq |
|---|
| [1969] | 23 | |
|---|
| [38] | 24 | IMPLICIT NONE |
|---|
| 25 | c======================================================================= |
|---|
| 26 | c subject: |
|---|
| 27 | c -------- |
|---|
| 28 | c Subroutine designed to update the aerosol size distribution used by |
|---|
| 29 | c the radiative transfer scheme. This size distribution is assumed |
|---|
| 30 | c to be a log-normal distribution, with effective radius "reffrad" and |
|---|
| 31 | c variance "nueffrad". |
|---|
| 32 | c At firstcall, "rice" and "nuice" are not known, because |
|---|
| 33 | c the H2O ice microphysical scheme is called after the radiative |
|---|
| 34 | c transfer in physiq.F. That's why we assess the size of the |
|---|
| 35 | c water-ice particles at firstcall (see part 1.2 below). |
|---|
| 36 | c |
|---|
| 37 | c author: |
|---|
| 38 | c ------ |
|---|
| 39 | c J.-B. Madeleine (2009-2010) |
|---|
| 40 | c |
|---|
| 41 | c======================================================================= |
|---|
| 42 | c |
|---|
| 43 | c Declarations : |
|---|
| 44 | c ------------- |
|---|
| 45 | c |
|---|
| [1969] | 46 | include "callkeys.h" |
|---|
| [38] | 47 | |
|---|
| 48 | c----------------------------------------------------------------------- |
|---|
| [1974] | 49 | c Inputs/outputs: |
|---|
| [38] | 50 | c ------ |
|---|
| 51 | |
|---|
| [1974] | 52 | INTEGER, INTENT(in) :: ngrid,nlayer |
|---|
| [38] | 53 | c Ice geometric mean radius (m) |
|---|
| [1974] | 54 | REAL, INTENT(out) :: rice(ngrid,nlayer) |
|---|
| [38] | 55 | c Estimated effective variance of the size distribution (n.u.) |
|---|
| [1974] | 56 | REAL, INTENT(out) :: nuice(ngrid,nlayer) |
|---|
| [38] | 57 | c Tracer mass mixing ratio (kg/kg) |
|---|
| [1974] | 58 | REAL, INTENT(in) :: pq(ngrid,nlayer,nqmx) |
|---|
| 59 | REAL, INTENT(out) :: rdust(ngrid,nlayer) ! Dust geometric mean radius (m) |
|---|
| 60 | REAL, INTENT(out) :: rstormdust(ngrid,nlayer) ! Dust geometric mean radius (m) |
|---|
| 61 | REAL, INTENT(in) :: pplay(ngrid,nlayer) ! altitude at the middle of the layers |
|---|
| 62 | REAL, INTENT(in) :: tau(ngrid,naerkind) |
|---|
| [38] | 63 | c Aerosol effective radius used for radiative transfer (meter) |
|---|
| [1974] | 64 | REAL, INTENT(out) :: reffrad(ngrid,nlayer,naerkind) |
|---|
| [38] | 65 | c Aerosol effective variance used for radiative transfer (n.u.) |
|---|
| [1974] | 66 | REAL, INTENT(out) :: nueffrad(ngrid,nlayer,naerkind) |
|---|
| 67 | REAL, INTENT(in) :: tauscaling(ngrid) ! Convertion factor for qccn and Nccn |
|---|
| 68 | |
|---|
| [38] | 69 | c Local variables: |
|---|
| 70 | c --------------- |
|---|
| 71 | |
|---|
| 72 | INTEGER :: ig,l ! 3D grid indices |
|---|
| 73 | INTEGER :: iaer ! Aerosol index |
|---|
| 74 | |
|---|
| 75 | c Number of cloud condensation nuclei near the surface |
|---|
| 76 | c (only used at firstcall). This value is taken from |
|---|
| 77 | c Montmessin et al. 2004 JGR 109 E10004 p5 (2E6 part m-3), and |
|---|
| 78 | c converted to part kg-1 using a typical atmospheric density. |
|---|
| 79 | |
|---|
| 80 | REAL, PARAMETER :: ccn0 = 1.3E8 |
|---|
| [629] | 81 | |
|---|
| 82 | c For microphysics only: |
|---|
| 83 | REAL Mo,No ! Mass and number of ccn |
|---|
| [1047] | 84 | REAL rhocloud(ngrid,nlayer) ! Cloud density (kg.m-3) |
|---|
| [38] | 85 | |
|---|
| [1224] | 86 | LOGICAL,SAVE :: firstcall=.true. |
|---|
| [38] | 87 | |
|---|
| 88 | REAL CBRT |
|---|
| 89 | EXTERNAL CBRT |
|---|
| 90 | |
|---|
| 91 | c================================================================== |
|---|
| [629] | 92 | c 1. Update radius from fields from dynamics or initial state |
|---|
| 93 | c================================================================== |
|---|
| [38] | 94 | |
|---|
| [358] | 95 | c 1.1 Dust particles |
|---|
| 96 | c ------------------ |
|---|
| 97 | IF (doubleq.AND.active) THEN |
|---|
| 98 | DO l=1,nlayer |
|---|
| 99 | DO ig=1, ngrid |
|---|
| [744] | 100 | call updaterdust(pq(ig,l,igcm_dust_mass), |
|---|
| 101 | & pq(ig,l,igcm_dust_number),rdust(ig,l)) |
|---|
| [358] | 102 | nueffdust(ig,l) = exp(varian**2.)-1. |
|---|
| 103 | ENDDO |
|---|
| 104 | ENDDO |
|---|
| 105 | ELSE |
|---|
| 106 | DO l=1,nlayer |
|---|
| 107 | DO ig=1, ngrid |
|---|
| 108 | rdust(ig,l) = 0.8E-6 |
|---|
| 109 | nueffdust(ig,l) = 0.3 |
|---|
| 110 | ENDDO |
|---|
| [38] | 111 | ENDDO |
|---|
| [358] | 112 | ENDIF |
|---|
| [1974] | 113 | |
|---|
| 114 | ! updating radius of stormdust particles |
|---|
| 115 | IF (rdstorm.AND.active) THEN |
|---|
| 116 | DO l=1,nlayer |
|---|
| 117 | DO ig=1, ngrid |
|---|
| 118 | call updaterdust(pq(ig,l,igcm_stormdust_mass), |
|---|
| 119 | & pq(ig,l,igcm_stormdust_number),rstormdust(ig,l)) |
|---|
| 120 | nueffdust(ig,l) = exp(varian**2.)-1. |
|---|
| 121 | ENDDO |
|---|
| 122 | ENDDO |
|---|
| 123 | ENDIF |
|---|
| [629] | 124 | |
|---|
| [358] | 125 | c 1.2 Water-ice particles |
|---|
| 126 | c ----------------------- |
|---|
| [744] | 127 | |
|---|
| 128 | IF (water.AND.activice) THEN |
|---|
| 129 | IF (microphys) THEN |
|---|
| [1208] | 130 | |
|---|
| 131 | c At firstcall, the true number and true mass of cloud condensation nuclei are not known. |
|---|
| 132 | c Indeed it is scaled on the prescribed dust opacity via a 'tauscaling' coefficient |
|---|
| 133 | c computed after radiative transfer. If tauscaling is not in startfi, we make an assumption for its value. |
|---|
| 134 | |
|---|
| [744] | 135 | IF (firstcall) THEN |
|---|
| [1974] | 136 | !IF (minval(tauscaling).lt.0) tauscaling(:) = 1.e-3 ! default value when non-read in startfi is -1 |
|---|
| 137 | !IF (freedust) tauscaling(:) = 1. ! if freedust, enforce no rescaling at all |
|---|
| [1208] | 138 | firstcall = .false. |
|---|
| 139 | ENDIF |
|---|
| 140 | |
|---|
| 141 | DO l=1,nlayer |
|---|
| 142 | DO ig=1,ngrid |
|---|
| 143 | call updaterice_micro(pq(ig,l,igcm_h2o_ice), |
|---|
| 144 | & pq(ig,l,igcm_ccn_mass), |
|---|
| 145 | & pq(ig,l,igcm_ccn_number), |
|---|
| 146 | & tauscaling(ig),rice(ig,l), |
|---|
| 147 | & rhocloud(ig,l)) |
|---|
| 148 | nuice(ig,l) = nuice_ref |
|---|
| [358] | 149 | ENDDO |
|---|
| [1208] | 150 | ENDDO |
|---|
| [744] | 151 | |
|---|
| 152 | ELSE ! if not microphys |
|---|
| 153 | |
|---|
| 154 | DO l=1,nlayer |
|---|
| 155 | DO ig=1,ngrid |
|---|
| 156 | call updaterice_typ(pq(ig,l,igcm_h2o_ice), |
|---|
| 157 | & tau(ig,1),pplay(ig,l),rice(ig,l)) |
|---|
| 158 | nuice(ig,l) = nuice_ref |
|---|
| [629] | 159 | ENDDO |
|---|
| [744] | 160 | ENDDO |
|---|
| 161 | |
|---|
| 162 | ENDIF ! of if microphys |
|---|
| 163 | ENDIF ! of if (water.AND.activice) |
|---|
| [38] | 164 | |
|---|
| 165 | c================================================================== |
|---|
| 166 | c 2. Radius used in the radiative transfer code (reffrad) |
|---|
| 167 | c================================================================== |
|---|
| 168 | |
|---|
| 169 | DO iaer = 1, naerkind ! Loop on aerosol kind |
|---|
| 170 | aerkind: SELECT CASE (name_iaer(iaer)) |
|---|
| 171 | c================================================================== |
|---|
| 172 | CASE("dust_conrath") aerkind ! Typical dust profile |
|---|
| 173 | c================================================================== |
|---|
| 174 | DO l=1,nlayer |
|---|
| 175 | DO ig=1,ngrid |
|---|
| [358] | 176 | reffrad(ig,l,iaer) = rdust(ig,l) * |
|---|
| 177 | & (1.e0 + nueffdust(ig,l))**2.5 |
|---|
| [38] | 178 | nueffrad(ig,l,iaer) = nueffdust(ig,l) |
|---|
| 179 | ENDDO |
|---|
| 180 | ENDDO |
|---|
| 181 | c================================================================== |
|---|
| 182 | CASE("dust_doubleq") aerkind! Two-moment scheme for dust |
|---|
| 183 | c================================================================== |
|---|
| 184 | DO l=1,nlayer |
|---|
| 185 | DO ig=1,ngrid |
|---|
| [358] | 186 | reffrad(ig,l,iaer) = rdust(ig,l) * ref_r0 |
|---|
| [38] | 187 | nueffrad(ig,l,iaer) = nueffdust(ig,l) |
|---|
| 188 | ENDDO |
|---|
| 189 | ENDDO |
|---|
| 190 | c================================================================== |
|---|
| 191 | CASE("dust_submicron") aerkind ! Small dust population |
|---|
| 192 | c================================================================== |
|---|
| 193 | DO l=1,nlayer |
|---|
| 194 | DO ig=1,ngrid |
|---|
| 195 | reffrad(ig,l,iaer)=radius(igcm_dust_submicron) |
|---|
| 196 | nueffrad(ig,l,iaer)=0.03 |
|---|
| 197 | ENDDO |
|---|
| 198 | ENDDO |
|---|
| 199 | c================================================================== |
|---|
| 200 | CASE("h2o_ice") aerkind ! Water ice crystals |
|---|
| 201 | c================================================================== |
|---|
| 202 | DO l=1,nlayer |
|---|
| 203 | DO ig=1,ngrid |
|---|
| [358] | 204 | c About reffice, do not confuse the mass mean radius |
|---|
| 205 | c (rayon moyen massique) and the number median radius |
|---|
| 206 | c (or geometric mean radius, rayon moyen géométrique). |
|---|
| 207 | c rice is a mass mean radius, whereas rdust |
|---|
| 208 | c is a geometric mean radius: |
|---|
| 209 | c number median rad = mass mean rad x exp(-1.5 sigma0^2) |
|---|
| 210 | c (Montmessin et al. 2004 paragraph 30). Therefore: |
|---|
| [38] | 211 | reffrad(ig,l,iaer)=rice(ig,l)*(1.+nuice_ref) |
|---|
| 212 | nueffrad(ig,l,iaer)=nuice_ref |
|---|
| 213 | ENDDO |
|---|
| 214 | ENDDO |
|---|
| 215 | c================================================================== |
|---|
| [1974] | 216 | CASE("stormdust_doubleq") aerkind! Two-moment scheme for |
|---|
| 217 | c stormdust; same distribution than normal dust |
|---|
| 218 | c================================================================== |
|---|
| 219 | DO l=1,nlayer |
|---|
| 220 | DO ig=1,ngrid |
|---|
| 221 | reffrad(ig,l,iaer) = rstormdust(ig,l) * ref_r0 |
|---|
| 222 | nueffrad(ig,l,iaer) = nueffdust(ig,l) |
|---|
| 223 | ENDDO |
|---|
| 224 | ENDDO |
|---|
| 225 | c================================================================== |
|---|
| [38] | 226 | END SELECT aerkind |
|---|
| 227 | ENDDO ! iaer (loop on aerosol kind) |
|---|
| 228 | |
|---|
| [1969] | 229 | END SUBROUTINE updatereffrad |
|---|
| 230 | |
|---|
| 231 | END MODULE updatereffrad_mod |
|---|
