Changeset 1336 for trunk/UTIL


Ignore:
Timestamp:
Aug 26, 2014, 6:06:26 PM (10 years ago)
Author:
aslmd
Message:

planetoplot. corrections to two examples scripts.

Location:
trunk/UTIL/PYTHON/planetoplot_v2/examples/ppclass_additional
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/UTIL/PYTHON/planetoplot_v2/examples/ppclass_additional/dynanalysis/eddyflux.py

    r1279 r1336  
    77#############################
    88ff = "diagfi.nc"
     9ff = "diagfi739139.nc"
    910tt = 240.
     11tt = 100000.
    1012#############################
    1113
    … …  
    1618
    1719# coordinates
    18 press = pp(var="p",file="diagfi.nc",t=tt,x=0,y=0).getf()
     20press = pp(var="p",file=ff,t=tt,x=0,y=0).getf()
    1921lat = np.linspace(-90.,90.,up.shape[1])
    2022
    … …  
    2931p.x = lat
    3032p.y = press/100.
    31 p.title = r"$\overline{u^\prime v^\prime}$"
     33#p.title = r'$\overline{u^\prime v^\prime}$'
     34p.title = ""
    3235p.xlabel = "Latitude"
    3336p.ylabel = "Pressure (mb)"
    34 p.units = r"$m^{2}s^{-2}$"
     37p.units = r'$m^{2}s^{-2}$'
    3538p.invert = True
    3639p.logy = True
    … …  
    4043p.vmax = 7
    4144p.div = 28
    42 p.makesave(mode="jpg",filename="eddyflux")
     45p.makesave(mode="png",filename="eddyflux")

  • trunk/UTIL/PYTHON/planetoplot_v2/examples/ppclass_additional/froude_jump.py

    r1192 r1336  
    88fi = 'testisaac/wrfout_d04_2024-03-46_06:00:00_zabg'
    99fi = 'testisaac/wrfout_d03_2024-03-46_06:00:00_zabg'
     10fi = 'wrfout_d05_2024-03-46_06:00:00_zabg'
     11fi = 'wrfout_d05_2024-03-57_06:00:00_zabg'
    1012ti = 1
    1113ti = 10
     14ti = 7
    1215zi = 25. # in meters
     16zi = 130. #ok vert velo
     17
    1318deltaT = 3.
     19deltaT = 4.
    1420grav = 3.72
    1521tpottest = 195. #test potential temperature for finding height of katabatic layer
     22#tpottest = 188.
     23tpottest = 200.
     24#tpottest = 205.
    1625
    1726### LOAD
    … …  
    4655p2d.title = "Topography"
    4756p2d.colorbar = "gist_rainbow"
    48 p2d.makeshow()
     57p2d.xlabel = "longitude"
     58p2d.ylabel = "latitude"
     59p2d.makesave(filename="topo")
    4960
    5061p2d = plot2d()
    … …  
    5566p2d.title = "Froude number"
    5667p2d.vmin = 0.
    57 p2d.vmax = 4.
    58 p2d.div = 40
    59 p2d.colorbar = "gist_stern"
     68p2d.vmax = 2. #4.
     69p2d.div = 20 #40
     70p2d.colorbar = "RdBu_r" #"gist_stern"
    6071p2d.fmt = "%.1f"
    61 p2d.makeshow()
     72p2d.xlabel = "longitude"
     73p2d.ylabel = "latitude"
     74p2d.makesave(filename="froude")
    6275
    6376p2d = plot2d()
    … …  
    7285p2d.colorbar = "spectral"
    7386p2d.fmt = "%.0f"
    74 p2d.makeshow()
     87p2d.xlabel = "longitude"
     88p2d.ylabel = "latitude"
     89p2d.makesave(filename="height")
    7590
    7691p2d = plot2d()
    … …  
    85100p2d.colorbar = "jet"
    86101p2d.fmt = "%.0f"
    87 p2d.makeshow()
     102p2d.xlabel = "longitude"
     103p2d.ylabel = "latitude"
     104p2d.makesave(filename="hwind")
    88105
    89106p2d = plot2d()
    … …  
    98115p2d.colorbar = "RdBu_r"
    99116p2d.fmt = "%.1f"
    100 p2d.makeshow()
     117p2d.xlabel = "longitude"
     118p2d.ylabel = "latitude"
     119p2d.makesave(filename="vwind")
    101120
Note: See TracChangeset for help on using the changeset viewer.