[ncl-talk] Error while handling multiple WRF output files
swati shewale
shewaleswati16 at gmail.com
Thu Nov 2 22:35:37 MDT 2017
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,plot2,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,plot2,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,plot5,plot0,plot1,plot2/),(/2,3/
>>>>>>>>>>>>>>> ),pres)
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> On Mon, Oct 9, 2017 at 7:10 AM, swati shewale <
>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> above attached plot is just to show the scale (rightmost)
>>>>>>>>>>>>>>>> and left side titles that I want. I know the resource for left side is
>>>>>>>>>>>>>>>> tiYAxisString / tiXAxisString. but unable to make plotting from top to
>>>>>>>>>>>>>>>> bottom
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> On Mon, Oct 9, 2017 at 6:25 PM, swati shewale <
>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Dear Mary,
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> The resources that you set in your demo script like left
>>>>>>>>>>>>>>>>> and right space for panel (first panel of six plot vertically and another
>>>>>>>>>>>>>>>>> panel of three plot vertically) as below
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> ;---Define space for leftmost column of plots, and their
>>>>>>>>>>>>>>>>> labelbar width
>>>>>>>>>>>>>>>>> pres at gsnPanelLeft = 0.
>>>>>>>>>>>>>>>>> pres at gsnPanelRight = 2./3
>>>>>>>>>>>>>>>>> pres at pmLabelBarWidthF = 0.6
>>>>>>>>>>>>>>>>> gsn_panel(wks,(/plot1,plot2,pl
>>>>>>>>>>>>>>>>> ot3,plot4,plot5,plot6/),(/3,2/),pres)
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> ;---Define space for righmost column of plots, and their
>>>>>>>>>>>>>>>>> labelbar width
>>>>>>>>>>>>>>>>> delta = 0.028
>>>>>>>>>>>>>>>>> pres at gsnPanelLeft = pres at gsnPanelRight - delta
>>>>>>>>>>>>>>>>> ; Left plots don't stretch all the way to
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> ; 2/3, so move right plots a little to left.
>>>>>>>>>>>>>>>>> pres at gsnPanelRight = 1.0 - delta
>>>>>>>>>>>>>>>>> pres at pmLabelBarWidthF = pres at pmLabelBarWidthF * 0.5
>>>>>>>>>>>>>>>>> ; half the width of other labelbar
>>>>>>>>>>>>>>>>> gsn_panel(wks,(/plot7,plot8,plot9/),(/3,1/),pres)
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> If I wish to do plots from top to bottom instead of left
>>>>>>>>>>>>>>>>> to right, along with both scales vertical (which should have position to
>>>>>>>>>>>>>>>>> the extreme right).
>>>>>>>>>>>>>>>>> I just did change in your script as follow
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> ;---Set some panel resources common to both sets of plots
>>>>>>>>>>>>>>>>> pres = True
>>>>>>>>>>>>>>>>> pres at txFontHeightF = .010
>>>>>>>>>>>>>>>>> pres at txString = "
>>>>>>>>>>>>>>>>> WUCM
>>>>>>>>>>>>>>>>> NOUCM WUCM - NOUCM"
>>>>>>>>>>>>>>>>> pres at gsnPanelLabelBar = True
>>>>>>>>>>>>>>>>> pres at gsnMaximize = True
>>>>>>>>>>>>>>>>> pres at pmLabelBarHeightF = 0.08
>>>>>>>>>>>>>>>>> pres at lbLabelFontHeightF = 0.01
>>>>>>>>>>>>>>>>> pres at gsnFrame = False ; Needed for
>>>>>>>>>>>>>>>>> paneling two sets of plots on same frame.
>>>>>>>>>>>>>>>>> pres at gsnPanelYWhiteSpacePercent = 4
>>>>>>>>>>>>>>>>> pres at gsnPanelXWhiteSpacePercent = 4
>>>>>>>>>>>>>>>>> pres at gsnPanelBoxes = True ; debugging tool
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> ;---Define space for leftmost column of plots, and their
>>>>>>>>>>>>>>>>> labelbar width
>>>>>>>>>>>>>>>>> pres at gsnPanelTop = 0.
>>>>>>>>>>>>>>>>> pres at gsnPanelBottom = 2./3
>>>>>>>>>>>>>>>>> pres at pmLabelBarWidthF = 0.6
>>>>>>>>>>>>>>>>> gsn_panel(wks,(/plot0,plot1,pl
>>>>>>>>>>>>>>>>> ot2,plot3,plot4,plot5/),(/2,3/),pres)
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> ;---Define space for righmost column of plots, and their
>>>>>>>>>>>>>>>>> labelbar width
>>>>>>>>>>>>>>>>> delta = 0.005
>>>>>>>>>>>>>>>>> pres at gsnPanelTop = pres at gsnPanelBottom - delta ;
>>>>>>>>>>>>>>>>> Left plots don't stretch all the way to
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> ; 2/3, so move right plots a little to left.
>>>>>>>>>>>>>>>>> pres at gsnPanelBottom = 1.0 - delta
>>>>>>>>>>>>>>>>> pres at pmLabelBarWidthF = pres at pmLabelBarWidthF * 0.5
>>>>>>>>>>>>>>>>> ; half the width of other labelbar
>>>>>>>>>>>>>>>>> gsn_panel(wks,(/plot6,plot7,plot8/),(/3,1/),pres)
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> txres = True
>>>>>>>>>>>>>>>>> txres at txFontHeightF = 0.02
>>>>>>>>>>>>>>>>> gsn_text_ndc(wks,"Temperature at 2m (degC) over Patna
>>>>>>>>>>>>>>>>> city",0.5,0.92,txres)
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> gives me error as
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> (0) Warning: gsn_panel: attribute gsnPanelTop must be >
>>>>>>>>>>>>>>>>> 0.0 and <= 1.0
>>>>>>>>>>>>>>>>> (0) Defaulting to 1.
>>>>>>>>>>>>>>>>> (0) Warning: gsn_panel: font height (0.01) of main
>>>>>>>>>>>>>>>>> string is too large to fit in space provided. Either decrease font size or
>>>>>>>>>>>>>>>>> set gsnPanelTop.
>>>>>>>>>>>>>>>>> (0) Warning: gsn_panel: font height (0.02) of main
>>>>>>>>>>>>>>>>> string is too large to fit in space provided. Either decrease font size or
>>>>>>>>>>>>>>>>> set gsnPanelTop.
>>>>>>>>>>>>>>>>> (0) Error: gsn_panel: attribute gsnPanelTop (0.661667)
>>>>>>>>>>>>>>>>> must be greater
>>>>>>>>>>>>>>>>> (0) than gsnPanelBottom (0.995).
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Also Please tell me, how should I fix scale vertically
>>>>>>>>>>>>>>>>> that to extreme right. Actually I want plot from top to bottom because I
>>>>>>>>>>>>>>>>> want to provide resources at left side (for example purpose attached a
>>>>>>>>>>>>>>>>> plot).
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> I have model runs two types ....named WUCM and NOUCM
>>>>>>>>>>>>>>>>> I want it like......first row contain 3 plot (WUCM)
>>>>>>>>>>>>>>>>> second row contain 3 plot (NOWUCM)
>>>>>>>>>>>>>>>>> above two rows have same scale
>>>>>>>>>>>>>>>>> third row contain 3 plot of above
>>>>>>>>>>>>>>>>> difference (WUCM - NOWUCM)
>>>>>>>>>>>>>>>>> this one have different scale (as it is difference between
>>>>>>>>>>>>>>>>> WUCM - WUCM)
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> I hope, I tried better to convey things I want.
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> On Mon, Oct 9, 2017 at 11:30 AM, swati shewale <
>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> Thanks a lot Mary for giving detail info along with demo
>>>>>>>>>>>>>>>>>> script.
>>>>>>>>>>>>>>>>>> I have mentioned special lat2d/lon2d attributes in gsn
>>>>>>>>>>>>>>>>>> script (named convention in my a script is lat/lon). Anyway I will go
>>>>>>>>>>>>>>>>>> through the script and link. Thanks again.
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> On Sat, Oct 7, 2017 at 12:46 AM, Mary Haley <
>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Hello,
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Thanks for providing clean scripts and explanations.
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> 1. This warning is because when you use
>>>>>>>>>>>>>>>>>>> gsn_csm_contour_map to plot curvilinear data, you must associate the
>>>>>>>>>>>>>>>>>>> lat/lon arrays with the data being plotted. Otherwise, it does't know where
>>>>>>>>>>>>>>>>>>> on the map the data is supposed to go.
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> If you looked at the wrf_gsn_1.ncl and wrf_gsn_2.ncl
>>>>>>>>>>>>>>>>>>> scripts, you'll see they are both using special lat2d/lon2d attributes to
>>>>>>>>>>>>>>>>>>> accomplish this:
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> hgt = wrf_user_getvar(a,"HGT",it) ; Terrain
>>>>>>>>>>>>>>>>>>> elevation
>>>>>>>>>>>>>>>>>>> hgt at lat2d = wrf_user_getvar(a,"XLAT",it) ;
>>>>>>>>>>>>>>>>>>> latitude/longitude
>>>>>>>>>>>>>>>>>>> hgt at lon2d = wrf_user_getvar(a,"XLONG",it) ; required
>>>>>>>>>>>>>>>>>>> for plotting
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> However, you would have to do this for every variable
>>>>>>>>>>>>>>>>>>> you plot.
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Instead, you can set these two resources before any
>>>>>>>>>>>>>>>>>>> calls to gsn_csm_contour_map:
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> res at sfXArray = lon
>>>>>>>>>>>>>>>>>>> res at sfYArray = lat
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> 2.
>>>>>>>>>>>>>>>>>>> I think the issue with cnConstFEnableFill
>>>>>>>>>>>>>>>>>>> is that you are setting that resource *after* plots
>>>>>>>>>>>>>>>>>>> #0-5 have been created. You need to make sure this line:
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> res at cnConstFEnableFill= True ; allow constant
>>>>>>>>>>>>>>>>>>> values (1=sig) to fill
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> is *before* all of your gsn_csm_contour_map calls.
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> 3. In order to panel two sets of plots on the same
>>>>>>>>>>>>>>>>>>> page, you first need to set:
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> pnlres at gsnFrame = False
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> so that the frame doesn't get advanced after each one.
>>>>>>>>>>>>>>>>>>> You will then need to advance the frame yourself when done:
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> frame(wks)
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> The other issue with paneling two sets of plots, is that
>>>>>>>>>>>>>>>>>>> you need to make sure they don't get drawn on top of each other. There are
>>>>>>>>>>>>>>>>>>> many examples on our "Panels
>>>>>>>>>>>>>>>>>>> <http://www.ncl.ucar.edu/Applications/panel.shtml>"
>>>>>>>>>>>>>>>>>>> example page that show how to do various configurations of plots and
>>>>>>>>>>>>>>>>>>> labelbars. Here are some that might help:
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> http://www.ncl.ucar.edu/Applications/panel.shtml#ex12
>>>>>>>>>>>>>>>>>>> http://www.ncl.ucar.edu/Applications/panel.shtml#ex22
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> The important thing is that you need to set
>>>>>>>>>>>>>>>>>>> gsnPanelRight / gsnPanelLeft to control which part of the frame each set of
>>>>>>>>>>>>>>>>>>> plots is drawn in. Since the first set is 3 x 2 plots and the second set
>>>>>>>>>>>>>>>>>>> is 3 x 1 plots, this means the 3 x 2 plots should take up 2/3 of the width,
>>>>>>>>>>>>>>>>>>> and the 3 x 1 plots will take up 1/3 of the width.
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> See attached script which uses dummy data.
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> On Fri, Oct 6, 2017 at 2:04 AM, swati shewale <
>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Hello Mary,
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> As per your suggestion I went through the links
>>>>>>>>>>>>>>>>>>>> provided by you. I will explain one by one the issues I came across while
>>>>>>>>>>>>>>>>>>>> doing plot using inbuilt wrf function as well as gsn function separately. I
>>>>>>>>>>>>>>>>>>>> have attached scripts along with plot. the kind of titles (as well as left
>>>>>>>>>>>>>>>>>>>> and side string), I want is reflected in gsn plot which am unable to do in
>>>>>>>>>>>>>>>>>>>> wrf plot. I have given suitable names to scripts as well as plot for your
>>>>>>>>>>>>>>>>>>>> convenience. Hope you wont get confuse.
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> In first script (1_plot_using_wrf_function)
>>>>>>>>>>>>>>>>>>>> 1. after mentioning re at cnConstFEnableFill = True,
>>>>>>>>>>>>>>>>>>>> getting warning as
>>>>>>>>>>>>>>>>>>>> warning:ContourPlotInitialize: scalar field is
>>>>>>>>>>>>>>>>>>>> constant; no contour lines will appear; use cnConstFEnableFill to enable
>>>>>>>>>>>>>>>>>>>> fill
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> actually using wrf in built function gives nice plot
>>>>>>>>>>>>>>>>>>>> but unable to make self explanatory (giving strings and titles) plot so
>>>>>>>>>>>>>>>>>>>> thought of ploting using gsn function
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> In second script (2_plot_using_gsn_function), issues are
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> 1. even after mentioning lat/lon getting warning as
>>>>>>>>>>>>>>>>>>>> follow
>>>>>>>>>>>>>>>>>>>> (0) check_for_y_lat_coord: Warning: Data either does
>>>>>>>>>>>>>>>>>>>> not contain a valid latitude coordinate array or doesn't contain one at all.
>>>>>>>>>>>>>>>>>>>> (0) A valid latitude coordinate array should have a
>>>>>>>>>>>>>>>>>>>> 'units' attribute equal to one of the following values:
>>>>>>>>>>>>>>>>>>>> (0) 'degrees_north' 'degrees-north'
>>>>>>>>>>>>>>>>>>>> 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N'
>>>>>>>>>>>>>>>>>>>> 'degreeN' 'degreesN' 'deg north'
>>>>>>>>>>>>>>>>>>>> (0) check_for_lon_coord: Warning: Data either does
>>>>>>>>>>>>>>>>>>>> not contain a valid longitude coordinate array or doesn't contain one at
>>>>>>>>>>>>>>>>>>>> all.
>>>>>>>>>>>>>>>>>>>> (0) A valid longitude coordinate array should have a
>>>>>>>>>>>>>>>>>>>> 'units' attribute equal to one of the following values:
>>>>>>>>>>>>>>>>>>>> (0) 'degrees_east' 'degrees-east' 'degree_east'
>>>>>>>>>>>>>>>>>>>> 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE'
>>>>>>>>>>>>>>>>>>>> 'deg east'
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> 2. after mentioning re at cnConstFEnableFill = True,
>>>>>>>>>>>>>>>>>>>> getting warning as
>>>>>>>>>>>>>>>>>>>> warning:ContourPlotInitialize: scalar field is
>>>>>>>>>>>>>>>>>>>> constant; no contour lines will appear; use cnConstFEnableFill to enable
>>>>>>>>>>>>>>>>>>>> fill
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> 3. I want all 9 plots on single panel with scales
>>>>>>>>>>>>>>>>>>>> vertical, instead I am getting it on two different pages
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> I think, I have to manage or add more resources in
>>>>>>>>>>>>>>>>>>>> pres1/pres2 such as (*vpXF
>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpXF>*
>>>>>>>>>>>>>>>>>>>> and *vpYF
>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpYF>*)
>>>>>>>>>>>>>>>>>>>> (*vpHeightF
>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpHeightF>*
>>>>>>>>>>>>>>>>>>>> and *vpWidthF
>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpWidthF>*)
>>>>>>>>>>>>>>>>>>>> , so that I may get 9 plots on single panel. If you know any link which
>>>>>>>>>>>>>>>>>>>> resembles like the plot i want,Please let me know. side by I will again go
>>>>>>>>>>>>>>>>>>>> through links given by you.
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> On Sat, Sep 23, 2017 at 9:30 AM, swati shewale <
>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> Thank you so much Mary....surely I will go through the
>>>>>>>>>>>>>>>>>>>>> link provided by you....If come across any problem, will let you know
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> On Sat, Sep 23, 2017 at 12:36 AM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> Please see this page:
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> http://www.ncl.ucar.edu/Applications/wrfgsn.shtml
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> and look at the first example (wrf_gsn_1.ncl and
>>>>>>>>>>>>>>>>>>>>>> wrf_nogsn_1.ncl) and the second example (wrf_gsn_2.ncl).
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> These show how to use gsn to plot data and it
>>>>>>>>>>>>>>>>>>>>>> compares the plot with using WRF functions.
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> In the wrf_gsn_1.ncl plot, you can see that the
>>>>>>>>>>>>>>>>>>>>>> gsnLeft/RightString resources are being set automatically, because the data
>>>>>>>>>>>>>>>>>>>>>> has a description and units attribute.
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> Try this script as a start on your own data, and also
>>>>>>>>>>>>>>>>>>>>>> look at wrf_gsn_2.ncl to see how to further customize your plot.
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> If you are unable to get the plot to look as nice as
>>>>>>>>>>>>>>>>>>>>>> you want, email ncl-talk at ucar.edu again with your
>>>>>>>>>>>>>>>>>>>>>> new script and what you want it to look like, and we can help further.
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> Good luck,
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> On Fri, Sep 22, 2017 at 3:21 AM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> Thank you Mary... its wrf file, I just renamed it
>>>>>>>>>>>>>>>>>>>>>>> for my convenience. Actually I am able to plot it using wrf function but
>>>>>>>>>>>>>>>>>>>>>>> unable to put left/right/centre strings, as those are related to gsn
>>>>>>>>>>>>>>>>>>>>>>> functions. so to overcome the difficulty of left/right/centre, I thought to
>>>>>>>>>>>>>>>>>>>>>>> make a plot using gsn functions injstead of wrf functions. but inbuilt wrf
>>>>>>>>>>>>>>>>>>>>>>> functions gives nicer plots than gsn one.
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> On Fri, Sep 22, 2017 at 3:48 AM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> It's hard for somebody to simply look at a script
>>>>>>>>>>>>>>>>>>>>>>>> to see if it's correct or not. I do see some errors, but even if I correct
>>>>>>>>>>>>>>>>>>>>>>>> the errors, I can't be sure it's 100% correct without having access to your
>>>>>>>>>>>>>>>>>>>>>>>> data.
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> The problems are:
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> [1]
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> Your data file is called "patna_data.nc". You are
>>>>>>>>>>>>>>>>>>>>>>>> using wrf_xxxx functions to plot this data. Since WRF-ARW output files are
>>>>>>>>>>>>>>>>>>>>>>>> normally called something like "wrfout_d02_2005-08-29_00:00:00",
>>>>>>>>>>>>>>>>>>>>>>>> I can't be sure that you are actually plotting WRF-ARW data. You can only
>>>>>>>>>>>>>>>>>>>>>>>> use wrf_xxx plotting functions on WRF-ARW data.
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> [2]
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> You have:
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> t = a->T2(0:7,:,:)
>>>>>>>>>>>>>>>>>>>>>>>> ter = dim_avg_n( ter1, 0 )
>>>>>>>>>>>>>>>>>>>>>>>> ter = ter - 273
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> I don't know what your objective, so I don't know
>>>>>>>>>>>>>>>>>>>>>>>> if the above is correct. However, you are calling "dim_avg_n" on "ter1"
>>>>>>>>>>>>>>>>>>>>>>>> which doesn't exist. Did you mean:
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> ter = dim_avg_n( t, 0 )
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> The above will take an average of the first 8 times
>>>>>>>>>>>>>>>>>>>>>>>> teps of your data, assuming that the leftmost dimension of T2 represents
>>>>>>>>>>>>>>>>>>>>>>>> time.
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> Further, you should use the _Wrap version of this
>>>>>>>>>>>>>>>>>>>>>>>> function, so metadata is maintained:
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> ter = dim_avg_n_Wrap( t, 0 )
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> [3]
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> You have:
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> ter = ter - 273
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> Since you are changing the units, you should update
>>>>>>>>>>>>>>>>>>>>>>>> the "units" attribute to something like this:
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> ter at units = "degC"
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> The rest of the script looks okay, but again, the
>>>>>>>>>>>>>>>>>>>>>>>> wrf_contour and wrf_map_overlays procedures will only work if you are
>>>>>>>>>>>>>>>>>>>>>>>> actually plotting WRF-ARW data.
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> On Thu, Sep 7, 2017 at 10:02 PM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> Thank you. I will try with your script. yesterday
>>>>>>>>>>>>>>>>>>>>>>>>> I tried it with below script. Please have a look and let me know, whether
>>>>>>>>>>>>>>>>>>>>>>>>> it is correct or not.
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Sep 8, 2017 at 12:54 AM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> If I understand you correctly, then you need to
>>>>>>>>>>>>>>>>>>>>>>>>>> reshape your X array so that it is NTIM x 8 x NLAT x NLON, and then you can
>>>>>>>>>>>>>>>>>>>>>>>>>> average across the second dimension from the left.
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> I've attached a modified version of your script,
>>>>>>>>>>>>>>>>>>>>>>>>>> which is UNTESTED. You may need to modify it to do the right thing with
>>>>>>>>>>>>>>>>>>>>>>>>>> regard to copying over the required metadata. I did this a little with
>>>>>>>>>>>>>>>>>>>>>>>>>> copy_VarAtts, and copy_coords.
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> On Thu, Sep 7, 2017 at 3:42 AM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> Hiii,
>>>>>>>>>>>>>>>>>>>>>>>>>>> using ncrcat I have single file containing 24
>>>>>>>>>>>>>>>>>>>>>>>>>>> time step (3hrly, 3 days). Now I want to plot temperature variable (average
>>>>>>>>>>>>>>>>>>>>>>>>>>> for each day, spatial, panel plot). I am unable to take out time step for
>>>>>>>>>>>>>>>>>>>>>>>>>>> 1st day (0:7) and then do average of it...then second day (8:15)
>>>>>>>>>>>>>>>>>>>>>>>>>>> ...likewise for third day (16:23). I wrote a script, in that it takes
>>>>>>>>>>>>>>>>>>>>>>>>>>> average over all time steps. what changes should I make so that I can get
>>>>>>>>>>>>>>>>>>>>>>>>>>> average over specific time step? attached herewith script.
>>>>>>>>>>>>>>>>>>>>>>>>>>> summary of variable temperature is as follow.
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> Variable: X
>>>>>>>>>>>>>>>>>>>>>>>>>>> Type: float
>>>>>>>>>>>>>>>>>>>>>>>>>>> Total Size: 5658552 bytes
>>>>>>>>>>>>>>>>>>>>>>>>>>> 1414638 values
>>>>>>>>>>>>>>>>>>>>>>>>>>> Number of Dimensions: 3
>>>>>>>>>>>>>>>>>>>>>>>>>>> Dimensions and sizes: [Time | 23] x
>>>>>>>>>>>>>>>>>>>>>>>>>>> [south_north | 201] x [west_east | 306]
>>>>>>>>>>>>>>>>>>>>>>>>>>> Coordinates:
>>>>>>>>>>>>>>>>>>>>>>>>>>> Number Of Attributes: 6
>>>>>>>>>>>>>>>>>>>>>>>>>>> FieldType : 104
>>>>>>>>>>>>>>>>>>>>>>>>>>> MemoryOrder : XY
>>>>>>>>>>>>>>>>>>>>>>>>>>> description : TEMP at 2 M
>>>>>>>>>>>>>>>>>>>>>>>>>>> units : K
>>>>>>>>>>>>>>>>>>>>>>>>>>> stagger :
>>>>>>>>>>>>>>>>>>>>>>>>>>> coordinates : XLONG XLAT XTIME
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> On Thu, Sep 7, 2017 at 3:04 PM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>> Thank you. I will try with it.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Wed, Sep 6, 2017 at 8:39 PM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> The error message is telling you that there's
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> something wrong with argument 0 passed to wrf_contour. Argument 0 refers
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> to the first argument in the list, which is "a(itime)" in this case.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> The issue is that when you use addfiles to
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> open a series of files, you need to use square brackets, [], to access a
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> single file, and not curved brackets, (). This is because addfiles returns
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> a "List" object, and lists in NCL need to be access with square brackets.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> You have:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> contour = wrf_contour
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (a(itime),wks,t,t_res)
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> and this should be:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> contour = wrf_contour
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (a[itime],wks,t,t_res)
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Mon, Sep 4, 2017 at 4:40 AM, swati shewale
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Hiiii
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> I want to do panel plotting using WRF output.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> I have total 8 files, each contains single time step (one day, 3 hrly
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> datafiles so total 8 files). out of eight times I want to plot 6 time step
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (panel plot) for variable slp and temperature.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> It shows below error.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> fatal:Argument type mismatch on argument (0)
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> of (wrf_contour) can not coerce
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> fatal:["Execute.c":8567]:Execute: Error
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> occurred at or near line 44 in file temp.ncl
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (0) Error: gsn_panel: all of the plots
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> passed to gsn_panel appear to be invalid
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Attached herewith script.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Thank you in advance. Any help will be
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> appreciated
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ______________________________
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> _________________
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ncl-talk mailing list
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ncl-talk at ucar.edu
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> List instructions, subscriber options,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> unsubscribe:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> http://mailman.ucar.edu/mailma
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> n/listinfo/ncl-talk
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>
>>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
>
> --
> Barry H. Lynn, Ph.D
> Senior Lecturer,
> The Institute of the Earth Science,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> C.E.O, Weather It Is, LTD
> Weather and Climate Focus
> http://weather-it-is.com
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
> US 914 432 3108
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171103/1ca91c68/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wind_at_850_using_gsn_function.ncl
Type: application/octet-stream
Size: 3942 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171103/1ca91c68/attachment.obj>
More information about the ncl-talk
mailing list