Changeset 3833 for trunk/LMDZ.PLUTO/util/script_figures/map_alb_ini.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/map_alb_ini.py (modified) (3 diffs)

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

@@ -1 +1 @@
 #! /usr/bin/env python
-from ppclass import pp
 from    netCDF4               import    Dataset
 from    numpy                 import    *
@@ -10 +9 @@
 import matplotlib.colors as mcolors
 from mpl_toolkits.basemap import Basemap
+from FV3_utils import *
 
 ############################
-filename="diagfi2015.nc"
+filename=name+".nc"
 nb_dataset=1  #change vect as well 
-tint=["0,0.001"] #Time must be as written in the input file
+tint=[0,0.001] #Time must be as written in the input file
 #zint=["0"] #alt in km
 #xarea="-180,179"
 #yarea="-90,90"
 #step=5  #step between each bin
-var="albedo" #variable
+var="ALB" #variable
 var2="phisinit"
 ############################
 
 nc=Dataset(filename)
-lat=nc.variables["lat"][:]
-lon=nc.variables["lon"][:]
+lat=getvar(nc,"latitude")
+lon=getvar(nc,"longitude")
+tim=getvar(nc,"Time")
 
-myvar = pp(file=filename,var=var,t=tint,compute="mean").getf()  # get data to be changed according to selected variable
-myvar2 = pp(file=filename,var=var2,t=tint,compute="mean").getf()  # get data to be changed according to selected variable
+myvar = getvar(nc,var,tint,t_mean=True)
+myvar2 = getvar(nc,var2)
 print(('shape var1 =',shape(myvar)))
 print(('shape lat = ',shape(lat)))
@@ -90 +91 @@
            myvar2bis[i,j]=myvar2[i,j-int(vec[1]/2)]
 
-lon=lon+180.
+if lon[0]<0:
+    lon=lon+180.