Changeset 515 for trunk/MESOSCALE/LMD_MM_MARS/SRC

Timestamp:

Feb 4, 2012, 10:29:01 AM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: new GCM surface fields. new test case: crater. a flag for the Titus cap.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Files:

• mars_lmd/libf/phymars/physiq.F (modified) (2 diffs)
• phys/module_lmd_driver.F (modified) (1 diff)

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

@@ -330 +330 @@
       integer igmin, lmin
       logical tdiag
+      logical tituscap
 
       REAL zplev(ngrid,nlayermx+1),zplay(ngrid,nlayermx)
@@ -1063 +1064 @@
 c   -------------------------------------------
 
-      !!! get the actual co2 seasonal cap from Titus observations
-      CALL geticecover( ngrid, 180.*zls/pi, 
+      tituscap = .false.
+      IF (tituscap) THEN
+         PRINT*,'TITUS CAP !!!'
+         !!! get the actual co2 seasonal cap from Titus observations
+         CALL geticecover( ngrid, 180.*zls/pi,
      .                  180.*long/pi, 180.*lati/pi, co2ice )
-      co2ice = co2ice * 10000.
+         co2ice = co2ice * 10000.
+      ENDIF
 
       IF(callcond) THEN

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

@@ -1211 +1211 @@
 subs = (j-jps)*(ipe-ips+1)+(i-ips+1)
 
-!!! SEMBLE T IL PROBLEME AVEC NEWPHYS .... qui marche tres bien sans !
-#ifndef NEWPHYS           
-            !!! sparadrap pour regler un probleme avec mpiifort en LES
-            !!! -- HFX apparaissait nul aux interfaces des tuiles
-            if ( output_tab2d(subs,ind_HFX) .le. 1.e-8 ) then
-               !print *, 'HFX is zero !!! ', i, j
-               !print *, 'I substituted the value right next to it ', output_tab2d(subs+1,ind_HFX) 
-               output_tab2d(subs,ind_HFX) = output_tab2d(subs+1,ind_HFX)
-            endif
-#endif
+!!!! SEMBLE T IL PROBLEME AVEC NEWPHYS .... qui marche tres bien sans !
+!#ifndef NEWPHYS           
+!            !!! sparadrap pour regler un probleme avec mpiifort en LES
+!            !!! -- HFX apparaissait nul aux interfaces des tuiles
+!            if ( output_tab2d(subs,ind_HFX) .le. 1.e-8 ) then
+!               !print *, 'HFX is zero !!! ', i, j
+!               !print *, 'I substituted the value right next to it ', output_tab2d(subs+1,ind_HFX) 
+!               output_tab2d(subs,ind_HFX) = output_tab2d(subs+1,ind_HFX)
+!            endif
+!#endif
 
 #include "module_lmd_driver_output3.inc"