How to setup wind forcing?

[dhtjtks]

Hi everyone!

I just want to put wind stress. It is same as wind-basin test case. So I turned off every other forcing terms except wind stress term. And I gave 1dyne/cm2 (0.1 Pa) in 'ana_smflux.h'. I needed to do unit conversion [dyne/cm2] to [m2/s2] unit. So I put 0.0001 m2/s2 which is equal to 1 dyne/cm2. Then I ran ROMS. The results show me more than 3 m/s currents at surface. I just used simple rectangle basin. I think the result is wrong. When I did this problem with 2D it worked well but with 3D it blew up. What do you think about this? Do you have any suggestion? Below is my cpp options.

#define SOLVE3D
#define ANA_GRID
#define ANA_INITIAL
#define ANA_SMFLUX
#define UV_QDRAG
#define ANA_STFLUX
#define ANA_BTFLUX
#define SPLINES

#define EASTERN_WALL
#define WESTERN_WALL
#define NORTHERN_WALL
#define SOUTHERN_WALL

#define AVERAGES

[kate]

You will probably need to add some form of vertical mixing to mix the momentum downwards. In the big basin case, we were applying the winds over the upper hundreds of meters as a body force, before the days of KPP or GLS or the like.