[Mitgcm-support] obcs in r1b1

[mitgcm-support at dev.mitgcm.org]

Hi,

try to use obcs in r1b1 and encountered a problem.
When I run verification experiment 4 with the following data.obcs 
everythings works fine:

# Open-boundaries
  &OBCS_PARM01
  OB_Jnorth=80*-1,
  OB_Jsouth=80*1,
  OB_Ieast=42*-1,
  OB_Iwest=42*1,
#useOrlanskiNorth=.TRUE.,
#useOrlanskiSouth=.TRUE.,
#useOrlanskiEast=.TRUE.,
#useOrlanskiWest=.FALSE.,
  &

But switching to pure Orlanski radiation conditions, i.e.:

# Open-boundaries
  &OBCS_PARM01
  OB_Jnorth=80*-1,
  OB_Jsouth=80*1,
  OB_Ieast=42*-1,
  OB_Iwest=42*1,
  useOrlanskiNorth=.TRUE.,
  useOrlanskiSouth=.TRUE.,
  useOrlanskiEast=.TRUE.,
  useOrlanskiWest=.FALSE.,


causes a core dump (floating divide by zero).
Does that mean Orlanski radiation conditions are not working? Or am I wrong?

-Arne



-- 

Dr. Arne Biastoch

Institute for Marine Research         phone: ++49 431 600-4013
FB1 Ocean Circulation and Climate     fax  : ++49 431 600-1515
Duesternbrooker Weg 20                email: abiastoch at ifm.uni-kiel.de
24105 Kiel, Germany

http://www.ifm.uni-kiel.de/fb/fb1/tm/data/pers/abiastoch.html