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
Access to additional scalar variables.
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
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
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.
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
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu
Re: Access to additional scalar variables.
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
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
Re: Access to additional scalar variables.
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
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