Changeset 332

Timestamp:

Oct 27, 2011, 12:57:37 AM (14 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

Files:

• LMDZ.MARS/README (modified) (1 diff)
• LMDZ.MARS/libf/phymars/geticecover.F90 (added)
• LMDZ.MARS/libf/phymars/physiq.F (modified) (1 diff)
• MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/geticecover.F90 (added)
• MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/physiq.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -1096 +1096 @@
 == 21/10/2011 == AS
 - Added possibility for CH4 tracer in tracer.h and initracer.F
+
+== 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.

trunk/LMDZ.MARS/libf/phymars/physiq.F

@@ -862 +862 @@
 c   -------------------------------------------
 
+#ifdef MESOSCALE
+      !!!! get the actual co2 seasonal cap from Titus observations
+      CALL geticecover( ngrid, zls, long, lati, co2ice )
+      co2ice = co2ice * 10000.  !! a whole bunch of ice
+#endif
+
       IF (callcond) THEN
          CALL newcondens(ngrid,nlayer,nq,ptimestep,

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