Skip to content

Instantly share code, notes, and snippets.

@brews
Created October 12, 2014 05:13
Show Gist options
  • Save brews/e846e6f5dd47cd576b3a to your computer and use it in GitHub Desktop.
Save brews/e846e6f5dd47cd576b3a to your computer and use it in GitHub Desktop.
GEMPAK vorticity advection plotting example
#! /bin/csh -f
# 2014-10-11
# S. B. Malevich <[email protected]
# Example of vorticity advection plotting script for weather analysis:
# You might want to change the fill plotting colors.
set IN_FILE = 'waf_assignment2_data.gem'
set OUT_FILE = 'totally_super_awesome_vadvec_map.gif'
set MAP_AREA = 'us'
rm $OUT_FILE
gpcolor << EOF
COLORS = 101=255:255:255
DEVICE = gif | $OUT_FILE | 1280 ; 960
run
exit
EOF
gdcntr << EOF
GDATTIM = f00
GLEVEL = 500
GVCORD = pres
GFUNC = hght
GDFILE = $IN_FILE
CINT = 60
LINE = 32//2
FLINE = 14-28
MAP = 32
TITLE = 32
DEVICE = gif | $OUT_FILE | 1280 ; 960
PROJ = STR/90;-100;0
GAREA = $MAP_AREA
IJSKIP = 1
CLEAR = 1
PANEL = 0
TEXT = 1/2//sw
CLRBAR = 32////0.75
CONTUR = 0
SKIP = 0
HILO = 32;32/X;N
HLSYM = 2///2
run
SCALE = 0
GFUNC = vor(wnd)
FINT = 2e-5
CTYPE = f
run
SCALE = 0
GFUNC = hght
CINT = 60
CTYPE = c
run
SCALE = 0
GFUNC = adv(avor(wnd),wnd)
LINE = 32/3
CTYPE = c
run
GFUNC = LOWS(VOR(wnd))
run
GFUNC = HIGH(VOR(wnd))
run
exit
EOF
gpend
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment