Changeset 310 for trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/2dplot.py

Timestamp:

Sep 29, 2011, 3:59:01 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE:
-- Code in storm scenario corresponds to 'OMEGA' reference case [Julien Faure]
-- Easier settings for dust lifting without the need to recompile [see below]
-- A few modifications to plot and dust-devil detection PYTHON routines

29/09/11 == AS

--> To easily explore sensitivity to lifting thresholds: in dustlift.F, ustar_seuil=sqrt(stress/rho)

and alpha_lift[dust_mass] can be prescribed through an external stress.def parameter file.
--- alpha_lift[dust_number] is computed from alpha_lift[dust_mass] as in initracer.F
--- ustar_seuil is more user-friendly than stress because direct comparison with ustar from model

--> For the moment this is MESOSCALE only, but potentially useful to everyone

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/2dplot.py (modified) (3 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/2dplot.py

@@ -26 +26 @@
 #title = "Atmospheric temperature"
 
+name = "wrfout_d01_2024-03-22_01:00:00_z"
+itime = 12
+ndiv = 10
+zey = 120
+var = "VMR_ICE"
+title = "Volume mixing ratio of water ice [ppm]"
+vmin = -0.5
+vmax = 400.
+
 nc = Dataset(name)
 
 what_I_plot, error = reducefield( getfield(nc,var), d4=itime, d2=zey )
+
 
 y = nc.variables["vert"][:]
@@ -34 +44 @@
 horinp = len(what_I_plot[0,:])
 x = np.linspace(0.,horinp*500.,horinp) / 1000.
+xlabel("Horizontal coordinate (km)")
+
+horinp = len(what_I_plot[0,:])
+x = np.linspace(0.,horinp,horinp) 
+xlabel("x grid points")
+
 
 zevmin, zevmax = calculate_bounds(what_I_plot,vmin=vmin,vmax=vmax)
@@ -46 +62 @@
                        extend='neither',spacing='proportional')
 ptitle(title)
-xlabel("Horizontal coordinate (km)")
 ylabel("Altitude (m)")
 makeplotres(var+str(itime),res=200.,disp=False)