[MITgcm-devel] open boundaries for sea ice

[Dimitris Menemenlis]

Martin, I am getting more and more pressure from one of your former AWI 
colleagues to implement OBCS for pkg/seaice.  In seaice_advection.F I notice:

> CML#ifdef ALLOW_OBCS
> CMLC-     Apply open boundary conditions
> CML          IF ( useOBCS ) THEN
> CML           IF (tracerIdentity.EQ.GAD_HEFF) THEN
> CML            CALL OBCS_APPLY_HEFF( bi, bj, k, localTij, myThid )
> CML           ELSEIF (tracerIdentity.EQ.GAD_AREA) THEN
> CML            CALL OBCS_APPLY_AREA( bi, bj, k, localTij, myThid )
> CML           ENDIF
> CML          ENDIF
> CML#endif /* ALLOW_OBCS */

Two questions:

1. Why do you need to do above twice?  Wouldn't once, after both the X and Y 
direction loops have completed suffice?

2. Do you already have stubs available for OBCS_APPLY_HEFF and OBCS_APPLY_AREA?

D.

-- 
Dimitris Menemenlis <menemenlis at sbcglobal.net>
5056 Oakwood Ave, La Canada, CA 91011-2450
tel/fax: 818-790-6735;   cell: 818-625-6498