Changeset 653

Timestamp:

May 10, 2012, 12:49:00 AM (13 years ago)

Author:

aslmd

Message:

UTIL PYTHON : inimeso.py : a routine to generate input profiles for idealized mesoscale modeling.

Location:

trunk

Files:

• MESOSCALE_DEV/MANUAL/SRC/advance.tex (modified) (1 diff)
• UTIL/PYTHON/mcd/inimeso.py (added)
• UTIL/PYTHON/mcd/input_coord (added)
• UTIL/PYTHON/mcd/mcd.py (modified) (2 diffs)

trunk/MESOSCALE_DEV/MANUAL/SRC/advance.tex

@@ -104 +104 @@
 \paragraph{Inputs/outputs} Large-Eddy Simulations need four input files \ttt{input\_coord}, \ttt{input\_sounding}, \ttt{input\_more}, \ttt{input\_therm} which define initial pressure, temperature, density, winds profiles at the location/season for which simulations are run, along with information about this location/season. Typical files are available upon request, or you might simply build your own profiles using the Mars Climate Database (see the sample \ttt{scilab} script \ttt{wrf\_sounding.sci} in \ttt{\$MMM/SIMU/RUN}). Examples for \ttt{input\_*} files are provided in \ttt{\$MMM/SRC/LES/modif\_mars/DEF} and correspond to the cases run in the study by \textit{Spiga et al.} [2010].
 
+%% now python inimeso.py in UTIL/PYTHON
+
 \begin{citemize}
 \item \ttt{input\_coord} contains longitude, latitude, $L_s$ and local time;

trunk/UTIL/PYTHON/mcd/mcd.py

@@ -238 +238 @@
       for i in range(nd): self.lat = self.xcoord[i] ; self.update() ; self.put1d(i)
 
-    def profile(self,nd=20,start=0.,end=100000.):
+    def profile(self,nd=20,start=0.,end=100000.,tabperso=None):
     ### retrieve an altitude slice (profile)
       self.xlabel = "Altitude (m)"
-      self.prepare(ndx=nd) ; self.xcoord = np.linspace(start,end,nd)
+      if tabperso is not None: nd = len(tabperso)
+      self.prepare(ndx=nd) 
+      if tabperso is None:  self.xcoord = np.linspace(start,end,nd)
+      else:                 self.xcoord = tabperso
       for i in range(nd): self.xz = self.xcoord[i] ; self.update() ; self.put1d(i)
 
@@ -293 +296 @@
 
     ### TODO: makeplot2d, plot2d, passer plot settings, vecteurs, plot loct pas fixe
+