begin a = addfile("NCEP_06_climo.nc","r") ; Read in first time step and u = a->U(0,:,{60},:) ; lon=60 of atmospheric data. wks = gsn_open_wks("png", "preshgt1a") ; Open png file "preshgt1a.png" res = True ; No plot mods desired plot = gsn_csm_pres_hgt(wks,u,res) ; Call gsn_csm function for ; drawing a pressure/height plot. end