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

[Jean-Michel Campin]

Hi Chichi,

For temporally constant wind stress:
a) without pkg/exf: 
   wind-stress is supposed to be provided on C-grid ;
   and with periodicExternalForcing=F (i.e., the default), it's constant in time.
b) with pkg/exf: 
   need to specify grid location (either: readStressOnAgrid=T or readStressOnCgrid=T,
   in data.exf, EXF_NML_01);
   and setting "ustressperiod = 0." will result in constant in time wind-stress forcing.
Note that for this case, if using spacial interpolation, the code for constant wind-stress 
is a bit more complicated than for other exf forcing fields (see comments in exf_init_varia.F,
lines 65-69).
 
So, in your case, I would try to change ustressperiod setting to zero and check if the 
comparison improves.

Cheers,
Jean-Michel

On Fri, Jul 12, 2024 at 09:24:09PM +0000, Chichi Zhou 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
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support