Search found 12 matches

by crowley
Thu Jul 02, 2015 3:22 pm
Forum: ROMS Problems
Topic: Changing the range would blow up
Replies: 4
Views: 1584

Re: Changing the range would blow up

I suspect the current velocities through that pass are just too high. Try reducing your timestep a little and see if that helps.
by crowley
Tue Jun 26, 2012 4:44 pm
Forum: ROMS Discussion
Topic: pressure gradient in the diagnostic output
Replies: 4
Views: 3184

Re: pressure gradient in the diagnostic output

Most of the work calculating the 2D diagnostics happens in step2d_LF_AM3.h. I believe the forcing terms for the coupling are computed (at least in part) in rhs3d.F.
by crowley
Fri Jun 24, 2011 4:24 pm
Forum: ROMS Problems
Topic: East-west tidal currents issue
Replies: 7
Views: 3107

Re: East-west tidal currents issue

You may want to try adding the N2 constituent. For the elevation it is at least as large as the S2 and larger than both the K1 and O1, within Long Island Sound anyway. Is the East River resolved in your grid? The tidal forcing through the East River and the western end of LIS is an important compone...
by crowley
Thu Dec 09, 2010 5:43 pm
Forum: ROMS Ecosystem
Topic: how to model COD (chemical oxygen demand) ?
Replies: 3
Views: 3620

Re: how to model COD (chemical oxygen demand) ?

2) If can, How to do add the decay rate in the passive tracer ?
I don't believe that ROMS has this capability, but you can modify the code to do it. I once applied a simple exponential decay to some passive tracers. I think I added code to step3d_t.F.
by crowley
Sat Jan 16, 2010 2:13 am
Forum: ROMS Source
Topic: Adding variables to history and averages files
Replies: 15
Views: 10335

Re: Adding variables to history and averages files

Do you have UV_VIS2 or UV_VIS4 defined? If not, then NDM2d is off by 1 (it was never incremented) You will need to adjust mod_scalars.F. I would try this: # ifdef DIAGNOSTICS_UV ! ! Indices for 2D momentum diagnostic variables. ! Order is important! ! M2pgrd=1 ic=1 # if defined UV_COR M2fcor=ic+1 ic...
by crowley
Mon May 11, 2009 4:07 pm
Forum: ROMS Information
Topic: Lateral Boundary Conditions?
Replies: 3
Views: 3823

Re: Lateral Boundary Conditions?

Check the bibliography link on the ROMS homepage for a paper (with pdf link) on boundary conditions by Marchesiello et al. This paper also has references to many other helpful sources (Orlanski, Flather, Chapman, etc.) I recall finding the Palma & Matano papers quite helpful when first learning abou...
by crowley
Thu Mar 19, 2009 5:01 pm
Forum: ROMS Information
Topic: Chemical Tracers
Replies: 4
Views: 2958

Re: Chemical Tracers

I did this for pesticide concentrations some time back. They were just passive tracers and I added a bit of code to step3d_t.F to handle the decay. My situation was a little more complicated than Radium would be (increased decay rate in daylight), but it was actually pretty straightforward.
by crowley
Thu Dec 04, 2008 2:10 am
Forum: ROMS Discussion
Topic: about tidal forcing
Replies: 5
Views: 2158

Re: about tidal forcing

In addition to the boundary condition problem, you should think about the behavior of tides in a semi-enclosed basin, as Kate said. Changing the depth of the basin will change the speed of the wave, hence the transit time through your basin. Think about what that can do to the relative phases of the...
by crowley
Sat Nov 08, 2008 1:10 am
Forum: ROMS Discussion
Topic: timestep musings
Replies: 8
Views: 3541

Re: timestep musings

This sounds familiar, but I think I was seeing it from the other side. Things would blow up after a restart, but then i would try restarting from a different time--sometimes an hour before or after, sometimes only a few timesteps away--and all would be well. I knew I was running really close to the ...
by crowley
Fri Jan 11, 2008 1:28 am
Forum: ROMS Information
Topic: What is the output of DIAGNOSTICS_TS/UV?
Replies: 6
Views: 5145

Re: What is the output of DIAGNOSTICS_TS/UV?

Way back when I submitted the code to output the diagnostic terms, I included a page or two of documentation. Not very wordy, but it does list the different diagnostic variables, etc.

It can be found by clicking Documentation on the ROMS homepage and choosing Diagnostics.
by crowley
Wed Jan 31, 2007 10:37 am
Topic: Bottom stress in ROMS
Replies: 4
Views: 10640

For a 3D application, the bottom stress is computed in set_vbc.F and is not computed in step2d.F at all. From step2d.F: # ifndef SOLVE3D ! !----------------------------------------------------------------------- ! Add in bottom stress. !---------------------------------------------------------------...
by crowley
Wed Nov 19, 2003 9:35 pm
Forum: ROMS Bugs
Topic: Array bound problem in get_initial
Replies: 1
Views: 2250

Array bound problem in get_initial

I am (finally) porting my application to version 2.0, but I have run into a problem. The array bound check is barfing on get initial: Run-Time Error 455: Array section bounds inconsistent with parent array In Procedure: get_initial Diagnostics Entered From Subroutine get_initial Line 249 Entered Fro...