[ncl-talk] problems with regriding

Vanúcia Schumacher vanucia-schumacher at hotmail.com
Fri Aug 22 14:50:37 MDT 2014


Hi Mary,
I tried with the ESM and the following error appears:

Variable: sitType: doubleTotal Size: 78099840 bytes            9762480 valuesNumber of Dimensions: 3Dimensions and sizes:	[TIME | 360] x [J | 149] x [I | 182]Coordinates:             TIME: [76.5..11002.5]            J: [1..149]            I: [1..182]Number Of Attributes: 4  missing_value :	-9.999999999999999e+33  _FillValue :	-9.999999999999999e+33  long_name :	TOS[L=3:122]  history :	From tos_Omon_CMCC-CM_decadal1980_r1i1p1_198011-199012warning:Non-monotonic coordinate array generated -- check validity of the aggregated (leftmost) dimensionwarning:Non-monotonic coordinate array generated -- check validity of the aggregated (leftmost) dimension
Variable: sitType: doubleTotal Size: 78099840 bytes            9762480 valuesNumber of Dimensions: 3Dimensions and sizes:	[time | 360] x [lat | 149] x [lon | 182]Coordinates:             time: [76.5..11002.5]            lat: [1..149]            lon: [1..182]Number Of Attributes: 6  lon2d :	<ARRAY of 546 elements>  lat2d :	<ARRAY of 447 elements>  missing_value :	-9.999999999999999e+33  _FillValue :	-9.999999999999999e+33  long_name :	TOS[L=3:122]  history :	From tos_Omon_CMCC-CM_decadal1980_r1i1p1_198011-199012
warning:Non-monotonic coordinate array generated -- check validity of the aggregated (leftmost) dimensionwarning:Non-monotonic coordinate array generated -- check validity of the aggregated (leftmost) dimension(0)	get_src_grid_info: source lat dims = (447)(0)	get_src_grid_info: source lon dims = (546)(0)	get_src_grid_info: SrcGridType and/or SrcGridLat/SrcGridLon were not set.(0)	          Cannot determine the source grid type.


---Vanúcia SchumacherMestranda em Meteorologia - UFVMeteorologista -UFPel
Departamento de Meteorologia Agrícola - DEACel: (31) 9978 2522 DEA: (31) 3899 1890

Date: Fri, 22 Aug 2014 14:29:44 -0600
Subject: Re: [ncl-talk] problems with regriding
From: haley at ucar.edu
To: vanucia-schumacher at hotmail.com

D, 
My guess is she's missing:

tos at lat2d = f->lattos at lon2d = f->lon

See:
http://www.ncl.ucar.edu/Applications/Scripts/ESMF_regrid_7.ncl


--M




On Fri, Aug 22, 2014 at 6:33 AM, Vanúcia Schumacher <vanucia-schumacher at hotmail.com> wrote:




Hi Dennis,
I still do not understand what can I do to solve this problem, because when I change the resolution, I choose the option of lat and lon (0,358),(-87.86, 87.86), ncview and ferret programs  I can plot the data without problems (global), but when I try plot the data without changing the resolution they are confused on the map (attached).

print(var&lon)print(var&lat)

Variable: lat (coordinate)Type: integer
Total Size: 596 bytes            149 valuesNumber of Dimensions: 1Dimensions and sizes:	[lat | 149]
Coordinates: Number Of Attributes: 4  units :	1
  long_name :	cell index along second dimension  point_spacing :	even
  axis :	Y
ncdump -h 

	time = UNLIMITED ; // (122 currently)
	j = 149 ;	i = 182 ;	bnds = 2 ;
	vertices = 4 ;variables:	double time(time) ;
		time:bounds = "time_bnds" ;		time:units = "days since 1980-11-01" ;
		time:calendar = "standard" ;		time:axis = "T" ;
		time:long_name = "time" ;		time:standard_name = "time" ;
	double time_bnds(time, bnds) ;	int j(j) ;
		j:units = "1" ;		j:long_name = "cell index along second dimension" ;
	int i(i) ;		i:units = "1" ;
		i:long_name = "cell index along first dimension" ;	float lat(j, i) ;
		lat:standard_name = "latitude" ;		lat:long_name = "latitude coordinate" ;
		lat:units = "degrees_north" ;		lat:bounds = "lat_vertices" ;
	float lon(j, i) ;		lon:standard_name = "longitude" ;
		lon:long_name = "longitude coordinate" ;		lon:units = "degrees_east" ;
		lon:bounds = "lon_vertices" ;	float lat_vertices(j, i, vertices) ;
		lat_vertices:units = "degrees_north" ;	float lon_vertices(j, i, vertices) ;
		lon_vertices:units = "degrees_east" ;	float tos(time, j, i) ;
		tos:standard_name = "sea_surface_temperature" ;		tos:long_name = "Sea Surface Temperature" ;
		tos:comment = "\"this may differ from \"\"surface temperature\"\" in regions of sea ice.\"" ;
		tos:units = "K" ;		tos:original_name = "sosstsst" ;
		tos:original_units = "degC" ;		tos:history = "2012-07-24T13:46:07Z altered by CMOR: Converted units from \'degC\' to \'K\'." ;
		tos:cell_methods = "time: mean (interval: 1 month)" ;		tos:cell_measures = "area: areacello" ;
		tos:missing_value = 1.e+20f ;		tos:_FillValue = 1.e+20f ;
		tos:associated_files = "baseURL: http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile: gridspec_ocean_fx_CMCC-CM_decadal1980_r0i0p0.nc areacello: areacello_fx_CMCC-CM_decadal1980_r0i0p0.nc" ;
		tos:coordinates = "lat lon" ;
// global attributes:
		:institution = "CMCC - Centro Euro-Mediterraneo per i Cambiamenti Climatici, Bologna, Italy" ;		:institute_id = "CMCC" ;
		:experiment_id = "decadal1980" ;		:source = "CMCC-CM" ;
		:model_id = "CMCC-CM" ;		:forcing = "Nat,Ant,GHG,SA,TO,Sl" ;
		:parent_experiment_id = "N/A" ;		:parent_experiment_rip = "N/A" ;
		:branch_time = 0. ;		:contact = "Silvio Gualdi (gualdi at bo.ingv.it)" ;
		:history = "Model output postprocessed with  CDO (https://code.zmaw.de/projects) 2012-07-24T13:46:07Z CMOR rewrote data to comply with CF standards and CMIP5 requirements." ;
		:comment = "Decadal hindcasts/predictions, some extended to 30 years" ;		:references = "model described in the documentation at http://www.cmcc.it/data-models/models" ;
		:initialization_method = 1 ;		:physics_version = 1 ;
		:tracking_id = "105250cc-615d-478c-b242-b72f4735fb74" ;		:product = "output" ;
		:experiment = "10- or 30-year run initialized in year 1980" ;		:frequency = "mon" ;
		:creation_date = "2012-07-24T13:46:07Z" ;		:Conventions = "CF-1.4" ;
		:project_id = "CMIP5" ;		:table_id = "Table Omon (27 April 2011) 340eddd4fd838d90fa9ffe1345ecbd73" ;
		:title = "CMCC-CM model output prepared for CMIP5 10- or 30-year run initialized in year 1980" ;		:parent_experiment = "N/A" ;
		:modeling_realm = "ocean" ;		:realization = 1 ;
		:cmor_version = "2.7.1" ;


---
Vanúcia SchumacherMestranda em Meteorologia - UFV
Meteorologista -UFPel
Departamento de Meteorologia Agrícola - DEA
Cel: (31) 9978 2522 DEA: (31) 3899 1890


Date: Thu, 21 Aug 2014 15:12:20 -0600
Subject: Re: [ncl-talk] problems with regriding
From: shea at ucar.edu
To: vanucia-schumacher at hotmail.com

CC: ncl-talk at ucar.edu

[1]
I would say it is highly unlikely that "NCL is not correctly plotting the data"!


NCL is plotting what you sent it. 

[2]
Did you examine the output from


printVarSummary(var)

[snip]



Dimensions and sizes:	[TIME | 360] x [lat | 149] x [lon | 182]Coordinates:             TIME: [76.5..3636.5]

            lat: [1..149]                  <================
            lon: [1..182]                <====================


The coordinates associated with 'var' are numbers spanning 1-149 and 1-182
These are *not* latitudes or longitudes. So, when you  use

new_var = linint2_Wrap(var&lon,var&lat,var,True,new_lon,new_lat,0)


You could do you own 'debugging by printing the input

print(var&lon)
print(var&lat)



Punch line ... you are inputting number 1 to 149 and 1 to 182
NCL is interpolating to those values.



Use 

ncl_filedump

or



ncdump  -h

on your source file. That should have some lat.lon info







On Thu, Aug 21, 2014 at 2:12 PM, Vanúcia Schumacher <vanucia-schumacher at hotmail.com> wrote:





Hi users,
I do regridding data for resolution of 2 degrees and the NCL is not correctly plotting the data, it is cutting the data to the new resolution as opposed to a new lat and lon.


>begin     

>var = f[:]->SST
>var!1="lat">var!2="lon"

>j=var!1>i=var!2   
newfile = "tos_90x180.nc"


file_create = addfile(newfile, "c")    
new_lon = fspan(0,358, 180)     

new_lon at units = "degrees_east"new_lon!0 = "lon"  
new_lat = fspan(-87.86, 87.86, 90)

new_lat at units = "degrees_north"new_lat!0 = "lat"        


new_var = linint2_Wrap(var&lon,var&lat,var,True,new_lon,new_lat,0)
file_create->tos = new_var



Variable: var

Type: doubleTotal Size: 78099840 bytes            9762480 valuesNumber of Dimensions: 3

Dimensions and sizes:	[TIME | 360] x [lat | 149] x [lon | 182]Coordinates:             TIME: [76.5..3636.5]

            lat: [1..149]            lon: [1..182]Number Of Attributes: 4  missing_value :	-9.999999999999999e+33

  _FillValue :	-9.999999999999999e+33  long_name :	TOS[L=3:122]

  history :	From tos_Omon_CMCC-CM_decadal1980_r1i1p1_198011-199012
Variable: new_var

Type: doubleTotal Size: 46656000 bytes            5832000 valuesNumber of Dimensions: 3

Dimensions and sizes:	[TIME | 360] x [lat | 90] x [lon | 180]Coordinates:             TIME: [76.5..3636.5]

            lat: [-87.86..87.86]            lon: [ 0..358]Number Of Attributes: 4

  missing_value :	-9.999999999999999e+33  _FillValue :	-9.999999999999999e+33

  long_name :	TOS[L=3:122]  history :	From tos_Omon_CMCC-CM_decadal1980_r1i1p1_198011-199012




---Vanúcia Schumacher

Mestranda em Meteorologia - UFVMeteorologista -UFPel

Departamento de Meteorologia Agrícola - DEA

Cel: (31) 9978 2522 DEA: (31) 3899 1890

 		 	   		  

_______________________________________________

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



 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140822/45274076/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: esm.regridding.ncl
Type: application/octet-stream
Size: 1733 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140822/45274076/attachment.obj 


More information about the ncl-talk mailing list