begin a = addfile("atmos2.nc","r") ; Read in first time step t = a->TS(0,:,:) ; of atmospheric data. wks = gsn_open_wks("png", "contour3c") ; Open png file "contour3c.png" res = True ; Plot mods desired res@gsnPolar = "NH" ; Specify northern hemisphere ; (which is the default). res@cnFillOn = True ; Turn on contour fill res@cnFillPalette = "wgne15" res@mpMinLatF = 60 ; Default is 0. ; ; This is ice data which has a gap in the tropics. So, rather than ; contouring the whole field, when we are zooming in on the area ; from latitude=60 to latitude=90, we can use coordinate subscripting to ; only pass in the data from latitude=60 to latitude=90. ; plot = gsn_csm_contour_map_polar(wks,t({60:90},:),res) end