Last modified:

Thermospheric General Circulation Models (TGCM)

zp/Height vs. Time

Master-Script Examples || Tutorial Examples



master script: zput_master.ncl

This script contours selected fields over a polar stereographic map projection at selected Zp pressure levels. User settable options include color, curly-vector overlay, and postscript file output.
1) user choice: black and white + lat/lon ut(hours)

     
     comment out colormap
     ; colormap = "BlWhRe"
     type       = "lon" 
     lon2plt    =  90 
     lat2plt    = -55.
     
2) user choice: color + lat/slt ut(hours)

      
     colormap = "BlWhRe" 
     type     = "slt" 
     slt      =  5 
     lat2plt  = -55.
     
     
3) user choice: color + lat/zonal mean ut(days)

     
     colormap = "BlWhRe" 
     type     = "zmean" 
     lat2plt  = -55.
     
4) user choice: color + lat/zonal mean ut(months)

     
     colormap = "BlAqGrYeOrRevi200" 
     type     = "zmean" 
     lat2plt  = 40.
     
     
5) user choice: color + lat/zonal mean + limit axis ut(months)

     
     colormap = "BlWhRe" 
     type     = "slt" 
     slt      =  5
     lat2plt  = -55.
     ymin     = -6 
     
6) user choice: color + lat/zonal mean + limit axis + manually select contour levels ut(months)

     
     colormap = "BlWhRe"
     type     = "slt" 
     slt      =  5 
     lat2plt  = -55.
     ymin     = -6 
     manlev   =  True
     minlev   = -160
     maxlev   =  160
     int      =  20
     
7) user choice: color + lat/lon + interpolate to height

     
     colormap = "BlAqGrYeOrRevi200"
     type     = "lon" 
     lon2plt  =  -112
     lat2plt  = -55.
     height   = True
     top      = 500
     bott     = 30
     dz       = 40 
     
8) user choice: color + lat/slt + density unit conversion

     
     colormap      = "BlAqGrYeOrRevi200"
     type          = "slt"
     slt           = 5	 
     lat2plt       = -55.
     ymin          = -9
     ymax          = -5
     density_units = "CM3"