[ncl-talk] Error while handling multiple WRF output files
Barry Lynn
barry.h.lynn at gmail.com
Fri Nov 3 07:50:09 MDT 2017
Hi Swati:
The errors tell you that you have a array size mismatch.
If you plotVarSummary and look at the dimension sizes, you will probably be
able to fix your error.
Barry
On Fri, Nov 3, 2017 at 2:48 PM, swati shewale <shewaleswati16 at gmail.com>
wrote:
> 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 <(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/5cd4651f/attachment.html>
More information about the ncl-talk
mailing list