Re: 2D varying z-axis

[Arne Biastoch <abiastoch@xxxxxxxx>]

Hi Ferreters,

sorry, I have to bother you again. Steves suggestion did not solved by
case. But I am sure that this is a general problem that all modelers
with variable bottom and free surface formulations should also have. 

I hope that there might be solutions around (Steve might already have 3
lines in mind that do this...).

Here's the case:

I have a standard netCDF file with one depth axis. Here I would like to
integrate not over the full water column but only down to a certain
depth given by a 2D field (the bottom is not at the edges of the grid
cells).

standard case      shallower            deeper
   +---+              +---+              +---+
   |   |              |   |              |   |
   | T |              | T |              | T |
   |   |              |   |              |   |
   +---+              +---+              +---+
   |   |              |   |              |   |
   | T |              | T |              | T |
   |   |              |   |              |   |
   +---+              +---+              +---+
   |   |              |   |              |   |
   |   |              |   |              |   |
   | T |              | T |              | T |
   |   |              |   |              |   |
   |   |              |...|  <---        |   |
   +---+  <---        +---+              +---+
   |...|              |...|              |   | 
   |...|              |...|              |...| <---
   |...|              |...|              | T |
   |...|              |...|              |...|
   |...|              |...|              |...|
   +---+              +---+              +---+

T are the values that are NOT missing values, "<--" marks the bottom,
given by a 2D field

The most intuitive solution would be:

shade T[z=0:depth] whereby depth[i,j]

but as I understand FERRET it is not possible to use variable
integration limits.

I appreciate any help.

Ciao,
	Arne


Steve Hankin wrote:
> 
> Hi Arne,
> 
> A good starting point might be the FAQ
> 
>     "How do I handle sigma coordinate output in Ferret?"
> 
> at
> http://www.ferret.noaa.gov/Ferret/FAQ/data_management/sigma_coordinate_demo.html
> 
>     Happy Ferreting - steve
> 
-- 

Dr. Arne Biastoch
Scripps Institution of Oceanography        phone: +1-858-822-3787
Physical Oceanography Research Division    fax  : +1-858-534-9820
MS: 0230                                   email: abiastoch@ucsd.edu
8605 La Jolla Shores Dr.            
La Jolla, CA 92093-0230, U.S.A.   http://www.ecco.ucsd.edu/~biastoch