Changeset 723 for trunk/UTIL

Timestamp:

Jul 16, 2012, 5:49:10 PM (12 years ago)

Author:

aslmd

Message:

UTIL PYTHON : added seasonal and wind capability to MCD python interface. updated READMEs for LMD insiders.

Location:

trunk/UTIL/PYTHON

Files:

• README.INSTALL (modified) (2 diffs)
• README.PP (modified) (1 diff)
• mcd/mcd.py (modified) (2 diffs)
• mcd/proto/index.html (modified) (1 diff)

trunk/UTIL/PYTHON/README.INSTALL

@@ +1 @@
+###########################
 ############ README.INSTALL
+###########################
+
 This README describes the required steps to install the python distribution and packages necessary to run the graphical routines.
 
@@ -42 +45 @@
 Step 4/- Read the README.PP ;)
 ------------------------------
+
+NB: To use the movie capability, you will have to install mencoder.
+NB: To use advance interpolation methods for GCM and mesoscale files, you will have to compile fortran routines with f2py (or compile those independently).
+
+###########################
+############ README.INSTALL.LMD
+###########################
+
+If you work on LMD or IPSL computers, you don't necessarily need to install EPD.
+
+[NB: You don't necessarily need to download pp.py sources either but in case
+you would like to try modifications of the code please use the command
+   svn co http://svn.lmd.jussieu.fr/Planeto/trunk/UTIL/PYTHON
+and change PYTHONPATH accordingly in the lines given in what follows]
+
+-------------------------
+account on lmd.jussieu.fr
+-------------------------
+
+A. Spiga installed on his account a version of EPD.
+--> Add the following lines to your env_pers.ksh, then source the file. You're OK! 
+--> (ask for an installation of mencoder if you'd like to use movie capabilities)
+
+#####################################
+if [ `uname` = 'Linux' ] ;then
+  if [ `uname -m` = 'x86_64' ] ;then
+  PATH="/donnees/aslmd/PY/epd-7.0-2-rh3-x86_64/bin/:$PATH"
+  else
+  PATH="/donnees/aslmd/PY/epd-7.0-2-rh5-x86/bin/:$PATH"
+  fi
+fi
+PYTHONPATH=/donnees/aslmd/MODELES/UTIL/PYTHON/
+export PYTHONPATH
+PATH=$PYTHONPATH:$PATH
+#####################################
+
+---------------------------------
+account on ciclad.ipsl.jussieu.fr
+---------------------------------
+
+P. Weill installed on the system a version of EPD (and mencoder for movie capabilities)
+--> Just add the following lines to your .bashrc, then source the file. You're OK!
+
+######################################
+PATH="/opt/epd-7.1-1/bin/:$PATH"
+PYTHONPATH=/home/spiga/MODELES/UTIL/PYTHON/
+export PYTHONPATH
+PATH=$PYTHONPATH:$PATH
+######################################
+
+-----------------------------------------------------------
+-----------------------------------------------------------
+-----------------------------------------------------------
+-----------------------------------------------------------
+Authors : A. Colaitis and A. Spiga
+

trunk/UTIL/PYTHON/README.PP

@@ -11 +11 @@
 **************************************
 **************************************
+
+The program pp.py works on virtually any NETCDF file. 
+If not, please contact us, it is usually an easy fix to make it work.
 
 ************************************

trunk/UTIL/PYTHON/mcd/mcd.py

@@ -247 +247 @@
       for i in range(nd): self.xz = self.xcoord[i] ; self.update() ; self.put1d(i)
 
+    def seasonal(self,nd=12,start=0.,end=360.):
+    ### retrieve a seasonal slice
+      self.xlabel = "Areocentric longitude (degrees)"
+      self.prepare(ndx=nd) ; self.xcoord = np.linspace(start,end,nd)
+      for i in range(nd): self.xdate = self.xcoord[i] ; self.update() ; self.put1d(i)
+
     def latlon(self,ndx=37,startx=-180.,endx=180.,ndy=19,starty=-90.,endy=90.):
     ### retrieve a latitude/longitude slice
@@ -283 +289 @@
       self.extvartab[i,j,1:100] = self.extvar[0:99] ## note: var numbering according to MCD manual is kept
 
-    def makemap2d(self,choice):
+    def makemap2d(self,choice,incwind=False):
     ### one 2D map is created for the user-defined variable in choice.
+      self.latlon() ## a map is implicitely a lat-lon plot. otherwise it is a plot (cf. makeplot2d)
       (field, fieldlab) = self.definefield(choice)
-      myplot.maplatlon(self.xcoord,self.ycoord,field,title=fieldlab,proj="moll")
-
-    def map2d(self,tabtodo):
+      if incwind:
+          (windx, fieldlabwx) = self.definefield("u")
+          (windy, fieldlabwy) = self.definefield("v")
+          myplot.maplatlon(self.xcoord,self.ycoord,field,title=fieldlab,proj="cyl",vecx=windx,vecy=windy)
+      else:
+          myplot.maplatlon(self.xcoord,self.ycoord,field,title=fieldlab,proj="moll")
+
+    def map2d(self,tabtodo,incwind=False):
     ### complete 2D figure with possible multiplots
       if isinstance(tabtodo,np.str): tabtodo=[tabtodo] ## so that asking one element without [] is possible.
       if isinstance(tabtodo,np.int): tabtodo=[tabtodo] ## so that asking one element without [] is possible.
       fig = mpl.figure() ; subv,subh = myplot.definesubplot( len(tabtodo) , fig ) 
-      for i in range(len(tabtodo)): mpl.subplot(subv,subh,i+1) ; self.makemap2d(tabtodo[i])
+      for i in range(len(tabtodo)): mpl.subplot(subv,subh,i+1) ; self.makemap2d(tabtodo[i],incwind=incwind)
 
     ### TODO: makeplot2d, plot2d, passer plot settings, vecteurs, plot loct pas fixe

trunk/UTIL/PYTHON/mcd/proto/index.html

@@ -51 +51 @@
 <li>Altitude   
 <input type="radio" name="isaltfree" value="0" checked> Value: <input type="text" size="2" name="altitude" value="10.">
-<input type="radio" name="isaltfree" value="1"> All [0 ; 100000] km
+<input type="radio" name="isaltfree" value="1"> All [0 ; 100000] m
 </li>
 </ul>