Re: overlay plot

Hi Oswald,
I don't see anything wrong with your code or the printVarSummary output.
I believe Dennis has written you offline asking you to copy your files
to our ftp site so we can run your code and examine the data. That would
be the next step.
Best regards,
Adam

Oswald Jason Lobo wrote:
> Hi,
>
> I'm trying to plot a ttest of pressure/height vs latitude for
> temperature. The problem i'm facing is that only half the plot gets
> hatched. Is there any way that i can get it to hatch over the entire
> plot. I have attached a picture of my plot and also my code and var
> summary. I also get an error which states
>
>
> warning:ScalarFieldSetValues: coordinate array sfXArray requires
> 149206043 elements:defaulting
> (0) get_lat_values: Warning: Your latitude values do not fall
> between -90 and 90 inclusive.
> (0) You will not get 'nice' latitude labels.
>
>
> My code
>
> load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>
> begin
>
>
>
> plot = new(3,graphic)
> header =
> (/"2xCO2-1xCO2~C~Temperature(K)","2xCO2_CCN-2xCO2~C~Temperature(K)","2xCO2_CCN-1xCO2~C~Temperature(K)"/)
>
> f1=addfile("ttest-ht-vs-lat-tCO2_atm800-CO2_400.nc","r")
> f2=addfile("ttest-ht-vs-lat-tgeo-erad115-CO2_atm800.nc","r")
> f3=addfile("ttest-ht-vs-lat-tgeo-erad115-CO2_400.nc","r")
>
> diff1=f1->diff
> t1 = f1->t
>
> diff2=f2->diff
> t2 = f2->t
>
> diff3=f3->diff
> t3 = f3->t
>
> printMinMax(diff3,True)
> printVarSummary(diff3)
> printMinMax(t3,True)
> printVarSummary(t3)
>
>
> wks = gsn_open_wks ("eps", "temp_ht_ttest" )
> gsn_define_colormap(wks,"BlueDarkRed18")
>
>
> res = True
> res_at_gsnDraw = False
> res_at_gsnFrame = False
> res_at_tmXBMode = "Explicit"
> res_at_tmXBValues=(/-90,-60,-30,0,30,60,90/)
> res_at_tmXBLabels=(/"SP","60S","30S","EQ","30N","60N","NP"/)
> res_at_tmYLMode = "Explicit"
> res_at_tmYLValues=(/1000,750,250,100,50,25,10/)
> res_at_tmYLLabels=(/1000,750,250,100,50,25,10/)
> res_at_tmXBLabelFontHeightF = 0.04
> res_at_tmXTOn = False
> res_at_tmYLOn = False
> res_at_tiYAxisString =""
> res_at_tiXAxisString =""
> res_at_cnInfoLabelOn=False
> res_at_tiYAxisFontHeightF = 0.05
> res_at_gsnLeftStringFontHeightF = 0.04
> res_at_cnLineLabelsOn = False
> res_at_cnFillOn = True
> res_at_lbLabelBarOn= False
> res_at_cnLevelSelectionMode = "ExplicitLevels"
> res_at_cnLevels = (/-8,-7,-6,-5,-4,-3,-2,-1,-.5,0.5,1,2,3,4,5/)
> res_at_cnFillColors =(/2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,18,19/)
> res_at_vpWidthF =0.7
> res_at_vpHeightF =0.7
>
>
>
> res2 = True
> res2_at_gsnDraw = False
> res2_at_gsnFrame = False
> res2_at_gsnAddCyclic = True
> res2_at_tiYAxisFontHeightF = 0.05
> res2_at_tmXBLabelsOn = False
> res2_at_cnFillOn = True
> res2_at_cnMonoFillColor= True
> res2_at_cnMonoFillPattern = True
> res2_at_cnFillPattern =3
> res2_at_cnLinesOn = False
> res2_at_cnLineLabelsOn = False
> res2_at_lbLabelBarOn = False
> res2_at_cnInfoLabelOn = False
> res2_at_cnMonoFillScale = True
> res2_at_cnFillScaleF=0.4
> res2_at_vpWidthF =0.7
> res2_at_vpHeightF =0.7
>
>
> res_at_gsnLeftString = header(0)
> plot(0)=gsn_csm_pres_hgt(wks,diff1,res)
> plott=gsn_csm_pres_hgt(wks,t1,res2)
> overlay(plot(0),plott)
> delete(res_at_gsnLeftString)
> delete(plott)
>
> res_at_gsnLeftString = header(1)
> plot(1)=gsn_csm_pres_hgt(wks,diff2,res)
> plott=gsn_csm_pres_hgt(wks,t2,res2)
> overlay(plot(1),plott)
> delete(res_at_gsnLeftString)
> delete(plott)
>
> res_at_gsnLeftString = header(2)
> plot(2)=gsn_csm_pres_hgt(wks,diff3,res)
> plott=gsn_csm_pres_hgt(wks,t3,res2)
> overlay(plot(2),plott)
> delete(res_at_gsnLeftString)
> delete(plott)
>
> res_pan = True
> res_pan_at_gsnFrame = False
> res_pan_at_gsnPanelRowSpec=True
> res_pan_at_gsnPanelXWhiteSpacePercent = 0
> res_pan_at_gsnPanelYWhiteSpacePercent=0
> res_pan_at_gsnPanelLabelBar = True
> res_pan_at_lbLabelFontHeightF = 0.02
> res_pan_at_lbLabelAutoStride = True
> res_pan_at_pmLabelBarWidthF =1
> res_pan_at_gsnPanelDebug = True
> res_pan_at_gsnPanelXF=(/0.05,0.5,0.3/)
> res_pan_at_gsnPanelYF=(/-1,-1,0.47/)
> gsn_panel(wks,plot,(/2,1/),res_pan)
>
>
>
>
> frame(wks)
> end
>
>
>
> Var summary
>
> (0)
> (0) min=-9.20668 max=0.762985
>
>
> Variable: diff3
> Type: float
> Total Size: 9984 bytes
> 2496 values
> Number of Dimensions: 2
> Dimensions and sizes: [lev | 26] x [lat | 96]
> Coordinates:
> lev: [3.54463800000001..992.5560999999998]
> lat: [ -90.. 90]
> Number Of Attributes: 2
> _FillValue : -999
> average_op_ncl : dim_avg over dimension(s): time
> (0)
> (0) min=-9.20668 max=0.762985
>
>
> Variable: t3
> Type: float
> Total Size: 9984 bytes
> 2496 values
> Number of Dimensions: 2
> Dimensions and sizes: [lev | 26] x [lat | 96]
> Coordinates:
> lev: [3.54463800000001..992.5560999999998]
> lat: [ -90.. 90]
> Number Of Attributes: 2
> average_op_ncl : dim_avg over dimension(s): time
> _FillValue : -999
>
> Thanks.
>
> ------------------------------------------------------------------------
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
--------------------------------------------------------------
Adam Phillips			             asphilli_at_ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
ESSL/CGD/CAS                               fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000	  http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk