Changeset 3833 for trunk/LMDZ.PLUTO/util/script_figures/movie_winds

Timestamp:

Jul 7, 2025, 1:39:13 PM (4 weeks ago)

Author:

afalco

Message:

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

Location:

trunk/LMDZ.PLUTO/util/script_figures/movie_winds

Files:

• anim.sh (added)
• input.py (added)
• tsurf_winds.py (modified) (5 diffs)

trunk/LMDZ.PLUTO/util/script_figures/movie_winds/tsurf_winds.py

@@ -16 +16 @@
 filename1="../"+name+"_A.nc"
 filename2="../"+name+".nc"
-filename3="../../phisinit.nc"
+filename3="../"+name+".nc"
+# filename3="../../phisinit.nc"
 var="tsurf" #variable
 phisinit="phisinit" #variable
@@ -87 +88 @@
     pal=get_cmap(name="jet")
     lev=np.linspace(37,51,15)
-    # newlon=lon+180
-    newlon=lon 
+    newlon=lon+180
     CF=mpl.contourf(newlon, lat, myvarbis,lev,cmap=pal,extend='both')
     yticks=[-90,-60,-30,0,30,60,90]
@@ -98 +98 @@
         t.set_fontsize(font)
 
-    c=mpl.contour(newlon, lat, phisinit2bis, 10,levels=np.linspace(-4,4,8), colors = 'k', linewidths = 3.5)
+    c=mpl.contour(newlon, lat, myvar2bis, 10,levels=np.linspace(-4,4,8), colors = 'k', linewidths = 3.5)
     mpl.clabel(c, fmt='%2.1f',inline=1, colors='k', fontsize=23,inline_spacing=1)
     #mpl.title('Local Time at Sputnik Planum='+str(i*3)+'H00',fontsize=font)
@@ -107 +107 @@
     getwinds(newlon,lat,u,v)
 
-# def getnumalt(choicealt,alt):
-#     numalt=np.where(abs(alt-choicealt)==min(abs(alt-choicealt)))
-#     numalt=numalt[0][0]
-#     return numalt
+def getnumalt(choicealt,alt):
+    numalt=np.where(abs(alt-choicealt)==min(abs(alt-choicealt)))
+    numalt=numalt[0][0]
+    return numalt
 
 #######################
-# numalt=getnumalt(30,alt)
-numalt=np.searchsorted(alt[...],30)
+numalt=getnumalt(30,alt)
 print(('numalt =',numalt,'altitude=',alt[numalt]))
 uini=getvar(nc1,vari,tint,tim)[:,numalt]
 vini=getvar(nc1,varj,tint,tim)[:,numalt]
 myvar=getvar(nc2,var,tint,tim)
-phisinit2=getvar(nc3,phisinit) # phisinitmyvar2=phisinit2/0.6169/1000.  # altitude km
+myvar2=getvar(nc3,phisinit) # phisinitmyvar2=myvar2/0.6169/1000.  # altitude km
 nbfig=uini.shape[0]
 print(("nbfig=",nbfig))
-phisinit2bis=np.copy(phisinit2)
+myvar2bis=switchlon(myvar2)
 
 for i in range(nbfig):
@@ -129 +128 @@
    v2=vini[i,:,:]
    myv=myvar[i,:,:]
-   u=np.copy(u2)
-   v=np.copy(v2)
-   myvarbis=np.copy(myv)
-   print(i)
+   u=switchlon(u2)
+   v=switchlon(v2)
+   myvarbis=switchlon(myv)
+   print(i,"/",nbfig)
    getfigvar(i)
    mpl.savefig('mapwinds'+str('{0:03}'.format(i))+'.eps',dpi=200)