Changeset 2562 for trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

Timestamp:

Sep 21, 2021, 9:02:49 AM (3 years ago)

Author:

cmathe

Message:

GCM MARS: CO2 clouds microphysics improvements

File:

• trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

@@ -8 +8 @@
 subroutine phyetat0 (fichnom,tab0,Lmodif,nsoil,ngrid,nlay,nq, &
                      day_ini,time0,tsurf,tsoil,albedo,emis,q2,qsurf,co2ice, &
-                     tauscaling,totcloudfrac,wstar,mem_Mccn_co2,mem_Nccn_co2,&
-                     mem_Mh2o_co2,watercap)
+                     tauscaling,totcloudfrac,wstar,watercap)
 
   use tracer_mod, only: noms ! tracer names
@@ -28 +27 @@
 !======================================================================
 ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
-!  Adaptation à Mars : Yann Wanherdrick 
+!  Adaptation � Mars : Yann Wanherdrick 
 ! Objet: Lecture de l etat initial pour la physique
 ! Modifs: Aug.2010 EM : use NetCDF90 to load variables (enables using
@@ -66 +65 @@
   real,intent(out) :: totcloudfrac(ngrid) ! total cloud fraction
   real,intent(out) :: wstar(ngrid) ! Max vertical velocity in thermals (m/s)
-  real,intent(out) :: mem_Mccn_co2(ngrid,nlay) ! Memory of CCN mass of H2O and dust used by CO2
-  real,intent(out) :: mem_Nccn_co2(ngrid,nlay) ! Memory of CCN number of H2O and dust used by CO2
-  real,intent(out) :: mem_Mh2o_co2(ngrid,nlay) ! Memory of H2O mass integred into CO2 crystal
   real,intent(out) :: watercap(ngrid) ! h2o_ice_cover
 !======================================================================
@@ -363 +359 @@
             minval(co2ice), maxval(co2ice)
 
-! Memory of the origin of the co2 particles
-if (startphy_file) then
-   call get_field("mem_Mccn_co2",mem_Mccn_co2,found,indextime)
-   if (.not.found) then
-     write(*,*) "phyetat0: <mem_Mccn_co2> not in file"
-     mem_Mccn_co2(:,:)=0
-   endif
-else
-   mem_Mccn_co2(:,:)=0
-endif !if (startphy_file)
-write(*,*) "phyetat0: Memory of CCN mass of H2O and dust used by CO2"
-write(*,*) " <mem_Mccn_co2> range:", &
-             minval(mem_Mccn_co2), maxval(mem_Mccn_co2)
-
-if (startphy_file) then
-   call get_field("mem_Nccn_co2",mem_Nccn_co2,found,indextime)
-   if (.not.found) then
-     write(*,*) "phyetat0: <mem_Nccn_co2> not in file"
-     mem_Nccn_co2(:,:)=0
-   endif
-else
-  mem_Nccn_co2(:,:)=0
-endif ! if (startphy_file)
-write(*,*) "phyetat0: Memory of CCN number of H2O and dust used by CO2"
-write(*,*) " <mem_Nccn_co2> range:", &
-             minval(mem_Nccn_co2), maxval(mem_Nccn_co2)
-
-if (startphy_file) then
-   call get_field("mem_Mh2o_co2",mem_Mh2o_co2,found,indextime)
-   if (.not.found) then
-     write(*,*) "phyetat0: <mem_Mh2o_co2> not in file"
-     mem_Mh2o_co2(:,:)=0
-   endif
-else
-   mem_Mh2o_co2(:,:)=0
-endif ! if (startphy_file)
-write(*,*) "phyetat0: Memory of H2O mass integred into CO2 crystal"
-write(*,*) " <mem_Mh2o_co2> range:", &
-             minval(mem_Mh2o_co2), maxval(mem_Mh2o_co2)
-
 ! Dust conversion factor
 if (startphy_file) then