Hi
Your listing looks like the output of SHOW DATA/ATTRIBUTE. (XCOASTS_C), which is listed in parentheses, is a coordinate axis. It won't be changed unless you use SET AXIS or DEFINE AXIS commands to explicitly change it.
To define new variables using the original variable names, and rewrite your file with the new definitions, you first need to rename the variables in the input file. Try a SHOW DATA after the SET VAR/NAME commands to see what you've got.
yes? use mydata.ncIf you wanted to keep some variables from the original dataset and replace others, you would still need to SAVE both the ones you are keeping and the new ones. /CLOBBER deletes the old file and replaces it with a new one.
yes? set var/name=lon_in lon_coasts_c
yes? set var/name=lat_in lat_coasts_c
! Now define the new variables in terms of the previous ones.
! Use /UNITS= and /TITLE= to give the new variables the correct attributes.
yes? let/units="degrees_east" lon_coasts_c = 2* lon_in
! Or, you can get the units and other information from the original variables
! using the RETURN= keyword.
yes? let/units="`lat_in`" lat_coasts_c = 2* lat_in
yes? save/clobber/file=mydata.nc lon_coasts_c, lat_coasts_c
Ansley
On 2/16/2011 9:20 AM, Szymon Roziewski wrote:Dear Ferreters!
I have some data in a netcdf file. I want to do some transformations to these data and then change existing in the netcdf file by freshly computed.
And also I would like to know what is this (XCOASTS_C), because this is not a variable - what is it?
Variable VarType(in dset) AttributeName AttType Size OutFlag AttValue
------------------------------------------------------------------------------------------
. history CHAR 24 T FERRET V5.41 4-Dec-07
(XCOASTS_C) DOUBLE point_spacing CHAR 4 T even
orig_file_axnameCHAR 9 F XCOASTS_C
LON_COASTS_C FLOAT missing_value FLOAT 1 T -1.000000E+34
_FillValue FLOAT 1 T -1.000000E+34
long_name CHAR 12 T lon_coasts_c
history CHAR 16 T From coast_c.dat
units CHAR 12 T degrees_east
LAT_COASTS_C FLOAT missing_value FLOAT 1 T -1.000000E+34
_FillValue FLOAT 1 T -1.000000E+34
long_name CHAR 12 T lat_coasts_c
history CHAR 16 T From coast_c.dat
units CHAR 13 T degrees_north
So I wanted to change e.g. existing LON_COASTS_C, LAT_COASTS_C by new ones. Simultaneously with no changing (XCOASTS_C).
Variable attributes could be rewritten I suppose.
How can I deal with that?
Many greetings!
Szymon Roziewski
Attachment:
exampleBaltic.gif
Description: GIF image