I'm working on getting our Princeton Ocean Model data into OpenDX. If
anyone has been successful with this already I'd be very interested in
talking with you.
The data is in a netCDF file and I'm starting with a simple set that
boils down to:
netcdf temp-no-time {
dimensions:
ypos = 120 ;
xpos = 180 ;
sigma = 22 ;
variables:
float depth(ypos, xpos) ;
depth:valid_range = 0.f, 5000.f ;
depth:long_name = "Bathymetry" ;
depth:units = "meters" ;
depth:missing_value = -99999.f ;
short temp(sigma, ypos, xpos) ;
temp:valid_range = -5.f, 30.f ;
temp:long_name = "Temperature" ;
temp:units = "Celsius" ;
temp:scale_factor = 0.0005340739f ;
temp:add_offset = 12.5f ;
temp:missing_value = -99999.f ;
}
I can open it and visualize the bathymetry fine. I use the Rubbersheet
module for that. I also am able to visualize Temperature data using
Isosurface. However, the depth is wrong.
The depth values in the Princeton Ocean Model (POM) are values from -1.0
to 0.0 and they are a scale factor such that:
z(i,j,k)=sigma(k)*depth(i,j)
I haven't succeeded in retrieving the sigma data yet so for now I'm just
faking it a bit. What I've got so far is Mark->Compute->Unmark and
feeding the Temperature data into the Compute module. In Mark, I've set
the "Name" input to "positions". The Expression in Compute is:
[a.x, a.y, a.z/22]
This just scales the depth back. What I really want is to feed Depth
into the Compute module as input B and use the Expression:
[a.x, a.y, a.z * [b.x, b.y] / 22]
When I do this, I get the following in the Message Window:
Error: Compute: Bad parameter: Inputs must be of same length / Input 1
not matching the master (input 0).
I'm unclear of the format for specifying the B part of the expression.
I have used the Print module to see that it is a 2D set but maybe that
is the problem. Do I have to somehow make that a 3D set in order to
make Input 1 match Input 0? Is that what the message means?
Am I going about this the wrong way?
I've put the network, cfg file and netCDF file at:
http://rocky.umeoce.maine.edu/cousins/pom.tgz
The Bathymetry data is dealt with in the left part of the network and
the Temperature data is dealt with in the right part of the network.
Thanks very much,
Steve
______________________________________________________________________
Steve Cousins, Ocean Modeling Group Email: [EMAIL PROTECTED]
Marine Sciences, 208 Libby Hall http://rocky.umeoce.maine.edu
Univ. of Maine, Orono, ME 04469 Phone: (207) 581-4302
