Changeset 332 for trunk/LMDZ.MARS/libf/phymars

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/LMDZ.MARS/libf/phymars

Files:

• geticecover.F90 (added)
• physiq.F (modified) (1 diff)

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,