Changeset 2252
- Timestamp:
- Mar 6, 2020, 5:46:50 PM (5 years ago)
- Location:
- trunk/LMDZ.MARS
- Files:
-
- 6 edited
-
README (modified) (1 diff)
-
libf/phymars/aeropacity_mod.F (modified) (3 diffs)
-
libf/phymars/callradite_mod.F (modified) (3 diffs)
-
libf/phymars/physiq_mod.F (modified) (2 diffs)
-
libf/phymars/rocketduststorm_mod.F90 (modified) (1 diff)
-
libf/phymars/topmons_mod.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/README
r2248 r2252 2872 2872 == 04/03/2020 == AB 2873 2873 Bug fixes in the slope winds parametrization and its interactions with rocketduststorm : now, it is considered that the sub-grid dust storm is decorrelated from the dust entrainment above the sub-grid scale topography, so that : 1) the background dust available for slope winds is updated after the rocketduststorm parametrization ; 2) the extra-heating on top of the mountains is calculated regarding topdust and background dust tracers' optical depths only (and not stormdust optical depth). 2874 2875 == 06/03/2020 == AB 2876 Bug fix following r2248 in aeropacity_mod and topmons_mod : since dsodust, dsords and dsotop are diagnostic physiq_mod variables, we don't want them to be reinitialized at each call of aeropacity_mod and topmons_mod, but we initialize them once and for all at the beginning of physiq_mod instead. -
trunk/LMDZ.MARS/libf/phymars/aeropacity_mod.F
r2246 r2252 89 89 REAL, INTENT(OUT) :: tau(ngrid,naerkind) 90 90 REAL, INTENT(OUT) :: aerosol(ngrid,nlayer,naerkind) 91 REAL, INTENT( OUT) :: dsodust(ngrid,nlayer)92 REAL, INTENT( OUT) :: dsords(ngrid,nlayer) !dso of stormdust93 REAL, INTENT( OUT) :: dsotop(ngrid,nlayer) !dso of topdust91 REAL, INTENT(INOUT) :: dsodust(ngrid,nlayer) 92 REAL, INTENT(INOUT) :: dsords(ngrid,nlayer) !dso of stormdust 93 REAL, INTENT(INOUT) :: dsotop(ngrid,nlayer) !dso of topdust 94 94 REAL, INTENT(INOUT) :: reffrad(ngrid,nlayer,naerkind) 95 95 REAL, INTENT(IN) :: QREFvis3d(ngrid,nlayer,naerkind) … … 169 169 INTEGER,SAVE :: naerdust ! number of dust scatterers 170 170 171 ! initializations 171 172 tau(1:ngrid,1:naerkind)=0 172 dsords(:,:)=0. !CW17: initialize dsords173 dsodust(:,:)=0.174 dsotop(:,:)=0.175 173 176 174 ! identify tracers … … 243 241 #endif 244 242 245 END IF 243 END IF ! end of if firstcall 246 244 247 245 c Vertical column optical depth at "odpref" Pa -
trunk/LMDZ.MARS/libf/phymars/callradite_mod.F
r2246 r2252 193 193 ! (direct comparison with TES) 194 194 REAL,INTENT(OUT) :: aerosol(ngrid,nlayer,naerkind) 195 REAL,INTENT( OUT) :: dsodust(ngrid,nlayer)195 REAL,INTENT(INOUT) :: dsodust(ngrid,nlayer) 196 196 REAL,INTENT(OUT) :: rdust(ngrid,nlayer) ! Dust geometric mean radius (m) 197 197 REAL,INTENT(OUT) :: rice(ngrid,nlayer) ! Ice geometric mean radius (m) … … 203 203 REAL,INTENT(IN) :: totstormfract(ngrid) ! dust storm mesh fraction 204 204 REAL,INTENT(OUT) :: rstormdust(ngrid,nlayer) ! Storm dust geometric mean radius (m) 205 REAL,INTENT( OUT) :: dsords(ngrid,nlayer) ! density scaled opacity for rocket dust storm dust205 REAL,INTENT(INOUT) :: dsords(ngrid,nlayer) ! density scaled opacity for rocket dust storm dust 206 206 207 207 c entrainment by slope wind … … 209 209 REAL, INTENT(IN) :: alpha_hmons(ngrid) ! sub-grid scale topography mesh fraction 210 210 REAL,INTENT(OUT) :: rtopdust(ngrid,nlayer) ! Topdust geometric mean radius (m) 211 REAL,INTENT( OUT) :: dsotop(ngrid,nlayer) ! density scaled opacity for topmons dust211 REAL,INTENT(INOUT) :: dsotop(ngrid,nlayer) ! density scaled opacity for topmons dust 212 212 213 213 c sub-grid scale water ice clouds -
trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
r2248 r2252 672 672 zdtsurf(:)=0 673 673 dqsurf(:,:)=0 674 674 dsodust(:,:)=0. 675 dsords(:,:)=0. 676 dsotop(:,:)=0. 677 675 678 #ifdef DUSTSTORM 676 679 pq_tmp(:,:,:)=0 … … 1124 1127 hsummit(:)=14000. 1125 1128 endif 1126 clearatm=.true. 1129 clearatm=.true. ! stormdust is not accounted in the extra heating on top of the mountains 1127 1130 nohmons=.false. 1128 1131 pdqtop(:,:,:)=0. -
trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90
r2246 r2252 92 92 REAL, INTENT(OUT) :: pdqrds(ngrid,nlayer,nq) ! tendancy field for dust when detraining 93 93 REAL, INTENT(OUT) :: wrad(ngrid,nlayer+1) ! vertical speed within the rocket dust storm 94 REAL, INTENT( OUT) :: dsodust(ngrid,nlayer) ! density scaled opacity of env. dust95 REAL, INTENT( OUT) :: dsords(ngrid,nlayer) ! density scaled opacity of storm dust96 REAL, INTENT( OUT) :: dsotop(ngrid,nlayer) ! density scaled opacity of topmons dust94 REAL, INTENT(INOUT) :: dsodust(ngrid,nlayer) ! density scaled opacity of env. dust 95 REAL, INTENT(INOUT) :: dsords(ngrid,nlayer) ! density scaled opacity of storm dust 96 REAL, INTENT(INOUT) :: dsotop(ngrid,nlayer) ! density scaled opacity of topmons dust 97 97 REAL, INTENT(OUT) :: tauref(ngrid) 98 98 -
trunk/LMDZ.MARS/libf/phymars/topmons_mod.F90
r2248 r2252 93 93 REAL, INTENT(OUT) :: wfin(ngrid,nlayer+1) ! final vertical wind velocity: combination of both wup and wrad 94 94 ! wfin < 0 means upward wind (in pressure levels/s) 95 REAL, INTENT(OUT) :: dsodust(ngrid,nlayer)96 95 97 96 ! output for second radiative transfer 98 97 REAL, INTENT(OUT) :: aerosol(ngrid,nlayer,naerkind) 99 REAL, INTENT(OUT) :: dsords(ngrid,nlayer) 100 REAL, INTENT(OUT) :: dsotop(ngrid,nlayer) 98 REAL, INTENT(INOUT) :: dsodust(ngrid,nlayer) 99 REAL, INTENT(INOUT) :: dsords(ngrid,nlayer) 100 REAL, INTENT(INOUT) :: dsotop(ngrid,nlayer) 101 101 REAL, INTENT(OUT) :: tauref(ngrid) 102 102 … … 210 210 ! ********************************************************************** 211 211 aerosol(:,:,:)=0. 212 dsodust(:,:)=0.213 dsords(:,:)=0.214 dsotop(:,:)=0.215 212 pdqtop(:,:,:) = 0. 216 213 dqvl_topdust_mass(:,:)=0.
Note: See TracChangeset
for help on using the changeset viewer.
