Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Sun Jul 21, 2019 2:25 am




Post new topic Reply to topic  [ 5 posts ] 

All times are UTC

Author Message
PostPosted: Sun Feb 10, 2019 11:39 pm 
Offline

Joined: Mon Feb 16, 2004 4:54 pm
Posts: 11
Location: CSIRO
Hi all,

just wondering if anyone has some code that calculates the soundspeed from ROMS output files ?

Cheers
Frank


Top
 Profile  
Reply with quote  
PostPosted: Mon Feb 11, 2019 5:12 am 
Offline
Site Admin
User avatar

Joined: Wed Feb 26, 2003 4:41 pm
Posts: 1078
Location: IMCS, Rutgers University
It has been in ROMS plotting package that uses the NCAR's GKS library since the beginning. Check the file src/Utility/svel.F in the plot repository. It uses the Che and Millero (1976) formula.

Chen, C.T. and F.J. Millero, 1976: re-evaluation of Willson sound-speed measurements in pure water. J. Acoust. Soc. of Amer., 60(6), 1270-1273

I am not familiar with the literature of sound speed computation. I bet that there are newer formulations out there.

Also, it will be very trivial to compute from ROMS nonlinear equation of state since it has all the variables that it needs:

Code:
svel2 = ABS(10000.0./(den * gamma))
svel = SQRT(svel2);

where gamma is the adiabatic and isentropic compressibility coefficient (1/Pa).


Top
 Profile  
Reply with quote  
PostPosted: Mon Feb 11, 2019 2:31 pm 
Offline
User avatar

Joined: Mon Apr 28, 2003 5:44 pm
Posts: 491
Location: Rutgers University
There is also code in the CSIRO seawater toolbox for Matlab (function sw_svel) though seawater has been superseded by the Gibbs SeaWater (GSW) toolbox for TEOS-10
http://www.teos-10.org/pubs/gsw/html/gsw_contents.html
and in there it is function gsw_sound_speed

_________________
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu


Top
 Profile  
Reply with quote  
PostPosted: Mon Feb 11, 2019 10:33 pm 
Offline

Joined: Mon Feb 16, 2004 4:54 pm
Posts: 11
Location: CSIRO
Thanks guys !

Hernan - OK - I remember having problems with compiling the plotting package in the past due to ncarg library. I'll give it another go.

Generally I think it would be nice to decouple the soundspeed calculation from the actual plotting routines and be able to write it into the average or history output files given T and S.

I look into this.


Top
 Profile  
Reply with quote  
PostPosted: Mon Feb 11, 2019 10:34 pm 
Offline
Site Admin
User avatar

Joined: Wed Feb 26, 2003 4:41 pm
Posts: 1078
Location: IMCS, Rutgers University
Just do it in Matlab.


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 5 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 4 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
cron
Powered by phpBB® Forum Software © phpBB Group