bustr and bvstr for BODYFORCE

[antoinemarie]

Hi Roms users!

I try to use BODYFORCE, but I have some questions:

- In witch file can I define my force bustr and bvstr?
- When I define a force in the file 'ana_smflux', is it sustr and svstr for each layer, or is it the total of all layers?

Thank you in advance for your reply,

Antoine

[kate]

Actually, get_data.F has code for reading sustr and svstr, but not for reading bustr or bvstr. It is assumed that you will be applying linear, quadratic, or a log layer bottom stress. If that's not the case, you can modify get_data.F and set_data.F, copying the sustr and svstr code. At least bottom stresses are known by varinfo.dat (idUbms and idVbms). If you do this, put the fields in the forcing file(s).

You want to provide the total values, not the per-level value. You can see the code in rhs3d.F where the thicknesses are summed, then divided out when applied via BODYFORCE.

[antoinemarie]

Thank you Kate for your reply.

I have another question about BODYFORCE. I use it, but it doesn't repect my numeber of level: LEVSFRC == 2, the force is applied on the entire water column. In my header file, I just define BODYFORCE. Is there another cpp to define?

Thanks,

Antoine

[kate]

Setting this "LEVSFRC == 2" means apply the surface stress to every level except the bottom. The surface is applied from levsfrc to N, while the bottom is applied from 1 to levbfrc. Does that clear up what you are getting?

[antoinemarie]

ok, it works. The problem cames from my understanding of the numbering of levels.

Thank you Kate