load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl" begin atilde = "a~H-13V2F35~D~FV-2H3~" ; ã iacute = "i~H-08V2F35~B~FV-2~" ; í eacute = "e~H-13V2F35~B~FV-2H3~" ; é cedil = "c~H-13F35~K~FH2~" ; ç aacute = "a~H-13V2F35~B~FV-2H3~" ; á ecirc = "e~H-13V2F35~C~FV-2H3~" ; ê eacute = "e~H-13V2F35~B~FV-2H3~" ; é oacute = "o~H-13V2F35~B~FV-2H3~" ; ó f1 = addfile("zonal.mcga.temp.z.1996.2005.nc","r") f2 = addfile("zonal.ncep.air.1986.2005.nc","r") f3 = addfile("temp.dif.nc","r") u1 = f1->temp(0,{1000:70},:,0) u2 = f1->temp(2,{1000:70},:,0) u3 = short2flt( f2->air(0,{1000:70},:,0) ) u4 = short2flt( f2->air(2,{1000:70},:,0) ) u5 = f3->temp(0,{1000:70},:,0) u6 = f3->temp(2,{1000:70},:,0) wks = gsn_open_wks ("eps", "teste") plot = new(6,graphic) gsn_define_colormap(wks,"gsdtol") ;shaded res = True res@gsnDraw = False res@gsnFrame = False res@tiMainString = "Temperatura do Ar (~S~o~N~C) - DJF" res@gsnLeftString = "MCGA" res@cnLevels = (/10,20,30/) res@cnFillColors = (/0,28,18,14/) res@gsnCenterString = "" res@gsnRightString = "" res@tiYAxisString = "Press"+atilde+"o (hPa)" res@tiXAxisString = "Latitude" res@cnLevelSelectionMode = "ExplicitLevels" res@cnFillDrawOrder = "PreDraw" res@cnFillOn = True res@cnLineLabelPlacementMode = "constant" res@cnLineLabelsOn = True ;contorno res1 = True res1@gsnDraw = False res1@gsnFrame = False res1@cnLevelSelectionMode = "ManualLevels" res1@gsnContourZeroLineThicknessF = 3. res1@gsnContourNegLineDashPattern = 1 res1@lbLabelStride = 1 res1@cnLineLabelsOn = True res1@gsnLeftString = " " res1@gsnCenterString = " " res1@gsnRightString = " " res1@tiYAxisString = " " res1@tiXAxisString = " " res1@cnLineLabelPlacementMode = "constant" res1@cnInfoLabelOn = False res1@cnLevelSpacingF = 10 res1@cnMinLevelValF = -70 ;res1@cnMaxLevelValF = 70 ; HERE IS WRONG!!! plot(0) = gsn_csm_pres_hgt(wks,u1,res) plot(0) = gsn_csm_pres_hgt(wks,u1,res1) plot(1) = gsn_csm_pres_hgt(wks,u2,res) plot(1) = gsn_csm_pres_hgt(wks,u2,res1) plot(2) = gsn_csm_pres_hgt(wks,u3,res) plot(2) = gsn_csm_pres_hgt(wks,u3,res1) plot(3) = gsn_csm_pres_hgt(wks,u4,res) plot(3) = gsn_csm_pres_hgt(wks,u4,res1) plot(4) = gsn_csm_pres_hgt(wks,u5,res) plot(4) = gsn_csm_pres_hgt(wks,u5,res1) plot(5) = gsn_csm_pres_hgt(wks,u6,res) plot(5) = gsn_csm_pres_hgt(wks,u6,res1) resP = True ; modify the panel plot resP@gsnFrame = False ; don't advance panel plot resP@gsnPanelBottom = 0.05 ; add space at bottom resP@gsnPanelFigureStrings = (/"a)","b)","c)","d)"/) ; add strings to panel resP@gsnPanelFigureStringsFontHeightF = 0.013 resP@amJust = "TopRight" resP@gsnPanelFigureStringsPerimOn = True gsn_panel(wks,plot,(/3,2/),resP) ; now draw as one plot frame(wks) end