Changeset 3833 for trunk/LMDZ.PLUTO/util/script_figures/temp_time_zoom.py

Timestamp:

Jul 7, 2025, 1:39:13 PM (24 hours ago)

Author:

afalco

Message:

Pluto: updated plots scripts.
Fixed some issues with reading XIOS, etc.
Included display_netcdf.py tool from Mars PCM.
AF

File:

• trunk/LMDZ.PLUTO/util/script_figures/temp_time_zoom.py (modified) (4 diffs)

trunk/LMDZ.PLUTO/util/script_figures/temp_time_zoom.py

@@ -1 +1 @@
 #! /usr/bin/env python
-from ppclass import pp
 from    netCDF4               import    Dataset
 from    numpy                 import    *
@@ -11 +10 @@
 import datetime
 from mpl_toolkits.basemap import Basemap, shiftgrid
+from FV3_utils import *
 
 ############################
-filename1="diagfi2015_S.nc"
+filename1=name+"_S.nc"
 var="temperature" #variable
-xarea="-169,-165"
-yarea="-19,-15"
+longitude=[-169,-165]
+# longitude=-165
+latitude=[-19,-15]
+# latitude=-19
 
 # local time de la longitude consideree a t=0
 loct=12
-sol0=30
+sol0=12
 t0=1./24*loct
-t1=t0+1
-tint=[str(sol0+t0)+','+str(sol0+t1)] #Time must be as written in the input file
+t1=t0+12
+tint=[sol0+t0,sol0+t1] #Time must be as written in the input file
 print(tint)
 nc1=Dataset(filename1)
 
-lat=nc1.variables["lat"][:]
-lon=nc1.variables["lon"][:]
-alt=nc1.variables["altitude"][:]
+lat=getvar(nc1,"latitude")
+lon=getvar(nc1,"longitude")
+alt=getvar(nc1,"altitude")
+tim=getvar(nc1,"Time",times=tint)
 ############################
-
-def getvar(filename,var,tint,xarea,yarea):
-    myvar = pp(file=filename,var=var,t=tint,x=xarea,y=yarea,compute="nothing").getf()
-    print(('shape myvar = ',shape(myvar)))
-    return myvar
-
 
 mpl.figure(figsize=(18, 10))
 
-
-myvar=getvar(filename1,var,tint,xarea,yarea)[:,:,0,0]
+myvar=getvar(nc1,var,times=tint,longitudes=longitude,latitudes=latitude)[:,:,0,0]
 font=23
-tim=np.linspace(0,24,9)
+# tim=np.linspace(0,24,9)
 print(("tim=",tim))
 print(('on prend les premiers indice, shape (tmps, alt, var) =',shape(tim), shape(alt), shape(myvar)))
@@ -49 +45 @@
 pal=get_cmap(name="Spectral_r")
 lev=np.linspace(40,50,10)
-xticks=[0,2,4,6,8,10,12,14,16,18,20,22,24] 
+# xticks=[0,2,4,6,8,10,12,14,16,18,20,22,24] 
 print(('hello:',np.linspace(0,24,13)))
 #yticks=np.linspace(0,240,9)
-alt=alt/1000.
 
 
@@ -72 +67 @@
 mpl.xlabel('Local Time (h)',labelpad=10,fontsize=font)
 mpl.ylabel('Altitude (km)',labelpad=10, fontsize=font)
-mpl.xticks(xticks,fontsize=font)
-#mpl.xticks(fontsize=font)
+# mpl.xticks(xticks,fontsize=font)
+mpl.xticks(fontsize=font)
 #mpl.yticks(yticks,fontsize=font)
 mpl.yticks(fontsize=font)