cdo

estrazione serie temporali su punti

per questo scopo si puo` usare la funzione di interpolazione (remapdis, ma non remapbil) su un grigliato di tipo cell, cioe` non strutturato, senza indicare il numero dei verici e la definizione dei boundaries delle celle.

Il descritore di griglia è: 

gridtype = cell
gridsize = 5
xvals = -4 -1 0 1 4
yvals = -4 -1 0 1 4
dove (xvals, yvals) contengono le coppie di coordinate dei punti su cui si vuole interpolare.

Il comando per estrarre la serie è: 

cdo remapdis,grid inputfile.nc outputfile.nc

L'intestazione dell'output file è 

netcdf outputfile {
dimensions:
	cell = 5 ;
	lev = 33 ;
	time = UNLIMITED ; // (24 currently)
variables:
	double lon(cell) ;
		lon:units = "degrees_east" ;
		lon:standard_name = "longitude" ;
	double lat(cell) ;
		lat:units = "degrees_north" ;
		lat:standard_name = "latitude" ;
	float lev(lev) ;
		lev:long_name = "sigma_level" ;
		lev:units = "sigma" ;
	double time(time) ;
		time:units = "day as %Y%m%d.%f" ;
	float var(time, lev, cell) ;
		bol:standard_name = "bol" ;
		bol:units = "micro_mol/m^3" ;
		bol:grid_type = "cell" ;
		bol:coordinates = "lon lat" ;


// global attributes:
		:CDI = "Climate Data Interface version 1.3.0" ;
		:Conventions = "CF-1.0" ;
		:history = "Wed Mar 11 11:34:43 2009: cdo remapdis,grid interp.nc qq.nc" ;
		:CDO = "Climate Data Operators version 1.3.0 (http://www.mpimet.mpg.de/cdo)" ;
}