Point sources of passive tracer

Suggest improvements/optimizations to the ROMS code.

Moderators: arango, robertson

Post Reply
Message
Author
staalstrom
Posts: 31
Joined: Mon Feb 04, 2008 3:43 pm
Location: NIVA, OSLO, NORWAY

Point sources of passive tracer

#1 Unread post by staalstrom »

Hello

I run roms 2.2.
The model is working with passive tracer dye_01 initiated in the ini NetCDF file. But I want the dye to be put into the model with a river.

TS_MPDATA, TS_DIF2, TS_PSOURCE and UV_PSOURCE is defined inn cppdefs.h.

I added this to varinfo.dat

'river_dye_' ! Input
'river runoff dye concentration'
'kilogram meter-3' ! [kg/m3]
'river_dye, scalar, series'
'river_time'
'idRtrc(inert)'
'nulvar'
1.0d0

And I made a river forcing file with river_dye_01 with concentration=1 in the top layer. The model run without any problem, but no dye concentration was added to dye_01 in the history file.

The river_flag is set to 3 in the river forcing file.

Andre Staalstrom

zhang
Posts: 27
Joined: Fri Mar 11, 2005 9:17 pm
Location: Woods Hole Oceanographic Institution

#2 Unread post by zhang »

I've never used 2.2. But for the latest version of ROMS we can get using svn, there are more than one place we need to change to add point source for a passive tracer.
Some changes in Nonlinear/get_idata.F is required to turn on SOURCES(ng)%Lsrc and SOURCES(ng)%Ltracer for the added passive tracer. in Modules/mod_ncparam.F, you need to add a case for the index variable name you added in varinfo.dat. There might be more than these. I did it more than a year ago. These are what I can think of now.

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

#3 Unread post by kate »

In modern ROMS, you need to check the writing of passive tracers in the input file as well:

Hout(inert) == T ! inert passive tracers

Check that you have this, that inp_par is reading this, and that def_his and wrt_his know about it as well.

gianni

Re: Point sources of passive tracer

#4 Unread post by gianni »

Hello,
I already use passive tracers in the initial condition successfully but now I want the dye to be put into the model with a river. Before starting modifying roms code I would ask to Andre Staalstrom if he finally have done it well and if I need to know more than what is written in this topic.
Thank you in advance,
Gianni.

Post Reply