Changeset 515


Ignore:
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:
5 added
4 edited
4 copied

Legend:

Unmodified
Added
Removed
  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars11_radact_therm

    r389 r515  
    1010#Directory where external input files are
    1111datadir = /san/home/spiga/DATAFILE/
     12datadir = /d5/aslmd/datafile/
    1213
    1314#Run with or without tracer transport ?
  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/namelist.input

    r346 r515  
    2020
    2121 &domains
    22  time_step   = 15            !! Dynamical timestep 
     22 time_step   = 20            !! Dynamical timestep 
    2323 dx = 30000,                 !! (p2) Horizontal resolution
    2424 dy = 30000,                 !! (p2) Horizontal resolution (should be equal to dx)
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/physiq.F

    r341 r515  
    330330      integer igmin, lmin
    331331      logical tdiag
     332      logical tituscap
    332333
    333334      REAL zplev(ngrid,nlayermx+1),zplay(ngrid,nlayermx)
     
    10631064c   -------------------------------------------
    10641065
    1065       !!! get the actual co2 seasonal cap from Titus observations
    1066       CALL geticecover( ngrid, 180.*zls/pi,
     1066      tituscap = .false.
     1067      IF (tituscap) THEN
     1068         PRINT*,'TITUS CAP !!!'
     1069         !!! get the actual co2 seasonal cap from Titus observations
     1070         CALL geticecover( ngrid, 180.*zls/pi,
    10671071     .                  180.*long/pi, 180.*lati/pi, co2ice )
    1068       co2ice = co2ice * 10000.
     1072         co2ice = co2ice * 10000.
     1073      ENDIF
    10691074
    10701075      IF(callcond) THEN
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

    r503 r515  
    12111211subs = (j-jps)*(ipe-ips+1)+(i-ips+1)
    12121212
    1213 !!! SEMBLE T IL PROBLEME AVEC NEWPHYS .... qui marche tres bien sans !
    1214 #ifndef NEWPHYS           
    1215             !!! sparadrap pour regler un probleme avec mpiifort en LES
    1216             !!! -- HFX apparaissait nul aux interfaces des tuiles
    1217             if ( output_tab2d(subs,ind_HFX) .le. 1.e-8 ) then
    1218                !print *, 'HFX is zero !!! ', i, j
    1219                !print *, 'I substituted the value right next to it ', output_tab2d(subs+1,ind_HFX)
    1220                output_tab2d(subs,ind_HFX) = output_tab2d(subs+1,ind_HFX)
    1221             endif
    1222 #endif
     1213!!!! SEMBLE T IL PROBLEME AVEC NEWPHYS .... qui marche tres bien sans !
     1214!#ifndef NEWPHYS           
     1215!            !!! sparadrap pour regler un probleme avec mpiifort en LES
     1216!            !!! -- HFX apparaissait nul aux interfaces des tuiles
     1217!            if ( output_tab2d(subs,ind_HFX) .le. 1.e-8 ) then
     1218!               !print *, 'HFX is zero !!! ', i, j
     1219!               !print *, 'I substituted the value right next to it ', output_tab2d(subs+1,ind_HFX)
     1220!               output_tab2d(subs,ind_HFX) = output_tab2d(subs+1,ind_HFX)
     1221!            endif
     1222!#endif
    12231223
    12241224#include "module_lmd_driver_output3.inc" 
Note: See TracChangeset for help on using the changeset viewer.