<div dir="ltr"><div dir="ltr"><div><font color="#000000" face="Calibri, sans-serif"><span style="font-size:21px">Forgot to cc ncl-talk. </span></font></div><div><font color="#000000" face="Calibri, sans-serif"><span style="font-size:21px"><br></span></font></div><div><font color="#000000" face="Calibri, sans-serif"><span style="font-size:21px">Your data are not on the same grid, so when you subtract you subtracted opposite sides of the globe from one another, the basic math statements don't pay attention to coordinates. See the extracted lines from your printVarSummary output. </span></font></div><div><font color="#000000" face="Calibri, sans-serif"><span style="font-size:21px">The addition of <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/lonFlip.shtml" target="_blank">LonFlip</a> will result in media2 being -180:178 like media. The rest should then fall into place. </span></font></div><div><font color="#000000" face="Calibri, sans-serif"><span style="font-size:21px"><br></span></font></div><div><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:21px">media2= LonFlip( dim_avg_n_Wrap(tnew, 0) )</span><br></div><div><br></div><div><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:21px"><br></span></div><div><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:21px">Variable: media</span><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:21px"><br></span></div><div><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:21px"> LON: [-180..178]</span><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:21px"><br></span></div><div><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:21px"><br></span></div><div><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:21px">Variable: media2</span><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:21px"><br></span></div><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:21px"> lon: [ 0..359.7000122070312]</span><br></div><div class="gmail_extra"><br style="font-family:arial,sans-serif;font-size:13px"><div class="gmail_quote" style="font-family:arial,sans-serif;font-size:13px"><div class="adm"><div id="q_148f647a8bc3610a_1" class="h4"></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 9, 2014 at 2:25 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"><div><span style="font-size:21px">Dear all,</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">My graphic output (attached) is having problems with the projections of the continents, can someone tell me how to solve this problem, when I plot each file data separately, this problem does not occur, only when calculating the anomaly between them.</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">My script:</span></div><div><span style="font-size:21px"><br></span></div><div><div><span style="font-size:21px">;========== open model datas ====================</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">fils =("<a href="http://tos_90x180.podac.nc" target="_blank">tos_90x180.podac.nc</a>")</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">f = addfile(fils, "r") </span></div><div><span style="font-size:21px"> </span></div><div><span style="font-size:21px">t =f->bsst</span></div><div><span style="font-size:21px">t1 = t + 273.15</span></div><div><span style="font-size:21px">copy_VarCoords(t, t1)</span></div><div><br></div><div><span style="font-size:21px">media= dim_avg_n_Wrap(t1, 0) </span></div><div><br></div><div><span style="font-size:21px">;========================= open reanalises datas =====================</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">fils2 =("<a href="http://tos_90x180.cfsr.nc" target="_blank">tos_90x180.cfsr.nc</a>")</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">f1 = addfile(fils2, "r") </span></div><div><span style="font-size:21px"> </span></div><div><span style="font-size:21px">t2 =f1->tos</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">t2@_FillValue = -9.999999999999999e+33</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">t2=where(t2.ge.271, t2, t2@_FillValue)</span></div><div><br></div><div><span style="font-size:21px">TIME= f1->TIME</span></div><div><span style="font-size:21px">YYYY= cd_calendar(TIME,-1)/100</span></div><div><br></div><div><span style="font-size:21px">time_start=ut_inv_calendar(1985,01,01,00,00,00,"minutes since 1980-01-01 00:00:00",0)</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">time_end=ut_inv_calendar(2009,12,31,00,00,00,"minutes since 1980-01-01 00:00:00",0)</span></div><div><br></div><div><span style="font-size:21px">tnew=t2({time_start:time_end},:,:)</span></div><div><br></div><div><span style="font-size:21px">media2= dim_avg_n_Wrap(tnew, 0) </span></div><div><br></div><div><span style="font-size:21px">;=============================== anomalia ===========================</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">anomalia= media - media2</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">copy_VarCoords(media2, anomalia)</span></div><div><br></div><div>....</div><div><br></div><div><span style="font-size:21px"> plot = gsn_csm_contour_map_ce(wks,anomalia,res)</span></div><div><span style="font-size:21px"> end</span></div></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">Variable: t1</span></div><div><span style="font-size:21px">Type: float</span></div><div><span style="font-size:21px">Total Size: 19440000 bytes</span></div><div><span style="font-size:21px"> 4860000 values</span></div><div><span style="font-size:21px">Number of Dimensions: 3</span></div><div><span style="font-size:21px">Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[time | 300] x [LAT | 90] x [LON | 180]</span></div><div><span style="font-size:21px">Coordinates: </span></div><div><span style="font-size:21px"> time: [ 0.. 0]</span></div><div><span style="font-size:21px"> LAT: [-89..89]</span></div><div><span style="font-size:21px"> LON: [-180..178]</span></div><div><span style="font-size:21px">Number Of Attributes: 1</span></div><div><span style="font-size:21px"> _FillValue :<span style="white-space:pre-wrap"> </span>9.96921e+36</span></div><div><span style="font-size:21px">(0)<span style="white-space:pre-wrap"> </span>t1: min=270.15 max=308.121</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">Variable: media</span></div><div><span style="font-size:21px">Type: float</span></div><div><span style="font-size:21px">Total Size: 64800 bytes</span></div><div><span style="font-size:21px"> 16200 values</span></div><div><span style="font-size:21px">Number of Dimensions: 2</span></div><div><span style="font-size:21px">Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[LAT | 90] x [LON | 180]</span></div><div><span style="font-size:21px">Coordinates: </span></div><div><span style="font-size:21px"> LAT: [-89..89]</span></div><div><span style="font-size:21px"> LON: [-180..178]</span></div><div><span style="font-size:21px">Number Of Attributes: 2</span></div><div><span style="font-size:21px"> _FillValue :<span style="white-space:pre-wrap"> </span>9.96921e+36</span></div><div><span style="font-size:21px"> average_op_ncl :<span style="white-space:pre-wrap"> </span>dim_avg_n over dimension(s): time</span></div><div><span style="font-size:21px">(0)<span style="white-space:pre-wrap"> </span>media: min=270.15 max=302.783</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">Variable: tnew</span></div><div><span style="font-size:21px">Type: double</span></div><div><span style="font-size:21px">Total Size: 38880000 bytes</span></div><div><span style="font-size:21px"> 4860000 values</span></div><div><span style="font-size:21px">Number of Dimensions: 3</span></div><div><span style="font-size:21px">Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[TIME | 300] x [lat | 90] x [lon | 180]</span></div><div><span style="font-size:21px">Coordinates: </span></div><div><span style="font-size:21px"> TIME: [2630880..15734880]</span></div><div><span style="font-size:21px"> lat: [-89.69999694824219..89.69999694824219]</span></div><div><span style="font-size:21px"> lon: [ 0..359.7000122070312]</span></div><div><span style="font-size:21px">Number Of Attributes: 5</span></div><div><span style="font-size:21px"> history :<span style="white-space:pre-wrap"> </span>From sst</span></div><div><span style="font-size:21px"> long_name :<span style="white-space:pre-wrap"> </span>SST[GX=X2DEG,GY=Y2DEG]</span></div><div><span style="font-size:21px"> _FillValue :<span style="white-space:pre-wrap"> </span>-9.999999790214768e+33</span></div><div><span style="font-size:21px"> missing_value :<span style="white-space:pre-wrap"> </span>-9.999999999999999e+33</span></div><div><span style="font-size:21px"> remap :<span style="white-space:pre-wrap"> </span>remapped via ESMF_regrid_with_weights: Bilinear remapping</span></div><div><span style="font-size:21px">(0)<span style="white-space:pre-wrap"> </span>tnew: min=271 max=304.748</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">Variable: media2</span></div><div><span style="font-size:21px">Type: double</span></div><div><span style="font-size:21px">Total Size: 129600 bytes</span></div><div><span style="font-size:21px"> 16200 values</span></div><div><span style="font-size:21px">Number of Dimensions: 2</span></div><div><span style="font-size:21px">Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[lat | 90] x [lon | 180]</span></div><div><span style="font-size:21px">Coordinates: </span></div><div><span style="font-size:21px"> lat: [-89.69999694824219..89.69999694824219]</span></div><div><span style="font-size:21px"> lon: [ 0..359.7000122070312]</span></div><div><span style="font-size:21px">Number Of Attributes: 6</span></div><div><span style="font-size:21px"> history :<span style="white-space:pre-wrap"> </span>From sst</span></div><div><span style="font-size:21px"> long_name :<span style="white-space:pre-wrap"> </span>SST[GX=X2DEG,GY=Y2DEG]</span></div><div><span style="font-size:21px"> _FillValue :<span style="white-space:pre-wrap"> </span>-9.999999790214768e+33</span></div><div><span style="font-size:21px"> missing_value :<span style="white-space:pre-wrap"> </span>-9.999999999999999e+33</span></div><div><span style="font-size:21px"> remap :<span style="white-space:pre-wrap"> </span>remapped via ESMF_regrid_with_weights: Bilinear remapping</span></div><div><span style="font-size:21px"> average_op_ncl :<span style="white-space:pre-wrap"> </span>dim_avg_n over dimension(s): TIME</span></div><div><span style="font-size:21px">(0)<span style="white-space:pre-wrap"> </span>media2: min=271.224 max=302.745</span></div><div><span style="font-size:21px"><br></span></div><div><span style="font-size:21px">Variable: anomalia</span></div><div><span style="font-size:21px">Type: double</span></div><div><span style="font-size:21px">Total Size: 129600 bytes</span></div><div><span style="font-size:21px"> 16200 values</span></div><div><span style="font-size:21px">Number of Dimensions: 2</span></div><div><span style="font-size:21px">Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[lat | 90] x [lon | 180]</span></div><div><span style="font-size:21px">Coordinates: </span></div><div><span style="font-size:21px"> lat: [-89.69999694824219..89.69999694824219]</span></div><div><span style="font-size:21px"> lon: [ 0..359.7000122070312]</span></div><div><span style="font-size:21px">Number Of Attributes: 1</span></div><div><span style="font-size:21px"> _FillValue :<span style="white-space:pre-wrap"> </span>9.969209968386869e+36</span></div><div><span style="font-size:21px">(0)<span style="white-space:pre-wrap"> </span>anomalia: min=-32.37 max=10.5374</span></div><div style="font-size:16pt"><br></div></font> </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>