Ocean Modeling Discussion

You can also use ROMSPath, which uses the same internal logic as floats in ROMS and some of the same formatting and code structure as LTRANS. (In the interest of disclosure, I was a low-level assistant to Eli Hunter in development.) Code: https://github.com/imcslatte/ROMSPath Hunter et al, 2022 ...

Tong, Look at equations 32 and 34 in the documentation here: https://www.myroms.org/wiki/Vertical_Mixing_Parameterizations The meaning of the gls variable that gets computed and saved is (I believe) determined by the p, m, and n coefficients that you choose. You can recover the dissipation from the ...

You'll need to install XCode (or at least the command line tools; I haven't kept up to date with the continual XCode changes) Then, as John said, install Macports. Parker's old notes (if that's what you're referring to) aren't completely obsolete, so they're a useful guide. I haven't installed ROMS ...

We invite applications for a Postdoctoral Associate to investigate how waves, turbulence, and flow-induced larval behavior affect larval transport in the Middle Atlantic Bight and adjacent estuaries. Waves and turbulence induce behavioral responses that alter larval vertical distributions and may ...

I am just starting a new LTRANS setup and have encountered the same problem. I haven't dug too far into it yet. Did you identify the cause?
--Greg

Using revision 467 I found that the diagnostics were not writing the vertical velocity terms. It looks like this is due to missing code in inp_par.F. I have added lines at 2767 that seem to have fixed this. ELSE IF (TRIM(KeyWord).eq.'Dout(M3vvis)') THEN IF (M3vvis.le.0) THEN IF (Master) WRITE (out ...

I would guess that most of us who use GLS use the Kantha and Clayson stability functions. The two-equation models in GLS require two choices: 1) choose a length-scale model (the numbers in ocean.in). 2) choose a stability function (choose in cppdefs). John Warner's 2005 paper (Ocean Modelling, 8, 81 ...

Alex, I've had a nearly identical problem recently. I would run a model for two weeks. Then use perfect restart to run for a few days more. I was only using PGI. When I optimized at level 3 (or 2) it blew up (sometimes it would go through a couple time-steps, I believe). When I optimize at level 1 ...

Peter,

I have had trouble with k-epsilon being a bit less stable than the other closure schemes. Reducing the time step has helped fix this. Not an entirely satisfactory solution, I know.

Greg

The mixing schemes are in gls_corstep.F and my25corstep.F, which are in the Nonlinear directory of your source code. You can get the essential features of Mellor-Yamada by choosing kkl in GLS.

Kate--

What settings are you using for GLS? I have found k-epsilon to be less stable than the other versions. Did you try k-omega or gen?

Thanks for checking that.
I get the same error with version 297.
I'll dig deeper into my files.

greg

Thanks.
FLT_TEST compiles without SOLVE3D, but not with SOLVE3D.
Does it work both ways for you? It doesn't seem to depend on the version I am using.

I am trying to run a simulation with the cpp option FLOATS, and I get a compilation error. Has anyone else encountered this? This is my first time with FLOATS. Are there other cpp options that I should be using? I am using revision 291. The error messages are: PGF90-S-0038-Symbol, flt_lagran, has ...

I found the same thing as jacopo. Reducing the time step made it run fine.
Have you tested it with all three stability functions (Kantha-Clayson, Canuto A and Canuto B?)