[ncl-talk] problem with continents projections

[Alan Brammer]

Forgot to cc ncl-talk.

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.
The addition of LonFlip
<https://www.ncl.ucar.edu/Document/Functions/Contributed/lonFlip.shtml> will
result in media2 being -180:178 like media.  The rest should then fall into
place.

media2= LonFlip( dim_avg_n_Wrap(tnew, 0)  )


Variable: media
 LON: [-180..178]

Variable: media2
 lon: [   0..359.7000122070312]


On Thu, Oct 9, 2014 at 2:25 PM, Vanúcia Schumacher <
vanucia-schumacher at hotmail.com> wrote:

> Dear all,
>
> 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.
>
> My script:
>
> ;========== open model datas ====================
>
> fils =("tos_90x180.podac.nc")
>
> f  = addfile(fils, "r")
>
> t =f->bsst
> t1 = t + 273.15
> copy_VarCoords(t, t1)
>
> media= dim_avg_n_Wrap(t1, 0)
>
> ;========================= open reanalises datas =====================
>
> fils2 =("tos_90x180.cfsr.nc")
>
> f1  = addfile(fils2, "r")
>
> t2 =f1->tos
>
> t2 at _FillValue = -9.999999999999999e+33
>
> t2=where(t2.ge.271, t2, t2 at _FillValue)
>
> TIME= f1->TIME
> YYYY= cd_calendar(TIME,-1)/100
>
> time_start=ut_inv_calendar(1985,01,01,00,00,00,"minutes since 1980-01-01
> 00:00:00",0)
>
> time_end=ut_inv_calendar(2009,12,31,00,00,00,"minutes since 1980-01-01
> 00:00:00",0)
>
> tnew=t2({time_start:time_end},:,:)
>
> media2= dim_avg_n_Wrap(tnew, 0)
>
> ;=============================== anomalia ===========================
>
> anomalia= media - media2
>
> copy_VarCoords(media2, anomalia)
>
> ....
>
>   plot = gsn_csm_contour_map_ce(wks,anomalia,res)
>   end
>
> Variable: t1
> Type: float
> Total Size: 19440000 bytes
>             4860000 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 300] x [LAT | 90] x [LON | 180]
> Coordinates:
>             time: [   0..   0]
>             LAT: [-89..89]
>             LON: [-180..178]
> Number Of Attributes: 1
>   _FillValue : 9.96921e+36
> (0) t1: min=270.15 max=308.121
>
> Variable: media
> Type: float
> Total Size: 64800 bytes
>             16200 values
> Number of Dimensions: 2
> Dimensions and sizes: [LAT | 90] x [LON | 180]
> Coordinates:
>             LAT: [-89..89]
>             LON: [-180..178]
> Number Of Attributes: 2
>   _FillValue : 9.96921e+36
>   average_op_ncl : dim_avg_n over dimension(s): time
> (0) media: min=270.15 max=302.783
>
> Variable: tnew
> Type: double
> Total Size: 38880000 bytes
>             4860000 values
> Number of Dimensions: 3
> Dimensions and sizes: [TIME | 300] x [lat | 90] x [lon | 180]
> Coordinates:
>             TIME: [2630880..15734880]
>             lat: [-89.69999694824219..89.69999694824219]
>             lon: [   0..359.7000122070312]
> Number Of Attributes: 5
>   history : From sst
>   long_name : SST[GX=X2DEG,GY=Y2DEG]
>   _FillValue : -9.999999790214768e+33
>   missing_value : -9.999999999999999e+33
>   remap : remapped via ESMF_regrid_with_weights: Bilinear remapping
> (0) tnew: min=271 max=304.748
>
> Variable: media2
> Type: double
> Total Size: 129600 bytes
>             16200 values
> Number of Dimensions: 2
> Dimensions and sizes: [lat | 90] x [lon | 180]
> Coordinates:
>             lat: [-89.69999694824219..89.69999694824219]
>             lon: [   0..359.7000122070312]
> Number Of Attributes: 6
>   history : From sst
>   long_name : SST[GX=X2DEG,GY=Y2DEG]
>   _FillValue : -9.999999790214768e+33
>   missing_value : -9.999999999999999e+33
>   remap : remapped via ESMF_regrid_with_weights: Bilinear remapping
>   average_op_ncl : dim_avg_n over dimension(s): TIME
> (0) media2: min=271.224 max=302.745
>
> Variable: anomalia
> Type: double
> Total Size: 129600 bytes
>             16200 values
> Number of Dimensions: 2
> Dimensions and sizes: [lat | 90] x [lon | 180]
> Coordinates:
>             lat: [-89.69999694824219..89.69999694824219]
>             lon: [   0..359.7000122070312]
> Number Of Attributes: 1
>   _FillValue : 9.969209968386869e+36
> (0) anomalia: min=-32.37 max=10.5374
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141009/38c1d6ca/attachment-0001.html