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

Timestamp:

Mar 19, 2020, 12:30:06 PM (5 years ago)

Author:

jnaar

Message:

[MARS GCM] Adding watercap in start2archive. When absent in the start files, it is initialized at 0,
and qsurf(igcm_h2o_ice) can no longer be negative. The negative fraction of qsurf(h2o_ice),
if it exists, is set in watercap.
JN

File:

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

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

@@ -64 +64 @@
       REAL tsoil(ngridmx,nsoilmx) ! Soil temperature
       REAL co2ice(ngridmx)        ! CO2 ice layer
+      REAL watercap(ngridmx)        ! h2o ice layer
       REAL tauscaling(ngridmx) ! dust conversion factor
       REAL totcloudfrac(ngridmx) ! sub-grid cloud fraction
@@ -85 +86 @@
       REAL ithS(ip1jmp1,nsoilmx) ! Soil Thermal Inertia
       REAL co2iceS(ip1jmp1)
+      REAL watercapS(ip1jmp1)
       REAL tauscalingS(ip1jmp1)
       REAL totcloudfracS(ip1jmp1)
@@ -165 +167 @@
      &      day_ini_fi,timefi,tsurf,tsoil,albedo,emis,q2,qsurf,co2ice,
      &      tauscaling,totcloudfrac,wstar,mem_Mccn_co2,mem_Nccn_co2,
-     &               mem_Mh2o_co2)
+     &               mem_Mh2o_co2,watercap)
 
        ierr = NF_OPEN (fichnom, NF_NOWRITE,nid1)
@@ -259 +261 @@
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,tsurf,tsurfS)
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,co2ice,co2iceS)
+      call gr_fi_dyn(1,ngridmx,iip1,jjp1,watercap,watercapS)
       call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,tsoil,tsoilS)
       ! Note: thermal inertia "inertiedat" is in comsoil.h
@@ -355 +358 @@
       call write_archive(nid,ntime,'co2ice','couche de glace co2',
      &  'kg/m2',2,co2iceS)
+      call write_archive(nid,ntime,'watercap','couche de glace h2o',
+     &  'kg/m2',2,watercapS)
       call write_archive(nid,ntime,'tauscaling',
      &  'dust conversion factor',' ',2,tauscalingS)