Changeset 242 for trunk/MESOSCALE/LMD_MM_MARS/SIMU/runmeso

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.

File:

• trunk/MESOSCALE/LMD_MM_MARS/SIMU/runmeso (modified) (5 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SIMU/runmeso

@@ -48 +48 @@
 echo '**** LMD Martian Mesoscale Model launcher ****'
 echo '**********************************************'
-#echo '** Hello '$(whoami)' on '$(hostname)' !'
-#echo '** ['$(uname -a)']'
 echo '** User                  >>> '$(whoami)
 echo '** Hostname              >>> '$(hostname)
@@ -146 +144 @@
 dx="$(more ${simu_folder}/namelist.input | grep -e '[^!]*dx *=' | awk '{print $3}' | sed s+','+''+g)"
 dy="$(more ${simu_folder}/namelist.input | grep -e '[^!]*dy *=' | awk '{print $3}' | sed s+','+''+g)"
-
 
 #
@@ -201 +198 @@
 echo '** Horizontal resolution >>> '${dx}
 
-#echo "OK ? ('no' to exit, any key if OK)" ; read wellok
-#  if [[ "${wellok}" = "no" ]] 
-#  then 
-#     exit
-#  fi
-
 ###########
 # ACTIONS #
@@ -324 +315 @@
              echo 'Please define the domain (not for nests for the moment) : '
              echo " - 1/4 Projection: polar lambert mercator ? previous is "${proj} ; read proj
-             echo " - 2/4 Data source (multiple if nests): 'gcm' '32ppd' '64ppd' '64ppd_noHRti' '32ppd_HRalb' ? previous is "${data_source} ; read data_source
+             echo " - 2/4 Data source (multiple if nests): gcm 32ppd 64ppd 64ppd_noHRti 32ppd_HRalb ? previous is "${data_source} ; read data_source
              echo " - 3/4 North Latitude of domain centre ? previous is "${lat} ; read lat
              echo " - 4/4 East  Longitude of domain centre ? previous is "${lon} ; read lon
@@ -330 +321 @@
              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
              ####
+             if [[ ${nest} -gt 1 ]]
+             then
+                 echo '**********************************************'
+                 echo "USEFUL RULES for NESTS :"
+                 echo "* parent_grid_ratio divides (e_we-s_we+1)"
+                 echo "* nproc divides e_we-1"
+                 echo "* nproc divides e_we-1"
+             fi
+             if [[ ${nest} -ge 4 ]]
+             then 
+                 echo more than 3 nests, not supported
+                 echo please edit the namelist.wps template in WPS to allow this
+                 exit
+             fi
+             ####
              ln -sf ${simu_folder}/namelist.wps ${work_folder}/WPS/namelist.wps
-             #vi       ${work_folder}/WPS/namelist.wps
              echo '**********************************************'
              echo OK, thank you. Please wait... ncview will pop up soon !
              cd ${work_folder}/WPS/
              \rm geo_em.d0*.nc 2> /dev/null
-             ./geogrid.exe > /dev/null 2> /dev/null 
+             ./geogrid.exe #> /dev/null 2> /dev/null 
              ncview geo_em.d01.nc > /dev/null 2> /dev/null &
                 if [[ -f "geo_em.d02.nc" ]]