Ocean Modeling Discussion

Me again, So I took a bit of time to go through the upward/downward motions in red_tide.h and merged it into fennel.h. I found that the two codes are essentially identical up until the line starting wdt=-Wbio(iswim)*dtdays (red_tide.h) or cff=dtdays*ABS(Wbio(isink)) (fennel.h) I found that I had som...

Thanks for all your replies!

This is a little bit of a daunting task looking at both codes, since they a structured differently, and a lot of potential places to code this incorrectly.

If anybody has any specific guidance for this, it would be most helpful!

Best wishes,
Ash

Many thanks for the quick replies, this is very handy indeed!

Thanks again,
Ash

Hi there, I would like to incorporate a floating mechanism to the phytoplankton component of the Fennel model, to mimic the behaviour of Microcystis cells. Currently, the model has a constant sinking term and I would like to mirror this for floating, i.e. moving the phytoplankton upwards at some con...

Right I see - I tried forcing it to true, but still couldn't get past the error. It's easy enough to do the interpolation myself anyhow.

Thanks for your help!
Ash

Hi there, I am currently trying to set up a nested grid (1 refinement grid). The atmospheric forcing files I provide fit the parent grid exactly (and works). How do I get the forcing files to work directly with the refinement grid without providing a new file (i.e. so that ROMS does the interpolatio...

Works perfectly! Thanks for the quick fix, it's much appreciated!

Thanks! One more ...

initial.f90(174): error #6631: A non-optional actual argument must be present when invoking a procedure with an explicit interface. [INIREC]
CALL get_wetdry (ng, iNLM, IniRec(ng))
---------------^
compilation aborted for initial.f90 (code 1)

Hi there, I ran a svn update this morning and the get_wetdry.f90 code updated and now most of the errors are gone (rev 1067). I'm just left with the following: get_wetdry.f90(31): error #6404: This name does not have a type, and must have an explicit type. [GET_GRID] PUBLIC :: get_grid -------------...

Thanks for the quick reply!

They are all there, does the compilation interact with varinfo.dat?

All the best,
Ash

Hi there, I'm having trouble compiling with the WET_DRY option on the latest source code. I get the following error: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% get_wetdry.f90(142): error #5082: Syntax error, found IDENTIFIER 'VID' when expecting one of: ) ( , : <EN...

I see both variables throughout in the new source code.

All the best,
Ash

Hi there, I recently updated to the latest source code and one of my bottom boundary conditions has stopped working. In particular, a nutrient flux associated with the Fennel biological model. I've added the following to ana_btflux.h ELSE IF (itrc.eq.iNH4_) THEN DO j=JstrT,JendT DO i=IstrT,IendT btf...

Thank you for your help on this.

I made the changes as you recommended and it works perfectly!

Thanks again,
Ash

Hi there, Thank you very much for your quick responses! Apologies I think I didn't explain my problem very clearly. My problem is that when I add evaporation and precipitation, the amount of water that I have in the system doesn't change. Actually in my case I'm not solving for salinity, as everythi...

Hi all, I'm currently running a lake model where E and P are important in terms of removing and adding water to the domain. However I found that ROMS does not support this function. I wondered if anybody added this feature to their model code and if so, could they share their experience. As far as I...

This change works perfectly! Thank you for considering the problem for me, it is greatly appreciated! So the idea is that sources add water and change the concentration of the tracer depending on the value of the tracer in the river, say. A sink on the other hand with remove water, but it won't chan...

So I tried to turn off the tracer at particular source locations. I want the water volume to change, but I don't want to insert any tracer differences. My logic is that I want to somehow mimic the behavior of "LtracerSrc = F" for certain rivers. But perhaps I am not fully understanding. I ...

Thanks for your reply! Ha yes, I went over this page with a fine tooth comb - and I think I can understand how it wants me to locate the rivers for a sink. But the stability was throwing me off a little (making me think I've done something wrong). How would I go about setting the tracer values of th...

I see! I’m using using both, some rivers are sources and some rivers are sinks. Is there a tidy way to handle that?

Thanks for the quick response!

Dear ROMS users, I'm having an issue with sinks in my river files. My understanding is (assume a square enclosed lake) that If I wanted a sink on my west bank, then I need the following: river_direction = 0 river_transport < 0 river location on the land (u-point) with adjacent water point. However, ...