Ocean Modeling Discussion

Hello,
I am trying to calculate vertical gradient du/dz and dv/dz from ROMS output. Usually, we use differential equation shown in picture below. However, it may lead some non-conservative problems, and I dont know how to deal with the the upper half part of the layer 8 in the picture (Red circle ...

Hi, ROMS users

I met some problems handling with RadNud boudary. In ocean.in, TNUDG and M3NUDG mean outflow nudging time scales, I am not sure that what is exact meaning of the ' time scale' .

In my current understanding, inflow time scale is decided by given boundary conditon and outflow conditon ...

Hz could be calculated by modifying set_depth.m:
for k=1:N(ng)
Hz(i,j,k)=z_w(i,j,k)-z_w(i,j,k-1)

zzy

kate wrote: ↑Wed Jan 15, 2020 5:01 pm No, the .h files are included in the fortran files, compiled by the Fortran compiler. Linking to C code is a different beast. Not impossible, just different.

Kate, thanks for your suggestion, I changed ana.h file successfully and I can run ROMSnow.

zzy

The answer might be in ROMS/Functionals//ana_smflux.h, which is the file with the analytic surface stresses for some of the standard test problems. Some of them have things which vary in time as examples. If you want to change things by the timestep number, that is iic in the code.


Thanks a ...

It is called every time step. Changes in time happen in your code.


Thanks, Sorry I have a quesiton and maybe it's stupid.

If i input a vector as sustr series, one element in the vector indicate one step's sustr, how does ROMS read it? It can read the first element in the first time step, the ...

It sounded like you wanted a linear progression in time for your stress. If so, ROMS does a linear interpolation in time from one record to the next and could be doing that work for you. Or for a very simple pattern you can use the ana_xxx.h code which gets called every timestep.


Oh, I want to ...

Thanks, I read the code, it only reads 'days' and 'seconds', it cannot read 's'.

And I also find another problem: the smallest dt of uvstress ROMS can recoganize is 1s. 0.1s is not permitted.

But I need to input uvstress with a dt of 0.1, I am trying to find the related code.

Hi all,

I have some new puzzles with the time units of sustr and svstr.

I set dt=0.1s and Ntime=4000 in ocean.in file, so the whole run time is 400s.

Firstly, I set unit of sustr and svstr as 's' in stress.nc file and input a set of sustr and svstr(dt is 0.1s). For example sustr 0.1s 0.2s 0.3s ...

This is not a bug. It always has to be equal or greater than one by design. If no input forcing files are needed, ROMS will ignore those values on standard input. I don't understand what it is the Fuzz about it. There are a lot of parameters in ROMS standard input that will be only relevant for a ...

1. Undef BULK_FLUXES in .h file. and define ANA_STFLUX, ANA_SSFLUX, ANA_BSFLUX and ANA_BTFLUX (totally 4 ?) in .h file (I wonder how ana function set these 4 variables, just as zero automatically？).

2. Set NFFFILE=1 and input test_stress .nc file Yes, exactly. The way to tell what the default ...

Well, you could try submitting a bug report for this case of no forcing files, but most of us just set NFFILES to 1 and move on.


Kate

Thanks for your reply, you really help me out! So people usually input at least one forcing file, my case is not common.

By the way, where can I submit a bug ...

I wonder that if I just input sustrs and svstrs, must I set ANA_STFLUX, ANA_SSFLUX, ANA_BSFLUX and ANA_BTFLUX? Yes, if you provide the stresses, you need to provide the tracer fluxes too, and turn off BULK_FLUXES.


Kate,

Thanks. So the step is :

1. Undef BULK_FLUXES in .h file, and define ...

I wonder that if I just input sustrs and svstrs, must I set ANA_STFLUX, ANA_SSFLUX, ANA_BSFLUX and ANA_BTFLUX? Yes, if you provide the stresses, you need to provide the tracer fluxes too, and turn off BULK_FLUXES.


Kate,

Thanks. So the step is :

1. Undef BULK_FLUXES in .h file. and define ANA ...

What is on this line?
at /home/ubuntu/app/COAWST_V3.2/Build/get_data.f90:101
Oops, didn't see this:
Fortran runtime error: Index '1' of dimension 2 of array 'frc' above upper bound of 0
Did you tell it NFFILES before the list of forcing files?


Hi Kate,

Thanks for your reply.

I don't ...

What is on this line?
at /home/ubuntu/app/COAWST_V3.2/Build/get_data.f90:101
Oops, didn't see this:
Fortran runtime error: Index '1' of dimension 2 of array 'frc' above upper bound of 0
Did you tell it NFFILES before the list of forcing files?


Hi Kate,

Thanks for your reply.

I don't want ...

Hi all,

I have some problems about bulk_fluxes.

I used bulk_fluxes before in my old case and I input Uwind, Vwind, Pair, Tair..... at that time.

But now I try to run a new case. I just have sustrs and svstrs data.

I wonder that if I just input sustrs and svstrs, must I set ANA_STFLUX, ANA_SSFLUX ...

Model Input Parameters: ROMS/TOMS version 3.7
Friday - January 3, 2020 - 4:44:40 PM
-----------------------------------------------------------------------------

test

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/mpif90
Compiler flags ...

Hi all,

I tried to run a easy case and met a segmentation fault. So I turned on the USER_DEBUG, the results show belows:

I tried to figured out the problem but still have no idea yet.

Hope anyone have met this problem. Thanks lol

You want river inputs to be from land to water as in your left arrow. The right arrow appears to be all within the land and will therefore do nothing at all.

How much transport are you asking for? What is the corresponding flow rate to go with it? Are you violating the ROMS speed limit? I have ...

Are we looking at depth there? Where is the land mask? Are you sure you are putting the rivers where you think you are? Blows up how? Are you providing the river tracers (T and S)?
Dear Kate,

Thank you so much,

1) The depth here is all 10m (hmin).

2) The yellow part is water mask and the blue ...

Hi all,

I am a new user of ROMS and I am trying to run my test case.

There is a river flows into coastal ocean in my case.

And when I change the XI-position and ETA-position in rivers.nc, the result changes or even blows up.

When I set the postion where the left arrow points, it blows up. But ...