| 1 | #!/usr/bin/python |
|---|
| 2 | ###!/usr/bin/env python |
|---|
| 3 | ###!/home/aymeric/Software/epd-7.0-2-rh5-x86/bin/python |
|---|
| 4 | ####!/home/marshttp/EPD/epd-7.0-2-rh5-x86_64/bin/python |
|---|
| 5 | ### here the version used to f2py the MCD Fortran routines |
|---|
| 6 | |
|---|
| 7 | ################################################## |
|---|
| 8 | ### A Python CGI for the Mars Climate Database ### |
|---|
| 9 | ### ------------------------------------------ ### |
|---|
| 10 | ### Aymeric SPIGA 18-19/04/2012 ~ 11/08/2012 ### |
|---|
| 11 | ### ------------------------------------------ ### |
|---|
| 12 | ### (see mcdtest.py for examples of use) ### |
|---|
| 13 | ################################################## |
|---|
| 14 | ### ajouts et corrections par Franck Guyon 09/2012 |
|---|
| 15 | ### ajouts suite a brainstorm equipe AS 10/2012 |
|---|
| 16 | |
|---|
| 17 | import cgi, cgitb |
|---|
| 18 | import numpy as np |
|---|
| 19 | from modules import * |
|---|
| 20 | from modules import mcd |
|---|
| 21 | |
|---|
| 22 | import cStringIO |
|---|
| 23 | import os as daos |
|---|
| 24 | import matplotlib.pyplot as mpl |
|---|
| 25 | |
|---|
| 26 | ### a function to read HTML arguments for coordinates |
|---|
| 27 | def gethtmlcoord(userinput,defmin,defmax): |
|---|
| 28 | import string |
|---|
| 29 | # accepted separators. the symbol - should always be last. |
|---|
| 30 | separators = [":",";",",","/","_"] |
|---|
| 31 | # initial values |
|---|
| 32 | val = -9999. ; vals = None ; vale = None ; foundinterv = False |
|---|
| 33 | if userinput == None: userinput = "1" |
|---|
| 34 | # the main work. either all, or an interval, or a single value. |
|---|
| 35 | # ... all |
|---|
| 36 | if userinput == "all": isfree = 1 ; vals = defmin ; vale = defmax ; foundinterv = True |
|---|
| 37 | # ... an interval |
|---|
| 38 | else: |
|---|
| 39 | for sep in separators: |
|---|
| 40 | if not foundinterv: |
|---|
| 41 | isfree = 1 ; ind = userinput.find(sep) |
|---|
| 42 | if ind != -1: vals = float(userinput[:ind]) ; vale = float(userinput[ind+1:]) ; foundinterv = True |
|---|
| 43 | # ... a single value (or an error) |
|---|
| 44 | if not foundinterv: |
|---|
| 45 | # treat the difficult case of possible - separator |
|---|
| 46 | test = userinput[1:].find("-") # at this stage: |
|---|
| 47 | # * if - is found in the first position, it could not be a separator |
|---|
| 48 | # * if - found at positions > 0, it must be considered as a separator |
|---|
| 49 | if test != -1: |
|---|
| 50 | isfree = 1 ; ind=test+1 |
|---|
| 51 | vals = float(userinput[:ind]) ; vale = float(userinput[ind+1:]) ; foundinterv = True |
|---|
| 52 | else: |
|---|
| 53 | # check if input is valid (each character is numeric or -) |
|---|
| 54 | for char in userinput: |
|---|
| 55 | if char not in string.digits: |
|---|
| 56 | if char not in ["-","."]: userinput="yorgl" |
|---|
| 57 | # either we are OK. if we are not we set isfree to -1. |
|---|
| 58 | if userinput != "yorgl": isfree = 0 ; val = float(userinput) |
|---|
| 59 | else: isfree = -1 |
|---|
| 60 | # return values |
|---|
| 61 | return isfree, val, vals, vale |
|---|
| 62 | |
|---|
| 63 | # set an errormess variable which must stay to None for interface to proceed |
|---|
| 64 | errormess = "" |
|---|
| 65 | |
|---|
| 66 | # for debugging in web browser |
|---|
| 67 | cgitb.enable() |
|---|
| 68 | |
|---|
| 69 | # Create instance of FieldStorage |
|---|
| 70 | form = cgi.FieldStorage() |
|---|
| 71 | |
|---|
| 72 | # create a MCD object |
|---|
| 73 | query=mcd.mcd() |
|---|
| 74 | |
|---|
| 75 | # set MCD version changes if needed |
|---|
| 76 | try: betatest = form.getvalue("betatest") |
|---|
| 77 | except: betatest = "off" |
|---|
| 78 | if betatest == "on": query.toversion5() |
|---|
| 79 | |
|---|
| 80 | # Get the kind of vertical coordinates and choose default behavior for "all" |
|---|
| 81 | try: query.zkey = int(form.getvalue("zkey")) |
|---|
| 82 | except: query.zkey = int(3) |
|---|
| 83 | if query.zkey == 2: minxz = -5000. ; maxxz = 150000. |
|---|
| 84 | elif query.zkey == 3: minxz = 0. ; maxxz = 250000. |
|---|
| 85 | elif query.zkey == 5: minxz = -5000. ; maxxz = 150000. |
|---|
| 86 | elif query.zkey == 4: minxz = 1.e3 ; maxxz = 1.e-6 |
|---|
| 87 | elif query.zkey == 1: minxz = 3396000. ; maxxz = 3596000. |
|---|
| 88 | |
|---|
| 89 | # Get data from user-defined fields and define free dimensions |
|---|
| 90 | islatfree, query.lat, query.lats, query.late = gethtmlcoord( form.getvalue("latitude"), -90., 90. ) |
|---|
| 91 | islonfree, query.lon, query.lons, query.lone = gethtmlcoord( form.getvalue("longitude"), -180., 180. ) |
|---|
| 92 | isloctfree, query.loct, query.locts, query.locte = gethtmlcoord( form.getvalue("localtime"), 0., 24. ) |
|---|
| 93 | isaltfree, query.xz, query.xzs, query.xze = gethtmlcoord( form.getvalue("altitude"), minxz, maxxz) |
|---|
| 94 | if minxz < 0.1: minxz=0.1 # otherwise bug with values smaller than 0.1m |
|---|
| 95 | |
|---|
| 96 | try: query.datekey = int(form.getvalue("datekeyhtml")) |
|---|
| 97 | except: query.datekey = float(1) |
|---|
| 98 | badlschar = False |
|---|
| 99 | if query.datekey == 1: |
|---|
| 100 | try: query.xdate = float(form.getvalue("ls")) |
|---|
| 101 | except: query.xdate = float(1) ; badlschar = True # comment the second part if in debug command line mode |
|---|
| 102 | else: |
|---|
| 103 | try: query.xdate = float(form.getvalue("julian")) |
|---|
| 104 | except: query.xdate = float(1) |
|---|
| 105 | query.loct = 0. |
|---|
| 106 | try: query.dust = int(form.getvalue("dust")) |
|---|
| 107 | except: query.dust = int(1) |
|---|
| 108 | |
|---|
| 109 | # Prevent the user from doing bad |
|---|
| 110 | badinterv = (islatfree == -1) or (islonfree == -1) or (isloctfree == -1) or (isaltfree == -1) |
|---|
| 111 | if badinterv: |
|---|
| 112 | errormess = errormess+"<li>Bad syntax. Write a value (or) a range val1;val2 (or) 'all'. Separator shall be either ; : , / _ " |
|---|
| 113 | badls = (query.datekey == 1 and (query.xdate < 0. or query.xdate > 360.)) |
|---|
| 114 | if badls: |
|---|
| 115 | errormess = errormess+"<li>Solar longitude must be between 0 and 360." |
|---|
| 116 | if badlschar: |
|---|
| 117 | errormess = errormess+"<li>Solar longitude is in the wrong format. It should be a positive number between 0 and 360. Intervals of solar longitude are not allowed." |
|---|
| 118 | badloct = (isloctfree == 0 and query.loct > 24.) \ |
|---|
| 119 | or (isloctfree == 1 and (query.locts > 24. or query.locte > 24.)) \ |
|---|
| 120 | or (isloctfree == 0 and query.loct < 0.) \ |
|---|
| 121 | or (isloctfree == 1 and (query.locts < 0. or query.locte < 0.)) |
|---|
| 122 | if badloct: |
|---|
| 123 | errormess = errormess+"<li>Local time must be less than 24 martian hours (and not a negative number)." |
|---|
| 124 | badlat = (islatfree == 0 and abs(query.lat) > 90.) \ |
|---|
| 125 | or (islatfree == 1 and (abs(query.lats) > 90. or abs(query.late) > 90.)) |
|---|
| 126 | if badlat: |
|---|
| 127 | errormess = errormess+"<li>Latitude coordinates must be between -90 and 90." |
|---|
| 128 | badlon = (islonfree == 0 and abs(query.lon) > 360.) \ |
|---|
| 129 | or (islonfree == 1 and (abs(query.lons) > 360. or abs(query.lone) > 360.)) |
|---|
| 130 | if badlon: |
|---|
| 131 | errormess = errormess+"<li>Longitude coordinates must be between -360 and 360." |
|---|
| 132 | badalt = (isaltfree == 0 and (query.zkey in [3]) and query.xz < 0.) \ |
|---|
| 133 | or (isaltfree == 1 and (query.zkey in [3]) and (query.xzs < 0. or query.xze < 0.)) |
|---|
| 134 | if badalt: |
|---|
| 135 | errormess = errormess+"<li>Vertical coordinates must be positive when requesting altitude above surface." |
|---|
| 136 | badalt2 = (isaltfree == 0 and (query.zkey in [1,4]) and query.xz <= 0.) \ |
|---|
| 137 | or (isaltfree == 1 and (query.zkey in [1,4]) and (query.xzs <= 0. or query.xze <= 0.)) |
|---|
| 138 | if badalt2: |
|---|
| 139 | errormess = errormess+"<li>Vertical coordinates must be <b>strictly</b> positive when requesting pressure levels or altitude above Mars center." |
|---|
| 140 | badalt3 = (isaltfree == 0 and query.zkey == 4 and query.xz > 1500.) \ |
|---|
| 141 | or (isaltfree == 1 and query.zkey == 4 and min(query.xzs,query.xze) > 1500.) |
|---|
| 142 | if badalt3: |
|---|
| 143 | errormess = errormess+"<li>Pressure values larger than 1500 Pa are unlikely to be encountered in the Martian atmosphere." |
|---|
| 144 | badrange = (isloctfree == 1 and query.locts == query.locte) \ |
|---|
| 145 | or (islatfree == 1 and query.lats == query.late) \ |
|---|
| 146 | or (islonfree == 1 and query.lons == query.lone) \ |
|---|
| 147 | or (isaltfree == 1 and query.xzs == query.xze) |
|---|
| 148 | if badrange: |
|---|
| 149 | errormess = errormess+"<li>One or several coordinate intervals are not... intervals. Set either a real range or an unique value." |
|---|
| 150 | stormls = ( (query.dust in [4,5,6]) and (query.datekey == 1 and query.xdate < 180.)) |
|---|
| 151 | if stormls: |
|---|
| 152 | errormess = errormess+"<li>When a dust storm scenario is selected, available dates must be within the dust storm season (180 < Ls < 360)." |
|---|
| 153 | if query.xdate == 666.: |
|---|
| 154 | errormess = "<li>CONGRATULATIONS! <br><img src='../surprise.jpg'><br> You reached secret mode.<br> You can <a href='http://www.youtube.com/watch?v=fTpQOZcNASw'>watch a nice video</a>." |
|---|
| 155 | |
|---|
| 156 | # Get how many free dimensions we have |
|---|
| 157 | sumfree = islatfree + islonfree + isloctfree + isaltfree |
|---|
| 158 | if sumfree >= 3: errormess = errormess + "<li>3 or more free dimensions are set... but only 1D and 2D plots are supported!" |
|---|
| 159 | |
|---|
| 160 | # Get additional parameters |
|---|
| 161 | try: query.hrkey = int(form.getvalue("hrkey")) |
|---|
| 162 | except: query.hrkey = int(1) |
|---|
| 163 | # self.perturkey = 0 #integer perturkey ! perturbation type (0: none) |
|---|
| 164 | # self.seedin = 1 #random number generator seed (unused if perturkey=0) |
|---|
| 165 | # self.gwlength = 0. #gravity Wave wavelength (unused if perturkey=0) |
|---|
| 166 | try: query.colorm = form.getvalue("colorm") |
|---|
| 167 | except: query.colorm = "jet" |
|---|
| 168 | |
|---|
| 169 | try: query.min2d = float(form.getvalue("minval")) |
|---|
| 170 | except: query.min2d = None |
|---|
| 171 | try: query.max2d = float(form.getvalue("maxval")) |
|---|
| 172 | except: query.max2d = None |
|---|
| 173 | |
|---|
| 174 | try: |
|---|
| 175 | query.dpi = form.getvalue("dpi") |
|---|
| 176 | if query.dpi == "eps": yeaheps = True ; query.dpi = 300. |
|---|
| 177 | else: yeaheps = False ; query.dpi = float(query.dpi) |
|---|
| 178 | except: |
|---|
| 179 | query.dpi = 80 |
|---|
| 180 | yeaheps = False |
|---|
| 181 | |
|---|
| 182 | # Get variables to plot |
|---|
| 183 | var1 = form.getvalue("var1") |
|---|
| 184 | var2 = form.getvalue("var2") |
|---|
| 185 | var3 = form.getvalue("var3") |
|---|
| 186 | var4 = form.getvalue("var4") |
|---|
| 187 | |
|---|
| 188 | # fg: init var as with form values |
|---|
| 189 | if var1 == None: var1="t" |
|---|
| 190 | |
|---|
| 191 | vartoplot = [] |
|---|
| 192 | if var1 != "none": vartoplot = np.append(vartoplot,var1) |
|---|
| 193 | if var2 != "none" and var2 != None: vartoplot = np.append(vartoplot,var2) |
|---|
| 194 | if var3 != "none" and var3 != None: vartoplot = np.append(vartoplot,var3) |
|---|
| 195 | if var4 != "none" and var4 != None: vartoplot = np.append(vartoplot,var4) |
|---|
| 196 | |
|---|
| 197 | iswind = form.getvalue("iswind") |
|---|
| 198 | if iswind == "on": iswindlog = True |
|---|
| 199 | else: iswindlog = False |
|---|
| 200 | isfixedlt = form.getvalue("isfixedlt") |
|---|
| 201 | if isfixedlt == "on": query.fixedlt=True |
|---|
| 202 | else: query.fixedlt=False |
|---|
| 203 | iszonmean = form.getvalue("zonmean") |
|---|
| 204 | if iszonmean == "on": query.zonmean=True |
|---|
| 205 | else: query.zonmean=False |
|---|
| 206 | |
|---|
| 207 | islog = form.getvalue("islog") |
|---|
| 208 | if islog == "on": query.islog=True |
|---|
| 209 | else: query.islog=False |
|---|
| 210 | |
|---|
| 211 | ### now, proceed... |
|---|
| 212 | if errormess == "": |
|---|
| 213 | |
|---|
| 214 | # reference name (to test which figures are already in the database) |
|---|
| 215 | try: reference = query.getnameset()+str(var1)+str(var2)+str(var3)+str(var4)+str(iswind)+str(isfixedlt)+str(iszonmean)+query.colorm+str(query.min2d)+str(query.max2d)+str(query.dpi)+str(islog) |
|---|
| 216 | except: reference = "test" |
|---|
| 217 | if yeaheps: figname = '../img/'+reference+'.eps' |
|---|
| 218 | else: figname = '../img/'+reference+'.png' |
|---|
| 219 | txtname = '../txt/'+reference+'.txt' |
|---|
| 220 | testexist = daos.path.isfile(figname) |
|---|
| 221 | |
|---|
| 222 | # extract data from MCD if needed |
|---|
| 223 | if not testexist: |
|---|
| 224 | |
|---|
| 225 | ### 1D plots |
|---|
| 226 | if sumfree == 1: |
|---|
| 227 | |
|---|
| 228 | ### getting data |
|---|
| 229 | if isloctfree == 1: query.diurnal(nd=25) |
|---|
| 230 | elif islonfree == 1: query.zonal(nd=64) |
|---|
| 231 | elif islatfree == 1: query.meridional(nd=48) |
|---|
| 232 | elif isaltfree == 1: query.profile(nd=35) |
|---|
| 233 | else: exit() |
|---|
| 234 | |
|---|
| 235 | ### generic building of figure |
|---|
| 236 | query.getascii(vartoplot,filename=txtname) |
|---|
| 237 | query.htmlplot1d(vartoplot,figname=figname) |
|---|
| 238 | |
|---|
| 239 | ### 2D plots |
|---|
| 240 | elif sumfree == 2: |
|---|
| 241 | |
|---|
| 242 | ### getting data |
|---|
| 243 | if islatfree == 1 and islonfree == 1: query.htmlmap2d(vartoplot,incwind=iswindlog,figname=figname) |
|---|
| 244 | else: query.htmlplot2d(vartoplot,figname=figname) |
|---|
| 245 | |
|---|
| 246 | #### NOW WRITE THE HTML PAGE TO USER |
|---|
| 247 | |
|---|
| 248 | ## This is quite common |
|---|
| 249 | print "Content-type:text/html\n" |
|---|
| 250 | print " " #Apache needs a space after content-type |
|---|
| 251 | |
|---|
| 252 | #entete="""<?xml version="1.0" encoding="UTF-8"?> |
|---|
| 253 | #<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.1 plus MathML 2.0//EN" "http://www.w3.org/Math/DTD/mathml2/xhtml-math11-f.dtd"> |
|---|
| 254 | #<html xmlns="http://www.w3.org/1999/xhtml"> """ |
|---|
| 255 | |
|---|
| 256 | header="""<html><head><title>Mars Climate Database: The Web Interface</title></head><body>""" |
|---|
| 257 | #if betatest == "on": |
|---|
| 258 | # print "<b>!!! THIS IS A BETA VERSION. RESULTS ARE NOT VALIDATED !!!</b>" |
|---|
| 259 | # if sumfree == 2: print "<br>" |
|---|
| 260 | |
|---|
| 261 | print header |
|---|
| 262 | #print query.printset() |
|---|
| 263 | #print "<br />" |
|---|
| 264 | |
|---|
| 265 | ## Now the part which differs |
|---|
| 266 | if errormess != "": |
|---|
| 267 | print "<h1>Ooops!</h1>" |
|---|
| 268 | print "Please correct the following problems before submitting again." |
|---|
| 269 | print "<ul>" |
|---|
| 270 | print errormess |
|---|
| 271 | print "</ul>" |
|---|
| 272 | else: |
|---|
| 273 | if sumfree == 0: query.update() ; query.htmlprinttabextvar(vartoplot) |
|---|
| 274 | elif sumfree == 2: |
|---|
| 275 | if yeaheps: print "<hr><a href='"+figname+"'>!!!! Click here to download the EPS figure file !!!!</a><br /><hr>" |
|---|
| 276 | else: print "<img src='"+figname+"'><br />" |
|---|
| 277 | elif sumfree == 1: |
|---|
| 278 | print "<a href='"+txtname+"'>Click here to download an ASCII file containing data</a><br />" |
|---|
| 279 | print "<hr>" |
|---|
| 280 | if yeaheps: print "<hr><a href='"+figname+"'>!!!! Click here to download the EPS figure file !!!!</a><br /><hr>" |
|---|
| 281 | else: print "<img src='"+figname+"'><br />" |
|---|
| 282 | |
|---|
| 283 | ## This is quite common |
|---|
| 284 | bottom = "</body></html>" |
|---|
| 285 | print bottom |
|---|
