[ncl-talk] Error while handling multiple WRF output files
swati shewale
shewaleswati16 at gmail.com
Fri Nov 3 00:10:17 MDT 2017
the provided program is quite long....going through it
On Fri, Nov 3, 2017 at 11:38 AM, Barry Lynn <barry.h.lynn at gmail.com> wrote:
> Hi:
>
> Did it work the way you want?
>
> On Fri, Nov 3, 2017 at 8:07 AM, swati shewale <shewaleswati16 at gmail.com>
> wrote:
>
>> Hello Barry,
>>
>> I have gone through the wrf function, you said (wrf_map_overlays).
>> Attached script as well as plot for the same. Actually I am prefering
>> combination of GSN and WRF function because while doing panel plot (3*3), I
>> need to use some of the gsn functions that helps to make self explanatory
>> plot (such as left/right/side strings, title, main title etc)
>>
>> On Fri, Nov 3, 2017 at 11:26 AM, Barry Lynn <barry.h.lynn at gmail.com>
>> wrote:
>>
>>> Swati:
>>>
>>> I assume you are aware that standard WRF programs can be used to make
>>> overlays, but you have chosen to use gsn instead.
>>> plot = wrf_map_overlays(a,wks,(/contour_tc,vector,contour_psl/),plt
>>> res,mpres)
>>>
>>> On Fri, Nov 3, 2017 at 7:22 AM, Barry Lynn <barry.h.lynn at gmail.com>
>>> wrote:
>>>
>>>> Hi:
>>>>
>>>> Just to check, from which data set are you trying to overlay? Is it
>>>> WRF or something else?
>>>>
>>>> The attached program works fine.
>>>>
>>>> Barry
>>>>
>>>> On Fri, Nov 3, 2017 at 7:05 AM, swati shewale <shewaleswati16 at gmail.com
>>>> > wrote:
>>>>
>>>>>
>>>>> I wish to overlay vector (u_plane, v_plane) over magnitude(spd i.e
>>>>> wind speed). so tried it another way as per your suggestion like
>>>>>
>>>>> u_plane!0 = "lat";
>>>>> u_plane!1 = "lon"
>>>>> u_new&lat= lat
>>>>> u_new&lon= lon
>>>>>
>>>>> v_plane!0 = "lat";
>>>>> v_plane!1 = "lon"
>>>>> v_new&lat= lat
>>>>> v_new&lon= lon
>>>>>
>>>>> spd_plane!0 = "lat";
>>>>> spd_plane!1 = "lon"
>>>>> spd_new&lat= lat
>>>>> spd_new&lon= lon
>>>>>
>>>>> still facing same issue
>>>>>
>>>>> On Fri, Nov 3, 2017 at 10:13 AM, swati shewale <
>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>
>>>>>> 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,pl
>>>>>>>>>>>>>>>>>>>>>>>>> ot9/),(/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,pl
>>>>>>>>>>>>>>>>>>>>>>>>> ot8/),(/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/Applic
>>>>>>>>>>>>>>>>>>>>>>>>>>> ations/panel.shtml#ex12
>>>>>>>>>>>>>>>>>>>>>>>>>>> http://www.ncl.ucar.edu/Applic
>>>>>>>>>>>>>>>>>>>>>>>>>>> ations/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/Applic
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ations/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 <(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 <(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 <(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/8627f053/attachment.html>
More information about the ncl-talk
mailing list