Changeset 782

Timestamp:

Sep 13, 2012, 12:29:15 PM (12 years ago)

Author:

acolaitis

Message:

UTIL PYTHON fixed mesoscale local time bugs introduced at revision 763. now working fine for LMD and MRAMS-generated wrfout files. plus polar domain now OK

Location:

trunk/UTIL/PYTHON

Files:

• myplot.py (modified) (2 diffs)
• planetoplot.py (modified) (2 diffs)

trunk/UTIL/PYTHON/myplot.py

@@ -30 +30 @@
 
 ## Author: AS + AC
-def localtime(time,lon,nc): # lon is the mean longitude of the plot, not of the domain. central lon of domain is taken from cen_lon
-    import numpy as np
-    dt_hour=1.
-    start=0.
+def localtime(time,lon,namefile): # lon is the mean longitude of the plot, not of the domain. central lon of domain is taken from cen_lon
+    import numpy as np
+    from netCDF4 import Dataset
+    ## THIS IS FOR MESOSCALE
+    nc = Dataset(namefile)
+    ## get start date and intervals
+    dt_hour=1. ; start=0.
+    if hasattr(nc,'TITLE'):
+        title=getattr(nc, 'TITLE')
+        if hasattr(nc,'DT') and hasattr(nc,'START_DATE') and 'MRAMS' in title: 
+            ## we must adapt what is done in getlschar to MRAMS (outputs from ic.py)
+            dt_hour=getattr(nc, 'DT')/60. 
+            start_date=getattr(nc, 'START_DATE')
+            start_hour=np.float(start_date[11:13])
+            start_minute=np.float(start_date[14:16])/60.
+            start=start_hour+start_minute # start is the local time of simu at longitude 0
+            #LMD MMM is 1 output/hour (and not 1 output/timestep)
+            #MRAMS is 1 output/timestep, unless stride is added in ic.py
+        elif 'WRF' in title: 
+            [dummy,start,dt_hour] = getlschar ( namefile ) # get start hour and interval hour
+    ## get longitude 
     if lon is not None:
        if lon[0,1]!=lon[0,0]: mean_lon_plot = 0.5*(lon[0,1]-lon[0,0])
@@ -39 +56 @@
     elif hasattr(nc, 'CEN_LON'): mean_lon_plot=getattr(nc, 'CEN_LON')
     else: mean_lon_plot=0.
-    if hasattr(nc,'TITLE'):
-       title=getattr(nc, 'TITLE')
-       if hasattr(nc,'DT') and 'MRAMS' in title: dt_hour=getattr(nc, 'DT')/60. #LMD MMM is 1 output/hour (and not 1 output/timestep)
-                        # MRAMS is 1 output/timestep, unless stride is added in ic.py
-       if hasattr(nc,'START_DATE'):
-          start_date=getattr(nc, 'START_DATE')
-          start_hour=np.float(start_date[11:13])
-          start_minute=np.float(start_date[14:16])/60.
-          start=start_hour+start_minute # start is the local time of simu at longitude 0
-          ltst = start + time*dt_hour + mean_lon_plot / 15.
-       else: ltst = time*dt_hour + mean_lon_plot / 15.
-    else: ltst = time*dt_hour + mean_lon_plot / 15.
+    ## calculate local time
+    ltst = start + time*dt_hour + mean_lon_plot / 15.
     ltst = int (ltst * 10) / 10.
     ltst = ltst % 24

trunk/UTIL/PYTHON/planetoplot.py

@@ -299 +299 @@
           if axtime in ["lt"]:
               ftime = np.zeros(len(time))
-              for i in range(len(time)): 
-                 ftime[i] = localtime ( time[i], slon , nc )
-              time=ftime
-              time=check_localtime(time)
+              for i in range(len(time)): ftime[i] = localtime ( time[i], slon , namefile )
+              time=ftime ; time=check_localtime(time)
               print "LOCAL TIMES.... ", time
           ### we define the vertical axis (or lack thereof) and cover possibilities for it to be altitude, pressure, geopotential. quite SPECIFIC.
@@ -463 +461 @@
        else:                                 indextime = getsindex(stime,((nplot-1)//(nlon*nlat*nvert))%ntime,time)
        ltst = None 
-       if typefile in ['meso'] and indextime is not None and len(indextime) < 2: ltst = localtime ( indextime, slon , nc)
+       if typefile in ['meso'] and indextime is not None and len(indextime) < 2: ltst = localtime ( indextime, slon , all_namefile[index_f])
        print "********** INDEX LON:",indexlon," LAT:",indexlat," VERT:",indexvert," TIME:",indextime
        ##var = nc.variables["phisinit"][:,:]