[ncl-talk] problems with rtest

[Vanúcia Schumacher]

I calculated the correlation between the variable sst to 2 different models, with same size lat, lon and time for both, and now I want to plot the test of significance, using the rtest. But my plot is very strange, does not plot anything (attached), someone can help?
I used this function:
... 
ccr    =   escorc(ts2,ts1)
  copy_VarAtts(ts1, ccr)                      copy_VarCoords_1(ts1,ccr)
....
;  statistical significance of ccr;*************************************************
  Nx    = num(.not.ismissing(ts1))
  prob  = rtest(ccr,Nx, 0) 	
  copy_VarAtts(ccr, prob)                      copy_VarCoords_1(ccr,prob)
prob!0="lat"prob!1="lon"prob&lat=ccr&latprob&lon=ccr&lon

Variable: tmp1Type: doubleTotal Size: 46656000 bytes            5832000 valuesNumber of Dimensions: 3Dimensions and sizes:	[TIME | 360] x [lat | 90] x [lon | 180]Coordinates:             TIME: [527040..16260480]            lat: [-89.69999694824219..89.69999694824219]            lon: [   0..359.7000122070312]Number Of Attributes: 5  remap :	remapped via ESMF_regrid_with_weights: Bilinear remapping  missing_value :	-9.999999999999999e+33  _FillValue :	-9.999999999999999e+33  long_name :	SST[GX=X2DEG,GY=Y2DEG]  history :	From sst
Variable: tmp2Type: floatTotal Size: 23328000 bytes            5832000 valuesNumber of Dimensions: 3Dimensions and sizes:	[time | 360] x [lat | 90] x [lon | 180]Coordinates:             time: [15.5..10934.5]            lat: [-87.86..87.86]            lon: [ 0..358]Number Of Attributes: 13  remap :	remapped via ESMF_regrid_with_weights: Bilinear remapping  associated_files :	baseURL: http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile: gridspec_ocean_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc areacello: areacello_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc  _FillValue :	1e+20  missing_value :	1e+20  cell_measures :	area: areacello  cell_methods :	time: mean (interval: 30 minutes)  history :	2011-11-09T14:11:22Z altered by CMOR: Converted units from 'degC' to 'K'. 2011-11-09T14:11:22Z altered by CMOR: replaced missing value flag (9.96921e+36) with standard missing value (1e+20).  original_units :	degC  original_name :	sosstsst  units :	K  comment :	"this may differ from ""surface temperature"" in regions of sea ice."  long_name :	Sea Surface Temperature  standard_name :	sea_surface_temperature
(0)	tmp2: min=3.88894 max=305.118(0)	tmp1: min=231.963 max=304.787
Variable: ts1Type: doubleTotal Size: 46656000 bytes            5832000 valuesNumber of Dimensions: 3Dimensions and sizes:	[lat | 90] x [lon | 180] x [TIME | 360]Coordinates:             lat: [-89.69999694824219..89.69999694824219]            lon: [   0..359.7000122070312]            TIME: [527040..16260480]Number Of Attributes: 5  history :	From sst  long_name :	SST[GX=X2DEG,GY=Y2DEG]  _FillValue :	-9.999999999999999e+33  missing_value :	-9.999999999999999e+33  remap :	remapped via ESMF_regrid_with_weights: Bilinear remapping
Variable: ts2Type: floatTotal Size: 23328000 bytes            5832000 valuesNumber of Dimensions: 3Dimensions and sizes:	[lat | 90] x [lon | 180] x [time | 360]Coordinates:             lat: [-87.86..87.86]            lon: [ 0..358]            time: [15.5..10934.5]Number Of Attributes: 13  standard_name :	sea_surface_temperature  long_name :	Sea Surface Temperature  comment :	"this may differ from ""surface temperature"" in regions of sea ice."  units :	K  original_name :	sosstsst  original_units :	degC  history :	2011-11-09T14:11:22Z altered by CMOR: Converted units from 'degC' to 'K'. 2011-11-09T14:11:22Z altered by CMOR: replaced missing value flag (9.96921e+36) with standard missing value (1e+20).  cell_methods :	time: mean (interval: 30 minutes)  cell_measures :	area: areacello  missing_value :	1e+20  _FillValue :	1e+20  associated_files :	baseURL: http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile: gridspec_ocean_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc areacello: areacello_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc  remap :	remapped via ESMF_regrid_with_weights: Bilinear remapping
Variable: ccrType: doubleTotal Size: 129600 bytes            16200 valuesNumber of Dimensions: 2Dimensions and sizes:	[lat | 90] x [lon | 180]Coordinates:             lat: [-89.69999694824219..89.69999694824219]            lon: [   0..359.7000122070312]Number Of Attributes: 5  remap :	remapped via ESMF_regrid_with_weights: Bilinear remapping  missing_value :	-9.999999999999999e+33  _FillValue :	-9.999999999999999e+33  long_name :	SST[GX=X2DEG,GY=Y2DEG]  history :	From sst(0)	ccr: min=-0.35971 max=1e+20
Variable: NxType: integerTotal Size: 4 bytes            1 valuesNumber of Dimensions: 1Dimensions and sizes:	[1]Coordinates: 
Variable: probType: doubleTotal Size: 129600 bytes            16200 valuesNumber of Dimensions: 2Dimensions and sizes:	[lat | 90] x [lon | 180]Coordinates:             lat: [-89.69999694824219..89.69999694824219]            lon: [   0..359.7000122070312]Number Of Attributes: 5  remap :	remapped via ESMF_regrid_with_weights: Bilinear remapping  missing_value :	-9.999999999999999e+33  _FillValue :	-9.999999999999999e+33  long_name :	SST[GX=X2DEG,GY=Y2DEG]  history :	From sst(0)	prob: min=0 max=0.862573
---Vanúcia SchumacherMestranda em Meteorologia - UFVMeteorologista -UFPel
Departamento de Meteorologia Agrícola - DEACel: (31) 9978 2522 DEA: (31) 3899 1890 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140909/c31463da/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cfsr-miroc.000002.png
Type: image/png
Size: 29406 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140909/c31463da/attachment.png