[MITgcm-support] Setting up external wind stress in ocean model

[Chichi Zhou]

Hello everyone,

I'm conducting a 2-D simulation in the x-z plane and would like to add spatially and temporally constant wind stress 0.3 N/m^2 at ocean surface. I have tried two methods:

(1) Specify the wind file in `input/data` and do not use the EXF package.

&PARM05
hydrogThetaFile='T.init',
bathyFile='flat_bottom',
uVelInitFile='U.init',
zonalWindFile='windx.bin',

Here, 'windx.bin' is a 1-D array with dimensions (nx, 1), where windx (1:nx, 1) = 0.3. ‘nx’ is the total number of grid points in the x direction.

(2) Set up the wind stress in `input/data.exf`.

&EXF_NML_02
ustressfile         = 'windx.bin',
ustressperiod       = 0.2,
&

Here  'windx.bin' is converted from a matrix  with dimensions (nx, 1, nt).  windx (1:nx, 1, 1:nt) = 0.3. `nx` is the grid points in the x direction, and ‘nt’ is the total time steps. ‘ustressperiod = 0.2’ because `dt = 0.2s`.

However, the output velocity and temperature fields from the two numerical simulations are different. I am confused about the reason of this and would appreciate any suggestions you might have. Thank you in advance!

Best,
Chichi

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20240712/bdb1d15d/attachment.html>