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

Timestamp:

Jul 7, 2025, 1:39:13 PM (2 days 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/ch4vmrsection.py (modified) (2 diffs)

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

@@ -8 +8 @@
 import matplotlib.colors as mcolors
 from FV3_utils import *
+from input import *
 
 ############################
-folder="../"
-filename1=folder+"diagfi_mean_A.nc"
 var="ch4_gas" #variable
 tint=[30,32] #Time must be as written in the input file
 xarea="-180,179"
 
-nc1=Dataset(filename1)
+nc1=Dataset(name+"_A.nc")
 
-lat=nc1.variables["lat"][:]
-lon=nc1.variables["lon"][:]
-alt=nc1.variables["altitude"][:]
-tim=nc1.variables["time_counter"][:]
+lat=getvar(nc1,"latitude")
+lon=getvar(nc1,"longitude")
+alt=getvar(nc1,"altitude")
+tim=getvar(nc1,"Time")
 ############################
 
-# def getvar(filename,var,tint,xarea):
-#     myvar = pp(file=filename,var=var,t=tint,x=xarea,compute="mean").getf()
-#     print((shape(myvar)))
-    # return myvar
-
-myvar=getvar(nc1,var,tint,tim)
+myvar=getvar(nc1,var,tint, l_mean=True, t_mean=True)
 myvar=myvar*28/16.*100.
 
@@ -42 +36 @@
 xticks=[-90,-60,-30,0,30,60,90]
 #yticks=np.linspace(0,240,9)
-alt=alt/1000.
 
 mymin=0.1