[Thread Prev][Thread Next][Index]

Re: Geostropy calcuklations?



Hello Sudheer,

First, use the script dynamic_height.jnl, which comes with Ferret, to
compute dynamic heights and geostrophic currents.  The script has
enough information to help you use it.  "Fdata dynamic_height" will
give you the location of the script.

Second, when you want to use an axis variable, like latitude (y), in
your journal file, you have to specify the grid from which to pick it.
For example, 

yes? sh grid temp
(Now suppose the answer is gfy1.)
Then, invoke y as follows.
yes? y[g=gfy1]

Regards, Shankar.

D. Shankar                           Ph. (O): [91](832) 22 6253 * 4312/4400
Physical Oceanography Division,                         22 1322 * 4312/4400
National Institute of Oceanography,      (R): [91](832) 23 8208      
Dona Paula,                          Fax    : [91](832) 22 3340/9102
Goa 403 004,                         email  : shankar@darya.nio.org
India.                                        shankar@csnio.ren.nic.in

On Mon, 11 Feb 2002, Dr.Sudheer Joseph wrote:

> 11/02/2002 12:40:13 PM
> 11/02/2002 11:32:43 AM
> Hellow,
>               I am a fresh user of ferret, Can any body help me to 
>               solve the following problem, I am trying to get 
>               Geostropic current vector plots usint the following 
>               script , 
> I am gettign the error message that
> 
> *** NOTE: Ambiguous coordinates on Y axis: G/(RHO*F)
> 
> How do I correct this?
> secondly,
> 
> When I change the set reion 
> set reg/x= 60:100/y= 0:30/z=1:20/l=6
> to
> set reg/x= 20:140/y= -50:30/z=1:20/l=6
> it is giving the following error message
> 
>  **ERROR: inconsistent sizes of data regions: Y axis
>           F has 30 points (J=01:30)
>           expression has 80 points (J=041:120)
> vector um, vm
> Command file, command group, or REPEAT execution aborted
> +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++====
> Please Help
> rgds
> sjo
> 
> used script
> 
> use levitus_climatology.cdf
> 
> let g = 9.81
> let pi = 3.14159
> let omega = 7.292e-5
> let f = 2*omega*sin(y*pi/180)
> 
> ! Set region to all space (use actual values for ni, nj and nk)
> ! set reg/i=1:ni/j=1:nj/k=1:nk
> set reg/x= 60:100/y= 0:30/z=1:20/l=6
>                        
> ! Compute density
> 
> let rho = rho_un(salt,temp,z)
> 
> ! Compute vertical shears
> let gamma = g/(rho*f)
> let u_z = gamma*rho[y=@DDC]
> let v_z = gamma*rho[x=@DDC]*(-1)
> 
> let u = u_z[k=@DIN]
> let v = v_z[k=@DIN]
> 
> ! The flow field at the base of a 100 m mixed layer assuming the
> ! depth of no motion is at 2000 m.
> !let um = u_z[z=2000:100@DIN]
> !let vm = v_z[z=2000:100@DIN]
> let um = u_z[z=2000:1@DIN]
> let vm = v_z[z=2000:1@DIN]
> vector um, vm
> go land green



[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement