[ncl-talk] problem with continents projections

[David Brown]

It's actually spelled lonFlip.
 -dave


On Thu, Oct 9, 2014 at 12:46 PM, Vanúcia Schumacher
<vanucia-schumacher at hotmail.com> wrote:
> How do I set this variable?
>
> Undefined identifier: (LonFlip) is undefined, can't continue
>
> Thanks,
>
>
>
> ________________________________
> Date: Thu, 9 Oct 2014 14:37:30 -0400
> Subject: Re: [ncl-talk] problem with continents projections
> From: abrammer at albany.edu
> To: vanucia-schumacher at hotmail.com
> CC: ncl-talk at ucar.edu
>
>
> 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 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
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>