Changeset 1944 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F

Timestamp:

Jun 5, 2018, 8:38:22 AM (7 years ago)

Author:

emillour

Message:

Mars GCM:
A first step towards 1+1=2 (for now only works without tracers):

• store and load "albedo" (surface albedo) and wstar (thermals' max vertical velocity) in physics (re)start file.
• turn phyetat0 into a module in the process.

EM

File:

• trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F

@@ -31 +31 @@
       USE temps_mod, ONLY: day_ini,hour_ini
       USE iniphysiq_mod, ONLY: iniphysiq
+      USE phyetat0_mod, ONLY: phyetat0
       implicit none
 
@@ -67 +68 @@
       REAL q2(ngridmx,llm+1)
       REAL emis(ngridmx)
+      REAL albedo(ngridmx,2)
+      REAL wstar(ngridmx)
+      DOUBLE PRECISION mem_Nccn_co2(ngridmx,llm)
+      DOUBLE PRECISION mem_Mccn_co2(ngridmx,llm)
+      DOUBLE PRECISION mem_Mh2o_co2(ngridmx,llm)
       INTEGER start,length
       PARAMETER (length = 100)
@@ -157 +163 @@
 
       CALL phyetat0 (fichnom,0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
-     &      day_ini_fi,timefi,tsurf,tsoil,emis,q2,qsurf,co2ice,
-     &      tauscaling,totcloudfrac)
+     &      day_ini_fi,timefi,tsurf,tsoil,albedo,emis,q2,qsurf,co2ice,
+     &      tauscaling,totcloudfrac,wstar,mem_Mccn_co2,mem_Nccn_co2,
+     &               mem_Mh2o_co2)
 
        ierr = NF_OPEN (fichnom, NF_NOWRITE,nid1)