- Timestamp:
- 2019-02-14T16:55:25+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/r8395_India_uncoupled/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90
r10684 r10685 167 167 ! Check and write out namelist parameters 168 168 ! ----------------------------------------- 169 170 169 ! IF( jperio /= 0 ) CALL ctl_stop( 'bdy_segs: Cyclic or symmetric,', & 170 ! & ' and general open boundary condition are not compatible' ) 171 171 172 172 IF( nb_bdy == 0 ) THEN … … 1161 1161 bdytmask(:,:) = ssmask(:,:) 1162 1162 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) 1176 1170 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 1194 1174 ! bdy masks are now set to zero on boundary points: 1195 1175 !
Note: See TracChangeset
for help on using the changeset viewer.