Hi Mike,
It's been my experience that whenever something like this happens,
something is wrong with the data that is being passed in to the plotting
function. I downloaded your file, and put the data thru a statistical
routine of mine, which showed that the min of the data = .1883 but the
max = infinity. You can see where the problem is if you do a print(t).
The first and last latitude points are all = inf...
Adam
brownmc_at_uci.edu wrote:
> Hi all,
> I am trying to plot a file I made using nco and when I run the script all
> the syntax is correct, but ncl sits and thinks forever and never plots the
> file. I am not sure what is wrong, I can view the file and its variables
> with ncview so I am pretty sure they exsist. I have posted the file
> online at dust.ess.uci.edu/brownmc/regrid/2x2_sresa1b_rms_all_models.nc
> and my script is below. Thanks for the help.
>
> Mike Brown
>
> 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/shea_util.ncl"
> begin
> a =
> addfile("/data/brownmc/sresa1b/regridded/2x2_sresa1b_rms_all_models.nc","r")
> t = a->tas(0,0,:,:)
> wks = gsn_open_wks("X11","rms_test")
> gsn_define_colormap(wks,"gui_default")
> res = True
> res@mpProjection = "Mollweide"
> res@mpGridAndLimbOn = True
> res@mpPerimOn = False
> res@mpFillOn = False
> res@cnFillOn = True
> res@cnLineLabelsOn = False
> res@txFontHeightF = 0.015
> res@vpXF = 0.1
> res@vpYF = 0.9
> res@vpWidthF = 0.8
> res@vpHeightF = 0.8
> res@gsnSpreadColors = True
> res@gsnSpreadColorStart = 2
> res@gsnSpreadColorEnd = 21
> res@lbLabelFontHeightF = 0.015
> res@tiMainString = "RMS test plot"
> res@tiMainFontHeightF = .018
> plot = gsn_csm_contour_map(wks,t,res)
> end
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> 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 ncl-talk_at_ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Wed Apr 26 2006 - 14:57:50 MDT
This archive was generated by hypermail 2.2.0 : Wed Apr 26 2006 - 17:00:19 MDT