Ignore:
Timestamp:
Jul 25, 2011, 12:24:59 PM (13 years ago)
Author:
aslmd
Message:

MESOSCALE: corrected runmeso and namelist.wps template so that nesting is supported by runmeso. and lowered smoothing coefficient in GEOGRID.TBL from 5 to 3 for topography. raised it from 5 to 10 for thermal inertia because of stripes. lastly: discovered that two-way nesting is now working! to be further investigated.

Location:
trunk/MESOSCALE/LMD_MM_MARS
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/namelist.wps

    r11 r242  
    11&share                                 
    22 wrf_core = 'ARW',                    !!   [do not modify: choice of dynamical core]
    3  __MAX_DOM__
     3__MAX_DOM__
    44 start_date = '0000-00-00_00:00:00','0000-00-00_00:00:00','0000-00-00_00:00:00'   !! YYYY-MM-DD_HH:mm:ss start date
    55 end_date   = '1111-11-11_11:11:11','0000-00-00_00:00:00','0000-00-00_00:00:00'   !! YYYY-MM-DD_HH:mm:ss end date
     
    1212
    1313&geogrid
    14  __PARENT_ID__
    15  __PARENT_GRID_RATIO__
    16  __I_PARENT_START__
    17  __J_PARENT_START__
    18  __E_WE__
    19  __E_SN__
    20  dx = __DX__
    21  dy = __DY__
    22  geog_data_res     = __DATA_SOURCE__,       !! choice of static data sources
    23                                             !! NB: possible: '64ppd', '32ppd', ...
    24                                             !! NB: please glance at geogrid/GEOGRID.TBL
     14__PARENT_ID__
     15__PARENT_GRID_RATIO__
     16__I_PARENT_START__
     17__J_PARENT_START__
     18__E_WE__
     19__E_SN__
     20 dx = __DX__,
     21 dy = __DY__,
     22 geog_data_res     = '__DATA_SOURCE__', '__DATA_SOURCE__', '__DATA_SOURCE__'       !! choice of static data sources
     23                                                                                   !! NB: possible: '64ppd', '32ppd', ...
     24                                                                                   !! NB: please glance at geogrid/GEOGRID.TBL
    2525 map_proj = '__PROJ__',           !! map projection: 'mercator', 'lambert' or 'polar'
    2626 ref_lat   =  __LAT__,            !! north latitude of the center of the domain
  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/runmeso

    r236 r242  
    4848echo '**** LMD Martian Mesoscale Model launcher ****'
    4949echo '**********************************************'
    50 #echo '** Hello '$(whoami)' on '$(hostname)' !'
    51 #echo '** ['$(uname -a)']'
    5250echo '** User                  >>> '$(whoami)
    5351echo '** Hostname              >>> '$(hostname)
     
    146144dx="$(more ${simu_folder}/namelist.input | grep -e '[^!]*dx *=' | awk '{print $3}' | sed s+','+''+g)"
    147145dy="$(more ${simu_folder}/namelist.input | grep -e '[^!]*dy *=' | awk '{print $3}' | sed s+','+''+g)"
    148 
    149146
    150147#
     
    201198echo '** Horizontal resolution >>> '${dx}
    202199
    203 #echo "OK ? ('no' to exit, any key if OK)" ; read wellok
    204 #  if [[ "${wellok}" = "no" ]]
    205 #  then
    206 #     exit
    207 #  fi
    208 
    209200###########
    210201# ACTIONS #
     
    324315             echo 'Please define the domain (not for nests for the moment) : '
    325316             echo " - 1/4 Projection: polar lambert mercator ? previous is "${proj} ; read proj
    326              echo " - 2/4 Data source (multiple if nests): 'gcm' '32ppd' '64ppd' '64ppd_noHRti' '32ppd_HRalb' ? previous is "${data_source} ; read data_source
     317             echo " - 2/4 Data source (multiple if nests): gcm 32ppd 64ppd 64ppd_noHRti 32ppd_HRalb ? previous is "${data_source} ; read data_source
    327318             echo " - 3/4 North Latitude of domain centre ? previous is "${lat} ; read lat
    328319             echo " - 4/4 East  Longitude of domain centre ? previous is "${lon} ; read lon
     
    330321             sed s+"__INTERVAL_SECONDS__"+"${interval_seconds}"+g ${meso_folder}/SIMU/namelist.wps | sed s+"__MAX_DOM__"+"${max_dom}"+g | sed s+"__PARENT_ID__"+"${parent_id}"+g |  sed s+"__PARENT_GRID_RATIO__"+"${parent_grid_ratio}"+g |  sed s+"__I_PARENT_START__"+"${i_parent_start}"+g | sed s+"__J_PARENT_START__"+"${j_parent_start}"+g | sed s+"__E_WE__"+"${e_we}"+g | sed s+"__E_SN__"+"${e_sn}"+g | sed s+"__DX__"+"${dx}"+g | sed s+"__DY__"+"${dy}"+g | sed s+"__PROJ__"+"${proj}"+g | sed s+"__DATA_SOURCE__"+"${data_source}"+g | sed s+"__LAT__"+"${lat}"+g | sed s+"__LON__"+"${lon}"+g | sed s+'0000-00-00_00:00:00'+${year}'-'${month}'-'${day}'_'${hour}':00:00'+g | sed s+'1111-11-11_11:11:11'+${eyear}'-'${emonth}'-'${eday}'_'${ehour}':00:00'+g > ${simu_folder}/namelist.wps
    331322             ####
     323             if [[ ${nest} -gt 1 ]]
     324             then
     325                 echo '**********************************************'
     326                 echo "USEFUL RULES for NESTS :"
     327                 echo "* parent_grid_ratio divides (e_we-s_we+1)"
     328                 echo "* nproc divides e_we-1"
     329                 echo "* nproc divides e_we-1"
     330             fi
     331             if [[ ${nest} -ge 4 ]]
     332             then
     333                 echo more than 3 nests, not supported
     334                 echo please edit the namelist.wps template in WPS to allow this
     335                 exit
     336             fi
     337             ####
    332338             ln -sf ${simu_folder}/namelist.wps ${work_folder}/WPS/namelist.wps
    333              #vi       ${work_folder}/WPS/namelist.wps
    334339             echo '**********************************************'
    335340             echo OK, thank you. Please wait... ncview will pop up soon !
    336341             cd ${work_folder}/WPS/
    337342             \rm geo_em.d0*.nc 2> /dev/null
    338              ./geogrid.exe > /dev/null 2> /dev/null
     343             ./geogrid.exe #> /dev/null 2> /dev/null
    339344             ncview geo_em.d01.nc > /dev/null 2> /dev/null &
    340345                if [[ -f "geo_em.d02.nc" ]]
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WPS/wps_mars/geogrid/GEOGRID.TBL.ARW_Mars

    r69 r242  
    88        df_dy=SLPY
    99        fill_missing=9000.
    10         smooth_option = 1-2-1; smooth_passes=5
     10        smooth_option = 1-2-1; smooth_passes=3 #5
    1111#        smooth_option = smth-desmth_special; smooth_passes=1
    1212        interp_option =         default:four_pt
     
    3333        dest_type = continuous
    3434        fill_missing=9000.
    35         smooth_option = 1-2-1; smooth_passes=5
     35        smooth_option = 1-2-1; smooth_passes=3 #5
    3636        interp_option =         default:four_pt
    3737        interp_option =         32ppd:four_pt
     
    5757        dest_type = continuous
    5858        fill_missing=9000.
    59         smooth_option = 1-2-1; smooth_passes=5
     59        smooth_option = 1-2-1; smooth_passes=3 #5
    6060        interp_option =         default:four_pt
    6161        interp_option =         32ppd:four_pt
     
    7979        priority = 1
    8080        dest_type = continuous
    81         smooth_option = 1-2-1; smooth_passes=5
     81        smooth_option = 1-2-1; smooth_passes=10 #5
    8282        interp_option =         32ppd:sixteen_pt+four_pt+average_4pt+average_16pt+search
    8383        rel_path=               32ppd:thermal_GCM/
Note: See TracChangeset for help on using the changeset viewer.