<div dir="ltr"><div class="gmail_default" style="font-size:small">It seems that nobody answered this question, so I'll take a stab at it.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Did you notice that "ccr" has a maximum value of 1e20? This indicates that something is likely wrong with a _FillValue attribute somewhere.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You included a bunch of "printVarSummary" output which is helpful, but it would also help to see the min/max of the ts1 and ts2 values. I noticed that each of these variables have different _FillValue attribute values, which makes me a little suspicious. It's almost like the escorc routine is using the same _FillValue for both sets of arrays.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">As a test, you could try setting the _FillValue of ts1 and ts2 to the same value, before calling escorc, to see if this helps. For example:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">ts1@_FillValue = 1e20</div><div class="gmail_default" style="font-size:small"><div class="gmail_default">ts2@_FillValue = 1e20</div><div><div style="font-family:arial,sans-serif;font-size:12.8000001907349px">ccr = escorc(ts2,ts1)</div></div><div><br></div><div>Also, before you do the above, can you include a printMinMax of both ts1 and ts2 to make sure the values look reasonable?</div><div><br></div><div>--Mary</div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 9, 2014 at 1:27 PM, Vanúcia Schumacher <span dir="ltr"><<a href="mailto:vanucia-schumacher@hotmail.com" target="_blank">vanucia-schumacher@hotmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr"><font face="Calibri,sans-serif" color="#000000"><span style="font-size:21px">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?</span><br></font><br>I used this function:<div><br></div><div>... </div><div><br></div><div><div>ccr = escorc(ts2,ts1)</div><div><br></div><div> copy_VarAtts(ts1, ccr) </div><div> copy_VarCoords_1(ts1,ccr)</div><div><br></div><div>....</div><div><br></div><div><div>; statistical significance of ccr</div><div>;*************************************************</div><div><br></div><div> Nx = num(.not.ismissing(ts1))</div><div><br></div><div> prob = rtest(ccr,Nx, 0) <span style="white-space:pre-wrap"> </span></div><div><br></div><div> copy_VarAtts(ccr, prob) </div><div> copy_VarCoords_1(ccr,prob)</div><div><br></div><div>prob!0="lat"</div><div>prob!1="lon"</div><div>prob&lat=ccr&lat</div><div>prob&lon=ccr&lon</div></div><div><br></div><div><br></div><div>Variable: tmp1</div><div>Type: double</div><div>Total Size: 46656000 bytes</div><div> 5832000 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[TIME | 360] x [lat | 90] x [lon | 180]</div><div>Coordinates: </div><div> TIME: [527040..16260480]</div><div> lat: [-89.69999694824219..89.69999694824219]</div><div> lon: [ 0..359.7000122070312]</div><div>Number Of Attributes: 5</div><div> remap :<span style="white-space:pre-wrap"> </span>remapped via ESMF_regrid_with_weights: Bilinear remapping</div><div> missing_value :<span style="white-space:pre-wrap"> </span>-9.999999999999999e+33</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>-9.999999999999999e+33</div><div> long_name :<span style="white-space:pre-wrap"> </span>SST[GX=X2DEG,GY=Y2DEG]</div><div> history :<span style="white-space:pre-wrap"> </span>From sst</div><div><br></div><div>Variable: tmp2</div><div>Type: float</div><div>Total Size: 23328000 bytes</div><div> 5832000 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[time | 360] x [lat | 90] x [lon | 180]</div><div>Coordinates: </div><div> time: [15.5..10934.5]</div><div> lat: [-87.86..87.86]</div><div> lon: [ 0..358]</div><div>Number Of Attributes: 13</div><div> remap :<span style="white-space:pre-wrap"> </span>remapped via ESMF_regrid_with_weights: Bilinear remapping</div><div> associated_files :<span style="white-space:pre-wrap"> </span>baseURL: <a href="http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation" target="_blank">http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation</a> gridspecFile: gridspec_ocean_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc areacello: areacello_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>1e+20</div><div> missing_value :<span style="white-space:pre-wrap"> </span>1e+20</div><div> cell_measures :<span style="white-space:pre-wrap"> </span>area: areacello</div><div> cell_methods :<span style="white-space:pre-wrap"> </span>time: mean (interval: 30 minutes)</div><div> history :<span style="white-space:pre-wrap"> </span>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).</div><div> original_units :<span style="white-space:pre-wrap"> </span>degC</div><div> original_name :<span style="white-space:pre-wrap"> </span>sosstsst</div><div> units :<span style="white-space:pre-wrap"> </span>K</div><div> comment :<span style="white-space:pre-wrap"> </span>"this may differ from ""surface temperature"" in regions of sea ice."</div><div> long_name :<span style="white-space:pre-wrap"> </span>Sea Surface Temperature</div><div> standard_name :<span style="white-space:pre-wrap"> </span>sea_surface_temperature</div><div><br></div><div>(0)<span style="white-space:pre-wrap"> </span><font color="#8c0095">tmp2: min=3.88894 max=305.118</font></div><div><font color="#8c0095">(0)<span style="white-space:pre-wrap"> </span>tmp1: min=231.963 max=304.787</font></div><div><br></div><div>Variable: ts1</div><div>Type: double</div><div>Total Size: 46656000 bytes</div><div> 5832000 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[lat | 90] x [lon | 180] x [TIME | 360]</div><div>Coordinates: </div><div> lat: [-89.69999694824219..89.69999694824219]</div><div> lon: [ 0..359.7000122070312]</div><div> TIME: [527040..16260480]</div><div>Number Of Attributes: 5</div><div> history :<span style="white-space:pre-wrap"> </span>From sst</div><div> long_name :<span style="white-space:pre-wrap"> </span>SST[GX=X2DEG,GY=Y2DEG]</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>-9.999999999999999e+33</div><div> missing_value :<span style="white-space:pre-wrap"> </span>-9.999999999999999e+33</div><div> remap :<span style="white-space:pre-wrap"> </span>remapped via ESMF_regrid_with_weights: Bilinear remapping</div><div><br></div><div>Variable: ts2</div><div>Type: float</div><div>Total Size: 23328000 bytes</div><div> 5832000 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[lat | 90] x [lon | 180] x [time | 360]</div><div>Coordinates: </div><div> lat: [-87.86..87.86]</div><div> lon: [ 0..358]</div><div> time: [15.5..10934.5]</div><div>Number Of Attributes: 13</div><div> standard_name :<span style="white-space:pre-wrap"> </span>sea_surface_temperature</div><div> long_name :<span style="white-space:pre-wrap"> </span>Sea Surface Temperature</div><div> comment :<span style="white-space:pre-wrap"> </span>"this may differ from ""surface temperature"" in regions of sea ice."</div><div> units :<span style="white-space:pre-wrap"> </span>K</div><div> original_name :<span style="white-space:pre-wrap"> </span>sosstsst</div><div> original_units :<span style="white-space:pre-wrap"> </span>degC</div><div> history :<span style="white-space:pre-wrap"> </span>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).</div><div> cell_methods :<span style="white-space:pre-wrap"> </span>time: mean (interval: 30 minutes)</div><div> cell_measures :<span style="white-space:pre-wrap"> </span>area: areacello</div><div> missing_value :<span style="white-space:pre-wrap"> </span>1e+20</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>1e+20</div><div> associated_files :<span style="white-space:pre-wrap"> </span>baseURL: <a href="http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation" target="_blank">http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation</a> gridspecFile: gridspec_ocean_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc areacello: areacello_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc</div><div> remap :<span style="white-space:pre-wrap"> </span>remapped via ESMF_regrid_with_weights: Bilinear remapping</div><div><br></div><div>Variable: ccr</div><div>Type: double</div><div>Total Size: 129600 bytes</div><div> 16200 values</div><div>Number of Dimensions: 2</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[lat | 90] x [lon | 180]</div><div>Coordinates: </div><div> lat: [-89.69999694824219..89.69999694824219]</div><div> lon: [ 0..359.7000122070312]</div><div>Number Of Attributes: 5</div><div> remap :<span style="white-space:pre-wrap"> </span>remapped via ESMF_regrid_with_weights: Bilinear remapping</div><div> missing_value :<span style="white-space:pre-wrap"> </span>-9.999999999999999e+33</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>-9.999999999999999e+33</div><div> long_name :<span style="white-space:pre-wrap"> </span>SST[GX=X2DEG,GY=Y2DEG]</div><div> history :<span style="white-space:pre-wrap"> </span>From sst</div><div><font color="#8c0095">(0)<span style="white-space:pre-wrap"> </span>ccr: min=-0.35971 max=1e+20</font></div><div><br></div><div>Variable: Nx</div><div>Type: integer</div><div>Total Size: 4 bytes</div><div> 1 values</div><div>Number of Dimensions: 1</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[1]</div><div>Coordinates: </div><div><br></div><div>Variable: prob</div><div>Type: double</div><div>Total Size: 129600 bytes</div><div> 16200 values</div><div>Number of Dimensions: 2</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[lat | 90] x [lon | 180]</div><div>Coordinates: </div><div> lat: [-89.69999694824219..89.69999694824219]</div><div> lon: [ 0..359.7000122070312]</div><div>Number Of Attributes: 5</div><div> remap :<span style="white-space:pre-wrap"> </span>remapped via ESMF_regrid_with_weights: Bilinear remapping</div><div> missing_value :<span style="white-space:pre-wrap"> </span>-9.999999999999999e+33</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>-9.999999999999999e+33</div><div> long_name :<span style="white-space:pre-wrap"> </span>SST[GX=X2DEG,GY=Y2DEG]</div><div> history :<span style="white-space:pre-wrap"> </span>From sst</div><div><font color="#8c0095">(0)<span style="white-space:pre-wrap"> </span>prob: min=0 max=0.862573</font></div><div><br></div><div><font face="Comic Sans MS"><span style="line-height:17px">---</span></font></div><span style="line-height:17px;background-color:rgb(255,255,255)"><font size="3" face="Comic Sans MS" style="font-size:12pt"><font>Vanúcia Schumacher</font></font></span><div><span style="line-height:17px;background-color:rgb(255,255,255)"><font face="Comic Sans MS">Mestranda em Meteorologia - UFV<font size="3" style="font-size:12pt"><font><br style="line-height:17px">Meteorologista -UFPel</font><br></font></font></span></div><div><span style="line-height:17px;background-color:rgb(255,255,255)"><font face="Comic Sans MS"><font size="3" style="font-size:12pt"><font>Departamento de Meteorologia Agrícola - DEA</font></font></font></span></div><div><span style="line-height:17px;background-color:rgb(255,255,255)"><font size="3" face="Comic Sans MS" style="font-size:12pt"><font>Cel: (31) 9978 2522 </font></font></span></div><div><span style="line-height:17px;background-color:rgb(255,255,255)"><font size="3" face="Comic Sans MS" style="font-size:12pt"><font>DEA: (31) 3899 1890</font></font></span></div></div> </div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>