[MITgcm-support] body force

[Dimitris Menemenlis]

Not sure if relevant to your problem — I sometimes answer too quickly and then
Jean-Michel tells me that I answered the wrong question — but have you looked
at MITgcm/pkg/obcs/obcs_add_tides.F ?

Examples how to use are here:
MITgcm/verification/seaice_obcs/input.tides

and here: MITgcm_contrib/tides

> On Nov 4, 2015, at 2:22 PM, guopu at scsio.ac.cn wrote:
> 
> Hello everyone,
> 
> I am trying to give tidal forcing to a 2D domain to generate internal tides. Beside prescribing currents on the open boundary, one way I found is to add forcing term to 'gU' in file 'external_forcing.F', which has been mentioned by past mails. So I add these codes into 'external_forcing.F' to generate M2 tide:
> 
>   Do j=0, sNy+1
> 
>     Do I=1, sNx+1
> 
>         gU(i, j, kLev, bi, bj) = gU(i, j, kLev, bi, bj)
> 
> &          +0.06*0.00014*cos(0.00014*myTime)
> 
>     ENDDO
> 
>   ENDDO
> 
> I use this file to force a rest, plat (2500 m depth and no walls at east and west boundaries) 2D domain. I thought the velocity all over the domain would vary sinusoidally with amplitude 0.06 m/s. However, the variation is really sinusoidal but the amplitude is varied.
> 
> Who can tell me what is my fault? By the way, in some past mails someone mentioned that RBCS package can give tidal forcing to the domain. Who know how to do this? or is it better or easier than body force?  I appreciate any help. Thanks.
> 
> Best,
> 
> Pu
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20151104/60d7984a/attachment.htm>