Changeset 8090 for branches/2015/dev_r5003_MERCATOR6_CRS
- Timestamp:
- 2017-05-30T09:59:39+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90
r7806 r8090 1390 1390 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zmask 1391 1391 INTEGER :: ni,nj 1392 !!---------------------------------------------------------------------- 1392 1393 1393 1394 ni=nlei-nldi+1 ; nj=nlej-nldj+1 … … 1404 1405 IF ( ln_mskland ) THEN 1405 1406 ! mask land points, keep values on coast line -> specific mask for U, V and W points 1406 SELECT CASE ( cdgrd ) 1407 CASE('T') ; zmask(:,:,:) = tmask(:,:,:) 1408 CASE('U') ; zmask(2:jpim1,:,:) = tmask(2:jpim1,:,:) + tmask(3:jpi,:,:) ; CALL lbc_lnk( zmask, 'U', 1. ) 1409 CASE('V') ; zmask(:,2:jpjm1,:) = tmask(:,2:jpjm1,:) + tmask(:,3:jpj,:) ; CALL lbc_lnk( zmask, 'V', 1. ) 1410 CASE('W') ; zmask(:,:,2:jpk ) = tmask(:,:,1:jpkm1) + tmask(:,:,2:jpk) ; zmask(:,:,1) = tmask(:,:,1) 1411 END SELECT 1407 1408 IF( lk_crs .AND. jpi == jpi_crs .AND. jpi .NE. 0 )THEN 1409 1410 SELECT CASE ( cdgrd ) 1411 CASE('T') ; zmask(:,:,:) = tmask_crs(:,:,:) 1412 CASE('U') ; zmask(2:jpim1,:,:) = tmask_crs(2:jpim1,:,:) + tmask_crs(3:jpi,:,:) ; CALL crs_lbc_lnk( zmask, 'U', 1. ) 1413 CASE('V') ; zmask(:,2:jpjm1,:) = tmask_crs(:,2:jpjm1,:) + tmask_crs(:,3:jpj,:) ; CALL crs_lbc_lnk( zmask, 'V', 1. ) 1414 CASE('W') ; zmask(:,:,2:jpk ) = tmask_crs(:,:,1:jpkm1) + tmask_crs(:,:,2:jpk) ; zmask(:,:,1) = tmask_crs(:,:,1) 1415 END SELECT 1416 1417 ELSE 1418 1419 SELECT CASE ( cdgrd ) 1420 CASE('T') ; zmask(:,:,:) = tmask(:,:,:) 1421 CASE('U') ; zmask(2:jpim1,:,:) = tmask(2:jpim1,:,:) + tmask(3:jpi,:,:) ; CALL lbc_lnk( zmask, 'U', 1. ) 1422 CASE('V') ; zmask(:,2:jpjm1,:) = tmask(:,2:jpjm1,:) + tmask(:,3:jpj,:) ; CALL lbc_lnk( zmask, 'V', 1. ) 1423 CASE('W') ; zmask(:,:,2:jpk ) = tmask(:,:,1:jpkm1) + tmask(:,:,2:jpk) ; zmask(:,:,1) = tmask(:,:,1) 1424 END SELECT 1425 1426 ENDIF 1412 1427 ! 1413 1428 #if ! defined key_xios2
Note: See TracChangeset
for help on using the changeset viewer.