[ncl-talk] Error while handling multiple WRF output files

swati shewale shewaleswati16 at gmail.com
Thu Nov 2 23:05:17 MDT 2017


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/Applications/panel.shtml#ex12
>>>>>>>>>>>>>>>>>>>>>> http://www.ncl.ucar.edu/Applications/panel.shtml#ex22
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> The important thing is that you need to set
>>>>>>>>>>>>>>>>>>>>>> gsnPanelRight / gsnPanelLeft to control which part of the frame each set of
>>>>>>>>>>>>>>>>>>>>>> plots is drawn in.  Since the first set is 3 x 2 plots and the second set
>>>>>>>>>>>>>>>>>>>>>> is 3 x 1 plots, this means the 3 x 2 plots should take up 2/3 of the width,
>>>>>>>>>>>>>>>>>>>>>> and the 3 x 1 plots will take up 1/3 of the width.
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> See attached script which uses dummy data.
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> On Fri, Oct 6, 2017 at 2:04 AM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> Hello Mary,
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> As per your suggestion I went through the links
>>>>>>>>>>>>>>>>>>>>>>> provided by you. I will explain one by one the issues I came across while
>>>>>>>>>>>>>>>>>>>>>>> doing plot using inbuilt wrf function as well as gsn function separately. I
>>>>>>>>>>>>>>>>>>>>>>> have attached scripts along with plot. the kind of titles (as well as left
>>>>>>>>>>>>>>>>>>>>>>> and side string), I want is reflected in gsn plot which am unable to do in
>>>>>>>>>>>>>>>>>>>>>>> wrf plot. I have given suitable names to scripts as well as plot for your
>>>>>>>>>>>>>>>>>>>>>>> convenience. Hope you wont get confuse.
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> In first script (1_plot_using_wrf_function)
>>>>>>>>>>>>>>>>>>>>>>> 1. after mentioning re at cnConstFEnableFill  = True,
>>>>>>>>>>>>>>>>>>>>>>> getting warning as
>>>>>>>>>>>>>>>>>>>>>>> warning:ContourPlotInitialize: scalar field is
>>>>>>>>>>>>>>>>>>>>>>> constant; no contour lines will appear; use cnConstFEnableFill to enable
>>>>>>>>>>>>>>>>>>>>>>> fill
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> actually using wrf in built function gives nice plot
>>>>>>>>>>>>>>>>>>>>>>> but unable to make self explanatory (giving strings and titles) plot so
>>>>>>>>>>>>>>>>>>>>>>> thought of ploting using gsn function
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> In second script (2_plot_using_gsn_function), issues
>>>>>>>>>>>>>>>>>>>>>>> are
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> 1. even after mentioning lat/lon getting warning as
>>>>>>>>>>>>>>>>>>>>>>> follow
>>>>>>>>>>>>>>>>>>>>>>> (0)    check_for_y_lat_coord: Warning: Data either
>>>>>>>>>>>>>>>>>>>>>>> does not contain a valid latitude coordinate array or doesn't contain one
>>>>>>>>>>>>>>>>>>>>>>> at all.
>>>>>>>>>>>>>>>>>>>>>>> (0)    A valid latitude coordinate array should have
>>>>>>>>>>>>>>>>>>>>>>> a 'units' attribute equal to one of the following values:
>>>>>>>>>>>>>>>>>>>>>>> (0)        'degrees_north' 'degrees-north'
>>>>>>>>>>>>>>>>>>>>>>> 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N'
>>>>>>>>>>>>>>>>>>>>>>> 'degreeN' 'degreesN' 'deg north'
>>>>>>>>>>>>>>>>>>>>>>> (0)    check_for_lon_coord: Warning: Data either
>>>>>>>>>>>>>>>>>>>>>>> does not contain a valid longitude coordinate array or doesn't contain one
>>>>>>>>>>>>>>>>>>>>>>> at all.
>>>>>>>>>>>>>>>>>>>>>>> (0)    A valid longitude coordinate array should
>>>>>>>>>>>>>>>>>>>>>>> have a 'units' attribute equal to one of the following values:
>>>>>>>>>>>>>>>>>>>>>>> (0)        'degrees_east' 'degrees-east'
>>>>>>>>>>>>>>>>>>>>>>> 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E'
>>>>>>>>>>>>>>>>>>>>>>> 'degreeE' 'degreesE' 'deg east'
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> 2. after mentioning re at cnConstFEnableFill  = True,
>>>>>>>>>>>>>>>>>>>>>>> getting warning as
>>>>>>>>>>>>>>>>>>>>>>> warning:ContourPlotInitialize: scalar field is
>>>>>>>>>>>>>>>>>>>>>>> constant; no contour lines will appear; use cnConstFEnableFill to enable
>>>>>>>>>>>>>>>>>>>>>>> fill
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> 3. I want all 9 plots on single panel with scales
>>>>>>>>>>>>>>>>>>>>>>> vertical, instead I am getting it on two different pages
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> I think, I have to manage or add more resources in
>>>>>>>>>>>>>>>>>>>>>>> pres1/pres2 such as (*vpXF
>>>>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpXF>*
>>>>>>>>>>>>>>>>>>>>>>> and *vpYF
>>>>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpYF>*)
>>>>>>>>>>>>>>>>>>>>>>> (*vpHeightF
>>>>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpHeightF>*
>>>>>>>>>>>>>>>>>>>>>>> and *vpWidthF
>>>>>>>>>>>>>>>>>>>>>>> <https://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml#vpWidthF>*)
>>>>>>>>>>>>>>>>>>>>>>> , so that I may get 9 plots on single panel. If you know any link which
>>>>>>>>>>>>>>>>>>>>>>> resembles like the plot i want,Please let me know. side by I will again go
>>>>>>>>>>>>>>>>>>>>>>> through links given by you.
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> On Sat, Sep 23, 2017 at 9:30 AM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> Thank you so much Mary....surely I will go through
>>>>>>>>>>>>>>>>>>>>>>>> the link provided by you....If come across any problem, will let you know
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> On Sat, Sep 23, 2017 at 12:36 AM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> Please see this page:
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> http://www.ncl.ucar.edu/Applications/wrfgsn.shtml
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> and look at the first example (wrf_gsn_1.ncl and
>>>>>>>>>>>>>>>>>>>>>>>>> wrf_nogsn_1.ncl) and the second example (wrf_gsn_2.ncl).
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> These show how to use gsn to plot data and it
>>>>>>>>>>>>>>>>>>>>>>>>> compares the plot with using WRF functions.
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> In the wrf_gsn_1.ncl plot, you can see that the
>>>>>>>>>>>>>>>>>>>>>>>>> gsnLeft/RightString resources are being set automatically, because the data
>>>>>>>>>>>>>>>>>>>>>>>>> has a description and units attribute.
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> Try this script as a start on your own data, and
>>>>>>>>>>>>>>>>>>>>>>>>> also look at wrf_gsn_2.ncl to see how to further customize your plot.
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> If you are unable to get the plot to look as nice
>>>>>>>>>>>>>>>>>>>>>>>>> as you want, email ncl-talk at ucar.edu again with
>>>>>>>>>>>>>>>>>>>>>>>>> your new script and what you want it to look like, and we can help further.
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> Good luck,
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Sep 22, 2017 at 3:21 AM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> Thank you Mary... its wrf file, I just renamed it
>>>>>>>>>>>>>>>>>>>>>>>>>> for my convenience. Actually I am able to plot it using wrf function but
>>>>>>>>>>>>>>>>>>>>>>>>>> unable to put left/right/centre strings, as those are related to gsn
>>>>>>>>>>>>>>>>>>>>>>>>>> functions. so to overcome the difficulty of left/right/centre, I thought to
>>>>>>>>>>>>>>>>>>>>>>>>>> make a plot using gsn functions injstead of wrf functions. but inbuilt wrf
>>>>>>>>>>>>>>>>>>>>>>>>>> functions gives nicer plots than gsn one.
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Sep 22, 2017 at 3:48 AM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> It's hard for somebody to simply look at a
>>>>>>>>>>>>>>>>>>>>>>>>>>> script to see if it's correct or not.  I do see some errors, but even if I
>>>>>>>>>>>>>>>>>>>>>>>>>>> correct the errors, I can't be sure it's 100% correct without having access
>>>>>>>>>>>>>>>>>>>>>>>>>>> to your data.
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> The problems are:
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> [1]
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> Your data file is called "patna_data.nc". You
>>>>>>>>>>>>>>>>>>>>>>>>>>> are using wrf_xxxx functions to plot this data. Since WRF-ARW output files
>>>>>>>>>>>>>>>>>>>>>>>>>>> are normally called something like "wrfout_d02_2005-08-29_00:00:00",
>>>>>>>>>>>>>>>>>>>>>>>>>>> I can't be sure that you are actually plotting WRF-ARW data. You can only
>>>>>>>>>>>>>>>>>>>>>>>>>>> use wrf_xxx plotting functions on WRF-ARW data.
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> [2]
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> You have:
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>  t   = a->T2(0:7,:,:)
>>>>>>>>>>>>>>>>>>>>>>>>>>>  ter = dim_avg_n( ter1, 0 )
>>>>>>>>>>>>>>>>>>>>>>>>>>>  ter = ter - 273
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> I don't know what your objective, so I don't
>>>>>>>>>>>>>>>>>>>>>>>>>>> know if the above is correct.  However, you are calling "dim_avg_n" on
>>>>>>>>>>>>>>>>>>>>>>>>>>> "ter1" which doesn't exist.  Did you mean:
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>  ter = dim_avg_n( t, 0 )
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> The above will take an average of the first 8
>>>>>>>>>>>>>>>>>>>>>>>>>>> times teps of your data, assuming that the leftmost dimension of T2
>>>>>>>>>>>>>>>>>>>>>>>>>>> represents time.
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> Further, you should use the _Wrap version of
>>>>>>>>>>>>>>>>>>>>>>>>>>> this function, so metadata is maintained:
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>  ter = dim_avg_n_Wrap( t, 0 )
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> [3]
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> You have:
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> ter = ter - 273
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> Since you are changing the units, you should
>>>>>>>>>>>>>>>>>>>>>>>>>>> update the "units" attribute to something like this:
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> ter at units = "degC"
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> The rest of the script looks okay, but again,
>>>>>>>>>>>>>>>>>>>>>>>>>>> the wrf_contour and wrf_map_overlays procedures will only work if you are
>>>>>>>>>>>>>>>>>>>>>>>>>>> actually plotting WRF-ARW data.
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> On Thu, Sep 7, 2017 at 10:02 PM, swati shewale <
>>>>>>>>>>>>>>>>>>>>>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>> Thank you. I will try with your script.
>>>>>>>>>>>>>>>>>>>>>>>>>>>> yesterday I tried it with below script. Please have a look and let me know,
>>>>>>>>>>>>>>>>>>>>>>>>>>>> whether it is correct or not.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Sep 8, 2017 at 12:54 AM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> If I understand you correctly, then you need
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> to reshape your X array so that it is NTIM x 8 x NLAT x NLON, and then you
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> can average across the second dimension from the left.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> I've attached a modified version of your
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> script, which is UNTESTED. You may need to modify it to do the right thing
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> with regard to copying over the required metadata.  I did this a little
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> with copy_VarAtts, and copy_coords.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Thu, Sep 7, 2017 at 3:42 AM, swati shewale
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Hiii,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> using ncrcat I have single file containing 24
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> time step (3hrly, 3 days). Now I want to plot temperature variable (average
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> for each day, spatial, panel plot). I am unable to take out time step for
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 1st day (0:7) and then do average of it...then second day (8:15)
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ...likewise for third day (16:23). I wrote a script, in that it takes
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> average over all time steps. what changes should I make so that I can get
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> average over specific time step? attached herewith script.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> summary of variable temperature is as follow.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Variable: X
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Type: float
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Total Size: 5658552 bytes
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>             1414638 values
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Number of Dimensions: 3
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Dimensions and sizes:    [Time | 23] x
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [south_north | 201] x [west_east | 306]
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Coordinates:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Number Of Attributes: 6
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>   FieldType :    104
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>   MemoryOrder :    XY
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>   description :    TEMP at 2 M
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>   units :    K
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>   stagger :
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>   coordinates :    XLONG XLAT XTIME
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Thu, Sep 7, 2017 at 3:04 PM, swati shewale
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Thank you. I will try with it.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Wed, Sep 6, 2017 at 8:39 PM, Mary Haley <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> haley at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> The error message is telling you that
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> there's something wrong with argument 0 passed to wrf_contour.  Argument 0
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> refers to the first argument in the list, which is "a(itime)" in this case.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> The issue is that when you use addfiles to
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> open a series of files, you need to use square brackets, [], to access a
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> single file, and not curved brackets, (). This is because addfiles returns
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> a "List" object, and lists in NCL need to be access with square brackets.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> You have:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>         contour = wrf_contour
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (a(itime),wks,t,t_res)
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> and this should be:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>         contour = wrf_contour
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (a[itime],wks,t,t_res)
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> --Mary
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Mon, Sep 4, 2017 at 4:40 AM, swati
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> shewale <shewaleswati16 at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Hiiii
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> I want to do panel plotting using WRF
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> output. I have total 8 files, each contains single time step (one day, 3
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> hrly datafiles so total 8 files). out of eight times I want to plot 6 time
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> step (panel plot) for variable slp and temperature.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> It shows below error.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> fatal:Argument type mismatch on argument
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (0) of (wrf_contour) can not coerce
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> fatal:["Execute.c":8567]:Execute: Error
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> occurred at or near line 44 in file temp.ncl
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> (0)    Error: gsn_panel: all of the plots
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> passed to gsn_panel appear to be invalid
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>  Attached herewith script.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Thank you in advance. Any help will be
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> appreciated
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ______________________________
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> _________________
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ncl-talk mailing list
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ncl-talk at ucar.edu
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> List instructions, subscriber options,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> unsubscribe:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> http://mailman.ucar.edu/mailma
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> n/listinfo/ncl-talk
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at ucar.edu
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>>
>>>>
>>>>
>>>> --
>>>> Barry H. Lynn, Ph.D
>>>> Senior Lecturer,
>>>> The Institute of the Earth Science,
>>>> The Hebrew University of Jerusalem,
>>>> Givat Ram, Jerusalem 91904, Israel
>>>> Tel: 972 547 231 170
>>>> Fax: (972)-25662581
>>>>
>>>> C.E.O, Weather It Is, LTD
>>>> Weather and Climate Focus
>>>> http://weather-it-is.com
>>>> Jerusalem, Israel
>>>> Local: 02 930 9525
>>>> Cell: 054 7 231 170
>>>> Int-IS: x972 2 930 9525
>>>> US 914 432 3108 <(914)%20432-3108>
>>>>
>>>
>>>
>>
>>
>> --
>> Barry H. Lynn, Ph.D
>> Senior Lecturer,
>> The Institute of the Earth Science,
>> The Hebrew University of Jerusalem,
>> Givat Ram, Jerusalem 91904, Israel
>> Tel: 972 547 231 170
>> Fax: (972)-25662581
>>
>> C.E.O, Weather It Is, LTD
>> Weather and Climate Focus
>> http://weather-it-is.com
>> Jerusalem, Israel
>> Local: 02 930 9525
>> Cell: 054 7 231 170
>> Int-IS: x972 2 930 9525
>> US 914 432 3108
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171103/9221ab64/attachment.html>


More information about the ncl-talk mailing list