Changeset 2324
- Timestamp:
- May 15, 2020, 3:24:17 PM (5 years ago)
- Location:
- trunk/LMDZ.MARS
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/README
r2323 r2324 3042 3042 == 14/05/2020 == MV 3043 3043 Follow-up of the last commit for the transport of the isotopic ratio: forgot the initialisation of zq0() in callsedim_mod.F ! 3044 3045 == 15/05/2020 == MV 3046 Follow-up of the last commit for the transport of the isotopic ratio: cleaning to help the user: 3047 -initracer.F: hdo must be defined as an isotope (meaning its "father" has to be precised in traceur.def) and must be placed at the end of the list of tracers in traceur.def. 3048 -/deftank/traceur.def.isotopes: writing of traceur.def in the case of hdo=true has been simplified to match both the 3D and the 1D simulations. 3049 -simpleclouds.F,hdo_surfex_mod.F,physiq_mod.F: thresholds of 1.e-16 have been replaced by the variable qperemin defined in tracer_mod.F. 3050 -
trunk/LMDZ.MARS/deftank/traceur.def.isotopes
r2296 r2324 1 1 7 2 10 10 co2 air 3 10 10 h2o_ice air 4 10 10 h2o_vap air 5 10 10 dust_number air6 10 10 dust_mass air 7 10 10hdo_ice h2o_ice8 10 10hdo_vap h2o_vap2 co2 3 h2o_ice 4 h2o_vap 5 dust_number 6 dust_mass 7 hdo_ice h2o_ice 8 hdo_vap h2o_vap -
trunk/LMDZ.MARS/libf/phymars/hdo_surfex_mod.F
r2316 r2324 11 11 12 12 use tracer_mod, only: igcm_h2o_vap, igcm_h2o_ice, 13 & igcm_hdo_vap, igcm_hdo_ice 14 13 & igcm_hdo_vap, igcm_hdo_ice, 14 & qperemin 15 15 use surfdat_h, only: watercaptag 16 16 use geometry_mod, only: longitude_deg,latitude_deg … … 70 70 if (h2oflux(ig).le.0.) then 71 71 72 if (pqsurf(ig,igcm_h2o_ice).gt. 1e-16) then72 if (pqsurf(ig,igcm_h2o_ice).gt.qperemin) then 73 73 pdqsdif(ig,igcm_hdo_ice) = 74 74 & pdqsdif(ig,igcm_h2o_ice)* … … 111 111 endif 112 112 113 if (old_h2o_vap(ig).gt. 1.e-16) then113 if (old_h2o_vap(ig).gt.qperemin) then 114 114 pdqsdif(ig,igcm_hdo_ice)= 115 115 & alpha_c(ig)*pdqsdif(ig,igcm_h2o_ice)* -
trunk/LMDZ.MARS/libf/phymars/initracer.F
r2321 r2324 433 433 write(*,*) "Tracer H2O must be initialized before HDO" 434 434 STOP 435 else if ((nqfils(igcm_h2o_ice).lt.1) 436 & .or. (nqfils(igcm_h2o_vap).lt.1)) then 437 write(*,*) "HDO must be transported as a son", 438 & " of H2O: complete traceur.def" 439 STOP 440 else if ((igcm_hdo_ice.lt.nq-2) 441 & .or. (igcm_hdo_vap.lt.nq-2)) then 442 write(*,*) "The isotopes (HDO) must be placed at", 443 & " the end of the list in traceur.def" 444 STOP 435 445 endif 436 446 endif -
trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
r2316 r2324 35 35 & igcm_he, igcm_stormdust_mass, 36 36 & igcm_stormdust_number, igcm_topdust_mass, 37 & igcm_topdust_number 37 & igcm_topdust_number, 38 & qperemin 38 39 use comsoil_h, only: inertiedat, ! soil thermal inertia 39 40 & tsoil, nsoilmx,!number of subsurface layers … … 2953 2954 do l=1,nlayer 2954 2955 do ig=1,ngrid 2955 if (zq(ig,l,igcm_h2o_vap).gt. 1e-16) then2956 if (zq(ig,l,igcm_h2o_vap).gt.qperemin) then 2956 2957 DoH_vap(ig,l) = ( zq(ig,l,igcm_hdo_vap)/ 2957 2958 & zq(ig,l,igcm_h2o_vap) )*1./(2.*155.76e-6) … … 2964 2965 do l=1,nlayer 2965 2966 do ig=1,ngrid 2966 if (zq(ig,l,igcm_h2o_ice).gt. 1e-16) then2967 if (zq(ig,l,igcm_h2o_ice).gt.qperemin) then 2967 2968 DoH_ice(ig,l) = ( zq(ig,l,igcm_hdo_ice)/ 2968 2969 & zq(ig,l,igcm_h2o_ice) )/(2.*155.76e-6) … … 3590 3591 C Calculation of the D/H ratio 3591 3592 do l=1,nlayer 3592 if (zq(1,l,igcm_h2o_vap).gt. 1e-16) then3593 if (zq(1,l,igcm_h2o_vap).gt.qperemin) then 3593 3594 DoH_vap(1,l) = 0.5*( zq(1,l,igcm_hdo_vap)/ 3594 3595 & zq(1,l,igcm_h2o_vap) )/155.76e-6 … … 3599 3600 3600 3601 do l=1,nlayer 3601 if (zq(1,l,igcm_h2o_ice).gt. 1e-16) then3602 if (zq(1,l,igcm_h2o_ice).gt.qperemin) then 3602 3603 DoH_ice(1,l) = 0.5*( zq(1,l,igcm_hdo_ice)/ 3603 3604 & zq(1,l,igcm_h2o_ice) )/155.76e-6 -
trunk/LMDZ.MARS/libf/phymars/simpleclouds.F
r2312 r2324 6 6 USE watersat_mod, ONLY: watersat 7 7 use tracer_mod, only: igcm_h2o_vap, igcm_h2o_ice, 8 & igcm_hdo_vap, igcm_hdo_ice 8 & igcm_hdo_vap, igcm_hdo_ice, 9 & qperemin 9 10 USE comcstfi_h 10 11 use dimradmars_mod, only: naerkind … … 177 178 endif 178 179 179 if (zq0(ig,l,igcm_h2o_vap).gt. 1.e-16) then180 if (zq0(ig,l,igcm_h2o_vap).gt.qperemin) then 180 181 pdqcloud(ig,l,igcm_hdo_ice)= 181 182 & pdqcloud(ig,l,igcm_h2o_ice)*alpha_c(ig,l)* … … 195 196 else ! sublimation 196 197 197 if (zq0(ig,l,igcm_h2o_ice).gt. 1.e-16) then198 if (zq0(ig,l,igcm_h2o_ice).gt.qperemin) then 198 199 pdqcloud(ig,l,igcm_hdo_ice)= 199 200 & pdqcloud(ig,l,igcm_h2o_ice)*
Note: See TracChangeset
for help on using the changeset viewer.