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

Timestamp:

Mar 13, 2020, 7:32:42 PM (5 years ago)

Author:

jnaar

Message:

[MARS GCM] Adding a new 2D variable "watercap", perennial ground water ice.
qsurf(ig,igcm_h2o_ice) can no longer be negative.
When watercaptag=true, ie when there is an infinite amount of water avalaible for sublimation,
and there is no more h2o frost, sublimation digs into watercap reservoir.
For now watercap can't grow, even on watercaptag, and has same albedo than qsurf(h2o_ice).
When using old start file, watercap is automatically initiated as 0.
Added watercap in newstart.F aswell.
JN

File:

• trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F (modified) (4 diffs)

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

@@ -28 +28 @@
       use surfdat_h, only: phisfi, z0, zmea, zstd, zsig, zgam, zthe,
      &                     albedodat, z0_default, qsurf, tsurf,
-     &                     co2ice, emis, hmons, summit, base
+     &                     co2ice, emis, hmons, summit, base, watercap
       use comsoil_h, only: inertiedat, layer, mlayer, nsoilmx, tsoil
       use control_mod, only: day_step, iphysiq, anneeref, planet_type
@@ -418 +418 @@
      &   date,tsurf,tsoil,emis,q2,
      &   t,ucov,vcov,ps,co2ice,teta,phisold_newgrid,q,qsurf,
-     &   tauscaling,totcloudfrac,surfith,nid)
+     &   tauscaling,totcloudfrac,surfith,nid,watercap)
         write(*,*) "OK, read start_archive file"
         ! copy soil thermal inertia
@@ -440 +440 @@
      &        day_ini,time,tsurf,tsoil,albedo,emis,
      &        q2,qsurf,co2ice,tauscaling,totcloudfrac,
-     &        wstar,mem_Mccn_co2,mem_Nccn_co2,mem_Mh2o_co2)
+     &        wstar,mem_Mccn_co2,mem_Nccn_co2,mem_Mh2o_co2,watercap)
         
         ! copy albedo and soil thermal inertia
@@ -1752 +1752 @@
      &              tsurf,tsoil,co2ice,albedo,emis,q2,qsurf,tauscaling,
      &              totcloudfrac,wstar,
-     &              mem_Mccn_co2,mem_Nccn_co2,mem_Mh2o_co2)
+     &              mem_Mccn_co2,mem_Nccn_co2,mem_Mh2o_co2,watercap)
 
 c=======================================================================