[ncl-talk] Error while handling multiple WRF output files
Barry Lynn
barry.h.lynn at gmail.com
Fri Nov 3 00:08:41 MDT 2017
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/),
>> pltres,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,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/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/94a06f7b/attachment.html>
More information about the ncl-talk
mailing list