[MITgcm-support] Tidal forcing with Orlanski boundary conditions.

[Jean-Michel Campin]

Hi Hanut,

Subroutines in external_forcing.F have been retired since Aug 2014
and replaced with the ones in apply_forcing.F
They are just kept in case someone has a customized version and want
to keep it (and in this case, need to compile with #define USE_OLD_EXTERNAL_FORCING
in CPP_OPTIONS.h). Are you in this case ?

Cheers,
Jean-Michel

On Tue, May 23, 2017 at 10:42:29PM +0530, Hanut Vemulapalli wrote:
> Hi all,
> 
> I am trying to simulate gravity waves produced by a gaussian topography in
> uniform stratificaiton. I have a 2D simulation with unit depth in the
> y-coordinate. I am giving tidal forcing to the model in
> code/external_forcing.F by adding the following lines .
> 
> DO j=0,sNy+1
>   DO i=1,sNx+1
>           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj) +
> 0.02*0.00014*cos(0.00014*myTime)
>   ENDDO
> ENDDO
> 
> 
> I am prescribing Orlanski boundary condition at the east and west
> boundaries.
> However, the velocity magnitude tends to very small values of the order
> 1E-5 times the tidal amplitude (i.e 1E-12) once I prescribe open boundary
> condition at the east and west boundaries.
> Can anyone tell me where I am going wrong.
> 
> Thanks in advance.
> 
> Best regards
> Hanut
> 
> 
> 
> 
>   <https://mailtrack.io/> Sent with Mailtrack
> <https://mailtrack.io/install?source=signature&lang=en&referral=hanut2@gmail.com&idSignature=22>

> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support