NCAR staff really rude

From: mark vogel <mdvogelii_at_nyahnyahspammersnyahnyah>
Date: Tue Jun 26 2012 - 16:15:56 MDT

Dear NCL talk
Sorry I repeat the message because I think the email did not send through
because I can not see the message.
I should see the message too because I am in ncl-talk list.
We are all need help here and use the ncl program that NCAR staff develop.
We are the tax payer that paid to the government.
I think it is really rude to say that they did not get pay to answer our
questions.
So they can just answer questions what ever they like.
So they can develop the software what ever people can use it or not.

Can anyone help me out beside NCAR staff that they did not get pay to
answer questions?
or anyone please recommend me other software that much easier to use to
draw nc or hdf files.

I wish my tax money does not go to NCAR people.

Mark

---------- Forwarded message ----------
From: Dennis Shea <shea@ucar.edu>
Date: Tue, Jun 26, 2012 at 5:45 PM
Subject: Re: still can not get the plot out -> want to do montly
average only the noon time value in each grid-
To: mark vogel <mdvogelii@gmail.com>

Offline ....

PLEASE do not repeat email messages.

Everyone sees them.

However, we have our own jobs to do.
We are not paid to answer questions.

Thank you

On 06/26/2012 03:37 PM, mark vogel wrote:

> Dear all
>
> I really stump here. Surprise why my code did not work.
> I still get the same results the whole picture turn to blue
>
> please help
>
> Mark
> On Mon, Jun 25, 2012 at 7:27 PM, mark vogel <mdvogelii@gmail.com
> <mailto:mdvogelii@gmail.com>> wrote:
>
> Dear all
> I already added
>
> lat2d@units = "degrees_north"
> lon2d@units = "degrees_east"
> xavg@lat2d = lat2d
> xavg@lon2d = lon2d
>
>
> in my code file but I still got the same results.
> The values are all lesser than -30 while in the text file, there
> are values > 100.
> I defined 1e-33 as the fill_value.
> What have I done wrong?
> Thanks
> Mark
>
> On Mon, Jun 25, 2012 at 6:52 PM, David Brown <dbrown@ucar.edu
> <mailto:dbrown@ucar.edu>> wrote:
>
> Hi Mark,
>
> You have the lines:
>
>> x@lat2d = lat2d
>> x@lon2d = lon2d
>> y@lat2d = lat2d
>> y@lon2d = lon2d
>>
>
> Have you added similar code for xavg and yavg? Note that
> these are not coordinate variables (because 2D coordinate
> variables are not supported). They are attributes.
> You would be better off doing:
> lat2d@units = "degrees_north"
> lon2d@units = "degrees_east"
> xavg@lat2d = lat2d
> xavg@lon2d = lon2d
>
> etc.
> -dave
>
>
>
> On Jun 25, 2012, at 4:30 PM, mark vogel wrote:
>
> Dear Ncl-talk
>> I try to plot the average out. The values look ok but the
>> plot give me very low values (blue color in screenshot)
>> I add xavg and yavg average in there but I am not sure
>> what unit is.
>> What a unit I should put it in? I put xavg&lat2d@unit =
>> "degree_north"
>> I got the error that said
>> fatal:(lat2d) is not a named dimension in variable (xavg).
>> fatal:Execute: Error occurred at or near line 64 in file
>> hfx_nov2.ncl
>>
>> After I comment it and I still get this error message.
>> (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'
>> (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'
>>
>>
>> On Wed, Jun 13, 2012 at 1:22 PM, David Brown
>> <dbrown@ucar.edu <mailto:dbrown@ucar.edu>> wrote:
>>
>> Hi Mark,
>> You are adding the coordinate variables to x and y,
>> but you are plotting xavg and yavg. You need to add
>> the coordinate variables to xavg and yavg. Also check
>> the coordinate variables
>> to make sure they have units that match what the
>> warning messages say they should be.
>> -dave
>>
>>
>> On Jun 13, 2012, at 10:52 AM, mark vogel wrote:
>>
>> Dear All
>>>
>>> I still can not generate the plot for my code, the
>>> average look correct but when I averaged it, it has
>>> an error message about coordinate. I think I specific
>>> something wrong in graphic part but could not find
>>> what is it.
>>> Please help me.
>>>
>>> (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'
>>> (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'
>>>
>>>
>>> The code is below
>>>
>>> xavg = dim_avg_n(x1,0)
>>> yavg = dim_avg_n(y1,0)
>>> xavg@units = " (a) (Default)"
>>> yavg@units = " (b) (GEM)"
>>> printVarSummary(yavg)
>>> printVarSummary(xavg)
>>> x@lat2d = lat2d
>>> x@lon2d = lon2d
>>> y@lat2d = lat2d
>>> y@lon2d = lon2d
>>>
>>>
>>> ;**************************************************
>>> ; The file should be examined via: ncdump -v
>>> grid_type static.wrsi
>>> ; This will print the print type. then enter below.
>>> ;**************************************************
>>> projection = "LambertConformal"
>>> ;**************************************************
>>> ; create plots
>>> ;**************************************************
>>> wks = gsn_open_wks("ps" ,"sh_Nov2005avg_noon")
>>> ; ps,pdf,x11,ncgm,eps
>>> gsn_define_colormap(wks ,"BlAqGrYeOrReVi200");
>>> choose colormap
>>>
>>> res = True ; plot
>>> mods desired
>>> res@gsnSpreadColors = True ; use
>>> full range of colormap
>>> res@cnFillOn = True ;
>>> color plot desired
>>> ; res@mpMinLonF = -120 ;
>>> set min lon
>>> ; res@mpMaxLonF = -85 ;
>>> set max lon
>>> ; res@mpMinLatF = 30. ;
>>> set min lat
>>> res@vpWidthF = 0.9 ; height and width
>>> of plot
>>> res@vpHeightF = 0.4
>>>
>>> res@cnLinesOn = False ; turn
>>> off contour lines
>>> res@cnLineLabelsOn = False ; turn
>>> off contour labels
>>> res@lbLabelBarOn = False ; turn
>>> off individual lb's
>>> res@cnFillMode = "RasterFill" ;
>>> activate raster mode
>>> res@cnLevelSelectionMode = "ManualLevels" ;
>>> set manual contour levels
>>> res@cnMinLevelValF = -50 ;
>>> set min contour level
>>> res@cnMaxLevelValF = 200 ; set max
>>> contour level
>>> res@cnLevelSpacingF = 20 ;
>>> set contour spacing
>>> ;**************************************************
>>> projection = "Mercator"
>>> ; projection = "LambertConformal"
>>>
>>>
>>> dimll = dimsizes(lat2d)
>>> nlat = dimll(0)
>>> mlon = dimll(1)
>>>
>>>
>>> res@mpProjection = projection
>>> res@mpLimitMode = "Corners"
>>> res@mpLeftCornerLatF = lat2d(0,0)
>>> res@mpLeftCornerLonF = lon2d(0,0)
>>> res@mpRightCornerLatF = lat2d(nlat-1,mlon-1)
>>> res@mpRightCornerLonF = lon2d(nlat-1,mlon-1)
>>>
>>> res@mpCenterLonF = 81.45 ; set
>>> center logitude
>>>
>>> if (projection.eq."**LambertConformal") then
>>> res@mpLambertParallel1F = 30
>>> res@mpLambertParallel2F = 60
>>> res@mpLambertMeridianF = -90
>>> end if
>>>
>>> res@mpFillOn = False ; turn
>>> off map fill
>>> res@mpOutlineDrawOrder = "PostDraw" ; draw
>>> continental outline last
>>> res@mpOutlineBoundarySets = "National" ; state
>>> boundaries
>>>
>>> res@tfDoNDCOverlay = False ;
>>> True for 'native' grid
>>> res@gsnAddCyclic = False ; data
>>> are not cyclic
>>> ; res@lbOrientation = "vertical"
>>>
>>> ;**************************************************
>>> ; allocate array for 3 plots
>>> ;**************************************************
>>> plts = new (2,"graphic")
>>>
>>> ;**************************************************
>>> ; Tell NCL not to draw or advance frame for
>>> individual plots
>>> res@gsnDraw = False ; (a)
>>> do not draw
>>> res@gsnFrame = False ; (b)
>>> do not advance 'frame'
>>>
>>> plts(0) =
>>> gsn_csm_contour_map(wks,xavg(:**,:),res)
>>> plts(1) =
>>> gsn_csm_contour_map(wks,yavg(:**,:),res)
>>> ;**************************************************
>>> ; create panel: panel plots have their own set of
>>> resources
>>> ;**************************************************
>>> resP = True ;
>>> modify the panel plot
>>> resP@txString = "Avg sensible heat flux (W
>>> m-2) Nov 2005"
>>> resP@gsnMaximize = True ;
>>> maximize panel area
>>> resP@gsnPanelLabelBar = True ;
>>> add common colorbar
>>> resP@lbLabelAutoStride= True ;
>>> let NCL figure lb stride
>>> resP@gsnPanelRowSpec = True ;
>>> specify 1 top, 2 lower level
>>> gsn_panel(wks,plts,(/1,2/),**resP) ;
>>> now draw as one plot
>>>
>>> end
>>>
>>>
>>>
>>> On Tue, Jun 12, 2012 at 5:03 PM, David Brown
>>> <dbrown@ucar.edu <mailto:dbrown@ucar.edu>> wrote:
>>>
>>> Mark,
>>> If your data is hourly and you want to start at
>>> the same time every day, the stride value should
>>> be 24. In other words,
>>>
>>> do n = 5, ntim, 24
>>> ; ... n will be positioned at the first hour
>>> you want
>>> end do
>>>
>>> On Jun 12, 2012, at 2:05 PM, mark vogel wrote:
>>>
>>> Dear Ncl _talk
>>>>
>>>> I did the average during afternoon time for
>>>> latent heat flux as the following code.
>>>> But when I plot it out the values are lesser
>>>> than -30
>>>>
>>>> It does not seem right to me.
>>>> Anyone can help me what wrong with my code.
>>>> Thanks
>>>>
>>>> dirs =
>>>> "/usr/rmt_share/tgdata/hrldas/**
>>>> workspace/jam2/india_out/"
>>>> ; change directory
>>>> filesn = systemfunc("ls " + dirs +
>>>> "200511*LDASOUT*.nc")
>>>> filesn1 = systemfunc("ls " + dirs +
>>>> "2006*LDASOUT*.nc")
>>>> f = addfiles(filesn, "r") ; change to
>>>> filesn1 if 2006 )
>>>> x =
>>>> f[:]->HFX(:,:,:) ;
>>>> (Time, south_north, west_east)
>>>> dirk =
>>>> "/usr/rmt_share/tgdata/hrldas/**
>>>> workspace/jam2/india_gemout/"
>>>> filesk = systemfunc("ls " + dirk +
>>>> "200511*LDASOUT*.nc")
>>>> filesk1 = systemfunc("ls " + dirs +
>>>> "2006*LDASOUT*.nc")
>>>>
>>>> g = addfiles(filesk, "r")
>>>> y = g[:]->HFX(:,:,:)
>>>> ntim = 700
>>>> nStrt = 5
>>>> nJump = 3
>>>> nSkip = 23
>>>>
>>>> do n = 5,ntim,nSkip
>>>> n1 = n
>>>> n2 = n1+nJump
>>>> print("n1="+n1+" n2="+n2)
>>>> x1 = x(n1:n2,:,:)
>>>> y1 = y(n1:n2,:,:)
>>>> xavg = dim_avg_n(x1,0)
>>>> yavg = dim_avg_n(y1,0)
>>>> end do
>>>> xavg@units = " (a) (Default)"
>>>> yavg@units = " (b) (GEM)"
>>>> printVarSummary(yavg)
>>>> printVarSummary(xavg)
>>>>
>>>> x@lat2d = lat2d
>>>> x@lon2d = lon2d
>>>> y@lat2d = lat2d
>>>> y@lon2d = lon2d
>>>>
>>>>
>>>> ;**************************************************
>>>> ; The file should be examined via: ncdump -v
>>>> grid_type static.wrsi
>>>> ; This will print the print type. then enter below.
>>>> ;**************************************************
>>>> projection = "LambertConformal"
>>>> ;**************************************************
>>>> ; create plots
>>>> ;**************************************************
>>>> wks = gsn_open_wks("ps"
>>>> ,"sh_Nov2005avg_noon") ; ps,pdf,x11,ncgm,eps
>>>> gsn_define_colormap(wks ,"BlAqGrYeOrReVi200");
>>>> choose colormap
>>>>
>>>> res = True ;
>>>> plot mods desired
>>>> res@gsnSpreadColors = True ;
>>>> use full range of colormap
>>>> res@cnFillOn = True ;
>>>> color plot desired
>>>> ; res@mpMinLonF =
>>>> -120 ; set min lon
>>>> ; res@mpMaxLonF =
>>>> -85 ; set max lon
>>>> ; res@mpMinLatF =
>>>> 30. ; set min lat
>>>> res@vpWidthF = 0.9 ; height and
>>>> width of plot
>>>> res@vpHeightF = 0.4
>>>>
>>>> res@cnLinesOn = False ;
>>>> turn off contour lines
>>>> res@cnLineLabelsOn = False ;
>>>> turn off contour labels
>>>> res@lbLabelBarOn = False ;
>>>> turn off individual lb's
>>>> res@cnFillMode = "RasterFill" ;
>>>> activate raster mode
>>>> res@cnLevelSelectionMode =
>>>> "ManualLevels" ; set manual contour levels
>>>> res@cnMinLevelValF = -50
>>>> ; set min contour level
>>>> res@cnMaxLevelValF = 200 ; set
>>>> max contour level
>>>> res@cnLevelSpacingF = 20
>>>> ; set contour spacing
>>>> ;**************************************************
>>>> projection = "Mercator"
>>>> res@mpProjection = projection
>>>> res@mpLimitMode = "Corners"
>>>> res@mpLeftCornerLatF = lat2d(0,0)
>>>> res@mpLeftCornerLonF = lon2d(0,0)
>>>> res@mpRightCornerLatF = lat2d(nlat-1,mlon-1)
>>>> res@mpRightCornerLonF = lon2d(nlat-1,mlon-1)
>>>>
>>>> res@mpCenterLonF = 81.45 ;
>>>> set center logitude
>>>>
>>>> if (projection.eq."**LambertConformal") then
>>>> res@mpLambertParallel1F = 30
>>>> res@mpLambertParallel2F = 60
>>>> res@mpLambertMeridianF = -90
>>>> end if
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> res@mpFillOn = False ;
>>>> turn off map fill
>>>> res@mpOutlineDrawOrder = "PostDraw" ;
>>>> draw continental outline last
>>>> res@mpOutlineBoundarySets = "National" ; state
>>>> boundaries
>>>>
>>>> res@tfDoNDCOverlay = False
>>>> ; True for 'native' grid
>>>> res@gsnAddCyclic = False ;
>>>> data are not cyclic
>>>> ; res@lbOrientation = "vertical"
>>>>
>>>> ;**************************************************
>>>> ; allocate array for 3 plots
>>>> ;**************************************************
>>>> plts = new (2,"graphic")
>>>>
>>>> ;**************************************************
>>>> ; Tell NCL not to draw or advance frame for
>>>> individual plots
>>>> ;**************************************************
>>>> res@gsnDraw = False ;
>>>> (a) do not draw
>>>> res@gsnFrame = False ;
>>>> (b) do not advance 'frame'
>>>>
>>>> plts(0) =
>>>> gsn_csm_contour_map(wks,xavg(:**,:),res)
>>>> plts(1) =
>>>> gsn_csm_contour_map(wks,yavg(:**,:),res)
>>>> ;**************************************************
>>>> ; create panel: panel plots have their own set
>>>> of resources
>>>> ;**************************************************
>>>> resP = True
>>>> ; modify the panel plot
>>>> resP@txString = "Avg sensible heat
>>>> flux (W m-2) Nov 2005"
>>>> resP@gsnMaximize = True
>>>> ; maximize panel area
>>>> resP@gsnPanelLabelBar = True
>>>> ; add common colorbar
>>>> resP@lbLabelAutoStride= True
>>>> ; let NCL figure lb stride
>>>> resP@gsnPanelRowSpec = True
>>>> ; specify 1 top, 2 lower level
>>>> gsn_panel(wks,plts,(/1,2/),**resP)
>>>> ; now draw as one plot
>>>>
>>>>
>>>>
>>>> On Mon, Jun 11, 2012 at 5:53 PM, Dennis Shea
>>>> <shea@ucar.edu <mailto:shea@ucar.edu>> wrote:
>>>>
>>>> Something like
>>>>
>>>> ntim = dimsizes(time)
>>>>
>>>> nStrt = 5
>>>> nJump = 6
>>>> nSkip = 19
>>>>
>>>> do n=5,ntim-1,nSkip
>>>> n1 = n
>>>> n2 = n1+nJump
>>>> print("n1="+n1+" n2="+n2)
>>>>
>>>> x = f->X(n1:n2,:,:)
>>>> xavg = dim_avg_n(x,0)
>>>> end do
>>>>
>>>>
>>>> On 6/11/12 1:40 PM, mark vogel wrote:
>>>>
>>>> Hi Dennis
>>>> I read it but I do not get it about the
>>>> stride syntax.
>>>> I want to average the time in the
>>>> afternoon for the whole month from 05
>>>> UTC to 11 UTC
>>>> so I start with array 5 to arry 11
>>>> (5:11) and then want to skip 19
>>>> values and start again at the array 30
>>>> to 36 and then skip 19 values
>>>> again until the end. So I get the values
>>>> only afternoon out before use
>>>> the dim_avg_n
>>>>
>>>> What syntax I should use?
>>>> Mark
>>>>
>>>> On Mon, Jun 11, 2012 at 3:24 PM, Dennis
>>>> Shea <shea@ucar.edu <mailto:shea@ucar.edu>
>>>> <mailto:shea@ucar.edu
>>>> <mailto:shea@ucar.edu>>> wrote:
>>>>
>>>> If you read the documentation for
>>>> dim_avg,
>>>>
>>>> http://www.ncl.ucar.edu/__**
>>>> Document/Functions/Built-in/__**dim_avg.shtml<http://www.ncl.ucar.edu/__Document/Functions/Built-in/__dim_avg.shtml>
>>>> <http://www.ncl.ucar.edu/**
>>>> Document/Functions/Built-in/**dim_avg.shtml<http://www.ncl.ucar.edu/Document/Functions/Built-in/dim_avg.shtml>
>>>> >
>>>>
>>>>
>>>>
>>>> you will see that it performs the
>>>> average of the *rightmost* dimension.
>>>>
>>>> (Time, south_north, west_east) ;
>>>> rightmost is 'west_east'
>>>>
>>>> It is suggested that you read and use
>>>>
>>>> http://www.ncl.ucar.edu/__**
>>>> Document/Functions/Built-in/__**dim_avg_n.shtml<http://www.ncl.ucar.edu/__Document/Functions/Built-in/__dim_avg_n.shtml>
>>>>
>>>>
>>>> <http://www.ncl.ucar.edu/**
>>>> Document/Functions/Built-in/**dim_avg_n.shtml<http://www.ncl.ucar.edu/Document/Functions/Built-in/dim_avg_n.shtml>
>>>> >
>>>>
>>>> qavg = dim_avg_n(q, 0)
>>>> printVarSummary(qavg)
>>>>
>>>> or with meta data
>>>>
>>>>
>>>> http://www.ncl.ucar.edu/__**
>>>> Document/Functions/__**Contributed/dim_avg_n_Wrap.__**shtml<http://www.ncl.ucar.edu/__Document/Functions/__Contributed/dim_avg_n_Wrap.__shtml>
>>>>
>>>>
>>>> <http://www.ncl.ucar.edu/**
>>>> Document/Functions/**Contributed/dim_avg_n_Wrap.**shtml<http://www.ncl.ucar.edu/Document/Functions/Contributed/dim_avg_n_Wrap.shtml>
>>>> >
>>>>
>>>> qavg = dim_avg_n_Wrap(q, 0)
>>>> printVarSummary(qavg)
>>>>
>>>> --
>>>> Also, this does not look correct
>>>>
>>>>
>>>> f = addfiles(filesn, "r") ;
>>>> change to filesn1 if 2006 )
>>>> x = f[:]->HFX(:,:,:)
>>>> ;
>>>> (Time, south_north, west_east)
>>>> x1 = x(5:10:19,:,:)
>>>>
>>>> The 5:10:19 makes no sense. It 'say'
>>>> to use index values 5 thru 10
>>>> in strides of 19.
>>>>
>>>> Please read about subscripting at:
>>>>
>>>>
>>>> http://www.ncl.ucar.edu/__**
>>>> Document/Manuals/Ref_Manual/__**NclVariables.shtml#Subscripts<http://www.ncl.ucar.edu/__Document/Manuals/Ref_Manual/__NclVariables.shtml#Subscripts>
>>>>
>>>>
>>>> <http://www.ncl.ucar.edu/**
>>>> Document/Manuals/Ref_Manual/**NclVariables.shtml#Subscripts<http://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclVariables.shtml#Subscripts>
>>>> >
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> On 6/11/12 1:02 PM, mark vogel wrote:
>>>>
>>>> Hi all
>>>> I used dim avg to do the average
>>>> for the afternoon time for each
>>>> lat lon.
>>>> But I am not sure dim-avg is
>>>> average for each grid.
>>>> It look like it averaged from
>>>> west to east domain.
>>>> What should I used to average the
>>>> time in each grid point?
>>>> Anyone answer me please
>>>> Mark
>>>> ---------- Forwarded message
>>>> ----------
>>>> From: *mark vogel*
>>>> <mdvogelii@gmail.com
>>>> <mailto:mdvogelii@gmail.com>
>>>> <mailto:mdvogelii@gmail.com
>>>> <mailto:mdvogelii@gmail.com>>
>>>> <mailto:mdvogelii@gmail.com
>>>> <mailto:mdvogelii@gmail.com>
>>>> <mailto:mdvogelii@gmail.com
>>>> <mailto:mdvogelii@gmail.com>>>**>
>>>> Date: Mon, Jun 11, 2012 at 1:43 PM
>>>> Subject: plot dim_avg
>>>> To: NCL USERS <ncl-talk@ucar.edu
>>>> <mailto:ncl-talk@ucar.edu>
>>>> <mailto:ncl-talk@ucar.edu
>>>> <mailto:ncl-talk@ucar.edu>>
>>>> <mailto:ncl-talk@ucar.edu
>>>> <mailto:ncl-talk@ucar.edu>
>>>> <mailto:ncl-talk@ucar.edu
>>>> <mailto:ncl-talk@ucar.edu>>>>
>>>>
>>>>
>>>> Dear all
>>>> I try to average data only for
>>>> daytime and plot it out.
>>>> After I used dim_avg, then
>>>> dimesion changed from 3 to 2
>>>> dimension and
>>>> when I used
>>>> gsm_contour_map to plot it out.
>>>>
>>>> Please help
>>>> Mark
>>>>
>>>>
>>>>
>>>> filesn1 = systemfunc("ls " +
>>>> dirs + "2006*LDASOUT*.nc")
>>>> f = addfiles(filesn, "r")
>>>> ; change to filesn1 if 2006 )
>>>> x = f[:]->HFX(:,:,:)
>>>> ; (Time,
>>>> south_north, west_east)
>>>> x1 = x(5:10:19,:,:)
>>>> x1@units = " (a) (Default)"
>>>> ; print(x)
>>>> dirk =
>>>> "/usr/rmt_share/tgdata/hrldas/**
>>>> __workspace/jam2/india_gemout/**"
>>>>
>>>>
>>>> filesk = systemfunc("ls
>>>> " + dirk + "200511*LDASOUT*.nc")
>>>> filesk1 = systemfunc("ls
>>>> " + dirs + "2006*LDASOUT*.nc")
>>>>
>>>> g = addfiles(filesk, "r")
>>>> y = g[:]->HFX(:,:,:)
>>>> y1 = y(5:10:19,:,:)
>>>>
>>>> y1@units = " (b) (GEM)"
>>>> xavg = dim_avg_Wrap(x1)
>>>> yavg = dim_avg_Wrap(y1)
>>>> ; print(x1)
>>>> print(yavg)
>>>> x@lat2d = lat2d
>>>> x@lon2d = lon2d
>>>> y@lat2d = lat2d
>>>> y@lon2d = lon2d
>>>>
>>>>
>>>> res@mpFillOn =
>>>> False ; turn off map fill
>>>> res@mpOutlineDrawOrder =
>>>> "PostDraw" ; draw continental
>>>> outline last
>>>> res@mpOutlineBoundarySets =
>>>> "National" ; state boundaries
>>>>
>>>> res@tfDoNDCOverlay =
>>>> False ; True for
>>>> 'native' grid
>>>> res@gsnAddCyclic =
>>>> False ; data are not
>>>> cyclic
>>>> ; res@lbOrientation =
>>>> "vertical"
>>>>
>>>>
>>>> ;*******************************
>>>> __*******************
>>>>
>>>>
>>>> ; allocate array for 3 plots
>>>>
>>>> ;*******************************
>>>> __*******************
>>>>
>>>>
>>>> plts =
>>>> new (2,"graphic")
>>>>
>>>>
>>>> ;*******************************
>>>> __*******************
>>>>
>>>>
>>>> ; Tell NCL not to draw or advance
>>>> frame for individual plots
>>>>
>>>> ;*******************************
>>>> __*******************
>>>>
>>>>
>>>> res@gsnDraw =
>>>> False ; (a) do not draw
>>>> res@gsnFrame =
>>>> False ; (b) do not
>>>> advance 'frame'
>>>>
>>>> plts(0) =
>>>>
>>>> gsn_csm_contour_map(wks,xavg(:**__,:),res)
>>>> plts(1) =
>>>>
>>>> gsn_csm_contour_map(wks,yavg(:**__,:),res)
>>>>
>>>> ;*******************************
>>>> __*******************
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> ______________________________**
>>>> ___________________
>>>>
>>>>
>>>> ncl-talk mailing list
>>>> List instructions, subscriber
>>>> options, unsubscribe:
>>>> http://mailman.ucar.edu/__**
>>>> mailman/listinfo/ncl-talk<http://mailman.ucar.edu/__mailman/listinfo/ncl-talk>
>>>> <http://mailman.ucar.edu/**
>>>> mailman/listinfo/ncl-talk<http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>>>> >
>>>>
>>>>
>>>>
>>>> ______________________________**_________________
>>>> ncl-talk mailing list
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>>>>
>>>
>>>
>>>
>>>
>>> On Tue, Jun 12, 2012 at 5:03 PM, David Brown
>>> <dbrown@ucar.edu <mailto:dbrown@ucar.edu>> wrote:
>>>
>>> Mark,
>>> If your data is hourly and you want to start at
>>> the same time every day, the stride value should
>>> be 24. In other words,
>>>
>>> do n = 5, ntim, 24
>>> ; ... n will be positioned at the first hour
>>> you want
>>> end do
>>>
>>> On Jun 12, 2012, at 2:05 PM, mark vogel wrote:
>>>
>>> Dear Ncl _talk
>>>>
>>>> I did the average during afternoon time for
>>>> latent heat flux as the following code.
>>>> But when I plot it out the values are lesser
>>>> than -30
>>>>
>>>> It does not seem right to me.
>>>> Anyone can help me what wrong with my code.
>>>> Thanks
>>>>
>>>> dirs =
>>>> "/usr/rmt_share/tgdata/hrldas/**
>>>> workspace/jam2/india_out/"
>>>> ; change directory
>>>> filesn = systemfunc("ls " + dirs +
>>>> "200511*LDASOUT*.nc")
>>>> filesn1 = systemfunc("ls " + dirs +
>>>> "2006*LDASOUT*.nc")
>>>> f = addfiles(filesn, "r") ; change to
>>>> filesn1 if 2006 )
>>>> x =
>>>> f[:]->HFX(:,:,:) ;
>>>> (Time, south_north, west_east)
>>>> dirk =
>>>> "/usr/rmt_share/tgdata/hrldas/**
>>>> workspace/jam2/india_gemout/"
>>>> filesk = systemfunc("ls " + dirk +
>>>> "200511*LDASOUT*.nc")
>>>> filesk1 = systemfunc("ls " + dirs +
>>>> "2006*LDASOUT*.nc")
>>>>
>>>> g = addfiles(filesk, "r")
>>>> y = g[:]->HFX(:,:,:)
>>>> ntim = 700
>>>> nStrt = 5
>>>> nJump = 3
>>>> nSkip = 23
>>>>
>>>> do n = 5,ntim,nSkip
>>>> n1 = n
>>>> n2 = n1+nJump
>>>> print("n1="+n1+" n2="+n2)
>>>> x1 = x(n1:n2,:,:)
>>>> y1 = y(n1:n2,:,:)
>>>> xavg = dim_avg_n(x1,0)
>>>> yavg = dim_avg_n(y1,0)
>>>> end do
>>>> xavg@units = " (a) (Default)"
>>>> yavg@units = " (b) (GEM)"
>>>> printVarSummary(yavg)
>>>> printVarSummary(xavg)
>>>>
>>>> x@lat2d = lat2d
>>>> x@lon2d = lon2d
>>>> y@lat2d = lat2d
>>>> y@lon2d = lon2d
>>>>
>>>>
>>>> ;**************************************************
>>>> ; The file should be examined via: ncdump -v
>>>> grid_type static.wrsi
>>>> ; This will print the print type. then enter below.
>>>> ;**************************************************
>>>> projection = "LambertConformal"
>>>> ;**************************************************
>>>> ; create plots
>>>> ;**************************************************
>>>> wks = gsn_open_wks("ps"
>>>> ,"sh_Nov2005avg_noon") ; ps,pdf,x11,ncgm,eps
>>>> gsn_define_colormap(wks ,"BlAqGrYeOrReVi200");
>>>> choose colormap
>>>>
>>>> res = True ;
>>>> plot mods desired
>>>> res@gsnSpreadColors = True ;
>>>> use full range of colormap
>>>> res@cnFillOn = True ;
>>>> color plot desired
>>>> ; res@mpMinLonF =
>>>> -120 ; set min lon
>>>> ; res@mpMaxLonF =
>>>> -85 ; set max lon
>>>> ; res@mpMinLatF =
>>>> 30. ; set min lat
>>>> res@vpWidthF = 0.9 ; height and
>>>> width of plot
>>>> res@vpHeightF = 0.4
>>>>
>>>> res@cnLinesOn = False ;
>>>> turn off contour lines
>>>> res@cnLineLabelsOn = False ;
>>>> turn off contour labels
>>>> res@lbLabelBarOn = False ;
>>>> turn off individual lb's
>>>> res@cnFillMode = "RasterFill" ;
>>>> activate raster mode
>>>> res@cnLevelSelectionMode =
>>>> "ManualLevels" ; set manual contour levels
>>>> res@cnMinLevelValF = -50
>>>> ; set min contour level
>>>> res@cnMaxLevelValF = 200 ; set
>>>> max contour level
>>>> res@cnLevelSpacingF = 20
>>>> ; set contour spacing
>>>> ;**************************************************
>>>> projection = "Mercator"
>>>> res@mpProjection = projection
>>>> res@mpLimitMode = "Corners"
>>>> res@mpLeftCornerLatF = lat2d(0,0)
>>>> res@mpLeftCornerLonF = lon2d(0,0)
>>>> res@mpRightCornerLatF = lat2d(nlat-1,mlon-1)
>>>> res@mpRightCornerLonF = lon2d(nlat-1,mlon-1)
>>>>
>>>> res@mpCenterLonF = 81.45 ;
>>>> set center logitude
>>>>
>>>> if (projection.eq."**LambertConformal") then
>>>> res@mpLambertParallel1F = 30
>>>> res@mpLambertParallel2F = 60
>>>> res@mpLambertMeridianF = -90
>>>> end if
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> res@mpFillOn = False ;
>>>> turn off map fill
>>>> res@mpOutlineDrawOrder = "PostDraw" ;
>>>> draw continental outline last
>>>> res@mpOutlineBoundarySets = "National" ; state
>>>> boundaries
>>>>
>>>> res@tfDoNDCOverlay = False
>>>> ; True for 'native' grid
>>>> res@gsnAddCyclic = False ;
>>>> data are not cyclic
>>>> ; res@lbOrientation = "vertical"
>>>>
>>>> ;**************************************************
>>>> ; allocate array for 3 plots
>>>> ;**************************************************
>>>> plts = new (2,"graphic")
>>>>
>>>> ;**************************************************
>>>> ; Tell NCL not to draw or advance frame for
>>>> individual plots
>>>> ;**************************************************
>>>> res@gsnDraw = False ;
>>>> (a) do not draw
>>>> res@gsnFrame = False ;
>>>> (b) do not advance 'frame'
>>>>
>>>> plts(0) =
>>>> gsn_csm_contour_map(wks,xavg(:**,:),res)
>>>> plts(1) =
>>>> gsn_csm_contour_map(wks,yavg(:**,:),res)
>>>> ;**************************************************
>>>> ; create panel: panel plots have their own set
>>>> of resources
>>>> ;**************************************************
>>>> resP = True
>>>> ; modify the panel plot
>>>> resP@txString = "Avg sensible heat
>>>> flux (W m-2) Nov 2005"
>>>> resP@gsnMaximize = True
>>>> ; maximize panel area
>>>> resP@gsnPanelLabelBar = True
>>>> ; add common colorbar
>>>> resP@lbLabelAutoStride= True
>>>> ; let NCL figure lb stride
>>>> resP@gsnPanelRowSpec = True
>>>> ; specify 1 top, 2 lower level
>>>> gsn_panel(wks,plts,(/1,2/),**resP)
>>>> ; now draw as one plot
>>>>
>>>>
>>>>
>>>> On Mon, Jun 11, 2012 at 5:53 PM, Dennis Shea
>>>> <shea@ucar.edu <mailto:shea@ucar.edu>> wrote:
>>>>
>>>> Something like
>>>>
>>>> ntim = dimsizes(time)
>>>>
>>>> nStrt = 5
>>>> nJump = 6
>>>> nSkip = 19
>>>>
>>>> do n=5,ntim-1,nSkip
>>>> n1 = n
>>>> n2 = n1+nJump
>>>> print("n1="+n1+" n2="+n2)
>>>>
>>>> x = f->X(n1:n2,:,:)
>>>> xavg = dim_avg_n(x,0)
>>>> end do
>>>>
>>>>
>>>> On 6/11/12 1:40 PM, mark vogel wrote:
>>>>
>>>> Hi Dennis
>>>> I read it but I do not get it about the
>>>> stride syntax.
>>>> I want to average the time in the
>>>> afternoon for the whole month from 05
>>>> UTC to 11 UTC
>>>> so I start with array 5 to arry 11
>>>> (5:11) and then want to skip 19
>>>> values and start again at the array 30
>>>> to 36 and then skip 19 values
>>>> again until the end. So I get the values
>>>> only afternoon out before use
>>>> the dim_avg_n
>>>>
>>>> What syntax I should use?
>>>> Mark
>>>>
>>>> On Mon, Jun 11, 2012 at 3:24 PM, Dennis
>>>> Shea <shea@ucar.edu <mailto:shea@ucar.edu>
>>>> <mailto:shea@ucar.edu
>>>> <mailto:shea@ucar.edu>>> wrote:
>>>>
>>>> If you read the documentation for
>>>> dim_avg,
>>>>
>>>> http://www.ncl.ucar.edu/____**
>>>> Document/Functions/Built-in/__**__dim_avg.shtml<http://www.ncl.ucar.edu/____Document/Functions/Built-in/____dim_avg.shtml>
>>>> <http://www.ncl.ucar.edu/__**
>>>> Document/Functions/Built-in/__**dim_avg.shtml<http://www.ncl.ucar.edu/__Document/Functions/Built-in/__dim_avg.shtml>
>>>> >
>>>> <http://www.ncl.ucar.edu/__**
>>>> Document/Functions/Built-in/__**dim_avg.shtml<http://www.ncl.ucar.edu/__Document/Functions/Built-in/__dim_avg.shtml>
>>>> <http://www.ncl.ucar.edu/**
>>>> Document/Functions/Built-in/**dim_avg.shtml<http://www.ncl.ucar.edu/Document/Functions/Built-in/dim_avg.shtml>
>>>> >>
>>>>
>>>>
>>>>
>>>> you will see that it performs the
>>>> average of the *rightmost* dimension.
>>>>
>>>> (Time, south_north, west_east) ;
>>>> rightmost is 'west_east'
>>>>
>>>> It is suggested that you read and use
>>>>
>>>> http://www.ncl.ucar.edu/____**
>>>> Document/Functions/Built-in/__**__dim_avg_n.shtml<http://www.ncl.ucar.edu/____Document/Functions/Built-in/____dim_avg_n.shtml>
>>>> <http://www.ncl.ucar.edu/__**
>>>> Document/Functions/Built-in/__**dim_avg_n.shtml<http://www.ncl.ucar.edu/__Document/Functions/Built-in/__dim_avg_n.shtml>
>>>> >
>>>>
>>>>
>>>> <http://www.ncl.ucar.edu/__**
>>>> Document/Functions/Built-in/__**dim_avg_n.shtml<http://www.ncl.ucar.edu/__Document/Functions/Built-in/__dim_avg_n.shtml>
>>>> <http://www.ncl.ucar.edu/**
>>>> Document/Functions/Built-in/**dim_avg_n.shtml<http://www.ncl.ucar.edu/Document/Functions/Built-in/dim_avg_n.shtml>
>>>> >>
>>>>
>>>> qavg = dim_avg_n(q, 0)
>>>> printVarSummary(qavg)
>>>>
>>>> or with meta data
>>>>
>>>>
>>>> http://www.ncl.ucar.edu/____**
>>>> Document/Functions/____**Contributed/dim_avg_n_Wrap.___**_shtml<http://www.ncl.ucar.edu/____Document/Functions/____Contributed/dim_avg_n_Wrap.____shtml>
>>>> <http://www.ncl.ucar.edu/__**
>>>> Document/Functions/__**Contributed/dim_avg_n_Wrap.__**shtml<http://www.ncl.ucar.edu/__Document/Functions/__Contributed/dim_avg_n_Wrap.__shtml>
>>>> >
>>>>
>>>>
>>>> <http://www.ncl.ucar.edu/__**
>>>> Document/Functions/__**Contributed/dim_avg_n_Wrap.__**shtml<http://www.ncl.ucar.edu/__Document/Functions/__Contributed/dim_avg_n_Wrap.__shtml>
>>>> <http://www.ncl.ucar.edu/**
>>>> Document/Functions/**Contributed/dim_avg_n_Wrap.**shtml<http://www.ncl.ucar.edu/Document/Functions/Contributed/dim_avg_n_Wrap.shtml>
>>>> >>
>>>>
>>>> qavg = dim_avg_n_Wrap(q, 0)
>>>> printVarSummary(qavg)
>>>>
>>>> --
>>>> Also, this does not look correct
>>>>
>>>>
>>>> f = addfiles(filesn, "r") ;
>>>> change to filesn1 if 2006 )
>>>> x = f[:]->HFX(:,:,:)
>>>> ;
>>>> (Time, south_north, west_east)
>>>> x1 = x(5:10:19,:,:)
>>>>
>>>> The 5:10:19 makes no sense. It 'say'
>>>> to use index values 5 thru 10
>>>> in strides of 19.
>>>>
>>>> Please read about subscripting at:
>>>>
>>>>
>>>> http://www.ncl.ucar.edu/____**
>>>> Document/Manuals/Ref_Manual/__**__NclVariables.shtml#**Subscripts<http://www.ncl.ucar.edu/____Document/Manuals/Ref_Manual/____NclVariables.shtml#Subscripts>
>>>> <http://www.ncl.ucar.edu/__**
>>>> Document/Manuals/Ref_Manual/__**NclVariables.shtml#Subscripts<http://www.ncl.ucar.edu/__Document/Manuals/Ref_Manual/__NclVariables.shtml#Subscripts>
>>>> >
>>>>
>>>>
>>>> <http://www.ncl.ucar.edu/__**
>>>> Document/Manuals/Ref_Manual/__**NclVariables.shtml#Subscripts<http://www.ncl.ucar.edu/__Document/Manuals/Ref_Manual/__NclVariables.shtml#Subscripts>
>>>> <http://www.ncl.ucar.edu/**
>>>> Document/Manuals/Ref_Manual/**NclVariables.shtml#Subscripts<http://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclVariables.shtml#Subscripts>
>>>> >**>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> On 6/11/12 1:02 PM, mark vogel wrote:
>>>>
>>>> Hi all
>>>> I used dim avg to do the average
>>>> for the afternoon time for each
>>>> lat lon.
>>>> But I am not sure dim-avg is
>>>> average for each grid.
>>>> It look like it averaged from
>>>> west to east domain.
>>>> What should I used to average the
>>>> time in each grid point?
>>>> Anyone answer me please
>>>> Mark
>>>> ---------- Forwarded message
>>>> ----------
>>>> From: *mark vogel*
>>>> <mdvogelii@gmail.com
>>>> <mailto:mdvogelii@gmail.com>
>>>> <mailto:mdvogelii@gmail.com
>>>> <mailto:mdvogelii@gmail.com>>
>>>> <mailto:mdvogelii@gmail.com
>>>> <mailto:mdvogelii@gmail.com>
>>>> <mailto:mdvogelii@gmail.com
>>>> <mailto:mdvogelii@gmail.com>>>**>
>>>> Date: Mon, Jun 11, 2012 at 1:43 PM
>>>> Subject: plot dim_avg
>>>> To: NCL USERS <ncl-talk@ucar.edu
>>>> <mailto:ncl-talk@ucar.edu>
>>>> <mailto:ncl-talk@ucar.edu
>>>> <mailto:ncl-talk@ucar.edu>>
>>>> <mailto:ncl-talk@ucar.edu
>>>> <mailto:ncl-talk@ucar.edu>
>>>> <mailto:ncl-talk@ucar.edu
>>>> <mailto:ncl-talk@ucar.edu>>>>
>>>>
>>>>
>>>> Dear all
>>>> I try to average data only for
>>>> daytime and plot it out.
>>>> After I used dim_avg, then
>>>> dimesion changed from 3 to 2
>>>> dimension and
>>>> when I used
>>>> gsm_contour_map to plot it out.
>>>>
>>>> Please help
>>>> Mark
>>>>
>>>>
>>>>
>>>> filesn1 = systemfunc("ls " +
>>>> dirs + "2006*LDASOUT*.nc")
>>>> f = addfiles(filesn, "r")
>>>> ; change to filesn1 if 2006 )
>>>> x = f[:]->HFX(:,:,:)
>>>> ; (Time,
>>>> south_north, west_east)
>>>> x1 = x(5:10:19,:,:)
>>>> x1@units = " (a) (Default)"
>>>> ; print(x)
>>>> dirk =
>>>> "/usr/rmt_share/tgdata/hrldas/**
>>>> ____workspace/jam2/india_**gemout/__"
>>>>
>>>>
>>>> filesk = systemfunc("ls
>>>> " + dirk + "200511*LDASOUT*.nc")
>>>> filesk1 = systemfunc("ls
>>>> " + dirs + "2006*LDASOUT*.nc")
>>>>
>>>> g = addfiles(filesk, "r")
>>>> y = g[:]->HFX(:,:,:)
>>>> y1 = y(5:10:19,:,:)
>>>>
>>>> y1@units = " (b) (GEM)"
>>>> xavg = dim_avg_Wrap(x1)
>>>> yavg = dim_avg_Wrap(y1)
>>>> ; print(x1)
>>>> print(yavg)
>>>> x@lat2d = lat2d
>>>> x@lon2d = lon2d
>>>> y@lat2d = lat2d
>>>> y@lon2d = lon2d
>>>>
>>>>
>>>> res@mpFillOn =
>>>> False ; turn off map fill
>>>> res@mpOutlineDrawOrder =
>>>> "PostDraw" ; draw continental
>>>> outline last
>>>> res@mpOutlineBoundarySets =
>>>> "National" ; state boundaries
>>>>
>>>> res@tfDoNDCOverlay =
>>>> False ; True for
>>>> 'native' grid
>>>> res@gsnAddCyclic =
>>>> False ; data are not
>>>> cyclic
>>>> ; res@lbOrientation =
>>>> "vertical"
>>>>
>>>>
>>>> ;*******************************
>>>> ____*******************
>>>>
>>>>
>>>> ; allocate array for 3 plots
>>>>
>>>> ;*******************************
>>>> ____*******************
>>>>
>>>>
>>>> plts =
>>>> new (2,"graphic")
>>>>
>>>>
>>>> ;*******************************
>>>> ____*******************
>>>>
>>>>
>>>> ; Tell NCL not to draw or advance
>>>> frame for individual plots
>>>>
>>>> ;*******************************
>>>> ____*******************
>>>>
>>>>
>>>> res@gsnDraw =
>>>> False ; (a) do not draw
>>>> res@gsnFrame =
>>>> False ; (b) do not
>>>> advance 'frame'
>>>>
>>>> plts(0) =
>>>>
>>>> gsn_csm_contour_map(wks,xavg(:**
>>>> ____,:),res)
>>>> plts(1) =
>>>>
>>>> gsn_csm_contour_map(wks,yavg(:**
>>>> ____,:),res)
>>>>
>>>> ;*******************************
>>>> ____*******************
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> ______________________________**
>>>> _____________________
>>>>
>>>>
>>>> ncl-talk mailing list
>>>> List instructions, subscriber
>>>> options, unsubscribe:
>>>> http://mailman.ucar.edu/____**
>>>> mailman/listinfo/ncl-talk<http://mailman.ucar.edu/____mailman/listinfo/ncl-talk>
>>>> <http://mailman.ucar.edu/__**
>>>> mailman/listinfo/ncl-talk<http://mailman.ucar.edu/__mailman/listinfo/ncl-talk>
>>>> >
>>>> <http://mailman.ucar.edu/__**
>>>> mailman/listinfo/ncl-talk<http://mailman.ucar.edu/__mailman/listinfo/ncl-talk>
>>>> <http://mailman.ucar.edu/**
>>>> mailman/listinfo/ncl-talk<http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>>>> >>
>>>>
>>>>
>>>>
>>>> ______________________________**_________________
>>>> ncl-talk mailing list
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>>>>
>>>
>>>
>>> ______________________________**_________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>>>
>>
>>
>> <hfx_nov2.ncl><hfx.png>
>>
>
>
>
>
>
>
> ______________________________**_________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jun 26 16:16:16 2012

This archive was generated by hypermail 2.1.8 : Fri Jun 29 2012 - 15:46:14 MDT