Changeset 3909 for trunk/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

Timestamp:

2013-06-11T13:01:03+02:00 (11 years ago)

Author:

davestorkey

Message:

Bug fixes for BDY module - see ticket #1106.

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90 (modified) (1 diff)

trunk/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

@@ -1053 +1053 @@
       
       bdytmask(:,:) = tmask(:,:,1)
+      IF( .not. ln_mask_file ) THEN
+         ! If .not. ln_mask_file then we need to derive mask on U and V grid 
+         ! from mask on T grid here.
+         bdyumask(:,:) = 0.e0
+         bdyvmask(:,:) = 0.e0
+         DO ij=1, jpjm1
+            DO ii=1, jpim1
+               bdyumask(ii,ij)=bdytmask(ii,ij)*bdytmask(ii+1, ij )
+               bdyvmask(ii,ij)=bdytmask(ii,ij)*bdytmask(ii  ,ij+1)  
+            END DO
+         END DO
+         CALL lbc_lnk( bdyumask(:,:), 'U', 1. )   ;   CALL lbc_lnk( bdyvmask(:,:), 'V', 1. )      ! Lateral boundary cond.
+      ENDIF
 
       ! bdy masks and bmask are now set to zero on boundary points: