Changeset 515 for trunk

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

Files:

• SIMU/DEF/TESTCRATER.tar.gz (added)
• SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars11_radact_therm (modified) (1 diff)
• SIMU/DEF/newphys_THARSIS_WATER/namelist.input (modified) (1 diff)
• SRC/WRFV2/mars_lmd/libf/phymars/physiq.F (modified) (2 diffs)
• SRC/WRFV2/phys/module_lmd_driver.F (modified) (1 diff)
• WPS_GEOG/albedo_GCM/00001-00180.00001-00180 (added)
• WPS_GEOG/albedo_GCM/00181-00360.00001-00180 (added)
• WPS_GEOG/albedo_GCM/old.00001-00180.00001-00180 (copied) (copied from trunk/MESOSCALE/LMD_MM_MARS/WPS_GEOG/albedo_GCM/00001-00180.00001-00180)
• WPS_GEOG/albedo_GCM/old.00181-00360.00001-00180 (copied) (copied from trunk/MESOSCALE/LMD_MM_MARS/WPS_GEOG/albedo_GCM/00181-00360.00001-00180)
• WPS_GEOG/thermal_GCM/00001-00180.00001-00180 (added)
• WPS_GEOG/thermal_GCM/00181-00360.00001-00180 (added)
• WPS_GEOG/thermal_GCM/old.00001-00180.00001-00180 (copied) (copied from trunk/MESOSCALE/LMD_MM_MARS/WPS_GEOG/thermal_GCM/00001-00180.00001-00180)
• WPS_GEOG/thermal_GCM/old.00181-00360.00001-00180 (copied) (copied from trunk/MESOSCALE/LMD_MM_MARS/WPS_GEOG/thermal_GCM/00181-00360.00001-00180)

trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars11_radact_therm

@@ -10 +10 @@
 #Directory where external input files are
 datadir = /san/home/spiga/DATAFILE/ 
+datadir = /d5/aslmd/datafile/
 
 #Run with or without tracer transport ?

trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/namelist.input

@@ -20 +20 @@
 
  &domains
- time_step   = 15            !! Dynamical timestep  
+ time_step   = 20            !! Dynamical timestep  
  dx = 30000,                 !! (p2) Horizontal resolution
  dy = 30000,                 !! (p2) Horizontal resolution (should be equal to dx)

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"