New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10685 for branches/UKMO/r8395_India_uncoupled/NEMOGCM/NEMO/OPA_SRC/BDY – NEMO

Ignore:
Timestamp:
2019-02-14T16:55:25+01:00 (5 years ago)
Author:
jcastill
Message:

Changes as Ash's files

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/r8395_India_uncoupled/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

    r10684 r10685  
    167167      ! Check and write out namelist parameters 
    168168      ! ----------------------------------------- 
    169       IF( jperio /= 0 )   CALL ctl_stop( 'bdy_segs: Cyclic or symmetric,',   & 
    170          &                               ' and general open boundary condition are not compatible' ) 
     169!     IF( jperio /= 0 )   CALL ctl_stop( 'bdy_segs: Cyclic or symmetric,',   & 
     170!        &                               ' and general open boundary condition are not compatible' ) 
    171171 
    172172      IF( nb_bdy == 0 ) THEN  
     
    11611161      bdytmask(:,:) = ssmask(:,:) 
    11621162 
    1163       IF( ln_mask_file ) THEN 
    1164          CALL iom_open( cn_mask_file, inum ) 
    1165          CALL iom_get ( inum, jpdom_data, 'bdy_msk', bdytmask(:,:) ) 
    1166          CALL iom_close( inum ) 
    1167  
    1168          ! Derive mask on U and V grid from mask on T grid 
    1169          bdyumask(:,:) = 0._wp 
    1170          bdyvmask(:,:) = 0._wp 
    1171          DO ij=1, jpjm1 
    1172             DO ii=1, jpim1 
    1173                bdyumask(ii,ij) = bdytmask(ii,ij) * bdytmask(ii+1, ij ) 
    1174                bdyvmask(ii,ij) = bdytmask(ii,ij) * bdytmask(ii  ,ij+1)   
    1175             END DO 
     1163      ! we need to derive mask on U and V grid from mask on T grid here. 
     1164      bdyumask(:,:) = 0._wp 
     1165      bdyvmask(:,:) = 0._wp 
     1166      DO ij = 1, jpjm1 
     1167         DO ii = 1, jpim1 
     1168            bdyumask(ii,ij) = bdytmask(ii,ij) * bdytmask(ii+1, ij ) 
     1169            bdyvmask(ii,ij) = bdytmask(ii,ij) * bdytmask(ii  ,ij+1)   
    11761170         END DO 
    1177          CALL lbc_lnk( bdyumask(:,:), 'U', 1. )   ;   CALL lbc_lnk( bdyvmask(:,:), 'V', 1. )      ! Lateral boundary cond. 
    1178  
    1179       ENDIF ! ln_mask_file=.TRUE. 
    1180        
    1181       IF( .NOT.ln_mask_file ) THEN 
    1182          ! If .not. ln_mask_file then we need to derive mask on U and V grid from mask on T grid here. 
    1183          bdyumask(:,:) = 0._wp 
    1184          bdyvmask(:,:) = 0._wp 
    1185          DO ij = 1, jpjm1 
    1186             DO ii = 1, jpim1 
    1187                bdyumask(ii,ij) = bdytmask(ii,ij) * bdytmask(ii+1, ij ) 
    1188                bdyvmask(ii,ij) = bdytmask(ii,ij) * bdytmask(ii  ,ij+1)   
    1189             END DO 
    1190          END DO 
    1191          CALL lbc_lnk( bdyumask(:,:), 'U', 1. )   ;   CALL lbc_lnk( bdyvmask(:,:), 'V', 1. )      ! Lateral boundary cond. 
    1192       ENDIF 
    1193  
     1171      END DO 
     1172      CALL lbc_lnk( bdyumask(:,:), 'U', 1. )   ;   CALL lbc_lnk( bdyvmask(:,:), 'V', 1. )      ! Lateral boundary cond. 
     1173      
    11941174      ! bdy masks are now set to zero on boundary points: 
    11951175      ! 
Note: See TracChangeset for help on using the changeset viewer.