[ncl-talk] problem with continents projections
Vanúcia Schumacher
vanucia-schumacher at hotmail.com
Thu Oct 9 12:25:50 MDT 2014
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141009/8c0b43c2/attachment.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: anomalia-cfsr-podac.png
Type: image/png
Size: 184255 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141009/8c0b43c2/attachment.png
More information about the ncl-talk
mailing list