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