Ferret script for checking bogus TC
> cancel data/all
> cancel view
> use $1 ### read in argument
> set data $1
> cancel mode logo
> set viewport ul
> show data
> let/title=”Surface pressure (hPa)” slp = psfc/100
> fill/level=(800)(920,1030,5)(1050) slp
> contour/overlay landmask
To plot the WPS/WRF domain setup
> cancel data/all
> cancel view
> cancel mode logo
> use geo_em.d02.nc ### read in argument
> list XLAT_M[X=1] or > list/file=lat.dat XLAT_M[X=1] ! store the output into the “lat.dat” file
> list XLONG_M[Y=1] or > list/file=lon.dat XLONG_M[Y=1]
> use geo_em.d03.nc ### read in argument
> list XLAT_M[X=1]
> list XLONG_M[Y=1]
> use geo_em.d01.nc ### read in argument
> list XLAT_M[X=1]
> list XLONG_M[Y=1]
> use etopo5
> fill ROSE[X=40:100,Y=0:40]
> go land_detail thick " " black black ! thick pen, continents, countries
> go box 21.31,64.51,11.03,26.85 ### box for domain 1
> go box xx.xx,xx.xx,xx.xx,xx.xx ### box for domain 2
> go box xx.xx,xx.xx,xx.xx,xx.xx ### box for domain 3
> contour/overlay
> frame/file=output.gif ! Saves the current graphics display image as a graphics file.
How to draw high resolution coastline
use the land_detail journal as follow
> go land_detail thick " " black black ! thick pen, continents, countries