[ncl-talk] Error while handling multiple WRF output files
swati shewale
shewaleswati16 at gmail.com
Thu Nov 2 22:43:48 MDT 2017
Hello Barry,
Thanks for quick response. I am trying the way you suggested. still can you
please share a program.
On Fri, Nov 3, 2017 at 10:09 AM, Barry Lynn <barry.h.lynn at gmail.com> wrote:
> Hi Swati:
>
> Please do the following: i) make sure you have copied the variables to all
> secondary variables you wish to overlay, and ii) add those other lines of
> code as well. I believe the other lines are what will get rid of the error.
>
> If this doesn't work, I can send you a program that does work and perhaps
> you can adapt it to your needs.
>
> Barry
>
> On Fri, Nov 3, 2017 at 6:35 AM, swati shewale <shewaleswati16 at gmail.com>
> wrote:
>
>> Hello Barry,
>> Thank you for suggestion. I tried it using function copy_VarCoords still
>> it shows warning below
>>
>> (0) check_for_y_lat_coord: Warning: Data either does not contain a
>> valid latitude coordinate array or doesn't contain one at all.
>> (0) A valid latitude coordinate array should have a 'units' attribute
>> equal to one of the following values:
>> (0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees north'
>> 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
>> (0) check_for_lon_coord: Warning: Data either does not contain a valid
>> longitude coordinate array or doesn't contain one at all.
>> (0) A valid longitude coordinate array should have a 'units' attribute
>> equal to one of the following values:
>> (0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
>> 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
>>
>> I have attached modified script. Please go through it.
>>
>>
>> On Thu, Nov 2, 2017 at 6:19 PM, Barry Lynn <barry.h.lynn at gmail.com>
>> wrote:
>>
>>> Dear Swati:
>>>
>>> The error you are getting occurs when the variable you are trying to map
>>> has no coordinate information.
>>>
>>> There are a few ways to provide this information. The first is just to
>>> copy it from one variable (that has the coordinates) to another,
>>>
>>> For example:
>>>
>>> copy_VarCoords(u,u_ave) ; copy coord
>>> vars to speed
>>>
>>> copy_VarCoords(v,v_ave) ; copy coord
>>> vars to speed
>>>
>>> copy_VarCoords(z,z_ave) ; copy coord
>>> vars to height
>>>
>>> copy_VarCoords(t,t_ave) ; copy coord
>>> vars to t
>>>
>>> copy_VarCoords(rh,rh_ave) ; copy coord
>>> vars to rh
>>>
>>>
>>> another way is like this:
>>>
>>>
>>> u_new!0 = "lat"
>>>
>>> u_new!1 = "lon"
>>>
>>> u_new&lat= lat
>>>
>>> u_new&lon= lon
>>>
>>> I suggest to try either approach and see if this fixes your problem.
>>>
>>> Barry
>>>
>>>
>>>
>>> On Thu, Nov 2, 2017 at 2:10 PM, swati shewale <shewaleswati16 at gmail.com>
>>> wrote:
>>>
>>>> Dear Mary,
>>>>
>>>> I wish to plot winds at 850 hpa (both magnitude as well as vector). I
>>>> wrote the script using combination of gsn and wrf function (so that to
>>>> avoid problem while doing panel plot). I have mentioned lat/ lon as
>>>> sfxarray/sfyarray, still facing below warning
>>>>
>>>> (0) check_for_y_lat_coord: Warning: Data either does not contain a
>>>> valid latitude coordinate array or doesn't contain one at all.
>>>> (0) A valid latitude coordinate array should have a 'units'
>>>> attribute equal to one of the following values:
>>>> (0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees
>>>> north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg
>>>> north'
>>>> (0) check_for_lon_coord: Warning: Data either does not contain a
>>>> valid longitude coordinate array or doesn't contain one at all.
>>>> (0) A valid longitude coordinate array should have a 'units'
>>>> attribute equal to one of the following values:
>>>> (0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
>>>> 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
>>>>
>>>> Also unable to see colour shade in plot (scale is there).I have
>>>> attached script as well as plot. Please let me know necessary steps to do
>>>> to avoid the issue. Thanks in advance.
>>>>
>>>>
>>>> On Fri, Oct 13, 2017 at 9:41 AM, swati shewale <
>>>> shewaleswati16 at gmail.com> wrote:
>>>>
>>>>> Thank you Mary
>>>>>
>>>>> On Thu, Oct 12, 2017 at 9:33 PM, Mary Haley <haley at ucar.edu> wrote:
>>>>>
>>>>>> BTW, I have an error in the comments of this line:
>>>>>>
>>>>>> pnl_lft = 0.08 ; Moves plots to left a little, 0 is
>>>>>> default
>>>>>> It should be:
>>>>>>
>>>>>> pnl_lft = 0.08 ; Moves plots to RIGHT a little, 0 is
>>>>>> default
>>>>>> It's the pnl_rgt variable, which sets the right side limit for the
>>>>>> panels (default =1) which moves plots to the left, if set to < 1.
>>>>>>
>>>>>>
>>>>>> On Thu, Oct 12, 2017 at 8:38 AM, Mary Haley <haley at ucar.edu> wrote:
>>>>>>
>>>>>>> By removing the tickmarks and the labels for the middle plots, this
>>>>>>> has made them even smaller. This means that the scale factor now won't be
>>>>>>> as large for resizing the plots, because it thinks the plots are smaller.
>>>>>>> This, in turn, means the scale factor applied to the left plots is not
>>>>>>> large enough.
>>>>>>>
>>>>>>> I think you simply need to increase the gsnPanelLeft value a little
>>>>>>> to give those left plots a little more room. In the script I sent you, it's
>>>>>>> this value you want to increase:
>>>>>>>
>>>>>>> pnl_lft = 0.08 ; Moves plots to left a little, 0 is
>>>>>>> default
>>>>>>>
>>>>>>> I think it was set to 0.04 before. I made it 0.08 and those left
>>>>>>> titles started showing up again, after turning off the other tickmarks.
>>>>>>>
>>>>>>> --Mary
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> On Thu, Oct 12, 2017 at 4:51 AM, swati shewale <
>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>
>>>>>>>> Dear Mary,
>>>>>>>>
>>>>>>>> I used your dummy script. It helped a lot.
>>>>>>>> When I made few changes like.. off the common axis labels/tick
>>>>>>>> marks, it caused disappearance of Y axis string.
>>>>>>>> Another issue is.....displaying plot in PDF, whether landscape or
>>>>>>>> portrait mode gives same output. Please find the attachment. Can you please
>>>>>>>> suggest me about the same so that I could make plot (in PDF format) which
>>>>>>>> should look similarly as it looks in PNG format. Thank you.
>>>>>>>>
>>>>>>>> On Wed, Oct 11, 2017 at 9:39 AM, swati shewale <
>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>
>>>>>>>>> Thanks a lot Mary. I will surely go through it.
>>>>>>>>>
>>>>>>>>> On Tue, Oct 10, 2017 at 10:01 PM, Mary Haley <haley at ucar.edu>
>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>>> The issue is that your left plots have a title on the Y axis,
>>>>>>>>>> while your middle and right plots do not. gsn_panel assumes all plots are
>>>>>>>>>> the same size, and uses the first plot in the list to determine the scale
>>>>>>>>>> factor for all the rest. Since the left plots are different in size, the
>>>>>>>>>> scale factor applied to these plots ends up being different than the scale
>>>>>>>>>> factor applied to the right plots.
>>>>>>>>>>
>>>>>>>>>> The key here is to set the special gsnPanelScalePlotIndex to the
>>>>>>>>>> index of one of the middle plots, so it uses one of these for the scaling
>>>>>>>>>> factor. Your middle plots are indexes 1, 3, or 5, so:
>>>>>>>>>>
>>>>>>>>>> pres at gsnPanelScalePlotIndex = 1 ; base resizing on plot1
>>>>>>>>>> (middle plot)
>>>>>>>>>>
>>>>>>>>>> I've added dummy data to your gsn_function_try.ncl script, and
>>>>>>>>>> cleaned it up for you. I had to leave a little bit of a margin on the left
>>>>>>>>>> side of the panels, so the Y axis strings didn't get cut off.
>>>>>>>>>>
>>>>>>>>>> If you don't want all that space between the X axis of each plot,
>>>>>>>>>> you can set the gsnPanelTop / gsnPanelBottom resources as desired to squish
>>>>>>>>>> them together. I also set the gsnPanelYWhiteSpacePercent to 0 to force them
>>>>>>>>>> closer.
>>>>>>>>>>
>>>>>>>>>> See attached.
>>>>>>>>>>
>>>>>>>>>> -Mary
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> On Tue, Oct 10, 2017 at 4:49 AM, swati shewale <
>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>
>>>>>>>>>>> Dear Mary,
>>>>>>>>>>> I followed your dummy script and plotted my dataset. It worked,
>>>>>>>>>>> now facing the issue of size of the plot (third vertical column). Please
>>>>>>>>>>> find the attachment and tell me suggestions for the same. Thanks a lot for
>>>>>>>>>>> your dummy script.
>>>>>>>>>>>
>>>>>>>>>>> On Tue, Oct 10, 2017 at 12:45 PM, swati shewale <
>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>
>>>>>>>>>>>> I think I can change the position of title itself....it
>>>>>>>>>>>> resolves the issue....I will try to write the script that way
>>>>>>>>>>>>
>>>>>>>>>>>> On Tue, Oct 10, 2017 at 11:53 AM, swati shewale <
>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>
>>>>>>>>>>>>> vertical scale for WUCM and NOUCM is common (single)
>>>>>>>>>>>>>
>>>>>>>>>>>>> On Tue, Oct 10, 2017 at 11:51 AM, swati shewale <
>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>
>>>>>>>>>>>>>> I am sorry raw idea about plot is correct....but mistaken
>>>>>>>>>>>>>> about giving title....please neglect title
>>>>>>>>>>>>>> it should be as
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> title ---> 13aug 14aug
>>>>>>>>>>>>>> 15aug
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> WUCM plot0
>>>>>>>>>>>>>> plot1 plot2 vertical scale for WUCM and
>>>>>>>>>>>>>> NOUCM vertical scale for WUCM-NOUCM
>>>>>>>>>>>>>> NOUCM plot3
>>>>>>>>>>>>>> plot4 plot5
>>>>>>>>>>>>>> (WUCM - NOUCM) plot6 plot7
>>>>>>>>>>>>>> plot8
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> WUCM, NOUCM and (WUCM - NOUCM) is left axis string....title
>>>>>>>>>>>>>> and string I can provide....about scale doubtful as well as arranging plots
>>>>>>>>>>>>>> from top to bottom
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> On Tue, Oct 10, 2017 at 11:44 AM, swati shewale <
>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> I understood gsn_panel(wks,(/plot0,plot1,pl
>>>>>>>>>>>>>>> ot2,plot3,plot4,plot5/),(/2,3/),pres) will give 3 plot in
>>>>>>>>>>>>>>> one row and another 3 plot in second row as follows
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> plot0 plot1 plot2
>>>>>>>>>>>>>>> plot3 plot4 plot5
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> but below resources
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> pres at gsnPanelLeft = 0.
>>>>>>>>>>>>>>> pres at gsnPanelRight = 2./3
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> tells arrange plot as
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> plot0 plot1
>>>>>>>>>>>>>>> plot2 plot3
>>>>>>>>>>>>>>> plot4 plot4
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> so using which resource should I replace it
>>>>>>>>>>>>>>> Also tell me about scale, that I want in extreme right
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> let me give you raw idea
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> WUCM = 3 plot (1st row =plot0, plot1, plot2)
>>>>>>>>>>>>>>> NOWUCM = 3 plot (1st row =plot3, plot4, plot5)
>>>>>>>>>>>>>>> WUCM-NOUCM = 3 plot (1st row =plot6, plot7, plot8)
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> title ---> WUCM NOUCM (WUCM-NOUCM)
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> plot0 plot1
>>>>>>>>>>>>>>> plot2 vertical scale for WUCM and NOUCM
>>>>>>>>>>>>>>> vertical scale for WUCM-NOUCM
>>>>>>>>>>>>>>> plot3 plot4 plot5
>>>>>>>>>>>>>>> plot6 plot7 plot8
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> On Tue, Oct 10, 2017 at 11:29 AM, swati shewale <
>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> but as We have mentioned that put first 6 plot in left
>>>>>>>>>>>>>>>> side ( pres at gsnPanelLeft = 0.
>>>>>>>>>>>>>>>> pres at gsnPanelRight = 2./3)....so even if we write
>>>>>>>>>>>>>>>> gsn_panel(wks,(/plot0,plot1,pl
>>>>>>>>>>>>>>>> ot2,plot3,plot4,plot5/),(/2,3/),pres)
>>>>>>>>>>>>>>>> as per my knowledge, it will arrange as follows
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> plot0 plot1
>>>>>>>>>>>>>>>> plot2 plot3
>>>>>>>>>>>>>>>> plot4 plot5
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> On Mon, Oct 9, 2017 at 11:41 PM, Mary Haley <haley at ucar.edu
>>>>>>>>>>>>>>>> > wrote:
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> There's no resource that changes the ordering of the plots.
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> To change the ordering of the plots, you simply need to
>>>>>>>>>>>>>>>>> rearrange them in the call to gsn_panel.
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> In order words, which ever plot you want in the upper left
>>>>>>>>>>>>>>>>> corner is the one that should be first listed in the call to gsn_panel:
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> gsn_panel(wks,(/plot0,plot1,pl
>>>>>>>>>>>>>>>>> ot2,plot3,plot4,plot5/),(/2,3/),pres)
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Right now, it's ordered:
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> plot0 plot1 plot2
>>>>>>>>>>>>>>>>> plot3 plot4 plot5
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> If you want:
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> plot3 plot4 plot5
>>>>>>>>>>>>>>>>> plot0 plot1 plot2
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> then use:
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> gsn_panel(wks,(/plot3,plot4,pl
>>>>>>>>>>>>>>>>> ot5,plot0,plot1,plot2/),(/2,3/),pres)
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> On Mon, Oct 9, 2017 at 7:10 AM, swati shewale <
>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> above attached plot is just to show the scale (rightmost)
>>>>>>>>>>>>>>>>>> and left side titles that I want. I know the resource for left side is
>>>>>>>>>>>>>>>>>> tiYAxisString / tiXAxisString. but unable to make plotting from top to
>>>>>>>>>>>>>>>>>> bottom
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> On Mon, Oct 9, 2017 at 6:25 PM, swati shewale <
>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Dear Mary,
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> The resources that you set in your demo script like left
>>>>>>>>>>>>>>>>>>> and right space for panel (first panel of six plot vertically and another
>>>>>>>>>>>>>>>>>>> panel of three plot vertically) as below
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> ;---Define space for leftmost column of plots, and their
>>>>>>>>>>>>>>>>>>> labelbar width
>>>>>>>>>>>>>>>>>>> pres at gsnPanelLeft = 0.
>>>>>>>>>>>>>>>>>>> pres at gsnPanelRight = 2./3
>>>>>>>>>>>>>>>>>>> pres at pmLabelBarWidthF = 0.6
>>>>>>>>>>>>>>>>>>> gsn_panel(wks,(/plot1,plot2,pl
>>>>>>>>>>>>>>>>>>> ot3,plot4,plot5,plot6/),(/3,2/),pres)
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> ;---Define space for righmost column of plots, and their
>>>>>>>>>>>>>>>>>>> labelbar width
>>>>>>>>>>>>>>>>>>> delta = 0.028
>>>>>>>>>>>>>>>>>>> pres at gsnPanelLeft = pres at gsnPanelRight -
>>>>>>>>>>>>>>>>>>> delta ; Left plots don't stretch all the way to
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> ; 2/3, so move right plots a little to left.
>>>>>>>>>>>>>>>>>>> pres at gsnPanelRight = 1.0 - delta
>>>>>>>>>>>>>>>>>>> pres at pmLabelBarWidthF = pres at pmLabelBarWidthF *
>>>>>>>>>>>>>>>>>>> 0.5 ; half the width of other labelbar
>>>>>>>>>>>>>>>>>>> gsn_panel(wks,(/plot7,plot8,plot9/),(/3,1/),pres)
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> If I wish to do plots from top to bottom instead of left
>>>>>>>>>>>>>>>>>>> to right, along with both scales vertical (which should have position to
>>>>>>>>>>>>>>>>>>> the extreme right).
>>>>>>>>>>>>>>>>>>> I just did change in your script as follow
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> ;---Set some panel resources common to both sets of plots
>>>>>>>>>>>>>>>>>>> pres = True
>>>>>>>>>>>>>>>>>>> pres at txFontHeightF = .010
>>>>>>>>>>>>>>>>>>> pres at txString = "
>>>>>>>>>>>>>>>>>>> WUCM
>>>>>>>>>>>>>>>>>>> NOUCM WUCM - NOUCM"
>>>>>>>>>>>>>>>>>>> pres at gsnPanelLabelBar = True
>>>>>>>>>>>>>>>>>>> pres at gsnMaximize = True
>>>>>>>>>>>>>>>>>>> pres at pmLabelBarHeightF = 0.08
>>>>>>>>>>>>>>>>>>> pres at lbLabelFontHeightF = 0.01
>>>>>>>>>>>>>>>>>>> pres at gsnFrame = False ; Needed for
>>>>>>>>>>>>>>>>>>> paneling two sets of plots on same frame.
>>>>>>>>>>>>>>>>>>> pres at gsnPanelYWhiteSpacePercent = 4
>>>>>>>>>>>>>>>>>>> pres at gsnPanelXWhiteSpacePercent = 4
>>>>>>>>>>>>>>>>>>> pres at gsnPanelBoxes = True ; debugging tool
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> ;---Define space for leftmost column of plots, and their
>>>>>>>>>>>>>>>>>>> labelbar width
>>>>>>>>>>>>>>>>>>> pres at gsnPanelTop = 0.
>>>>>>>>>>>>>>>>>>> pres at gsnPanelBottom = 2./3
>>>>>>>>>>>>>>>>>>> pres at pmLabelBarWidthF = 0.6
>>>>>>>>>>>>>>>>>>> gsn_panel(wks,(/plot0,plot1,pl
>>>>>>>>>>>>>>>>>>> ot2,plot3,plot4,plot5/),(/2,3/),pres)
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> ;---Define space for righmost column of plots, and their
>>>>>>>>>>>>>>>>>>> labelbar width
>>>>>>>>>>>>>>>>>>> delta = 0.005
>>>>>>>>>>>>>>>>>>> pres at gsnPanelTop = pres at gsnPanelBottom - delta
>>>>>>>>>>>>>>>>>>> ; Left plots don't stretch all the way to
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> ; 2/3, so move right plots a little to left.
>>>>>>>>>>>>>>>>>>> pres at gsnPanelBottom = 1.0 - delta
>>>>>>>>>>>>>>>>>>> pres at pmLabelBarWidthF = pres at pmLabelBarWidthF *
>>>>>>>>>>>>>>>>>>> 0.5 ; half the width of other labelbar
>>>>>>>>>>>>>>>>>>> gsn_panel(wks,(/plot6,plot7,plot8/),(/3,1/),pres)
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> txres = True
>>>>>>>>>>>>>>>>>>> txres at txFontHeightF = 0.02
>>>>>>>>>>>>>>>>>>> gsn_text_ndc(wks,"Temperature at 2m (degC) over Patna
>>>>>>>>>>>>>>>>>>> city",0.5,0.92,txres)
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> gives me error as
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> (0) Warning: gsn_panel: attribute gsnPanelTop must be
>>>>>>>>>>>>>>>>>>> > 0.0 and <= 1.0
>>>>>>>>>>>>>>>>>>> (0) Defaulting to 1.
>>>>>>>>>>>>>>>>>>> (0) Warning: gsn_panel: font height (0.01) of main
>>>>>>>>>>>>>>>>>>> string is too large to fit in space provided. Either decrease font size or
>>>>>>>>>>>>>>>>>>> set gsnPanelTop.
>>>>>>>>>>>>>>>>>>> (0) Warning: gsn_panel: font height (0.02) of main
>>>>>>>>>>>>>>>>>>> string is too large to fit in space provided. Either decrease font size or
>>>>>>>>>>>>>>>>>>> set gsnPanelTop.
>>>>>>>>>>>>>>>>>>> (0) Error: gsn_panel: attribute gsnPanelTop
>>>>>>>>>>>>>>>>>>> (0.661667) must be greater
>>>>>>>>>>>>>>>>>>> (0) than gsnPanelBottom (0.995).
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Also Please tell me, how should I fix scale vertically
>>>>>>>>>>>>>>>>>>> that to extreme right. Actually I want plot from top to bottom because I
>>>>>>>>>>>>>>>>>>> want to provide resources at left side (for example purpose attached a
>>>>>>>>>>>>>>>>>>> plot).
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> I have model runs two types ....named WUCM and NOUCM
>>>>>>>>>>>>>>>>>>> I want it like......first row contain 3 plot (WUCM)
>>>>>>>>>>>>>>>>>>> second row contain 3 plot (NOWUCM)
>>>>>>>>>>>>>>>>>>> above two rows have same scale
>>>>>>>>>>>>>>>>>>> third row contain 3 plot of above
>>>>>>>>>>>>>>>>>>> difference (WUCM - NOWUCM)
>>>>>>>>>>>>>>>>>>> this one have different scale (as it is difference
>>>>>>>>>>>>>>>>>>> between WUCM - WUCM)
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> I hope, I tried better to convey things I want.
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> On Mon, Oct 9, 2017 at 11:30 AM, swati shewale <
>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Thanks a lot Mary for giving detail info along with
>>>>>>>>>>>>>>>>>>>> demo script.
>>>>>>>>>>>>>>>>>>>> I have mentioned special lat2d/lon2d attributes in gsn
>>>>>>>>>>>>>>>>>>>> script (named convention in my a script is lat/lon). Anyway I will go
>>>>>>>>>>>>>>>>>>>> through the script and link. Thanks again.
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> On Sat, Oct 7, 2017 at 12:46 AM, Mary Haley <
>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> Hello,
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> Thanks for providing clean scripts and explanations.
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> 1. This warning is because when you use
>>>>>>>>>>>>>>>>>>>>> gsn_csm_contour_map to plot curvilinear data, you must associate the
>>>>>>>>>>>>>>>>>>>>> lat/lon arrays with the data being plotted. Otherwise, it does't know where
>>>>>>>>>>>>>>>>>>>>> on the map the data is supposed to go.
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> If you looked at the wrf_gsn_1.ncl and wrf_gsn_2.ncl
>>>>>>>>>>>>>>>>>>>>> scripts, you'll see they are both using special lat2d/lon2d attributes to
>>>>>>>>>>>>>>>>>>>>> accomplish this:
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> hgt = wrf_user_getvar(a,"HGT",it) ; Terrain
>>>>>>>>>>>>>>>>>>>>> elevation
>>>>>>>>>>>>>>>>>>>>> hgt at lat2d = wrf_user_getvar(a,"XLAT",it) ;
>>>>>>>>>>>>>>>>>>>>> latitude/longitude
>>>>>>>>>>>>>>>>>>>>> hgt at lon2d = wrf_user_getvar(a,"XLONG",it) ;
>>>>>>>>>>>>>>>>>>>>> required for plotting
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> However, you would have to do this for every variable
>>>>>>>>>>>>>>>>>>>>> you plot.
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> Instead, you can set these two resources before any
>>>>>>>>>>>>>>>>>>>>> calls to gsn_csm_contour_map:
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> res at sfXArray = lon
>>>>>>>>>>>>>>>>>>>>> res at sfYArray = lat
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> 2.
>>>>>>>>>>>>>>>>>>>>> I think the issue with cnConstFEnableFill
>>>>>>>>>>>>>>>>>>>>> is that you are setting that resource *after* plots
>>>>>>>>>>>>>>>>>>>>> #0-5 have been created. You need to make sure this line:
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> res at cnConstFEnableFill= True ; allow constant
>>>>>>>>>>>>>>>>>>>>> values (1=sig) to fill
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> is *before* all of your gsn_csm_contour_map calls.
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> 3. In order to panel two sets of plots on the same
>>>>>>>>>>>>>>>>>>>>> page, you first need to set:
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> pnlres at gsnFrame = False
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> so that the frame doesn't get advanced after each one.
>>>>>>>>>>>>>>>>>>>>> You will then need to advance the frame yourself when done:
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> frame(wks)
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> The other issue with paneling two sets of plots, is
>>>>>>>>>>>>>>>>>>>>> that you need to make sure they don't get drawn on top of each other. There
>>>>>>>>>>>>>>>>>>>>> are many examples on our "Panels
>>>>>>>>>>>>>>>>>>>>> <http://www.ncl.ucar.edu/Applications/panel.shtml>"
>>>>>>>>>>>>>>>>>>>>> example page that show how to do various configurations of plots and
>>>>>>>>>>>>>>>>>>>>> labelbars. Here are some that might help:
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> http://www.ncl.ucar.edu/Applications/panel.shtml#ex12
>>>>>>>>>>>>>>>>>>>>> http://www.ncl.ucar.edu/Applications/panel.shtml#ex22
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> The important thing is that you need to set
>>>>>>>>>>>>>>>>>>>>> gsnPanelRight / gsnPanelLeft to control which part of the frame each set of
>>>>>>>>>>>>>>>>>>>>> plots is drawn in. Since the first set is 3 x 2 plots and the second set
>>>>>>>>>>>>>>>>>>>>> is 3 x 1 plots, this means the 3 x 2 plots should take up 2/3 of the width,
>>>>>>>>>>>>>>>>>>>>> and the 3 x 1 plots will take up 1/3 of the width.
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> See attached script which uses dummy data.
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> On Fri, Oct 6, 2017 at 2:04 AM, swati shewale <
>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> Hello Mary,
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> As per your suggestion I went through the links
>>>>>>>>>>>>>>>>>>>>>> provided by you. I will explain one by one the issues I came across while
>>>>>>>>>>>>>>>>>>>>>> doing plot using inbuilt wrf function as well as gsn function separately. I
>>>>>>>>>>>>>>>>>>>>>> have attached scripts along with plot. the kind of titles (as well as left
>>>>>>>>>>>>>>>>>>>>>> and side string), I want is reflected in gsn plot which am unable to do in
>>>>>>>>>>>>>>>>>>>>>> wrf plot. I have given suitable names to scripts as well as plot for your
>>>>>>>>>>>>>>>>>>>>>> convenience. Hope you wont get confuse.
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> In first script (1_plot_using_wrf_function)
>>>>>>>>>>>>>>>>>>>>>> 1. after mentioning re at cnConstFEnableFill = True,
>>>>>>>>>>>>>>>>>>>>>> getting warning as
>>>>>>>>>>>>>>>>>>>>>> warning:ContourPlotInitialize: scalar field is
>>>>>>>>>>>>>>>>>>>>>> constant; no contour lines will appear; use cnConstFEnableFill to enable
>>>>>>>>>>>>>>>>>>>>>> fill
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> actually using wrf in built function gives nice plot
>>>>>>>>>>>>>>>>>>>>>> but unable to make self explanatory (giving strings and titles) plot so
>>>>>>>>>>>>>>>>>>>>>> thought of ploting using gsn function
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> In second script (2_plot_using_gsn_function), issues
>>>>>>>>>>>>>>>>>>>>>> are
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> 1. even after mentioning lat/lon getting warning as
>>>>>>>>>>>>>>>>>>>>>> follow
>>>>>>>>>>>>>>>>>>>>>> (0) check_for_y_lat_coord: Warning: Data either
>>>>>>>>>>>>>>>>>>>>>> does not contain a valid latitude coordinate array or doesn't contain one
>>>>>>>>>>>>>>>>>>>>>> at all.
>>>>>>>>>>>>>>>>>>>>>> (0) A valid latitude coordinate array should have
>>>>>>>>>>>>>>>>>>>>>> a 'units' attribute equal to one of the following values:
>>>>>>>>>>>>>>>>>>>>>> (0) 'degrees_north' 'degrees-north'
>>>>>>>>>>>>>>>>>>>>>> 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N'
>>>>>>>>>>>>>>>>>>>>>> 'degreeN' 'degreesN' 'deg north'
>>>>>>>>>>>>>>>>>>>>>> (0) check_for_lon_coord: Warning: Data either does
>>>>>>>>>>>>>>>>>>>>>> not contain a valid longitude coordinate array or doesn't contain one at
>>>>>>>>>>>>>>>>>>>>>> all.
>>>>>>>>>>>>>>>>>>>>>> (0) A valid longitude coordinate array should have
>>>>>>>>>>>>>>>>>>>>>> a 'units' attribute equal to one of the following values:
>>>>>>>>>>>>>>>>>>>>>> (0) 'degrees_east' 'degrees-east'
>>>>>>>>>>>>>>>>>>>>>> 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E'
>>>>>>>>>>>>>>>>>>>>>> 'degreeE' 'degreesE' 'deg east'
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> 2. after mentioning re at cnConstFEnableFill = True,
>>>>>>>>>>>>>>>>>>>>>> getting warning as
>>>>>>>>>>>>>>>>>>>>>> warning:ContourPlotInitialize: scalar field is
>>>>>>>>>>>>>>>>>>>>>> constant; no contour lines will appear; use cnConstFEnableFill to enable
>>>>>>>>>>>>>>>>>>>>>> fill
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> 3. I want all 9 plots on single panel with scales
>>>>>>>>>>>>>>>>>>>>>> vertical, instead I am getting it on two different pages
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> I think, I have to manage or add more resources in
>>>>>>>>>>>>>>>>>>>>>> pres1/pres2 such as (*vpXF
>>>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpXF>*
>>>>>>>>>>>>>>>>>>>>>> and *vpYF
>>>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpYF>*)
>>>>>>>>>>>>>>>>>>>>>> (*vpHeightF
>>>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpHeightF>*
>>>>>>>>>>>>>>>>>>>>>> and *vpWidthF
>>>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpWidthF>*)
>>>>>>>>>>>>>>>>>>>>>> , so that I may get 9 plots on single panel. If you know any link which
>>>>>>>>>>>>>>>>>>>>>> resembles like the plot i want,Please let me know. side by I will again go
>>>>>>>>>>>>>>>>>>>>>> through links given by you.
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> On Sat, Sep 23, 2017 at 9:30 AM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> Thank you so much Mary....surely I will go through
>>>>>>>>>>>>>>>>>>>>>>> the link provided by you....If come across any problem, will let you know
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> On Sat, Sep 23, 2017 at 12:36 AM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> Please see this page:
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> http://www.ncl.ucar.edu/Applications/wrfgsn.shtml
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> and look at the first example (wrf_gsn_1.ncl and
>>>>>>>>>>>>>>>>>>>>>>>> wrf_nogsn_1.ncl) and the second example (wrf_gsn_2.ncl).
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> These show how to use gsn to plot data and it
>>>>>>>>>>>>>>>>>>>>>>>> compares the plot with using WRF functions.
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> In the wrf_gsn_1.ncl plot, you can see that the
>>>>>>>>>>>>>>>>>>>>>>>> gsnLeft/RightString resources are being set automatically, because the data
>>>>>>>>>>>>>>>>>>>>>>>> has a description and units attribute.
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> Try this script as a start on your own data, and
>>>>>>>>>>>>>>>>>>>>>>>> also look at wrf_gsn_2.ncl to see how to further customize your plot.
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> If you are unable to get the plot to look as nice
>>>>>>>>>>>>>>>>>>>>>>>> as you want, email ncl-talk at ucar.edu again with
>>>>>>>>>>>>>>>>>>>>>>>> your new script and what you want it to look like, and we can help further.
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> Good luck,
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Sep 22, 2017 at 3:21 AM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> Thank you Mary... its wrf file, I just renamed it
>>>>>>>>>>>>>>>>>>>>>>>>> for my convenience. Actually I am able to plot it using wrf function but
>>>>>>>>>>>>>>>>>>>>>>>>> unable to put left/right/centre strings, as those are related to gsn
>>>>>>>>>>>>>>>>>>>>>>>>> functions. so to overcome the difficulty of left/right/centre, I thought to
>>>>>>>>>>>>>>>>>>>>>>>>> make a plot using gsn functions injstead of wrf functions. but inbuilt wrf
>>>>>>>>>>>>>>>>>>>>>>>>> functions gives nicer plots than gsn one.
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Sep 22, 2017 at 3:48 AM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> It's hard for somebody to simply look at a script
>>>>>>>>>>>>>>>>>>>>>>>>>> to see if it's correct or not. I do see some errors, but even if I correct
>>>>>>>>>>>>>>>>>>>>>>>>>> the errors, I can't be sure it's 100% correct without having access to your
>>>>>>>>>>>>>>>>>>>>>>>>>> data.
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> The problems are:
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> [1]
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> Your data file is called "patna_data.nc". You
>>>>>>>>>>>>>>>>>>>>>>>>>> are using wrf_xxxx functions to plot this data. Since WRF-ARW output files
>>>>>>>>>>>>>>>>>>>>>>>>>> are normally called something like "wrfout_d02_2005-08-29_00:00:00",
>>>>>>>>>>>>>>>>>>>>>>>>>> I can't be sure that you are actually plotting WRF-ARW data. You can only
>>>>>>>>>>>>>>>>>>>>>>>>>> use wrf_xxx plotting functions on WRF-ARW data.
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> [2]
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> You have:
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> t = a->T2(0:7,:,:)
>>>>>>>>>>>>>>>>>>>>>>>>>> ter = dim_avg_n( ter1, 0 )
>>>>>>>>>>>>>>>>>>>>>>>>>> ter = ter - 273
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> I don't know what your objective, so I don't know
>>>>>>>>>>>>>>>>>>>>>>>>>> if the above is correct. However, you are calling "dim_avg_n" on "ter1"
>>>>>>>>>>>>>>>>>>>>>>>>>> which doesn't exist. Did you mean:
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> ter = dim_avg_n( t, 0 )
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> The above will take an average of the first 8
>>>>>>>>>>>>>>>>>>>>>>>>>> times teps of your data, assuming that the leftmost dimension of T2
>>>>>>>>>>>>>>>>>>>>>>>>>> represents time.
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> Further, you should use the _Wrap version of this
>>>>>>>>>>>>>>>>>>>>>>>>>> function, so metadata is maintained:
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> ter = dim_avg_n_Wrap( t, 0 )
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> [3]
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> You have:
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> ter = ter - 273
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> Since you are changing the units, you should
>>>>>>>>>>>>>>>>>>>>>>>>>> update the "units" attribute to something like this:
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> ter at units = "degC"
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> The rest of the script looks okay, but again, the
>>>>>>>>>>>>>>>>>>>>>>>>>> wrf_contour and wrf_map_overlays procedures will only work if you are
>>>>>>>>>>>>>>>>>>>>>>>>>> actually plotting WRF-ARW data.
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> On Thu, Sep 7, 2017 at 10:02 PM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> Thank you. I will try with your script.
>>>>>>>>>>>>>>>>>>>>>>>>>>> yesterday I tried it with below script. Please have a look and let me know,
>>>>>>>>>>>>>>>>>>>>>>>>>>> whether it is correct or not.
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Sep 8, 2017 at 12:54 AM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>> If I understand you correctly, then you need to
>>>>>>>>>>>>>>>>>>>>>>>>>>>> reshape your X array so that it is NTIM x 8 x NLAT x NLON, and then you can
>>>>>>>>>>>>>>>>>>>>>>>>>>>> average across the second dimension from the left.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>> I've attached a modified version of your
>>>>>>>>>>>>>>>>>>>>>>>>>>>> script, which is UNTESTED. You may need to modify it to do the right thing
>>>>>>>>>>>>>>>>>>>>>>>>>>>> with regard to copying over the required metadata. I did this a little
>>>>>>>>>>>>>>>>>>>>>>>>>>>> with copy_VarAtts, and copy_coords.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Thu, Sep 7, 2017 at 3:42 AM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Hiii,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> using ncrcat I have single file containing 24
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> time step (3hrly, 3 days). Now I want to plot temperature variable (average
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> for each day, spatial, panel plot). I am unable to take out time step for
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 1st day (0:7) and then do average of it...then second day (8:15)
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ...likewise for third day (16:23). I wrote a script, in that it takes
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> average over all time steps. what changes should I make so that I can get
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> average over specific time step? attached herewith script.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> summary of variable temperature is as follow.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Variable: X
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Type: float
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Total Size: 5658552 bytes
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 1414638 values
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Number of Dimensions: 3
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Dimensions and sizes: [Time | 23] x
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [south_north | 201] x [west_east | 306]
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Coordinates:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Number Of Attributes: 6
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> FieldType : 104
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> MemoryOrder : XY
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> description : TEMP at 2 M
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> units : K
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> stagger :
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> coordinates : XLONG XLAT XTIME
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Thu, Sep 7, 2017 at 3:04 PM, swati shewale
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Thank you. I will try with it.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Wed, Sep 6, 2017 at 8:39 PM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> The error message is telling you that
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> there's something wrong with argument 0 passed to wrf_contour. Argument 0
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> refers to the first argument in the list, which is "a(itime)" in this case.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> The issue is that when you use addfiles to
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> open a series of files, you need to use square brackets, [], to access a
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> single file, and not curved brackets, (). This is because addfiles returns
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> a "List" object, and lists in NCL need to be access with square brackets.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> You have:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> contour = wrf_contour
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (a(itime),wks,t,t_res)
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> and this should be:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> contour = wrf_contour
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (a[itime],wks,t,t_res)
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Mon, Sep 4, 2017 at 4:40 AM, swati
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> shewale <shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Hiiii
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> I want to do panel plotting using WRF
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> output. I have total 8 files, each contains single time step (one day, 3
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> hrly datafiles so total 8 files). out of eight times I want to plot 6 time
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> step (panel plot) for variable slp and temperature.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> It shows below error.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> fatal:Argument type mismatch on argument
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (0) of (wrf_contour) can not coerce
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> fatal:["Execute.c":8567]:Execute: Error
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> occurred at or near line 44 in file temp.ncl
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (0) Error: gsn_panel: all of the plots
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> passed to gsn_panel appear to be invalid
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Attached herewith script.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Thank you in advance. Any help will be
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> appreciated
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ______________________________
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> _________________
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ncl-talk mailing list
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ncl-talk at ucar.edu
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> List instructions, subscriber options,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> unsubscribe:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> http://mailman.ucar.edu/mailma
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> n/listinfo/ncl-talk
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>
>>>
>>> --
>>> Barry H. Lynn, Ph.D
>>> Senior Lecturer,
>>> The Institute of the Earth Science,
>>> The Hebrew University of Jerusalem,
>>> Givat Ram, Jerusalem 91904, Israel
>>> Tel: 972 547 231 170
>>> Fax: (972)-25662581
>>>
>>> C.E.O, Weather It Is, LTD
>>> Weather and Climate Focus
>>> http://weather-it-is.com
>>> Jerusalem, Israel
>>> Local: 02 930 9525
>>> Cell: 054 7 231 170
>>> Int-IS: x972 2 930 9525
>>> US 914 432 3108 <(914)%20432-3108>
>>>
>>
>>
>
>
> --
> Barry H. Lynn, Ph.D
> Senior Lecturer,
> The Institute of the Earth Science,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> C.E.O, Weather It Is, LTD
> Weather and Climate Focus
> http://weather-it-is.com
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
> US 914 432 3108
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171103/7a038eb3/attachment-0001.html>
More information about the ncl-talk
mailing list