Changeset 1212 for trunk/MESOSCALE/LMD_MM_MARS/makemeso

Timestamp:

Mar 24, 2014, 9:46:50 AM (11 years ago)

Author:

aslmd

Message:

LMDZ.MARS + MESOSCALE

A quite major commit, at least for MESOSCALE.
In a word: any ngrid deserves to be free.

• no need to recompile when changing number of horizontal grid points or number of processors
• latest version of LMDZ.MARS physics can be used
• WARNING! Nesting is still yet to be fixed (since r1027)

Also some small bug fixes to LMDZ.MARS.

Changes in LMDZ.MARS

---

--> fixed a potential bug in thermal plume model because zlmax was computed both in thermcell_main_mars and calltherm_interface... so made it an OUT argument of calltherm_interface. also: changed the name to limz. and added precompiling flags to avoid the use of planetwide in MESOSCALE. in MESOSCALE we just go high enough (nlayer-5) and do not care about computational cost (although we certainly gain from not using MAXVAL).
--> moved allocations upward in inifis. does not change anything for GCM, but make MESOSCALE modifications simpler, and overall make inifis better organized: first allocations, then reading callphys.def file.
--> added precompiling flags around lines that are both useless for MESOSCALE (notably I/O) and recently adapted to parallel computations in the GCM
--> tidied up what is MESOSCALE vs. GCM in surfini

Changes in MESOSCALE

---

--> changed makemeso to allow dynamically set nx ny nprocs
--> changed makemeso to remove links to Fortran code adapted to parallel GCM and useless for mesoscale
--> changed ngridmx to ngrid in inifis includes

File:

• trunk/MESOSCALE/LMD_MM_MARS/makemeso (modified) (6 diffs)

trunk/MESOSCALE/LMD_MM_MARS/makemeso

@@ -134 +134 @@
   echo Grid points in longitude ? ; read lon 
   echo Grid points in latitude ? ; read lat 
+  # ngrid: now dynamically set in newphys 03/2014
+   if [[ "${phys}" == "newphys_" ]]
+   then
+      lon=999
+      lat=999
+      numproc=999
+   fi
   echo Number of vertical levels ? ; read level 
   fi
@@ -174 +181 @@
          else
             echo Number of tracers ? ; read tra
-            # tracers: no dynamically set in newphys 09/2013
+            # tracers: now dynamically set in newphys 09/2013
             if [[ "${phys}" == "newphys_" ]]
             then
@@ -243 +250 @@
      #ln -sf meso_callkeys_newphys.h meso_callkeys.h
      cd Registry ; Registry.bash ; cd ..
+
+     # remove files adapted to parallel GCM and useless for mesoscale
+     cd mars_lmd
+     cd libf/phymars
+     \rm tabfi* soil_settings* phyetat0* phyredem* iostart* writediagfi* mkstat* writediagsoil* eofdump* wstats* inistats* mod_*
+     mv ../dyn3d/control_mod.F90 .
+     cd ../..
+     cd ..
+
     else
          if [[ "${config}" == "les" ]] 
@@ -594 +610 @@
   cd mars_lmd/
 
+ if [[ "${phys}" == "newphys_" ]]
+ then
+  divx=1
+  divy=1
+ else
   # required size
   #----------------
@@ -611 +632 @@
       *) echo not yet supported ; exit ;;
   esac
+ fi
   physx=$(expr ${lon} - 1)
   restex=$(expr ${physx} \% ${divx})   
@@ -629 +651 @@
   if [[ "${phys}" == "newphys_" ]]
   then
-     soilsize=18 ## nouvelle physique
+     soilsize=18 ## nouvelle physique. ne sert a rien ici. voir comsoil.
   else
      soilsize=10