begin a = addfile("sst8292a.nc","r") b = addfile("olr7991a.nc","r") ; ; Read in sea surface temperature anomalies. ; Read in olr anomalies. ; Determine the indexes corresponding to Dec 1982. ; sst = a->SSTA olr = b->OLRA date_sst = a->date date_olr = b->date ind_sst = ind(date_sst.eq.198212) ind_olr = ind(date_olr.eq.198212) wks = gsn_open_wks("png","contour9a") ; Open "contour9a.png" res1 = True ; Res list for 1st data array res2 = True ; Res list for 2nd data array plot = gsn_csm_contour_map_overlay(wks,sst(ind_sst,:,:),\ olr(ind_olr,:,:),res1,res2) end