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

Timestamp:

Jun 8, 2011, 6:01:00 PM (14 years ago)

Author:

aslmd

Message:

MESOSCALE: modified makemeso and runmeso to comply with changes performed in r148 about scatterers in newphys. a new question is now asked if you run makemeso -p or if you run runmeso with run.def in the folder. note that runmeso only asks the question for scenarios mars=1 or mars=11 where water cycle is on, otherwise by default number of scatterers is 1 (i.e. only dust).

File:

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

trunk/MESOSCALE/LMD_MM_MARS/SIMU/runmeso

@@ -27 +27 @@
 
 scenario=""
-#scenario="storm"
+scenario="storm"
 
 if [[ "${LMDMOD}" = "" ]]
@@ -167 +167 @@
         exit
   fi
-tracers=$(more ${simu_folder}/namelist.input | grep mars | awk '{print $3}' | sed s+','+''+g)
-  case ${tracers} in
+mars_scenario=$(more ${simu_folder}/namelist.input | grep mars | awk '{print $3}' | sed s+','+''+g)
+  case ${mars_scenario} in
       0) tracers='1' ;;
       1) tracers='2' ;;
@@ -178 +178 @@
 echo '** Domain(s)             >>> '${nest} 
 echo '** Tracer(s)             >>> '${tracers}
+
+   ############ SCATTERERS SCATTERERS ######
+   if [[ ${newphys} -eq 1 ]]
+   then
+     #### CHANGE THIS IF YOU ADD NEW SCENARIOS
+     case ${mars_scenario} in
+       0) scat=1 ;;
+       1) echo "!! How many scatterers are you using : 1,2,..." ; read scat ;;
+       2) scat=1 ;;
+       3) scat=1 ;;
+      11) echo "!! How many scatterers are you using : 1,2,..." ; read scat ;;
+      20) scat=1 ;;
+       *) scat=1 ;;
+     esac
+     echo '** Scatterers            >>> '${scat}
+   else 
+     #### USELESS here BUT COMPATIBILITY
+     scat=1
+   fi
+   #########################################
+
 echo '** Grid points           >>> '${nx}" x "${ny}" x "${nz}
 echo '** Horizontal resolution >>> '${dx}
+
 #echo "OK ? ('no' to exit, any key if OK)" ; read wellok
 #  if [[ "${wellok}" = "no" ]] 
@@ -326 +348 @@
  cd ${simu_folder}
      ## corresponding executables
-    wrfexe=${work_folder}/wrf_x${nx}_y${ny}_z${nz}_d${nest}_t${tracers}_p${numproc}.exe
-    realexe=${work_folder}/real_x${nx}_y${ny}_z${nz}_d${nest}_t${tracers}_p${numproc}.exe
+    wrfexe=${work_folder}/wrf_x${nx}_y${ny}_z${nz}_d${nest}_t${tracers}_p${numproc}_s${scat}.exe
+    realexe=${work_folder}/real_x${nx}_y${ny}_z${nz}_d${nest}_t${tracers}_p${numproc}_s${scat}.exe
  if [[ -f ${wrfexe} ]]
  then
@@ -360 +382 @@
    echo ${nest} >> ${simu_folder}/makemeso_input
    echo ${tracers} >> ${simu_folder}/makemeso_input
-   #echo 'y' >> ${simu_folder}/makemeso_input
+   if [[ ${newphys} -eq 1 ]]
+   then
+     echo 'n' >> ${simu_folder}/makemeso_input
+     echo ${scat} >> ${simu_folder}/makemeso_input
+   fi
    #
    # COMPILE THE MODEL