Changeset 793

Timestamp:

Sep 21, 2012, 5:11:37 PM (12 years ago)

Author:

aslmd

Message:

UTIL PYTHON : a working version of MCD python on an Apache server thanks to F. Guyon and L. Fairhead. plus added specific functions in mcd.py to make plots without using GUI, it was necessary to rewrite some stuff to avoid calling to matplotlib.pyplot.

Location:

trunk/UTIL/PYTHON/mcd

Files:

• gale.py (modified) (3 diffs)
• mcd.py (modified) (4 diffs)
• proto/cgi-bin/mcdcgi.py (modified) (7 diffs)
• proto/index.html (modified) (2 diffs)

trunk/UTIL/PYTHON/mcd/gale.py

@@ -1 +1 @@
 #! /usr/bin/env python
+
+#import sys
+#
+#class NullDevice():
+#    def write(self, s):
+#        pass
+#
+#original_stdout = sys.stdout  # keep a reference to STDOUT
+#
+#sys.stdout = NullDevice()  # redirect the real STDOUT
+#sys.stderr = NullDevice()  # redirect the real STDOUT
 
 from mcd import mcd
@@ -18 +29 @@
 #gale.printextvar("rho")
 
-#gale.seasonal()
-#gale.plot1d(["tsurf","u","v"])
+gale.seasonal()
+gale.plot1d(["tsurf","u","v"])
+
+import matplotlib.pyplot as mpl
+mpl.savefig("temp.png",dpi=85,bbox_inches='tight',pad_inches=0.25)
+
 
 #gale.xdate = 270.
@@ -25 +40 @@
 #gale.plot1d(["u","v"])
 
-#gale.latlon()
-#gale.map2d("tsurf")
+gale.latlon()
+gale.map2d("tsurf")
 
 import matplotlib.pyplot as mpl
-mpl.show()
+#mpl.show()
+mpl.savefig("img/temp.png",dpi=110,bbox_inches='tight',pad_inches=0.4)

trunk/UTIL/PYTHON/mcd/mcd.py

@@ -12 +12 @@
 import myplot
 
-class mcd:
-
+
+class mcd():
+ 
     def __repr__(self):
     # print out a help string when help is invoked on the object
@@ -27 +28 @@
         ## 0. general stuff
         self.name      = "MCD v4.3 output"
-        self.dset      = '/home/aymeric/Science/MCD_v4.3/data/'
+        #self.dset      = '/home/aymeric/Science/MCD_v4.3/data/'
+        self.dset      = '/home/marshttp/MCD_v4.3/data/'
         ## 1. spatio-temporal coordinates
         self.lat       = 0.
@@ -303 +305 @@
       for i in range(len(tabtodo)): mpl.subplot(subv,subh,i+1).grid(True, linestyle=':', color='grey') ; self.makeplot1d(tabtodo[i],vertplot)
 
+    def htmlplot1d(self,tabtodo,vertplot=0,figname="temp.png"):
+    ### complete 1D figure with possible multiplots
+    ### added in 09/2012 for online MCD
+    ### see http://www.dalkescientific.com/writings/diary/archive/2005/04/23/matplotlib_without_gui.html
+      from matplotlib.figure import Figure
+      from matplotlib.backends.backend_agg import FigureCanvasAgg
+      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 = Figure(figsize=(8,8)) ; subv,subh = myplot.definesubplot( len(tabtodo) , fig )
+      for i in range(len(tabtodo)):
+        yeah = fig.add_subplot(subv,subh,i+1) #.grid(True, linestyle=':', color='grey') 
+        choice = tabtodo[i]
+        (field, fieldlab) = self.definefield(choice)
+        if vertplot != 1:  absc = self.xcoord ; ordo = field ; ordolab = fieldlab ; absclab = self.xlabel
+        else:              ordo = self.xcoord ; absc = field ; absclab = fieldlab ; ordolab = self.xlabel
+        yeah.plot(absc,ordo,'-bo') #; mpl.xticks(query.xcoord)
+        ax = fig.gca() ; ax.set_ylabel(ordolab) ; ax.set_xlabel(absclab)
+      fig.text(0.5, 0.01, "Mars Climate Database (c) LMD/OU/IAA/ESA/CNES", ha='center')
+      canvas = FigureCanvasAgg(fig)
+      # The size * the dpi gives the final image size
+      #   a4"x4" image * 80 dpi ==> 320x320 pixel image
+      canvas.print_figure(figname, dpi=80)
+
 ###################
 ### 2D analysis ###
@@ -346 +371 @@
       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 ) 
+      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],incwind=incwind,fixedlt=fixedlt)
 
+    def htmlmap2d(self,tabtodo,incwind=False,fixedlt=False,figname="temp.png"):
+    ### complete 2D figure with possible multiplots
+    ### added in 09/2012 for online MCD
+    ### see http://www.dalkescientific.com/writings/diary/archive/2005/04/23/matplotlib_without_gui.html
+      from matplotlib.figure import Figure
+      from matplotlib.backends.backend_agg import FigureCanvasAgg
+      from matplotlib.cm import get_cmap
+      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 = Figure(figsize=(8,8)) ; subv,subh = myplot.definesubplot( len(tabtodo) , fig )
+
+      ### topocontours
+      fieldc = self.getextvar(self.convertlab("topo"))
+
+      for i in range(len(tabtodo)):
+        yeah = fig.add_subplot(subv,subh,i+1)
+        choice = tabtodo[i]
+        self.latlon(fixedlt=fixedlt) ## a map is implicitely a lat-lon plot. otherwise it is a plot (cf. makeplot2d)
+        (field, fieldlab) = self.definefield(choice)
+        if incwind: (windx, fieldlabwx) = self.definefield("u") ; (windy, fieldlabwy) = self.definefield("v")
+
+        proj="cyl" ; colorb="jet" ; ndiv=20 ; zeback="molabw" ; trans=1.0 #0.6
+        title="" ; vecx=None ; vecy=None ; stride=2
+        lon = self.xcoord
+        lat = self.ycoord
+
+        ### get lon and lat in 2D version. get lat/lon intervals
+        #numdim = len(np.array(lon).shape)
+        #if numdim == 2:     [lon2d,lat2d] = [lon,lat]
+        #elif numdim == 1:   [lon2d,lat2d] = np.meshgrid(lon,lat)
+        #else:               errormess("lon and lat arrays must be 1D or 2D")
+        #[wlon,wlat] = myplot.latinterv()
+        ### define projection and background. define x and y given the projection
+        #m = basemap.Basemap(projection='moll') marche pas
+        #m = myplot.define_proj(proj,wlon,wlat,back=zeback,blat=None,blon=None)
+        #x, y = m(lon2d, lat2d)
+        ### TEMP
+        x = lon ; y = lat
+        ## define field. bound field.
+        what_I_plot = np.transpose(field)
+        zevmin, zevmax = myplot.calculate_bounds(what_I_plot)  ## vmin=min(what_I_plot_frame), vmax=max(what_I_plot_frame))
+        what_I_plot = myplot.bounds(what_I_plot,zevmin,zevmax)
+        ## define contour field levels. define color palette
+        ticks = ndiv + 1
+        zelevels = np.linspace(zevmin,zevmax,ticks)
+        palette = get_cmap(name=colorb)
+        ## contours topo
+        zelevc = np.linspace(-8000.,20000.,20)
+        yeah.contour( x, y, np.transpose(fieldc), zelevc, colors='black',linewidths = 0.4)
+        # contour field
+        c = yeah.contourf( x, y, what_I_plot, zelevels, cmap = palette, alpha = trans )
+        Figure.colorbar(fig,c,orientation='vertical',format="%.1e")
+        ax = fig.gca() ; ax.set_title(fieldlab) ; ax.set_ylabel("Latitude") ; ax.set_xlabel("Longitude")
+        ax.set_xticks(np.arange(-180,181,45)) ; ax.set_xbound(lower=-180, upper=180)
+        ax.set_yticks(np.arange(-90,91,30)) ; ax.set_ybound(lower=-90, upper=90)
+        if incwind:
+          [x2d,y2d] = np.meshgrid(x,y)
+          yeah.quiver(x2d,y2d,np.transpose(windx),np.transpose(windy))
+      fig.text(0.5, 0.01, "Mars Climate Database (c) LMD/OU/IAA/ESA/CNES", ha='center')
+      canvas = FigureCanvasAgg(fig)
+      # The size * the dpi gives the final image size
+      #   a4"x4" image * 80 dpi ==> 320x320 pixel image
+      canvas.print_figure(figname, dpi=80)
+
+
     ### TODO: makeplot2d, plot2d, passer plot settings
 

trunk/UTIL/PYTHON/mcd/proto/cgi-bin/mcdcgi.py

@@ -1 @@
-#!/home/aymeric/Software/epd-7.0-2-rh5-x86/bin/python
+#!/usr/bin/python
+###!/home/aymeric/Software/epd-7.0-2-rh5-x86/bin/python
+####!/home/marshttp/EPD/epd-7.0-2-rh5-x86_64/bin/python
 ### here the version used to f2py the MCD Fortran routines
 
@@ -9 +11 @@
 ### (see mcdtest.py for examples of use)       ###
 ##################################################
+### ajouts et corrections par Franck Guyon 09/2012
 
 import cgi, cgitb 
 import numpy as np
-from mcd import mcd
+#from mcd import mcd
+from modules import *
+from modules import mcd
+
 import cStringIO
 import os as daos
 import matplotlib.pyplot as mpl
-import Image
+from PIL import Image
 
 # for debugging in web browser
@@ -25 +31 @@
 
 # create a MCD object
-query = mcd()
+#query = mcd()
+query=mcd.mcd() #FG: import from module mcd
 
 # Get data from user-defined fields and define free dimensions
+# FG: add tests if var==None to have values in local without forms ones
 getlat = form.getvalue("latitude")
+if getlat == None: getlat = 1
+
 if getlat == "all":  islatfree = 1 ; query.lat = -9999.
 else:                islatfree = 0 ; query.lat = float(getlat)
 getlon = form.getvalue("longitude")
+if getlon == None: getlon = 1
 if getlon == "all":  islonfree = 1 ; query.lon = -9999.
 else:                islonfree = 0 ; query.lon = float(getlon)
+
 getloct = form.getvalue("localtime")
+
+if getloct == None: getloct = 1
 if getloct == "all": isloctfree = 1 ; query.loct = -9999.
 else:                isloctfree = 0 ; query.loct = float(getloct)
+
 getalt = form.getvalue("altitude")
+if getalt == None: getalt = 1
 if getalt == "all":  isaltfree = 1 ; query.xz = -9999.
 else:                isaltfree = 0 ; query.xz = float(getalt)
+
 sumfree = islatfree + islonfree + isloctfree + isaltfree 
 if sumfree > 2: exit() ## only 1D or 2D plots for the moment
-query.xdate = float(form.getvalue("ls"))
-query.hrkey = int(form.getvalue("hrkey"))
-query.dust = int(form.getvalue("dust"))
+
+try: query.xdate = float(form.getvalue("ls"))
+except: query.xdate = float(1)
+try: query.hrkey = int(form.getvalue("hrkey"))
+except: query.hrkey = int(1)
+try: query.dust = int(form.getvalue("dust"))
+except: query.dust  = int(1)
 #        self.zkey      = 3  # specify that xz is the altitude above surface (m)
 #        self.perturkey = 0  #integer perturkey ! perturbation type (0: none)
@@ -55 +76 @@
 var3 = form.getvalue("var3")
 var4 = form.getvalue("var4")
-vartoplot = [var1]
+
+# fg: vartoplot is not None without form values
+# vartoplot = [var1]
+# fg: init var as with form values
+if var1 == None: var1="t"
+if var2 == None: var2="p"
+
+vartoplot = []
+if var1 != "none": vartoplot = np.append(vartoplot,var1)
 if var2 != "none": vartoplot = np.append(vartoplot,var2)
-if var3 != "none": vartoplot = np.append(vartoplot,var3)
-if var4 != "none": vartoplot = np.append(vartoplot,var4)
+if var3 != "none" and var3 != None: vartoplot = np.append(vartoplot,var3)
+if var4 != "none" and var4 != None: vartoplot = np.append(vartoplot,var4)
+
 iswind = form.getvalue("iswind")
 if iswind == "on": iswindlog = True
@@ -68 +98 @@
 # reference name (to test which figures are already in the database)
 reference = str(islatfree)+str(islonfree)+str(isloctfree)+str(isaltfree)+query.getnameset()+str(var1)+str(var2)+str(var3)+str(var4)+str(iswind)+str(isfixedlt)
-figname = 'img/'+reference+'.jpg'
+figname = '../img/'+reference+'.png'
+
 testexist = daos.path.isfile(figname)
 
@@ -85 +116 @@
 
     ### generic building of figure
-    #query.plot1d(["t","p","u","v"],vertplot=isaltfree)
-    query.plot1d(vartoplot,vertplot=isaltfree)
-    mpl.savefig("img/temp.png",dpi=85,bbox_inches='tight',pad_inches=0.25)
-    Image.open("img/temp.png").save(figname,'JPEG')
+    query.htmlplot1d(vartoplot,vertplot=isaltfree,figname=figname)
+    #mpl.savefig("img/temp.png",dpi=85,bbox_inches='tight',pad_inches=0.25)
+    #Image.open("../img/temp.png").save(figname,'JPEG')
 
   ### 2D plots
@@ -104 +134 @@
 
 ## This is quite common
-print "Content-type:text/html\r\n\r\n"
-print "<html>"
-print "<head>"
-print "<title>MCD. Simple Python interface</title>"
-print "</head>"
-print "<body>"
+print "Content-type:text/html\n"
+print "  "  #Apache needs a space after content-type
+
+#entete="""<?xml version="1.0" encoding="UTF-8"?> 
+#<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.1 plus MathML 2.0//EN" "http://www.w3.org/Math/DTD/mathml2/xhtml-math11-f.dtd">
+#<html xmlns="http://www.w3.org/1999/xhtml"> """
+
+header="""<html><head><title>Mars Climate Database: The Web Interface</title></head><body>"""
+
+print header
 
 ## Now the part which differs
 if sumfree == 0:        query.update() ; query.htmlprinttabextvar(vartoplot)  #query.printmeanvar()
-elif sumfree >= 1:      print "<img src='../"+figname+"'><br />"
-else:                   exit()
+elif sumfree >= 1:      print "<img src='"+figname+"'><br />"
+else:                   print "<h1>ERROR : sumfree is not or badly defined ...</h1></body></html>"
+
 
 ## This is quite common
-#print "Based on the <a href='http://www-mars.lmd.jussieu.fr'>Mars Climate Database</a> (c) LMD/OU/IAA/ESA/CNES.<br />"
-print "<hr>"
-print "<a href='../index.html'>Click here to start a new query</a>."
-#query.printset()
-print "<hr>"
-print "</body>"
-print "</html>"
+bottom = "<hr><a href='../index.html'>Click here to start a new query</a>.<hr></body></html>"
+print bottom
 
 ##write to file object

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

@@ -15 +15 @@
 <!-- aussi possible: get a la place de post. pour avoir un beau lien -->
 
+<!--<form name="calendar" action="/marscgi-bin/mcdcgi.py" method="post"> <!--target="_new">-->
 <form name="calendar" action="./cgi-bin/mcdcgi.py" method="post"> <!--target="_new">-->
 
@@ -315 +316 @@
 <td align="center">
 <input type="submit" value="SUBMIT" style="font-weight:bold"/>
+<input type="button" value="RESET" style="font-weight:bold" onClick="DefaultDateValues();DefaultTimeValues();Convert2Ls();PlaceValues(0.,0.);DefaultSpaceTime()"/><br />
 </td>
 <td align="center">