Changeset 332 for trunk/MESOSCALE/LMD_MM_MARS/SRC

Timestamp:

Oct 27, 2011, 12:57:37 AM (13 years ago)

Author:

aslmd

Message:

LMDZ.MARS and MESOSCALE: added evolutive CO2 seasonal caps according to Titus crocus lines

compiles fine. to be tested in mesoscale polar simulations.

27/10/11 == AS

--> Added geticecover.F90 which computes seasonal ice cover given ls, lati(ngrid), long(ngrid)

as proposed by T. Titus from TES observations [fitting functions for crocus line]
... output is icecover(ngrid) which value is 0 [no ice cover] or 1 [ice cover]
... no calculations are done for latitudes between -40 and +40 [ice cover is directly set to 0]

--> In physiq.F, co2ice is set to a dummy high value to simulate a CO2 cap

wherever icecover(ngrid) is 1. This is done at each timestep before newcondens is called.

--> For the moment this is MESOSCALE only, but potentially useful to everyone.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars

Files:

• geticecover.F90 (added)
• physiq.F (modified) (2 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/physiq.F

@@ -1072 +1072 @@
 c   -------------------------------------------
 
+      !!! get the actual co2 seasonal cap from Titus observations
+      CALL geticecover( ngrid, zls, long, lati, co2ice )
+      co2ice = co2ice * 10000.  !! a whole bunch of ice
+
       IF(callcond) THEN
          CALL newcondens(ngrid,nlayer,nq,ptimestep,
@@ -1079 +1083 @@
      $              zdtc,zdtsurfc,pdpsrf,zduc,zdvc,zdqc,
      $              fluxsurf_sw) 
+
 
          DO l=1,nlayer