[ncl-talk] problem with continents projections

[Vanúcia Schumacher]

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->bsstt1 = t + 273.15copy_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->TIMEYYYY= 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: t1Type: floatTotal Size: 19440000 bytes            4860000 valuesNumber of Dimensions: 3Dimensions 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: mediaType: floatTotal Size: 64800 bytes            16200 valuesNumber of Dimensions: 2Dimensions 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: tnewType: doubleTotal Size: 38880000 bytes            4860000 valuesNumber of Dimensions: 3Dimensions 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: media2Type: 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: 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: anomaliaType: 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: 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/e7c554ab/attachment.html