[MITgcm-support] orlanski open boundary conditions

[Alistair Adcroft]

Do you really mean "only North and East" has problems? Clearly a bug if so -
make sure your initial conditions are symmetric and beta=0. If the OBCS are
assymetric then lets find the bug.

A.
--
Dr Alistair Adcroft            http://www.mit.edu/~adcroft
MIT Climate Modeling Initiative        tel: (617) 253-5938
EAPS 54-1523,  77 Massachusetts Ave,  Cambridge,  MA,  USA

-----Original Message-----
From: mitgcm-support-bounces at mitgcm.org
[mailto:mitgcm-support-bounces at mitgcm.org] On Behalf Of Martin Losch
Sent: Tuesday, December 02, 2003 7:33 AM
To: MITgcm Support
Cc: Samar Khatiwala
Subject: [MITgcm-support] orlanski open boundary conditions


Hi,

is there any experience on how well the open boudary conditions with 
orlanski radiation  are working? I would like to study the 
relaxation/spin-down of a solitary geostrophic eddy and from my initial 
adjustment I get all sorts of odd processes near the open boundaries 
(only "north" and "east"). I had hoped that the gravity waves would be 
more or less radiated out of the domain, which they mostly do, but 
there are strong vertical velocities near the open boundaries and 
eventually, the (linear free) sea-surface height is affected (its mean 
goes down by meters) and the solution diverges (explodes).

Any experience with this? Any ideas?

Martin

PS. this is my data.obcs
  &OBCS_PARM01
  OB_Jnorth=100*-1,
  OB_Jsouth=100*1,
  OB_Ieast=100*-1,
  OB_Iwest=100*1,
  useOrlanskiNorth=.TRUE.,
  useOrlanskiSouth=.TRUE.,
  useOrlanskiEast =.TRUE.,
  useOrlanskiWest =.TRUE.,
  &

  &OBCS_PARM02
  cmax=0.45D0,
  cvelTimeScale=500.D0,
  &

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