Changeset 1813 for tags/nemo_v3_2/NEMO

Timestamp:

2010-03-25T12:22:36+01:00 (14 years ago)

Author:

rblod

Message:

Fixes for OBC, see ticket #648

Location:

tags/nemo_v3_2/NEMO/OPA_SRC/OBC

Files:

• obc_oce.F90 (modified) (1 diff)
• obcini.F90 (modified) (5 diffs)

tags/nemo_v3_2/NEMO/OPA_SRC/OBC/obc_oce.F90

@@ -27 +27 @@
    !
    !                                  !!* Namelist namobc: open boundary condition *
-   INTEGER           ::   nn_nbobc    = 2        !: number of open boundaries ( 1=< nbobc =< 4 ) 
    INTEGER           ::   nn_obcdta   = 0        !:  = 0 use the initial state as obc data
    !                                             !   = 1 read obc data in obcxxx.dta files

tags/nemo_v3_2/NEMO/OPA_SRC/OBC/obcini.F90

@@ -70 +70 @@
 
       ! convert DOCTOR namelist name into the OLD names
-      nbobc    = nn_nbobc
       nobc_dta = nn_obcdta
       cffile   = cn_obcdta
@@ -101 +100 @@
       IF(lwp) WRITE(numout,*) 'obc_init : initialization of open boundaries'
       IF(lwp) WRITE(numout,*) '~~~~~~~~'
-      IF(lwp) WRITE(numout,*) '   Number of open boundaries    nn_nbobc = ', nn_nbobc
+      IF(lwp) WRITE(numout,*) '   Number of open boundaries    nbobc = ', nbobc
       IF(lwp) WRITE(numout,*)
 
@@ -306 +305 @@
       IF( lp_obc_east ) THEN
          !... (jpjed,jpjefm1),jpieob
+         bmask(nie0p1:nie1p1,nje0:nje1m1) = 0.e0
 
          ! ... initilization to zero
@@ -341 +341 @@
       IF( lp_obc_north ) THEN
          ! ... jpjnob,(jpind,jpisfm1)
+         bmask(nin0:nin1m1,njn0p1:njn1p1) = 0.e0
 
          ! ... initilization to zero
@@ -440 +441 @@
             END DO
          END IF
-  
+
          IF( lp_obc_north ) THEN ! ... North open boundary lateral surface
             DO jj = njn0, njn1