Re: Simple line plot -need help

Hi Madeleine,

If you want multiple PDF files, you might try using the "newpdf" workstation type; you can have multiple workstations open simultaneously (provided they are named differently).

Rick

On Jan 27, 2011, at 12:19 PM, Madeleine Patterson wrote:

> Hi Again,
>
> I was able to graph the line plots I wanted, but even though I've had a look at various bits of the NCL site, I haven't been able to produce lots of different individual plots (both x11 and pdf) from the one script. I figured out that I can only have one pdf workstation at a time, but I was trying to get around it by naming them differently. (doesn't work, but would if I were trying to produce ps plots -why???)
>
> Is there any way round this - can I close a pdf workstation and open another one in the same script? I tried delete(wks) but my script still hangs after the first x11 plot pops up onto my screen (i.e):
>
> wks1a = gsn_open_wks("pdf","tsa-diff")
> wks1b= gsn_open_wks("x11","tsa-diff")
>
> res=True
> res@tiMainString = "monthly output TSA-difference over time"
>
> plot = gsn_csm_xy(wks1a,time,tsa_diff,res) <---- GETS STUCK HERE
> plot = gsn_csm_xy(wks1b,time,tsa_diff,res)
>
> delete(wks1a)
> delete(wks1b)
>
> ;;;
>
> wks2a = gsn_open_wks("pdf","fcev_diff")
> wks2b= gsn_open_wks("x11","fcev_diff")
>
> res=True
> res@tiMainString = "monthly output FCEV -difference over time"
>
> plot = gsn_csm_xy(wks2a,time,fcev_diff,res)
> plot = gsn_csm_xy(wks2b,time,fcev_diff,res)
>
> Any suggestions? Or will I have to be content to output ps files?
>
> Thanks
>
> Maddy
>
>
>
>
> On Tue, Jan 25, 2011 at 12:28 PM, Dennis Shea <shea@ucar.edu> wrote:
> Use wgt_areaave_Wrap rather than wgt_areaave.
> The _Wrap preserves the meta data
>
> http://www.ncl.ucar.edu/Document/Functions/Contributed/wgt_areaave_Wrap.shtml
>
> ---
>
> TsoiAve_clat = wgt_areaave(tsoi_all, clat, 1.0, 1)
> printVarSummary(TsoiAve_clat)
>
> you will see that there is no meta data
>
> Using
>
> TsoiAve_clat = wgt_areaave_Wrap(tsoi_all, clat, 1.0, 1)
> printVarSummary(TsoiAve_clat)
>
> ====
> The coordinates TsoiAve_clat will be (say) (time,depth)
>
> Do you wish the plot the global average at each level???
> or
> Do you wish the plot the global average of all levels
> where there is only one coordinate (time):
>
> If the latter ...
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/wgt_volave.shtml
>
> glAve = wgt_volave(tsoi_all, depth, clat, 1.0, 0) ; glAve(ktime)
> glAve!0 = "time"
> glAve&time = tsoi_all&time
> printVarSummary(glAve)
> =====
>
>
>
> Also, There is no need to explicitly read "time" to memory.
> time=fall[:]->time
>
> You can use the coordinate associated with the variable you are trying to plot eg ...
>
> plot = gsn_csm_xy(wks1,glAve&time,glAve,res)
>
>
> On 1/25/11 10:14 AM, Madeleine Patterson wrote:
> Hi,
>
> I am trying to plot weighted global area averaged soil temperature over
> time (10 layers of soil temps) over 30 years of monthly values.
>
> However I am getting a long list of TsoiAve_clat = 1e+36 at different
> coords when I run the script and a bunch of warnings and error messages
> at the end
>
> (0) gsn_csm_xy: Fatal: X and Y must have the same dimensions sizes,
> or one must be one-dimensional and both have the same rightmost dimension.
>
> warning:["CoordArrays.c":591]:No Valid values in Array, unable to
> compute Min or Max
> warning:LegendInitialize: lgItemTypes exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgDashIndexes exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgMarkerIndexes exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineLabelStrings exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineColors exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgMarkerColors exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineDashSegLens exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineThicknesses exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgMarkerThicknesses exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineLabelFontHeights exceeds maximum number
> of elements, 256: ignoring
> warning:LegendInitialize: lgMarkerSizes exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLabelStrings exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineLabelFontColors exceeds maximum number
> of elements, 256: ignoring
> warning:LegendInitialize: lgItemTypes exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgDashIndexes exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgMarkerIndexes exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineLabelStrings exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineColors exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgMarkerColors exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineDashSegLens exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineThicknesses exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgMarkerThicknesses exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineLabelFontHeights exceeds maximum number
> of elements, 256: ignoring
> warning:LegendInitialize: lgMarkerSizes exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLabelStrings exceeds maximum number of
> elements, 256: ignoring
> warning:LegendInitialize: lgLineLabelFontColors exceeds maximum number
> of elements, 256: ignoring
> warning:LegendSetValues: lgItemTypes exceeds maximum number of elements,
> 2::56: ignoring
> warning:LegendSetValues: lgDashIndexes exceeds maximum number of
> elements, 256: ignoring
> warning:LegendSetValues: lgMarkerIndexes exceeds maximum number of
> elements, 256: ignoring
> warning:LegendSetValues: lgLineLabelStrings exceeds maximum number of
> elements, 256: ignoring
> warning:LegendSetValues: lgLineColors exceeds maximum number of
> elements, 256: ignoring
> warning:LegendSetValues: lgMarkerColors exceeds maximum number of
> elements, 256: ignoring
> warning:LegendSetValues: lgLineDashSegLens exceeds maximum number of
> elements, 256: ignoring
> warning:LegendSetValues: lgLineThicknesses exceeds maximum number of
> elements, 256: ignoring
> warning:LegendSetValues: lgMarkerThicknesses exceeds maximum number of
> elements, 256: ignoring
> warning:LegendSetValues: lgLineLabelFontHeights exceeds maximum number
> of elements, 256: ignoring
> warning:LegendSetValues: lgMarkerSizes exceeds maximum number of
> elements, 256: ignoring
> warning:LegendSetValues: lgLabelStrings exceeds maximum number of
> elements, 256: ignoring
> warning:LegendSetValues: lgLineLabelFontColors exceeds maximum number of
> elements, 256: ignoring
> warning:LegendSetValues: lgItemTypes exceeds maximum number of elements,
> 256: ignoring
> warning:LegendSetValues: lgDashIndexes exceeds maximum number of
> elements, 256: ignoring
>
> My various print commands reveal that tsoi_all has 4 dims and TsoiAvg
> has 2, and I'm not really sure why this would be (or if its the cause of
> my problem).
>
> Can anyone point to where my script goes wrong?
>
> Thanks
>
> Maddy
>
> *************************************
>
> My script:
>
> begin
>
> 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"
>
> setfileoption("nc","SuppressClose",False)
>
> all_files = systemfunc("ls model-output2000.*")
>
> fall = addfiles(all_files,"r")
> ListSetType(fall,"cat")
> tsoi_all = fall[:]->TSOI
> printVarSummary(tsoi_all)
>
> lon = fall[1] ->lon
> lat = fall[1] ->lat
>
> rad = 4.0*atan(1.0)/180.0
>
> ; use the cosine of the latitudes to calc weighted global avg
>
> clat = cos(lat*rad)
>
> TsoiAve_clat = wgt_areaave(tsoi_all, clat, 1.0, 1)
>
> print(" TsoiAve_clat = " + TsoiAve_clat)
>
> TsoiAvg= TsoiAve_clat
>
> printVarSummary(TsoiAvg)
>
> time=fall[:]->time
>
> wks = gsn_open_wks("pdf","tsoi'")
> wks1= gsn_open_wks("x11","tsoi") ;open two types of workstations
>
> res=True
> res@tiMainString = "Soil Temperature (1970-2000)"
>
> plot = gsn_csm_xy(wks1,time,TsoiAvg,res)
>
> end
>
>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jan 27 13:38:32 2011