How to extract vertical grid from soda/history file.

Facts, news, and guidance about ROMS software

Moderators: arango, robertson

Post Reply
Message
Author
User avatar
p.heidary
Posts: 13
Joined: Fri Aug 10, 2012 4:11 pm
Location: khaje nasir unversity

How to extract vertical grid from soda/history file.

#1 Post by p.heidary » Wed Jul 30, 2014 4:19 pm

Hello all
I am working on Lagrangian particle tracking (tracpy) and i decided to use soda as input file. But I don't know how to find vertical grid parameters such as hc, Cs_r,
Cs_w, theta_b, theta_s, hmin, tcline.
I have the same problem with my old roms history file.

Thanx in advance.

User avatar
kate
Posts: 3678
Joined: Wed Jul 02, 2003 5:29 pm
Location: IMS/UAF, USA

Re: How to extract vertical grid from soda/history file.

#2 Post by kate » Wed Jul 30, 2014 4:46 pm

SODA is not a ROMS run and therefore has a z-coordinate vertical grid. Its grid file has:

Code: Select all

 DEPTH = 5.01, 15.07, 25.28, 35.76, 46.61, 57.98, 70.02, 82.92, 96.92, 
    112.32, 129.49, 148.96, 171.4, 197.79, 229.48, 268.46, 317.65, 381.39, 
    465.91, 579.31, 729.35, 918.37, 1139.15, 1378.57, 1625.7, 1875.11, 
    2125.01, 2375, 2625, 2875, 3125, 3375, 3625, 3875, 4125, 4375, 4625, 
    4875, 5125, 5375 ;
For an old ROMS run, those parameters should be in the history file. You can run "ncdump" on the file, pipe to more (or less) and see:

Code: Select all

        int Vtransform ;
                Vtransform:long_name = "vertical terrain-following transformation equation" ;
        int Vstretching ;
                Vstretching:long_name = "vertical terrain-following stretching function" ;
        double theta_s ;
                theta_s:long_name = "S-coordinate surface control parameter" ;
        double theta_b ;
                theta_b:long_name = "S-coordinate bottom control parameter" ;
        double Tcline ;
                Tcline:long_name = "S-coordinate surface/bottom layer width" ;
                Tcline:units = "meter" ;
        double hc ;
                hc:long_name = "S-coordinate parameter, critical depth" ;
                hc:units = "meter" ;
then further down:

Code: Select all

 Vtransform = 2 ;

 Vstretching = 4 ;

 theta_s = 7 ;

 theta_b = 2 ;

 Tcline = 250 ;

 hc = 250 ;

User avatar
p.heidary
Posts: 13
Joined: Fri Aug 10, 2012 4:11 pm
Location: khaje nasir unversity

Re: How to extract vertical grid from soda/history file.

#3 Post by p.heidary » Thu Jul 31, 2014 8:54 am

Great, I think I have to make some changes to tracpy to recieve z-coordinate.

Post Reply