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

[Menemenlis, Dimitris (US 329B)]

As a sanity check, you could write out the surface stress to make sure you are forcing the ocean surface correctly.  D.

On Jul 12, 2024, at 2:24 PM, Chichi Zhou <chichi.zhou at uwaterloo.ca> wrote:

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

_______________________________________________
MITgcm-support mailing list
MITgcm-support at mitgcm.org<mailto:MITgcm-support at mitgcm.org>
https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!PyliTId6w07S3NcXZwgcfmy00WCbEdUh9_IFZM-eoThtGDgtzRCyyWrMwv-lS5Ev6NAbgIXqalyXAYKe9XPr4yZnV22llxQ1TTdS5Q$

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