Access to additional scalar variables.

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
jtkirby
Posts: 2
Joined: Fri Jun 04, 2004 8:50 pm
Location: Center for Applied Coastal Research, U. of Delawar

Access to additional scalar variables.

#1 Post by jtkirby » Thu Nov 02, 2006 2:17 pm

I am new to any of ROMS' parts other than basic hydrodynamics. I am interested in adding a computation of a few scalar fields that probably fall outside the range of the existing variables in ECOSIM or BIO_FASHAM. I am wondering how to access additional scalar fields (where defined? mod_scalars.F, inert(NTP)? Where to go to modify governing equations for the particular tracers I'd like to model, etc.)

If someone could point me in the right direction, I'd appreciate it.

Jim Kirby

George
Posts: 9
Joined: Mon Jun 26, 2006 1:15 pm
Location: International Arctic Research Center

#2 Post by George » Fri Nov 03, 2006 5:50 pm

Hi Jim
There a a number of files that you will have to edit to add a biological tracer. I have included a list of the files and a short description of what they do below. Let me know if you have any questions.
George

mod_scalars.F-Tracer identification indices

varinfo.dat -Input/Output variables information

analytical.F -sets initial conditions for biological tracer fields

fasham.h/ecosim.h - computes the biological sources and sinks
- this is where the governing equations are.

mod_ncparam.F -input/output info for NetCDF files

mod_param.F -set number of bio tracers

User avatar
wilkin
Posts: 526
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

#3 Post by wilkin » Sat Nov 04, 2006 1:53 am

This is already done for you (but yes I know we haven't documented it...)

In Include/cppdefs.h ...
#define T_PASSIVE

this will activate "number of passive tracers" in Modules/mod_param.F which by default is
NPT = 1

If you wish to have more than one passive tracer edit mod_param.F accordingly
! NPT Number of extra passive tracer type variables to !
! advect and diffuse only (dyes, etc). !

You DO NOT need to do anything to varinfo.dat. By default the varinfo.dat configuration sets the naming convention for the passive tracers to the variable name stem "dye_" and the variables in the netcdf file will be named dye_01, dye_02, dye_03 ... depending on the value of NPT.

So to initialize your passive "dye" create a variable dye_01 in the initial conditions netcdf file. I make this file separately and merge the two files using the netcdf kitchen sink 'ncks' from the nco toolbox - but how you do it is up to you.

As is often the case these days, the hard work is all in the netcdf file manipulation not with ROMS itself.

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

User avatar
susonic
Posts: 160
Joined: Tue Aug 21, 2007 5:44 pm
Location: Jeju National University
Contact:

Re: Access to additional scalar variables.

#4 Post by susonic » Fri Sep 18, 2009 1:18 am

Hi Dr.Wilkins,

I wanted to try to use the passive tracer. And I was able to add dye_01 in my initail condition and also boundary
condition file. The model was run well. Now I want to run with passive tracer in my river forcing file.
I have been searching out the answer for this question in ROMS forum, but I couldn't get specific information.
So, would you help me out and give me specific information about 'add passive tracer in river forcing' ?

Thanks in advance,

-Peter

User avatar
susonic
Posts: 160
Joined: Tue Aug 21, 2007 5:44 pm
Location: Jeju National University
Contact:

Re: Access to additional scalar variables.

#5 Post by susonic » Mon Sep 21, 2009 9:39 pm

Dear all,

I wanted to try to use the passive tracer in river forcing as mentioned above.
Would anyone help me out and give specific informations to me about
'add passive tracer in river forcing' ?
Thanks in advance,

-Peter

Post Reply