Changeset 1922 for trunk/LMDZ.MARS/libf/phymars/phyetat0.F90

Timestamp:

Apr 19, 2018, 3:16:47 PM (7 years ago)

Author:

emillour

Message:

Mars GCM:
CO2 code updates:

• make co2cloud a module and save mem_* variables (initialized via phys_state_var_init)
• make improvedCO2cloud a module
• read/write mem_* variables in phyetat0.F and phyredem.F

DB

File:

• trunk/LMDZ.MARS/libf/phymars/phyetat0.F90 (modified) (3 diffs)

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

@@ -1 +1 @@
 subroutine phyetat0 (fichnom,tab0,Lmodif,nsoil,ngrid,nlay,nq, &
                      day_ini,time0,tsurf,tsoil,emis,q2,qsurf,co2ice, &
-                     tauscaling,totcloudfrac)
-!  use netcdf
+                     tauscaling,totcloudfrac,mem_Mccn_co2,mem_Nccn_co2,&
+                     mem_Mh2o_co2)
+
   use tracer_mod, only: noms ! tracer names
   use surfdat_h, only: phisfi, albedodat, z0, z0_default,&
@@ -47 +48 @@
   real,intent(out) :: tauscaling(ngrid) ! dust conversion factor
   real,intent(out) :: totcloudfrac(ngrid) ! total cloud fraction
+  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
 !======================================================================
 !  Local variables:
@@ -230 +234 @@
 endif
 
+! Memory of the origin of the co2 particles
+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
+else
+  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)
+endif
+
+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
+else
+  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)
+endif
+
+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
+else
+  write(*,*) "phyetat0: Memory of H2O mass integred into CO2 crystal"
+  write(*,*) " <mem_Mh2o_co2> range:", &
+             minval(mem_Mh2o_co2), maxval(mem_Mh2o_co2)
+endif
 
 ! Dust conversion factor