Hi Jean-
The HYCOM model has layer
thicknesses rather than layer depths, which ZAXREPLACE uses. The layer
thicknesses may vary with space and time.We've written a function
function ZAXREPLACE_ZLEV to convert this data to a fixed Z axis:
yes? show function ZAXREPLACE_ZLEV
ZAXREPLACE_ZLEV(V,THICKNESS,ZAX)
regrid V onto Z axis of ZAX based on layer thickness
V: Variable on native z axis
THICKNESS: Thickness corresponding to input Z axis layers in units
of output Z axis
ZAX: Variable with desired z (depth) axis points
Jean, this function will be in the next Ferret release. Or, if you let
me know what platform you're running Ferret on, I could send you a
shared object file with the function.
Ansley
jean oliveira wrote:
Hi friends,
I`m working with the HYCOM ocean model
and I`d like to manage the temperature in Z coordinates. I dont know if
a simple linear interpolation is the best method since a temperature
represent a entire isopycnal layer. For example, suppose it has a
isopycnal layer with thickness of
5 meters
and temperature of 20ºC and the below isopycnal
layer has a thickness of 20 meters and temperature
of 18ºC and it wants to represent de
temperature to each one meter. What I doing: I take the middle of the
layer as the position of the temperature and then make the linear
interpolation to eache meter. Mainly for those that use HYCOM and
FERRET, is this the best way to represent the hycom`s output in Z
coordinates? If not, how can I do a different method in the FERRET?
Thanks a lot for all
Jean
|