Ocean Modeling Discussion

Hi Kate,
Thanks so much for the quick response. (Yes, I turned off periodicity.) You mentioned stratification. I used ANA_INITIAL to set the initial T and S to constants, so there should be no stratification. Would you still expect the pressure gradient problem to create velocities? It appears so ...

At the risk of asking a question with an obvious answer, here goes...

Running the upwelling example with the following cpp defs, with T & S set to constants and sustr and svstr set to 0, produces zero velocities as expected. My question in why, when turning on NONLIN_EOS (so that density is a ...

I changed AKT_BAK in ocean.in to try (one at a time) all three of the following:

AKT_BAK == 1.0d-6 1.0d-6 ! m2/s
AKT_BAK == 5.0d-6 5.0d-6 ! m2/s
AKT_BAK == 0.5d-6 0.5d-6 ! m2/s

The model results for the three different runs appear to be exactly the same. In ocean_his.nc, the variable AKT_BAK ...

I am trying to use CORE daily forcing data with time origin 1/1/48, but something is setup incorrectly. I set TIME_REF= 19480101.0d0, but the model is not accessing the correct record (shows Rec=**** in the output file) and the initial time is = 21899 21:00 when Tmax = 21899.875, even though DSTART ...

Hi, I'd like to use Vslice (part of Rslice) to make figures of ROMS output. I haven't yet been able to figure out how to get the same exact location of the profile for Vslice when opening different netcdf his files. I think I need to modify the mfiles to temporarily hardcode in the beginning and ...

It is possible to output winds (and many other variables) to the ocean_his.nc history file by setting the following in ocean.in

Hout(idUair) == T
Hout(idVair) == T

But is it possible to output Tair?

Hout(idTair) == T

had no effect for me (and looking in wrt_his.nc, the code is not there to ...

(Min = 1.00000000E+35 Max = -1.00000000E+35)

I had the same problem with one of my netcdf forcing files from NCEP and (thanks to your post) I similarly fixed it by flipping my latitude values - and my data to match.

Yep, pn and pm numbers make sense.
Thanks for clarifying that Kate!

-Tasha

I have 2 grid files (one 150x100 and one 330x170) for approximately the same area lat/lon coverage. The domain is 62N to 80N and 170E to 250E. I created the grids using seagrid and then ran seagrid2roms on them. I was looking in my grid files at xl and el (and x_rho, etc) and they seem to be off by ...

I would like to set ssh in ana_ssh.h higher in one part of the domain and then observe the effect on ocean velocities in response to the raised ssh. In theory this should work, no?

I defined ANA_SSH and modified ana_ssh.h (and have it in an appropriate directory where the modified file should be ...

Yes, I was missing SOUTH_TNUDGING. Oops.

Thanks so much for your reply Kate!

I'm doing something wrong and your help in figuring out what is very much appreciated!

I am trying to set an open boundary condition for temperature along the southern boundary, but the warm waters aren't showing up in the model results. The model is running fine, not blowing up, just not having ...

Has anyone used ROMS for a pan-Arctic simulation? If so, how was the North Pole handled?

Thank you, Tasha

Thank you for the suggestion, Kate. I could edit the new "initial" file and run the model again, but I was hoping to run the model for 15 yrs in a row, which would not be ideal to edit the file again and again.

I ended up getting my code modification to work, I just had to fix one little bug ...

Hi all,
What is the most straightforward way to reset certain tracers back to their initial state (from the *ini*.nc file) once per year?

I tried copying and modifying get_state and calling it from nl_ocean.h. It calls get_state okay, but the reset values don't show up in ocean_his.nc. Also the ...