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

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

Location:

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Files:

• newstart.F (modified) (3 diffs)
• start2archive.F (modified) (3 diffs)

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

@@ -30 +30 @@
       use comsoil_h, only: inertiedat, layer, mlayer, nsoilmx, tsoil
       use control_mod, only: day_step, iphysiq, anneeref, planet_type
+      use phyetat0_mod, only: phyetat0
       use phyredem, only: physdem0, physdem1
       use iostart, only: open_startphy
-      use dimradmars_mod, only: tauscaling
-      use turb_mod, only: q2
+      use dimradmars_mod, only: tauscaling, albedo
+      use turb_mod, only: q2, wstar
+      use co2cloud_mod, only: mem_Mccn_co2, mem_Mh2o_co2,
+     &                        mem_Nccn_co2
       use filtreg_mod, only: inifilr
       USE mod_const_mpi, ONLY: COMM_LMDZ
@@ -418 +421 @@
         fichnom = 'start.nc'
         CALL dynetat0(fichnom,vcov,ucov,teta,q,masse,
-     .       ps,phis,time)
+     &       ps,phis,time)
 
         write(*,*) 'Reading file STARTFI'
         fichnom = 'startfi.nc'
         CALL phyetat0 (fichnom,tab0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
-     .        day_ini,time,
-     .        tsurf,tsoil,emis,q2,qsurf,co2ice,tauscaling,totcloudfrac)
+     &        day_ini,time,tsurf,tsoil,albedo,emis,
+     &        q2,qsurf,co2ice,tauscaling,totcloudfrac,
+     &        wstar,mem_Mccn_co2,mem_Nccn_co2,mem_Mh2o_co2)
         
         ! copy albedo and soil thermal inertia
@@ -1640 +1644 @@
      .              airefi,albfi,ithfi,zmea,zstd,zsig,zgam,zthe)
       call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,
-     .              dtphys,hour_ini,
-     .              tsurf,tsoil,co2ice,emis,q2,qsurf,tauscaling,
-     .              totcloudfrac)
+     &              dtphys,hour_ini,
+     &              tsurf,tsoil,co2ice,albedo,emis,q2,qsurf,tauscaling,
+     &              totcloudfrac,wstar,
+     &              mem_Mccn_co2,mem_Nccn_co2,mem_Mh2o_co2)
 
 c=======================================================================

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)