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

Timestamp:

Jul 7, 2025, 1:39:13 PM (36 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/proftempNH.py (modified) (5 diffs)

trunk/LMDZ.PLUTO/util/script_figures/proftempNH.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 *
 
 ############################
-zefile="diagfi2015_A.nc"
-d1="restart_ref/"
-
-f1=d1+zefile
-
-
 var="temperature" #variable
 
-nc1=Dataset(f1)
-lat=nc1.variables["lat"][:]
-lon=nc1.variables["lon"][:]
-alt=nc1.variables["altitude"][:]
-tim=nc1.variables["time_counter"][:]
+nc1=Dataset(name+"_A.nc")
+lat=getvar(nc1,"latitude")
+lon=getvar(nc1,"longitude")
+alt=getvar(nc1,"altitude")
+tim=getvar(nc1,"Time")
 
 print(('Time = ',tim))
@@ -43 +37 @@
     print(('Point =',p,' Time=',tim[indt1]))
     return indt1
- 
-def getvar(filename,var):
-    myvar = pp(file=filename,var=var,compute="nothing").getf()
-    return myvar
 
 def getindex(lat,lon,mylat,mylon):
@@ -57 +47 @@
     indt=findindextime(tini,lt0,ltchoice,p)
     indlat,indlon=getindex(lat,lon,p[1],p[0])
-    myvar=getvar(f1,var)[indt,:,indlat,indlon]
+    myvar=getvar(nc1,var)[indt,:,indlat,indlon]
     return myvar
 ############################
@@ -92 +82 @@
 font=23
 
-alt=alt/1000.
-
 mpl.plot(myvar1,alt,'r',label='Entry')
 mpl.plot(myvar2,alt,'g',label='Entry-like on edge')