# how to deal the negtive value?

From: yanling_1220 <yanling_1220_at_nyahnyahspammersnyahnyah>
Date: Thu, 10 Sep 2009 22:23:32 +0800 (CST)

Hi£¬

When I use "obj_anal_ic" to interpolate the observed station rainfall to grid£¬and find some grid with a negative value£¬and it's no possible obviously,and how do you solve this problem usually? If I want to using "0" to fill this value,how should modify the following main segment of my script,any help will be appreciate.

----------------------------------------------------------------------------------------------------------------------
fo =new((/nlat,nlon,ntim/),typeof(xObs),-999)
do nmo =0,ntim-1
rscan = (/5,3,2,1,0.5/)
opt = True
opt_at_zonal = True ; use zonal averages as a 1st guess
fo(:,:,nmo) = obj_anal_ic(lon1,lat1,xObs(nmo,:),lon,lat, rscan,opt) ;(nlat*mlon*ntim)
end do

fo!0 = "lat"
fo!1 = "lon"
fo!2 = "time"
pr_ratio =fo(time|:,lat|:,lon|:)
dim_pr_ratio = dimsizes( pr_ratio)
ntim = dim_pr_ratio(0)
nlat = dim_pr_ratio(1)
mlon = dim_pr_ratio(2)

pr_ratio_at_original_name = " observed pr_ratio "
pr_ratio_at_coordinates = " height"
pr_ratio@_FillValue = 1e+20 ;define the FillValue of the pr_ratio in the output
;nc file,it can't be "1e+020 " with the quotes ("...")
pr_ratio_at_units = " mm/d"
pr_ratio_at_long_name = "observed grid pr_ratio "
printVarSummary(pr_ratio)

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

2009-09-10£¬ncl-talk-request_at_ucar.edu £º
>Send ncl-talk mailing list submissions to
> ncl-talk_at_ucar.edu
>
>To subscribe or unsubscribe via the World Wide Web, visit
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>or, via email, send a message with subject or body 'help' to
> ncl-talk-request_at_ucar.edu
>
>You can reach the person managing the list at
> ncl-talk-owner_at_ucar.edu
>
>When replying, please edit your Subject line so it is more specific
>than "Re: Contents of ncl-talk digest..."
>
>
>Today's Topics:
>
> 1. Array Error (DJ Rasmussen)
> 2. Drawing order for polymarkers (Dave Allured)
> 3. Re: Array Error (Dave Allured)
> 4. Re: Array Error (David Brown)
> 5. Re: Array Error (Dennis Shea)
> 6. Correlation significant test.... (Sabeerali(sebi))
> 7. Re: Correlation significant test.... (Debasish PaiMazumder)
>
>
>----------------------------------------------------------------------
>
>Message: 1
>Date: Wed, 09 Sep 2009 17:57:45 -0500
>From: DJ Rasmussen <dj.rasmussen_at_ssec.wisc.edu>
>Subject: [ncl-talk] Array Error
>To: ncl-talk_at_ucar.edu
>Message-ID: <4AA832E9.80608_at_ssec.wisc.edu>
>Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
>Hi NCL users,
>
>I am attempting to make an NCL script that reads in multiple monthly
>NARR files with wind data and then averages them all and, finally
>outputs the averages to a netCDF file that I can use to make a contour
>plot with. To do the averages, however, I need to create an array. For
>some reason, NCL will not allow me to do it the way I am trying to: by
>multiplying lat and lon dimensions to get the number of data points for
>wind values at every lat/lon on the NARR grid.
>
>As always, your assistance is appreciated.
>
>Cheers
>
>
>
>
>
>Using NCL version 5.1.1
>;*************************************************
>; NARRAVG.ncl
>;
>;************************************************
>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"
>
>;************************************************
>begin
>
>;************************************************
>; read in multiple grib files
>;************************************************
>diri = "./" ; directory where files reside
>fils = systemfunc("ls "+diri+"narr*.grb")
>
>setfileoption("grb","SingleElementDimensions","initial_time")
>
>f = addfiles(fils,"r"); read in all files in the folder
>
>;*************************Get Variables to Average**********
>U_GRD_221 = addfiles_GetVar(f, fils, "U_GRD_221_HTGL_113") ;Get these
>variables from all files
>V_GRD_221 = addfiles_GetVar(f, fils, "V_GRD_221_HTGL_113")
>
>printVarSummary(U_GRD_221)
>printVarSummary(V_GRD_221)
>
>dimU_GRD = dimsizes(U_GRD_221)
>
>gridx_221 = dimU_GRD(2) ; get dimensions of lat
>gridy_221 = dimU_GRD(3)
>wind_data = dimU_GRD(2)*dimU_GRD(3) ; multiply lat and lon
>dimensions to get wind data points needed
>
>store = new((/wind_data,gridx_221,gridy_221/), "float") ; create a new
>array for new netCDF file.
>
>
>;*************************Calculate Average******
>
>U_GRD_221_AVG = dim_avg_n_Wrap(U_GRD_221(:,0,:,:),0) ; calculate
>average of zonal wind at level 0 for all grid points
>V_GRD_221_AVG = dim_avg_n_Wrap(V_GRD_221(:,0,:,:),0)
>
>
>wspd = sqrt( U_GRD_221_AVG2 + V_GRD_221_AVG2 ) ; calculate wind speed
>magnitude
>
>
>;*************create NetCDF output*******************
>store!0 = "wspd" ; assign named dimensions
>store!1 = "gridx_221"
>store!2 = "gridy_221"
>
>store_at_long_name = "Avg Wind Speed"
>store_at_units = "m/s"
>
>system("/bin/rm -f dailyAvg.WSPD10_NARR.nc") ; remove any pre-existing
>file
>ncdf = addfile("dailyAvg.WSPD10_NARR.nc","c")
>filedimdef(ncdf,"initial_time0_hours",-1,True)
>fileattdef (ncdf, f)
>
>ncdf->U_GRD_221_HTGL_113 = store
>
>end
>
>
>
>--
>
>DJ Rasmussen
>
>Antarctic Meteorological Research Center
>
>University of Wisconsin - Madison Space Science and Engineering Center
>
>1225 W. Dayton St., Rm 947, Madison, WI 53706
>
>Email: dj.rasmussen_at_ssec.wisc.edu
>
>http://amrc.ssec.wisc.edu
>
>
>
>
>------------------------------
>
>Message: 2
>Date: Wed, 09 Sep 2009 17:29:51 -0600
>From: Dave Allured <dave.allured_at_noaa.gov>
>Subject: [ncl-talk] Drawing order for polymarkers
>To: ncl-talk_at_ucar.edu
>Message-ID: <4AA83A6F.2020903_at_noaa.gov>
>Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
>NCL team,
>
>I need to control the drawing order of polylines and polymarkers.
>In a typical case, the base plot is made with
>gsn_csm_vector_scalar_map_ce, plus an overlay made with gsn_csm_contour.
>
>I find a number of DrawOrder resources, but none seem to apply to
>objects made with gsn_add_poly* functions. In general, what is the
>right way to control drawing order for these? I'm using NCL 5.1.1.
> Thanks for your help.
>
>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/
>
>
>------------------------------
>
>Message: 3
>Date: Wed, 09 Sep 2009 17:40:13 -0600
>From: Dave Allured <dave.allured_at_noaa.gov>
>Subject: Re: [ncl-talk] Array Error
>To: ncl-talk_at_ucar.edu
>Message-ID: <4AA83CDD.6000706_at_noaa.gov>
>Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
>store = new((/wind_data,gridx_221,gridy_221/), "float")
>
>This line makes a 3-D array of dimensions nx*ny, nx, ny. The first
>dimension does not seem right to me. Should this be the time
>dimension or the level dimension instead?
>
>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/
>
>DJ Rasmussen wrote:
>> Hi NCL users,
>>
>> I am attempting to make an NCL script that reads in multiple monthly
>> NARR files with wind data and then averages them all and, finally
>> outputs the averages to a netCDF file that I can use to make a contour
>> plot with. To do the averages, however, I need to create an array. For
>> some reason, NCL will not allow me to do it the way I am trying to: by
>> multiplying lat and lon dimensions to get the number of data points for
>> wind values at every lat/lon on the NARR grid.
>>
>> As always, your assistance is appreciated.
>>
>> Cheers
>>
>>
>>
>>
>>
>> Using NCL version 5.1.1
>> ;*************************************************
>> ; NARRAVG.ncl
>> ;
>> ;************************************************
>> 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"
>>
>> ;************************************************
>> begin
>>
>> ;************************************************
>> ; read in multiple grib files
>> ;************************************************
>> diri = "./" ; directory where files reside
>> fils = systemfunc("ls "+diri+"narr*.grb")
>>
>> setfileoption("grb","SingleElementDimensions","initial_time")
>>
>> f = addfiles(fils,"r"); read in all files in the folder
>>
>> ;*************************Get Variables to Average**********
>> U_GRD_221 = addfiles_GetVar(f, fils, "U_GRD_221_HTGL_113") ;Get these
>> variables from all files
>> V_GRD_221 = addfiles_GetVar(f, fils, "V_GRD_221_HTGL_113")
>>
>> printVarSummary(U_GRD_221)
>> printVarSummary(V_GRD_221)
>>
>> dimU_GRD = dimsizes(U_GRD_221)
>>
>> gridx_221 = dimU_GRD(2) ; get dimensions of lat
>> gridy_221 = dimU_GRD(3)
>> wind_data = dimU_GRD(2)*dimU_GRD(3) ; multiply lat and lon
>> dimensions to get wind data points needed
>>
>> store = new((/wind_data,gridx_221,gridy_221/), "float") ; create a new
>> array for new netCDF file.
>>
>>
>> ;*************************Calculate Average******
>>
>> U_GRD_221_AVG = dim_avg_n_Wrap(U_GRD_221(:,0,:,:),0) ; calculate
>> average of zonal wind at level 0 for all grid points
>> V_GRD_221_AVG = dim_avg_n_Wrap(V_GRD_221(:,0,:,:),0)
>>
>>
>> wspd = sqrt( U_GRD_221_AVG2 + V_GRD_221_AVG2 ) ; calculate wind speed
>> magnitude
>>
>>
>> ;*************create NetCDF output*******************
>> store!0 = "wspd" ; assign named dimensions
>> store!1 = "gridx_221"
>> store!2 = "gridy_221"
>>
>> store_at_long_name = "Avg Wind Speed"
>> store_at_units = "m/s"
>>
>> system("/bin/rm -f dailyAvg.WSPD10_NARR.nc") ; remove any pre-existing
>> file
>> ncdf = addfile("dailyAvg.WSPD10_NARR.nc","c")
>> filedimdef(ncdf,"initial_time0_hours",-1,True)
>> fileattdef (ncdf, f)
>>
>> ncdf->U_GRD_221_HTGL_113 = store
>>
>> end
>>
>>
>>
>
>
>------------------------------
>
>Message: 4
>Date: Wed, 9 Sep 2009 18:08:12 -0600
>From: David Brown <dbrown_at_ucar.edu>
>Subject: Re: [ncl-talk] Array Error
>To: ncl-talk_at_ucar.edu
>Message-ID: <66FA3C42-59C8-47C2-A817-F5AD9F1B9D00_at_ucar.edu>
>Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes
>
>
>I'll add that for GRIB grid 221, nx is 349 and ny is 277. This makes
>the memory you are trying to allocate for the "store" variable about
>37 Gigabytes -- more than all but a few supercomputer class systems
>could supply.
>
>
>On Sep 9, 2009, at 5:40 PM, Dave Allured wrote:
>
>> store = new((/wind_data,gridx_221,gridy_221/), "float")
>>
>> This line makes a 3-D array of dimensions nx*ny, nx, ny. The first
>> dimension does not seem right to me. Should this be the time
>> dimension or the level dimension instead?
>>
>> 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/
>>
>> DJ Rasmussen wrote:
>>> Hi NCL users,
>>>
>>> I am attempting to make an NCL script that reads in multiple monthly
>>> NARR files with wind data and then averages them all and, finally
>>> outputs the averages to a netCDF file that I can use to make a
>>> contour
>>> plot with. To do the averages, however, I need to create an array.
>>> For
>>> some reason, NCL will not allow me to do it the way I am trying to:
>>> by
>>> multiplying lat and lon dimensions to get the number of data points
>>> for
>>> wind values at every lat/lon on the NARR grid.
>>>
>>> As always, your assistance is appreciated.
>>>
>>> Cheers
>>>
>>>
>>>
>>>
>>>
>>> Using NCL version 5.1.1
>>> ;*************************************************
>>> ; NARRAVG.ncl
>>> ;
>>> ;************************************************
>>> 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"
>>>
>>> ;************************************************
>>> begin
>>>
>>> ;************************************************
>>> ; read in multiple grib files
>>> ;************************************************
>>> diri = "./" ; directory where files reside
>>> fils = systemfunc("ls "+diri+"narr*.grb")
>>>
>>> setfileoption("grb","SingleElementDimensions","initial_time")
>>>
>>> f = addfiles(fils,"r"); read in all files in the folder
>>>
>>> ;*************************Get Variables to Average**********
>>> U_GRD_221 = addfiles_GetVar(f, fils, "U_GRD_221_HTGL_113") ;Get these
>>> variables from all files
>>> V_GRD_221 = addfiles_GetVar(f, fils, "V_GRD_221_HTGL_113")
>>>
>>> printVarSummary(U_GRD_221)
>>> printVarSummary(V_GRD_221)
>>>
>>> dimU_GRD = dimsizes(U_GRD_221)
>>>
>>> gridx_221 = dimU_GRD(2) ; get dimensions of lat
>>> gridy_221 = dimU_GRD(3)
>>> wind_data = dimU_GRD(2)*dimU_GRD(3) ; multiply lat and lon
>>> dimensions to get wind data points needed
>>>
>>> store = new((/wind_data,gridx_221,gridy_221/), "float") ; create a
>>> new
>>> array for new netCDF file.
>>>
>>>
>>> ;*************************Calculate Average******
>>>
>>> U_GRD_221_AVG = dim_avg_n_Wrap(U_GRD_221(:,0,:,:),0) ; calculate
>>> average of zonal wind at level 0 for all grid points
>>> V_GRD_221_AVG = dim_avg_n_Wrap(V_GRD_221(:,0,:,:),0)
>>>
>>>
>>> wspd = sqrt( U_GRD_221_AVG2 + V_GRD_221_AVG2 ) ; calculate wind
>>> speed
>>> magnitude
>>>
>>>
>>> ;*************create NetCDF output*******************
>>> store!0 = "wspd" ; assign named dimensions
>>> store!1 = "gridx_221"
>>> store!2 = "gridy_221"
>>>
>>> store_at_long_name = "Avg Wind Speed"
>>> store_at_units = "m/s"
>>>
>>> system("/bin/rm -f dailyAvg.WSPD10_NARR.nc") ; remove any pre-
>>> existing
>>> file
>>> ncdf = addfile("dailyAvg.WSPD10_NARR.nc","c")
>>> filedimdef(ncdf,"initial_time0_hours",-1,True)
>>> fileattdef (ncdf, f)
>>>
>>> ncdf->U_GRD_221_HTGL_113 = store
>>>
>>> end
>>>
>>>
>>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
>------------------------------
>
>Message: 5
>Date: Wed, 09 Sep 2009 22:20:46 -0600
>From: Dennis Shea <shea_at_ucar.edu>
>Subject: Re: [ncl-talk] Array Error
>To: DJ Rasmussen <dj.rasmussen_at_ssec.wisc.edu>
>Cc: ncl-talk_at_ucar.edu
>Message-ID: <4AA87E9E.4070503_at_ucar.edu>
>Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
>As noted in the documentation, addfiles_GetVar has been deprecated:
>http://www.ncl.ucar.edu/Document/Functions/Contributed/addfiles_GetVar.shtml
>ie: don't use it
>
>addfiles is more verstile.
>----
>
>If you are trying to compute the temporal mean of the wind speed
>at each grid point:
>
> f = addfiles (fils, "r")
> U_GRD_221 = f[:]->U_GRD_221_HTGL_113(:,0,:,:) ; (time,lat,lon)
> V_GRD_221 = f[:]->V_GRD_221_HTGL_113(:,0,:,:)
> WSPD = sqrt(U_GRD_221^2 + V_GRD_221^2)
> copy_VarMeta(U_GRD_221, WSPD)
> WSPD_at_long_name = "wind speed"
> printVarSummary(WSPD)
>
> WSPD_AVG = dim_avg_n_Wrap(WSPD,0) ; (lat,lon)
> WSPD_AVG_at_long_name = "mean wind speed"
>
> printVarSummary(WSPD_AVG)
>
>
>
>
>
>Dave Allured wrote:
>> store = new((/wind_data,gridx_221,gridy_221/), "float")
>>
>> This line makes a 3-D array of dimensions nx*ny, nx, ny. The first
>> dimension does not seem right to me. Should this be the time
>> dimension or the level dimension instead?
>>
>> 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/
>>
>> DJ Rasmussen wrote:
>>> Hi NCL users,
>>>
>>> I am attempting to make an NCL script that reads in multiple monthly
>>> NARR files with wind data and then averages them all and, finally
>>> outputs the averages to a netCDF file that I can use to make a contour
>>> plot with. To do the averages, however, I need to create an array. For
>>> some reason, NCL will not allow me to do it the way I am trying to: by
>>> multiplying lat and lon dimensions to get the number of data points for
>>> wind values at every lat/lon on the NARR grid.
>>>
>>> As always, your assistance is appreciated.
>>>
>>> Cheers
>>>
>>>
>>>
>>>
>>>
>>> Using NCL version 5.1.1
>>> ;*************************************************
>>> ; NARRAVG.ncl
>>> ;
>>> ;************************************************
>>> 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"
>>>
>>> ;************************************************
>>> begin
>>>
>>> ;************************************************
>>> ; read in multiple grib files
>>> ;************************************************
>>> diri = "./" ; directory where files reside
>>> fils = systemfunc("ls "+diri+"narr*.grb")
>>>
>>> setfileoption("grb","SingleElementDimensions","initial_time")
>>>
>>> f = addfiles(fils,"r"); read in all files in the folder
>>>
>>> ;*************************Get Variables to Average**********
>>> U_GRD_221 = addfiles_GetVar(f, fils, "U_GRD_221_HTGL_113") ;Get these
>>> variables from all files
>>> V_GRD_221 = addfiles_GetVar(f, fils, "V_GRD_221_HTGL_113")
>>>
>>> printVarSummary(U_GRD_221)
>>> printVarSummary(V_GRD_221)
>>>
>>> dimU_GRD = dimsizes(U_GRD_221)
>>>
>>> gridx_221 = dimU_GRD(2) ; get dimensions of lat
>>> gridy_221 = dimU_GRD(3)
>>> wind_data = dimU_GRD(2)*dimU_GRD(3) ; multiply lat and lon
>>> dimensions to get wind data points needed
>>>
>>> store = new((/wind_data,gridx_221,gridy_221/), "float") ; create a new
>>> array for new netCDF file.
>>>
>>>
>>> ;*************************Calculate Average******
>>>
>>> U_GRD_221_AVG = dim_avg_n_Wrap(U_GRD_221(:,0,:,:),0) ; calculate
>>> average of zonal wind at level 0 for all grid points
>>> V_GRD_221_AVG = dim_avg_n_Wrap(V_GRD_221(:,0,:,:),0)
>>>
>>>
>>> wspd = sqrt( U_GRD_221_AVG2 + V_GRD_221_AVG2 ) ; calculate wind speed
>>> magnitude
>>>
>>>
>>> ;*************create NetCDF output*******************
>>> store!0 = "wspd" ; assign named dimensions
>>> store!1 = "gridx_221"
>>> store!2 = "gridy_221"
>>>
>>> store_at_long_name = "Avg Wind Speed"
>>> store_at_units = "m/s"
>>>
>>> system("/bin/rm -f dailyAvg.WSPD10_NARR.nc") ; remove any pre-existing
>>> file
>>> ncdf = addfile("dailyAvg.WSPD10_NARR.nc","c")
>>> filedimdef(ncdf,"initial_time0_hours",-1,True)
>>> fileattdef (ncdf, f)
>>>
>>> ncdf->U_GRD_221_HTGL_113 = store
>>>
>>> end
>>>
>>>
>>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>------------------------------
>
>Message: 6
>Date: Thu, 10 Sep 2009 13:22:47 +0530
>From: "Sabeerali(sebi)" <sabeerl_at_gmail.com>
>Subject: [ncl-talk] Correlation significant test....
>To: ncl-talk_at_ucar.edu
>Message-ID:
> <4792511f0909100052g353f443am7d66a9f56bb52520_at_mail.gmail.com>
>Content-Type: text/plain; charset="iso-8859-1"
>
>Hello,
>
> By using the escorc function I found the correlation between variables
>'a' and 'b'. I want to do its significance test and plot. Which function I
>can use and how I can use it? The dimension of 'a' is lat=73,lon=144,time=23
>and that of 'b' is time=23
>
>
>Thanks in advance
>-------------- next part --------------
>An HTML attachment was scrubbed...
>URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20090910/1d19dccf/attachment.html
>
>------------------------------
>
>Message: 7
>Date: Thu, 10 Sep 2009 09:54:10 -0400
>From: Debasish PaiMazumder <debasish_at_sca.uqam.ca>
>Subject: Re: [ncl-talk] Correlation significant test....
>To: "Sabeerali(sebi)" <sabeerl_at_gmail.com>
>Cc: ncl-talk_at_ucar.edu
>Message-ID: <9A282345-6E5F-479A-9366-8C4F39F9CF9C_at_sca.uqam.ca>
>Content-Type: text/plain; charset="iso-8859-1"
>
>Dear Sabeerali,
>
>If you already calculated correlation coefficient by using "escorc",
>you need to use "rtest" to determine the statistical significance of
>correlation coefficients.
>Please check the web page http://www.ncl.ucar.edu/Document/Functions/
>Built-in/rtest.shtml
>
>To plot significance please check the example conOnCon_4.ncl in web
>page http://www.ncl.ucar.edu/Applications/conOncon.shtml
>
>Good Luck
>
>Debasish
>
>On 09-09-10, at 03:52, Sabeerali(sebi) wrote:
>
>> Hello,
>>
>> By using the escorc function I found the correlation between
>> variables 'a' and 'b'. I want to do its significance test and plot.
>> Which function I can use and how I can use it? The dimension of 'a'
>> is lat=73,lon=144,time=23 and that of 'b' is time=23
>>
>>
>> Thanks in advance
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>Debasish PaiMazumder, PhD
>Hydrology & Climate Science
>Postdoctoral Researcher at UQ?M/Ouranos
>550, Sherbrooke west, 19th floor, west tower,
>Montr?al, QC, Canada H3A 1B9
>Tel: (514) 282-6464-343
> (514) 937-0042 (R)
>email: debasish_at_sca.uqam.ca
>
>
>
>
>-------------- next part --------------
>An HTML attachment was scrubbed...
>URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20090910/accc8fff/attachment.html
>
>------------------------------
>
>_______________________________________________
>ncl-talk mailing list
>ncl-talk_at_ucar.edu
>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>End of ncl-talk Digest, Vol 70, Issue 25
>****************************************

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Sep 10 2009 - 08:23:32 MDT

This archive was generated by hypermail 2.2.0 : Tue Sep 15 2009 - 16:03:21 MDT