Ocean Modeling Discussion

Oh yes - you are right! I initialized ROMS from a larger domain by downscaling, and zeta is around 0.15 m. That's exactly the missing difference!
Thanks, Reiner

In ROMS, I have launched various isobaric floats at Fz0=-1.00 m depth. In the floats output file, however, the variable depth is around -0.85 m all the time:

depth ...

I have inserted 67 Lagrangian floats in a ROMS model in shallow water
(max water depth around 40 m). I expected the floats to move along
(potential) density surfaces in order to conserve isopycnal potential
vorticity, but several floats penetrate the density surfaces.
Here are my questions:
(1) The ...

I have run 2 identical ROMS runs named ROMS_13 and ROMS_13a but the results are different!

and these are the specifications:
ROMS_13
compiled with NetCDF3
NHIS=2880 in ocean.in

ROMS_13a
export USE_NETCDF4=on in build.bash
export NC_CONFIG=nf-config in build.bash
NHIS=360 in ocean.in ...

Is the wind stress calculation in ROMS corrected with the water velocity at the sea surface?

Reiner

Sometimes, atmospheric models provide NEGATIVE relative humidities. It is highly recommended to remove such values because ROMS will immediately blow up because lwrad produces NaNs. It took me 2 days to locate the bug!
Therefore, when producing the forcing fields, I limit
(1) Qair to 0 <= Qair ...

I ran ROMS for 10 days and accidently, my forcing file covered only the first 2 days. Surprisingly, ROMS did not complain or blow up. What is ROMS doing in this case? Does it run without any atmospheric forcing? Does it lock the forcing to the last data in the forcing file? Or an average?

Reiner

Thanks a lot - I will try the link on GITHUB.

Cheers, have a nice weekend!

Reiner

Are there any tools available to calculate energetics in ROMS, i.e. eddy kinetic/potential energy, energy conversions?
Reiner

(1) What is the default scheme for horizontal eddy viscosity and diffusivity? In the options.h file, I did neither define UV_VIS2 nor UV_VIS4, TS_DIF2, or TS_DIF4. The CPU time of that ROMS run was much less compared to a run with UV_VIS4 and TS_DIF4 but the results were approximately the same.
(2 ...

I want to create an animation of SST and SSS from a ROMS history file. Therefore, I need only the uppermost layer in the output file. Is there any way to suppress the output of the other layers? This would be desirable to keep the history file small.

Who knows where I can find a version of "seagrid.m" which works with MATLAB R2014b on LINUX 64bit? Since I upgraded from the previous MATLAB version, seagrid crashes with the following error messages:
------------------------------------------------------
Undefined function 'any' for input arguments ...

I try to generate the mex files for SEAGRID with >>make_seagrid_mex but I did not succeed. I have Ubuntu 11.04 and MATLAB R2010b.
I have hacked the mexopts.sh in the following ways:
# change g77 with ifort
# installed gcc-4.4 and gfortran-4.4
# try different FFLAGS
and
# change links in .../MATLAB ...