Changeset 228 for trunk/MESOSCALE/LMD_MM_MARS

Timestamp:

Jul 15, 2011, 8:33:51 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: corrected the diagfi bug. updated test case and check those are still working OK. other changes are cosmetic

Location:

trunk/MESOSCALE/LMD_MM_MARS/SIMU

Files:

• DEF/newphys_test/callphys.def (modified) (6 diffs)
• DEF/newphys_test/namelist.input (modified) (2 diffs)
• DEF/newphys_test/namelist.wps (deleted)
• DEF/newphys_test/nougaro_facteur_0.75 (deleted)
• runmeso (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/callphys.def

@@ -1 +1 @@
 ##General options
 ##~~~~~~~~~~~~~~~
+#Directory where external input files are
+datadir=/u/forget/WWW/datagcm/datafile/
+
 #Run with or without tracer transport ?
 tracer=.false.
@@ -14 +17 @@
 
 #Save statistics in file "stats.nc" ?
-stats =.false.
+callstats =.false.
 
 #Save EOF profiles in file "profiles" for Climate Database?
@@ -26 +29 @@
 #  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
 #  =26 Mars Year 26 from TES assimilation
-iaervar = 24 
+iaervar = 1 
 # Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
 tauvis=0.2
@@ -41 +44 @@
 callrad   = .true.
 # call NLTE radiative schemes ?   matters only if callrad=T
-callnlte = .true.
+callnlte = .false.
 # call CO2 NIR absorption ?   matters only if callrad=T
 callnirco2 = .true.
@@ -58 +61 @@
 ## ~~~~~~~~~~~~~~~~~~~~~~~~~~
 # the rad.transfer is computed every "iradia" physical timestep
-iradia    = 1
+iradia    = 10
 # Output of the exchange coefficient mattrix ? for diagnostic only
 callg2d  = .false.
@@ -89 +92 @@
 #         is a cold trap)
 caps  = .false.  
-#caps  = .true.
 # PHOTOCHEMISTRY: include chemical species 
 photochem  = .false.

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

@@ -2 +2 @@
  start_year       = 2024,    !! (p1) Start Martian Year (20XX for MY XX)
  start_month      = 01,      !! (p1) Start Martian Month
- start_day        = 01,      !! (p1) Start Martian Day
- start_hour       = 06,      !! (p1) Start Martian Hour (at longitude 0)
+ start_day        = 17,      !! (p1) Start Martian Day
+ start_hour       = 02,      !! (p1) Start Martian Hour (at longitude 0)
  end_year         = 2024,    !! (p1) End Martian Year (20XX for MY XX)
  end_month        = 01,      !! (p1) End Martian Month
- end_day          = 04,      !! (p1) End Martian Day
- end_hour         = 06,      !! (p1) End Martian Hour (at longitude 0)
- history_interval    = 37,     
- frames_per_outfile  = 24, 
- restart          = .false.
- restart_interval = 8880,
- io_form_history  = 2
- io_form_restart  = 2
- io_form_input    = 2
- io_form_boundary = 2
- debug_level      = 0
+ end_day          = 17,      !! (p1) End Martian Day
+ end_hour         = 12,      !! (p1) End Martian Hour (at longitude 0)
+ history_interval    = 37,   !! Frequency of outputs (37 --> 3700s = 1 Martian hour)  
+ frames_per_outfile  = 24,   !! Size of time dimension in files 
+ restart          = .false.  !! (*) Output restart files ?
+ restart_interval = 8880     !! (*) Frequency of output restart files ?
+ io_form_history  = 2        !! (*) Choice of NETCDF for ouputs
+ io_form_restart  = 2        !! (*) Choice of NETCDF for ouputs     
+ io_form_input    = 2        !! (*) Choice of NETCDF for ouputs
+ io_form_boundary = 2        !! (*) Choice of NETCDF for ouputs
+ debug_level      = 0        !! (*) Verbose level
  /
 
+
  &domains
- time_step   = 25, 
- dx = 20000,
- dy = 20000,
- e_we   =  61,               
- e_sn   =  61,              
- e_vert =  61,           !! (r)(p2) Number of vertical levels
- p_top_requested = 5.    !! (p3) Chosen value of pressure at the top of the model
- num_metgrid_levels = 19
+ time_step   = 50            !! Dynamical timestep
+ dx = 20000,                 !! (p2) Horizontal resolution
+ dy = 20000,                 !! (p2) Horizontal resolution (should be equal to dx)
+ e_we   = 61,                !! (r)(p2) Number of longitude grid points
+ e_sn   = 61,                !! (r)(p2) Number of latitude grid points
+ e_vert = 61,                !! (r)(p2) Number of vertical levels
+ p_top_requested = 5         !! (p3) Chosen value of pressure at the top of the model
+ !!
+ !! OPTIONAL
+ !!
+ num_metgrid_levels = 33     !! (p1) number of vertical levels in GCM inputs (+1)
  /
 
  &physics
- !init_TI = 85.,           !! Define constant thermal inertia value
+ !!
+ !! OPTIONAL
+ !!
  num_soil_layers = 18 !!! NOUVELLE PHYSIQUE
+ init_TI = 100.,             !! Define constant thermal inertia value
+ init_AL = 0.27,             !! Define constant albedo value
+ /
+
+ &dynamics
+ !!
+ !! OPTIONAL
+ !!
+ non_hydrostatic = F,        !! Integrate in non-hydrostatic/hydrostatic mode
+ tiso = 100.
+ /
+
+ &bdy_control
+ !!
+ !! OPTIONAL
+ !!
+ spec_bdy_width = 3          !! (p3) Width of transition zone with specified=T
+                             !!      NB: spec_bdy_width must equal relax_zone+1
+ relax_zone = 2              !! (p3) Width of relaxation zone with specified=T
+ /
+
+ &grib2
  /
 
@@ -38 +66 @@
  /
 
- &dynamics
- tiso = 100.
- /
-
- &grib2
- /
-
- &bdy_control
- /
-
-!! DO NOT MODIFY
- &namelist_quilt
- nio_tasks_per_group = 0,
- nio_groups = 1,
- /
-!! DO NOT MODIFY
+ &namelist_quilt              !! (*)
+ nio_tasks_per_group = 0,     !! (*)
+ nio_groups = 1,              !! (*)
+ /                            !! (*)

trunk/MESOSCALE/LMD_MM_MARS/SIMU/runmeso

@@ -384 +384 @@
    if [[ ${newphys} -eq 1 ]]
    then
-     echo 'n' >> ${simu_folder}/makemeso_input
      echo ${scat} >> ${simu_folder}/makemeso_input
+     #echo 'n' >> ${simu_folder}/makemeso_input
    fi
    #