From: <debasish_at_nyahnyahspammersnyahnyah>

Date: Thu, 25 Jun 2009 12:03:40 -0800

Date: Thu, 25 Jun 2009 12:03:40 -0800

Dear Ibraheem

Please check the web page

http://www.ncl.ucar.edu/Document/Graphics/error_msg.shtml

It will show you how to get rid of this error massage.

Thanks

Debasish

Quoting Ibraheem Khan <khan18_at_atmos.uiuc.edu>:

*> printVarSummary(x1):
*

*>
*

*> Variable: x1
*

*> Type: float
*

*> Total Size: 585728 bytes
*

*> 146432 values
*

*> Number of Dimensions: 3
*

*> Dimensions and sizes: [time | 44] x [lat | 52] x [lon | 64]
*

*> Coordinates:
*

*> time: [ 0..516]
*

*> lat: [24.125..36.875]
*

*> lon: [61.125..76.875]
*

*> Number Of Attributes: 1
*

*>
*

*> and
*

*>
*

*> printVarSummary(c):
*

*>
*

*> Variable: c
*

*> Type: float
*

*> Total Size: 13312 bytes
*

*> 3328 values
*

*> Number of Dimensions: 2
*

*> Dimensions and sizes: [52] x [64]
*

*> Coordinates:
*

*> Number Of Attributes: 1
*

*> _FillValue : -999
*

*>
*

*>
*

*>
*

*> Please guide how to handle this data...as i have started over scripts of NCL.
*

*> Thanks
*

*>
*

*> -Ibraheem
*

*>
*

*> ----- Original Message -----
*

*> From: "Mary Haley" <haley_at_ucar.edu>
*

*> To: "Ibraheem Khan" <khan18_at_atmos.uiuc.edu>
*

*> Cc: ncl-talk_at_ucar.edu
*

*> Sent: Thursday, June 25, 2009 2:46:55 PM GMT -06:00 US/Canada Central
*

*> Subject: Re: question
*

*>
*

*>
*

*> In order to plot data over a map, your data must either have 1D
*

*> latitude/longitude coordinate arrays attached to it, or you must have
*

*> 2D lat/lon arrays the same size as your data array containing the
*

*> lat/lon locations of your data.
*

*>
*

*> It looks like you attempted to copy the coordinate information
*

*> from "x1" to "c", but that it may not have been successful.
*

*>
*

*> What does:
*

*>
*

*> printVarSummary(x1)
*

*> printVarSummary(c)
*

*>
*

*> report?
*

*>
*

*> --Mary
*

*>
*

*>
*

*> On Thu, 25 Jun 2009, Ibraheem Khan wrote:
*

*>
*

*>> Mary,
*

*>>
*

*>> Now its showing map but without data. And giving following:
*

*>>
*

*>> (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'
*

*>>
*

*>> How to handle this?
*

*>>
*

*>> -Ibraheem
*

*>>
*

*>>
*

*>>
*

*>>
*

*>>
*

*>> ----- Original Message -----
*

*>> From: "Mary Haley" <haley_at_ucar.edu>
*

*>> To: "Ibraheem Khan" <khan18_at_atmos.uiuc.edu>
*

*>> Cc: ncl-talk_at_ucar.edu
*

*>> Sent: Thursday, June 25, 2009 2:32:52 PM GMT -06:00 US/Canada Central
*

*>> Subject: Re: question
*

*>>
*

*>>
*

*>> You have
*

*>>
*

*>> res_at_gsnDraw = False
*

*>>
*

*>> This turns off drawing by the gsn_csm_xxxx function.
*

*>>
*

*>> Remove or comment this line out.
*

*>>
*

*>> --Mary
*

*>>
*

*>>
*

*>> On Thu, 25 Jun 2009, Ibraheem Khan wrote:
*

*>>
*

*>>> Hi David,
*

*>>>
*

*>>> I am using following script to calculate the point time series
*

*>>> correlation with the time series grid data. But no output is
*

*>>> turning up. Can you or anyone look at it to point out anything
*

*>>> missing.
*

*>>> Thanks
*

*>>>
*

*>>> Script:
*

*>>>
*

*>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
*

*>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
*

*>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
*

*>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
*

*>>>
*

*>>> begin
*

*>>> ;**********************************************
*

*>>> ; Reading Data and correlation function
*

*>>> ;**********************************************
*

*>>> f1 = addfile("6104_01.nc","r")
*

*>>> f2 = addfile("avg_6104_01.nc","r")
*

*>>> x1 = f1->a ;time series grid data
*

*>>> x2 = f2->a ;point time series
*

*>>> print(x1)
*

*>>> print(x2)
*

*>>> st = dim_standardize_Wrap(x2,1)
*

*>>> printVarSummary(st)
*

*>>> c = escorc(x1(lat|:,lon|:,time|:),x2)
*

*>>> copy_VarCoords(x1(0,:,:),c)
*

*>>> printVarSummary(c)
*

*>>>
*

*>>> ;**********************************************
*

*>>> ; create plot
*

*>>> ;*********************************************
*

*>>> wks = gsn_open_wks ("x11","cor_prcp_01") ; open ps file
*

*>>>
*

*>>> gsn_define_colormap(wks,"radar")
*

*>>> res = True
*

*>>> res_at_gsnDraw = False
*

*>>> res_at_cnFillOn = True
*

*>>> res_at_cnLinesOn = False
*

*>>>
*

*>>> res_at_gsnSpreadColors = True
*

*>>> res_at_lbAutoLabelStride = True
*

*>>>
*

*>>>
*

*>>> res_at_cnLevelselectionMode = "ManualLevels" ; manually set cn levels
*

*>>> res_at_cnMinLevelValF = -1. ; min level
*

*>>> res_at_cnMaxLevelValF = 1. ; max level
*

*>>> res_at_tiMainString = ":F25:Prcp Point Corr"
*

*>>>
*

*>>> res_at_gsnLeftString = " "
*

*>>> res_at_tiMainFontHeightF = 0.03
*

*>>> res_at_gsnAddCyclic = False
*

*>>> ;printVarSummary(c)
*

*>>> plot = gsn_csm_contour_map_ce(wks,c,res)
*

*>>> end
*

*>>>
*

*>>>
*

*>>>
*

*>>> -Ibraheem
*

*>>>
*

*>>>
*

*>>> ----- Original Message -----
*

*>>> From: "Dave Allured" <dave.allured_at_noaa.gov>
*

*>>> To: ncl-talk_at_ucar.edu
*

*>>> Sent: Wednesday, June 24, 2009 8:58:50 PM GMT -06:00 US/Canada Central
*

*>>> Subject: Re: question
*

*>>>
*

*>>> Ibraheem,
*

*>>>
*

*>>> Point to grid correlation is easy. At each grid point, just compute
*

*>>> the correlation between the single point time series and the time
*

*>>> series for the grid point. The output from this is a single grid of
*

*>>> correlation numbers. Just plot this field, using any generic
*

*>>> contour plot function from NCL.
*

*>>>
*

*>>> How do you define correlation between two unequal regions? There
*

*>>> are various definitions. First you need to have a good mathematical
*

*>>> description of what you want, independent from any programming
*

*>>> considerations.
*

*>>>
*

*>>> One definition is a simple cross correlation where single grid
*

*>>> points at matching coordinates are correlated. If the regions are
*

*>>> different sizes, then this makes a plot that is only the size of the
*

*>>> small region, and blank in the remainder of the large region. This
*

*>>> is the third example on the correlation example page:
*

*>>>
*

*>>> http://www.ncl.ucar.edu/Applications/corel.shtml
*

*>>>
*

*>>> Another definition is to reduce the small region to a single time
*

*>>> series by a method of your choice, such as weighted or unweighted
*

*>>> averaging. Then compute and plot the gridded correlation of the
*

*>>> average time series to the larger region, similar to the single
*

*>>> point case.
*

*>>>
*

*>>> Yet another definition is to start at one grid point in the larger
*

*>>> region, and compute independent correlations to each point in the
*

*>>> smaller region. Then combine the resulting correlation numbers to a
*

*>>> single number by your choice of method, such as taking the mean or
*

*>>> median. Repeat over all grid points of the larger region, and once
*

*>>> again you have a single grid of numbers to plot in the usual way.
*

*>>>
*

*>>> Once you have decided on the definition, NCL has good matrix
*

*>>> capability to help with the computations.
*

*>>>
*

*>>> Dave Allured
*

*>>> CU/CIRES Climate Diagnostics Center (CDC)
*

*>>> http://cires.colorado.edu/science/centers/cdc/
*

*>>> NOAA/ESRL/PSD, Climate Analysis Branch (CAB)
*

*>>> http://www.cdc.noaa.gov/psd1/
*

*>>>
*

*>>> Ibraheem Khan wrote:
*

*>>>> Hi,
*

*>>>>
*

*>>>> I am thinking about drawing a correlation map. The correlation
*

*>>>> must be between a point value and gridded data. Also a
*

*>>>> correlation between two unequal gridded datasets (one covering a
*

*>>>> small region and the other a larger domain). Is it possible to do
*

*>>>> both correlations?
*

*>>>> Thanks
*

*>>>>
*

*>>>> -Ibraheem
*

*>>>>
*

*>>>>
*

*>>>> _______________________________________________
*

*>>>> ncl-talk mailing list
*

*>>>> List instructions, subscriber options, unsubscribe:
*

*>>>> 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
*

*>>> _______________________________________________
*

*>>> ncl-talk mailing list
*

*>>> List instructions, subscriber options, unsubscribe:
*

*>>> 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
*

*>
*

----------------------------------------------------------------

This message was sent using IMP, the Internet Messaging Program.

_______________________________________________

ncl-talk mailing list

List instructions, subscriber options, unsubscribe:

http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Thu Jun 25 2009 - 14:03:40 MDT

*
This archive was generated by hypermail 2.2.0
: Wed Jul 01 2009 - 06:11:45 MDT
*