Ignore:
Timestamp:
Aug 5, 2011, 4:48:44 AM (13 years ago)
Author:
aslmd
Message:

MESOSCALE: major commit at an absurd hour. new version for user manual finished.

Location:
trunk/MESOSCALE/LMD_MM_MARS
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/calendar

    r11 r262  
    1       init --   GCM sol --   GCM ls --     WRF
     1      init --   GCM sol --   GCM ls --     LMD Martian Mesoscale Model
    22            0            0    0.000000     2024-01-01_00:00:00
    33            0            1   0.5125165     2024-01-02_00:00:00
  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/namelist.input_full

    r258 r262  
    2323
    2424 &domains
    25  time_step   = 50            !! Dynamical timestep
     25 time_step   = 30            !! Dynamical timestep
    2626 dx = 20000,                 !! (p2) Horizontal resolution
    2727 dy = 20000,                 !! (p2) Horizontal resolution (should be equal to dx)
     
    3333 time_step_fract_num = 0     !! Additional fraction to time_step: numerator
    3434 time_step_fract_den = 1     !! Additional fraction to time_step: denominator
    35  num_metgrid_levels = 26     !! (p1) number of vertical levels in GCM inputs (+1)
     35 num_metgrid_levels = 26     !! (p1) Number of vertical levels in GCM inputs (+1)
    3636 force_sfc_in_vinterp = 8    !! (p3) Number of levels hardwired in the PBL
    37                              !!         NB: decrease this parameter when low p_top_requested
     37                             !!       NB: decrease this parameter if low p_top_requested
    3838 max_dz = 1500.              !! (p3) Maximal interval (m) between vertical levels
    3939 eta_levels = -1.            !! (p3)  Specify a list of e_vert eta levels
     
    4545 parent_grid_ratio      = 1, !! (n)(p2) Ratio of horizontal resolution parent/nest
    4646 parent_time_step_ratio = 1, !! (n) Ratio of time step parent/nest
     47 feedback               = 0  !! (n) Define one-way nesting [0] or two-way nesting [1]
     48 smooth_option          = 2  !! (n) Smoothing option for parent domain (feedback=1)
     49                             !!      0: none; 1: 1-2-1 smooth; 2: smooth-desmooth
     50 blend_width            = 5  !! (n) Extent of blending zone for parent/nest boundaries
    4751 /
    4852
     
    5054 !!!!! OPTIONAL !!!!!!!!!!!!!!!
    5155 radt = 1,                   !! Ratio between physical and dynamical time step
     56                             !! ... 1 is recommended, >5 is not recommended
    5257 mars =  0,                  !! (r)(p2) Configuration of tracers:
    5358                             !!         0: no tracers, 1: water vapor + ice, 2: dust
     
    5863 init_WX = 0.,               !! (p3) \  Ini/bdy wind profile is everywhere equal to
    5964 init_WY = 0.,               !! (p3) /  the wind profile @ grid pt (init_WX,init_WY)
    60  init_MU = 0.,               !! (p3) Multiply ini & bdy zonal wind by init_U
    61  init_MV = 0.,               !! (p3) Multiply ini & bdy meridional wind by init_V
    62  init_LES = F,               !! (p3) Uniform domain initialization for LES
    63  num_soil_layers = 10.       !! (r)(p1) Number of layers in soil model [set to 18 for newphys]
     65 init_MU = 0.,               !! (p3) Multiply ini & bdy zonal wind by init_MU
     66 init_MV = 0.,               !! (p3) Multiply ini & bdy meridional wind by init_MV
     67 num_soil_layers = 10.       !! (r)(p1) Number of layers in soil model
    6468 /
    6569
     
    7276 tiso = 140.,                !! (p3) Temperature at which base profile is isotherm
    7377 !!
    74  diff_opt = 1                !! (*d) Diffusion option [set to 0 if LES or GCM]
     78 diff_opt = 1                !! (*d) Diffusion option [set to 0 if GCM-like]
    7579 km_opt = 4                  !! (*d) Eddy coefficient option           
    76  diff_6th_opt = 2,           !! (*d) Knievel numerical diffusion [set to 0 if LES]     
     80 diff_6th_opt = 2,           !! (*d) Knievel numerical diffusion
    7781 diff_6th_factor = 0.2,      !! (*d) Knievel numerical coeff. [set to 0.5 if GCM]       
    7882 smdiv = 0.1,                !! (*d) Divergence damping
     
    96100 open_ye = F,                !! (p3) Open boundary conditions @ northern boundary
    97101 spec_bdy_width = 5          !! (p3) Width of transition zone with specified=T
    98                              !!      NB: spec_bdy_width must equal relax_zone+1
     102                             !!      (spec_bdy_width must be equal to relax_zone+1)
    99103 relax_zone = 4              !! (p3) Width of relaxation zone with specified=T
    100104 /
     
    111115 nio_groups = 1,              !! (*)
    112116 /                            !! (*)
     117 !!!!! DO NOT MODIFY !!!!!!!!!!
  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/namelist.input_minim

    r223 r262  
    2020
    2121 &domains
    22  time_step   = 50         
     22 time_step   = 30         
    2323 dx = 20000,               
    2424 dy = 20000,               
  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/namelist.input_nests

    r258 r262  
    3636 e_sn                   =    117,   121,  121,  !! IMPORTANT: e_sn[parent domain] = e_sn[child domains] - 4
    3737  !! WARNING: number of processors must divide e_we minus 1 and e_sn minus 1
    38  e_vert                 =     61,    61,   61,  !! Vertical levels must be duplicates         
     38 e_vert                 =     61,    61,   61,  !! Vertical levels must be duplicated
    3939 dx                     =  63000, 21000, 7000,  !! Compute dx[child] with parent_grid_ratio and dx[parent]
    4040 dy                     =  63000, 21000, 7000,  !! .... (keep 10^-2 accuracy if result of division is not int)
  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/namelist.wps_nests

    r258 r262  
    2020 dy = 63000,
    2121 geog_data_res     = '64ppd', '64ppd', '64ppd'       
    22  map_proj = 'lambert',           
    23  ref_lat   =  -35.         
    24  ref_lon   =  50.,           
    25  truelat1  =  -35.,           
    26  truelat2  =  0.0,               
    27  stand_lon =  50.,           
    28  geog_data_path = './WPS_GEOG',   
     22 map_proj = 'lambert'
     23 ref_lat   =  -35.         
     24 ref_lon   =  50.           
     25 truelat1  =  -35.           
     26 truelat2  =  0.0               
     27 stand_lon =  50.           
     28 geog_data_path = './WPS_GEOG'
    2929/
    3030
  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/runmeso

    r261 r262  
    3030###########################################
    3131
    32 changeregis=0
    33 while getopts "r" options; do
     32fresh=0
     33while getopts "f" options; do
    3434  case $options in
    35    r ) changeregis=1;;       
     35   r ) fresh=1;;       
    3636  esac
    3737done   
     
    412412   #
    413413   cd ${meso_folder}
    414                       if [[ "${changeregis}" = "1" ]]
    415                       then
    416                           cd ${work_folder}/WRFV2/Registry
    417                           \rm Registry 2> /dev/null
    418                           ./Registry.bash
    419                           cd ../../..
    420                       fi
    421 
     414                      #if [[ "${changeregis}" = "1" ]]
     415                      #then
     416                      #    cd ${work_folder}/WRFV2/Registry
     417                      #    \rm Registry 2> /dev/null
     418                      #    ./Registry.bash
     419                      #    cd ../../..
     420                      #fi
     421   if [[ ${fresh} -eq 0 ]]
     422   then
    422423     ##******
    423424     if [[ ${newphys} -eq 0 ]]
     
    433434     fi
    434435     ##******
    435 
     436   else
     437     ##******
     438     if [[ ${newphys} -eq 0 ]]
     439     then
     440       ./makemeso -r -f < ${simu_folder}/makemeso_input
     441     else
     442         if [[ "${scenario}" == "" ]]
     443         then
     444          ./makemeso -r -p -f < ${simu_folder}/makemeso_input
     445         else
     446          ./makemeso -r -p -f -s ${scenario} < ${simu_folder}/makemeso_input
     447         fi
     448     fi
     449     ##******
     450   fi
    436451
    437452   if [[ -f ${wrfexe} ]]
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM

    r170 r262  
    10871087rconfig   integer parent_time_step_ratio  namelist,domains      max_domains    1       h     "parent_time_step_ratio"                ""      ""
    10881088rconfig   integer feedback                namelist,domains      1    0       h     "feedback"          ""      ""
    1089 rconfig   integer smooth_option           namelist,domains      1    0       h     "smooth_option"          ""      ""
     1089rconfig   integer smooth_option           namelist,domains      1    2       h     "smooth_option"          ""      ""
    10901090rconfig   integer blend_width             namelist,domains      1    5       h     "blend_width"  "width of cg fg terrain blended zone"      ""
    10911091rconfig   real    ztop                    namelist,domains      max_domains    35000.  h    "ztop"          ""      ""
  • trunk/MESOSCALE/LMD_MM_MARS/makemeso

    r247 r262  
    259259      clean -a > /dev/null 2> /dev/null
    260260      cd Registry ; mv Registry Registry.bak ; Registry.bash ; cd ..
    261       answer='y' ## a voir... pas si sur
    262     else
    263             if [[ "${phys}" != "nophys_" ]]
    264             then
    265       echo Did you modify anything in the Registry or clean ? y for yes, any key for no ; read answer
    266       case ${answer} in
    267         y) cd Registry ; mv Registry Registry.bak ; Registry.bash ; cd .. ;;
    268         *) answer='no' ;;
    269       esac
    270             fi
     261#     answer='y' ## a voir... pas si sur
     262#    else
     263#            if [[ "${phys}" != "nophys_" ]]
     264#            then
     265#      echo Did you modify anything in the Registry or clean ? y for yes, any key for no ; read answer
     266#      case ${answer} in
     267#        y) cd Registry ; mv Registry Registry.bak ; Registry.bash ; cd .. ;;
     268#        *) answer='no' ;;
     269#      esac
     270#            fi
    271271    fi
    272272  fi
     
    311311  echo ${tra} >> last
    312312  echo ${scat} >> last
    313   echo ${answer} >> last
     313  #echo ${answer} >> last
    314314
    315315#------------
Note: See TracChangeset for help on using the changeset viewer.