Changeset 2324


Ignore:
Timestamp:
May 15, 2020, 3:24:17 PM (5 years ago)
Author:
mvals
Message:

Mars GCM:
Follow-up of the last commit for the transport of the isotopic ratio: cleaning to help the user:
-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.
-/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.
-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.
MV

Location:
trunk/LMDZ.MARS
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/README

    r2323 r2324  
    30423042== 14/05/2020 == MV
    30433043Follow-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
     3046Follow-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  
    117
    2 10 10 co2 air
    3 10 10 h2o_ice air
    4 10 10 h2o_vap air
    5 10 10 dust_number air
    6 10 10 dust_mass air
    7 10 10 hdo_ice h2o_ice
    8 10 10 hdo_vap h2o_vap
     2co2
     3h2o_ice
     4h2o_vap
     5dust_number
     6dust_mass
     7hdo_ice h2o_ice
     8hdo_vap h2o_vap
  • trunk/LMDZ.MARS/libf/phymars/hdo_surfex_mod.F

    r2316 r2324  
    1111
    1212      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
    1515      use surfdat_h, only: watercaptag
    1616      use geometry_mod, only: longitude_deg,latitude_deg
     
    7070            if (h2oflux(ig).le.0.) then
    7171
    72                if (pqsurf(ig,igcm_h2o_ice).gt.1e-16) then
     72               if (pqsurf(ig,igcm_h2o_ice).gt.qperemin) then
    7373                pdqsdif(ig,igcm_hdo_ice) =
    7474     &            pdqsdif(ig,igcm_h2o_ice)*
     
    111111               endif
    112112
    113                if (old_h2o_vap(ig).gt.1.e-16) then
     113               if (old_h2o_vap(ig).gt.qperemin) then
    114114                         pdqsdif(ig,igcm_hdo_ice)=
    115115     &                      alpha_c(ig)*pdqsdif(ig,igcm_h2o_ice)*
  • trunk/LMDZ.MARS/libf/phymars/initracer.F

    r2321 r2324  
    433433           write(*,*) "Tracer H2O must be initialized before HDO"
    434434           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
    435445        endif
    436446      endif
  • trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

    r2316 r2324  
    3535     &                      igcm_he, igcm_stormdust_mass,
    3636     &                      igcm_stormdust_number, igcm_topdust_mass,
    37      &                      igcm_topdust_number
     37     &                      igcm_topdust_number,
     38     &                      qperemin
    3839      use comsoil_h, only: inertiedat, ! soil thermal inertia
    3940     &                     tsoil, nsoilmx,!number of subsurface layers
     
    29532954            do l=1,nlayer
    29542955                do ig=1,ngrid
    2955                 if (zq(ig,l,igcm_h2o_vap).gt.1e-16) then
     2956                if (zq(ig,l,igcm_h2o_vap).gt.qperemin) then
    29562957                    DoH_vap(ig,l) = ( zq(ig,l,igcm_hdo_vap)/
    29572958     &              zq(ig,l,igcm_h2o_vap) )*1./(2.*155.76e-6)
     
    29642965            do l=1,nlayer
    29652966                do ig=1,ngrid
    2966                 if (zq(ig,l,igcm_h2o_ice).gt.1e-16) then
     2967                if (zq(ig,l,igcm_h2o_ice).gt.qperemin) then
    29672968                    DoH_ice(ig,l) = ( zq(ig,l,igcm_hdo_ice)/
    29682969     &                  zq(ig,l,igcm_h2o_ice) )/(2.*155.76e-6)
     
    35903591C           Calculation of the D/H ratio
    35913592            do l=1,nlayer
    3592                 if (zq(1,l,igcm_h2o_vap).gt.1e-16) then
     3593                if (zq(1,l,igcm_h2o_vap).gt.qperemin) then
    35933594                    DoH_vap(1,l) = 0.5*( zq(1,l,igcm_hdo_vap)/
    35943595     &                zq(1,l,igcm_h2o_vap) )/155.76e-6
     
    35993600
    36003601            do l=1,nlayer
    3601                 if (zq(1,l,igcm_h2o_ice).gt.1e-16) then
     3602                if (zq(1,l,igcm_h2o_ice).gt.qperemin) then
    36023603                    DoH_ice(1,l) = 0.5*( zq(1,l,igcm_hdo_ice)/
    36033604     &                  zq(1,l,igcm_h2o_ice) )/155.76e-6
  • trunk/LMDZ.MARS/libf/phymars/simpleclouds.F

    r2312 r2324  
    66      USE watersat_mod, ONLY: watersat
    77      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
    910      USE comcstfi_h
    1011      use dimradmars_mod, only: naerkind
     
    177178                endif
    178179               
    179                 if (zq0(ig,l,igcm_h2o_vap).gt.1.e-16) then
     180                if (zq0(ig,l,igcm_h2o_vap).gt.qperemin) then
    180181                  pdqcloud(ig,l,igcm_hdo_ice)=
    181182     &              pdqcloud(ig,l,igcm_h2o_ice)*alpha_c(ig,l)*
     
    195196          else  ! sublimation
    196197
    197              if (zq0(ig,l,igcm_h2o_ice).gt.1.e-16) then
     198             if (zq0(ig,l,igcm_h2o_ice).gt.qperemin) then
    198199                pdqcloud(ig,l,igcm_hdo_ice)=
    199200     &               pdqcloud(ig,l,igcm_h2o_ice)*
Note: See TracChangeset for help on using the changeset viewer.