[Thread Prev][Thread Next][Index]
[ferret_users] How to grab z axis info?
- To: ferret_users@xxxxxxxx
- Subject: [ferret_users] How to grab z axis info?
- From: Szymon Roziewski <szymon.roziewski@xxxxxxxxx>
- Date: Wed, 20 Apr 2011 15:29:39 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:date:message-id:subject:from:to :content-type; bh=i0Uk4oANVE68nsIcVS/+MmUhtf2jaBcG82kWb+BF7Io=; b=bxVTuEjHMHsfQTbCjx+UD1LvOfFIjipboaD+dyVkppL53PQurCbxhswITgKLgVzq6s rZrWe4CsRYNHqBSBdKfXDmyNTlhme26cA8x+hjkkEwfiUWkJe8Lj5+eFB7jwwKzvmQpS NiVO/+84Z/kIQcOJfoRX75AxcRUK9lQNxtgVA=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:date:message-id:subject:from:to:content-type; b=cRA/3DSZME4C5Ay/E8NIqH+UyJn+ndZ0MBfD7FnbGPV0cYHsK6vRILTA3y1yfeBkjb SiTXInIgZET0QCktoYhAjRGn1wP/4lXgwBhf/spB28o4sP+8E9ZRjqggYCRL8FdcSCu4 Jrqo4fsYm01rM6r5/J7YVFiPe1Axu0WL9AaxM=
- Sender: owner-ferret_users@xxxxxxxx
Hello Ferreters,
I do an interpolation and as a result I get a variable without time axis and z axis information.
So I need to fill these gaps. With time axis I have no problem.
let tmp_t = temperature[k=10, l=1]
define axis/t="`t[g=tmp_t,l=1],return=tend`":"`t[g=tmp_t,l=1],return=tend`":1 time
let timestamp = 0* t[gt=time]
This above works.
The problem starts when I proceed to z axis info.
I have a problem with getting this info from origin variable.
I've been trying to get there by that
list "`Z[g=tmp_t,k=10],return=tend`"
but with failure.
**ERROR: command syntax: *
* may not start an _expression_
list "*"
When I renamed z variable into e.g. zs using ncrename I didn't help me.
Could you give me please any suggestions?
I also attach the picture in order to see what I want to get (red line).
I need to add that my netCDF has a following structure (ncdump):
dimensions:
time = UNLIMITED ; // (86 currently)
z = 10 ;
y = 145 ;
x = 132 ;
variables:
float temperature(time, z, y, x) ;
temperature:long_name = "potential temperature" ;
temperature:units = "C" ;
temperature:coordinates = "east_e north_e zz" ;
float east_e(y, x) ;
east_e:long_name = "easting of elevation points" ;
east_e:units = "metre" ;
float north_e(y, x) ;
north_e:long_name = "northing of elevation points" ;
north_e:units = "metre" ;
float zz(z) ;
zz:long_name = "sigma of cell centre" ;
zz:units = "sigma_level" ;
zz:standard_name = "ocean_sigma_coordinate" ;
zz:formula_terms = "sigma: zz eta: elb depth: h" ;
float z(z) ;
z:long_name = "sigma of cell face" ;
z:units = "sigma_level" ;
z:standard_name = "ocean_sigma_coordinate" ;
z:formula_terms = "sigma: z eta: elb depth: h" ;
and in ferret (show data/att)
(time) FLOAT long_name CHAR 4 T time
orig_file_axnameCHAR 4 F time
(z) FLOAT long_name CHAR 18 T sigma of cell face
units CHAR 11 T sigma_level
standard_name CHAR 22 F ocean_sigma_coordinate
formula_terms CHAR 26 F sigma: z eta: elb depth: h
orig_file_axnameCHAR 1 F z
temperature FLOAT long_name CHAR 21 T potential temperature
units CHAR 1 T C
coordinates CHAR 17 F east_e north_e zz
--
Z wyrazami szacunku,
Szymon Roziewski
Attachment:
temp.gif
Description: GIF image
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
Ferret
Privacy Policy | Disclaimer | Accessibility Statement