Changeset 235 for trunk/MESOSCALE/LMD_MM_MARS

Timestamp:

Jul 20, 2011, 3:07:02 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: corrected runmeso and gcm initialization so that it could be used on GNOME ifort cluster and virtually everywhere. changed choice numbering and added choices to initialize but not execute (good for job submission). Also corrected what was broken in the two previous commits on the version of the mesoscale model with old physics

Location:

trunk/MESOSCALE/LMD_MM_MARS

Files:

• SIMU/runmeso (modified) (6 diffs)
• SRC/PREP_MARS/compile_and_exec (modified) (1 diff)
• SRC/PREP_MARS/compile_ifort (modified) (1 diff)
• SRC/PREP_MARS/compile_pgf (moved) (moved from trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/compile)
• SRC/WRFV2/mars_lmd/libf/phymars/inifis.F (moved) (moved from trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/meso_inifis.F) (2 diffs)
• SRC/WRFV2/mars_lmd/libf/phymars/meso_testphys1d.F (deleted)
• SRC/WRFV2/mars_lmd/libf/phymars/physiq.F (moved) (moved from trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/meso_physiq.F) (4 diffs)
• SRC/WRFV2/mars_lmd/libf/phymars/testphys1d.F (added)
• makemeso (modified) (2 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SIMU/runmeso

@@ -210 +210 @@
 # ACTIONS #
 ###########
-#echo '***********************'
-#echo '****     Steps     ****'
-#echo '***********************'
 echo '**********************************************'
-echo "** This script will compile the model then :"
-echo "**      <1> Prepare GCM inputs          <-1> to skip compilation" 
-echo "**      <2> Perform WPS interpolations  <-2> to skip compilation" 
-echo "**      <3> Execute real.exe            <-3> to skip compilation" 
-echo "**      <4> Execute wrf.exe             <-4> to skip compilation"
-echo "**      <5> Exit"
-echo "**      <9> Oops ! Exit now please."
+echo "** The runmeso script does it for you !"
 echo '**********************************************'
-#echo "Start at which step ? (1-4)" ; read step
-echo "Where shall we start ? 1, 2, 3, 4, 5 or 9" ; read step
+echo "** Step 0: compile the model if not done yet"
+echo "** Step 1: run GCM and prepare inputs"
+echo "** Step 2: interpolate on domain horizontal grid"
+echo "** Step 3: interpolate on chosen vertical levels"
+echo "** Step 4: run the model"
+echo "************************* Please make a choice"
+echo " 1 --> Perform steps 0 1 2 3 4  [ -1: step 0 skipped]" 
+echo " 2 --> Perform steps 0 2 3 4    [ -2: step 0 skipped]"
+echo " 3 --> Perform steps 0 3 4      [ -3: step 0 skipped]" 
+echo " 4 --> Perform steps 0 4        [ -4: step 0 skipped]"
+echo "11 --> Perform steps 0 1 2 3    [-11: step 0 skipped]" 
+echo "12 --> Perform steps 0 2 3      [-12: step 0 skipped]" 
+echo "13 --> Perform steps 0 3        [-13: step 0 skipped]"
+echo " 0 --> Perform step  0 only"
+echo "99 --> Oops ! Exit now please."
+echo '**********************************************'
+echo "What is your choice? " ; read step
 if [[ ${step} -lt 0 ]]
 then
@@ -231 +237 @@
 fi
 case ${step} in
-        9) exit ;;
-        5) gcmstep=0 ; wpsstep=0 ; realstep=0 ; wrfstep=0 ;;
-        1) gcmstep=1 ; wpsstep=1 ; realstep=1 ; wrfstep=1 ;;
-        2) gcmstep=0 ; wpsstep=1 ; realstep=1 ; wrfstep=1 ;;
-        3) gcmstep=0 ; wpsstep=0 ; realstep=1 ; wrfstep=1 ;;
-        4) gcmstep=0 ; wpsstep=0 ; realstep=0 ; wrfstep=1 ;;
-        *) echo "** I guess you do not know what to do. Let's start from the beginning."  
-           gcmstep=1 ; wpsstep=1 ; realstep=1 ; wrfstep=1 ;;
+        99) exit ;;
+        0)  gcmstep=0 ; wpsstep=0 ; realstep=0 ; wrfstep=0 ;;
+        1)  gcmstep=1 ; wpsstep=1 ; realstep=1 ; wrfstep=1 ;;
+        2)  gcmstep=0 ; wpsstep=1 ; realstep=1 ; wrfstep=1 ;;
+        3)  gcmstep=0 ; wpsstep=0 ; realstep=1 ; wrfstep=1 ;;
+        4)  gcmstep=0 ; wpsstep=0 ; realstep=0 ; wrfstep=1 ;;
+        11) gcmstep=1 ; wpsstep=1 ; realstep=1 ; wrfstep=0 ;;
+        12) gcmstep=0 ; wpsstep=1 ; realstep=1 ; wrfstep=0 ;;
+        13) gcmstep=0 ; wpsstep=0 ; realstep=1 ; wrfstep=0 ;;
+        *)  echo "** I guess you do not know what to do. Let's start from the beginning."  
+            gcmstep=1 ; wpsstep=1 ; realstep=1 ; wrfstep=1 ;;
 esac
 
@@ -262 +271 @@
    if [[ ${compilo} -eq 3 || ${compilo} -eq 5 ]]
    then
-       echo How many processors ? 1, 2, 4, 6, 8, 12, 16, 20, 24, 32 ; read numproc
+       echo How many processors ? e.g. 1, 2, 4, 6, 8, 12, 16, 20, 24, 32, 64, 128 ; read numproc
    else
        numproc=1 
@@ -446 +455 @@
         if [[ ${newphys} -eq 0 ]]
         then
-          cd ${main_folder}'/LMDZ.MARS/myGCM' 
+          cd ${main_folder}'/LMDZ.MARS/' 
         else
-          cd ${main_folder}'/LMDZ.MARS.new/myGCM'
+          cd ${main_folder}'/LMDZ.MARS.new/'
         fi
-         #echo '(during the run, you might like to check domain settings in '${work_folder}'/WPS)'
+        if [[ ! ( -f "gcm.e" ) ]]
+        then
+          echo "Oops. There is no gcm.e, you shall compile the GCM first."
+          exit
+        fi
+        cd myGCM
         echo ${daygcm} | ./launch_gcm
         #
@@ -458 +472 @@
         if [[ -f ${work_folder}/PREP_MARS/create_readmeteo.exe && -f ${work_folder}/PREP_MARS/readmeteo.exe ]]
         then
-              cd ${work_folder}/PREP_MARS/  
-                ./compile_and_exec
+              ### This is now generic so PREP_MARS could be moved out of DIRCOMPs files
+              cd ${work_folder}/PREP_MARS/ 
+              if [[ ${compilo} -eq 1 || ${compilo} -eq 3 ]]
+              then
+                 ln -sf compile_pgf compile
+              fi
+              if [[ ${compilo} -eq 2 ]]
+              then
+                 ln -sf compile_g95 compile
+              fi
+              if [[ ${compilo} -eq 4 || ${compilo} -eq 5 ]]
+              then
+                 ln -sf compile_ifort compile
+              fi
+              if [[ ${compilo} -gt 5 ]]
+              then
+                 echo "please adapt compile in PREP_MARS"
+              fi
+              ./compile_and_exec
               #echo 1 | ./create_readmeteo.exe
               #./readmeteo.exe < readmeteo.def
@@ -521 +552 @@
         ln -sf ${simu_folder}/namelist.input .
         ln -sf ${simu_folder}/real.exe .
+        ln -sf ${simu_folder}/levels .
         rm -rf ${simu_folder}/met_em/wrfbdy* 2> /dev/null
         rm -rf ${simu_folder}/met_em/wrfinput* 2> /dev/null
-        if [[ ${compilo} -eq 3 ]]
+        if [[ ${compilo} -eq 3 || ${compilo} -eq 5 ]]
         then
                 rm -rf ${simu_folder}/met_em/rsl.* 2> /dev/null

trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/compile_and_exec

@@ -8 +8 @@
 echo 1 | create_readmeteo.exe 
 readmeteo.exe < readmeteo.def 
-

trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/compile_ifort

@@ -18 +18 @@
 \rm *.o 2> /dev/null
 
-echo I assume input_diagfi is correctly linked
-echo I assume WPSFEED folder exists or is linked
-
-echo 1 | create_readmeteo.exe
-readmeteo.exe < readmeteo.def
+#echo I assume input_diagfi is correctly linked
+#echo I assume WPSFEED folder exists or is linked
+#echo 1 | create_readmeteo.exe
+#readmeteo.exe < readmeteo.def

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/inifis.F

@@ -1 @@
-      SUBROUTINE meso_inifis(ngrid,nlayer,nq,
-     $           wdt,
+      SUBROUTINE inifis(ngrid,nlayer,
      $           wday_ini,wdaysec,
-     $           wappel_phys,wecri_phys,
+     $           wappel_phys,
      $           plat,plon,parea,
      $           prad,pg,pr,pcpp,
+     $           nq, wdt,
      $           womeg,wmugaz,
      $           wyear_day,wperiheli,waphelie,wperi_day,wobliquit, 
@@ -528 +528 @@
         dtphys=wdt*float(wappel_phys)
         print*,'Physical timestep (s) ',dtphys
-        ecri_phys=wecri_phys
+        ecri_phys=8.e18  !! a dummy low frequency
         print*,'Physical frequency for writing ',ecri_phys
 c

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/physiq.F

@@ -1 @@
-      SUBROUTINE meso_physiq(ngrid,nlayer,nq,
+      SUBROUTINE physiq(ngrid,nlayer,nq,
      $            firstcall,lastcall,
-     $            wday_ini,
      $            pday,ptime,ptimestep,
      $            pplev,pplay,pphi,
@@ -8 +7 @@
      $            pdu,pdv,pdt,pdq,pdpsrf,tracerdyn,
      $            wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,
-     $            wecritphys,
+     $            wday_ini,
      $            output_tab2d, output_tab3d,
      $            flag_LES)
-
 
 
@@ -174 +172 @@
       REAL wqsurf(ngridmx,nqmx)       
       REAL wq2(ngridmx,nlayermx+1) 
-      REAL wecritphys
       REAL output_tab2d(ngridmx,n2d)
       REAL output_tab3d(ngridmx,nlayer,n3d)
@@ -395 +392 @@
       iphysiq=ptimestep
 c
-      ecritphy=wecritphys
+      ecritphy=8.e18  !! a dummy low frequency
       PRINT*,'Write LMD physics each:',ecritphy,' seconds'
               !!PRINT*,ecri_phys  

trunk/MESOSCALE/LMD_MM_MARS/makemeso

@@ -234 +234 @@
      ./copy_model_tmp
      \rm copy_model_tmp
-     mv WRFV2 ${conf_wrf}/
+     mv zeWRFV2 ${conf_wrf}/WRFV2
      cd ${conf_wrf}/WRFV2
      cd Registry ; Registry.bash ; cd .. 
@@ -581 +581 @@
   cd libf/phymars
   \rm dimphys.h 2> /dev/null
+  if [[ "${phys}" == "newphys_" ]]
+  then
+    ln -sf meso_inc/meso_dimphys.h_ref .
+  fi
   physize=$(expr ${physx} \* ${physy})
-  sed s/--xsize--/${physx}/g meso_inc/meso_dimphys.h_ref | sed s/--ysize--/${physy}/g | sed s/--physize--/${physize}/g | sed s/--zsize--/${physz}/g | sed s/--soilsize--/${soilsize}/g > dimphys.h
+  sed s/--xsize--/${physx}/g meso_dimphys.h_ref | sed s/--ysize--/${physy}/g | sed s/--physize--/${physize}/g | sed s/--zsize--/${physz}/g | sed s/--soilsize--/${soilsize}/g > dimphys.h
   head -15 dimphys.h
      ### TEST new new phys