Changeset 8090


Ignore:
Timestamp:
2017-05-30T09:59:39+02:00 (3 years ago)
Author:
cbricaud
Message:

adapt maskland option to CRS

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90

    r7806 r8090  
    13901390      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zmask 
    13911391      INTEGER  :: ni,nj 
     1392      !!---------------------------------------------------------------------- 
    13921393       
    13931394      ni=nlei-nldi+1 ; nj=nlej-nldj+1 
     
    14041405      IF ( ln_mskland ) THEN 
    14051406         ! 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 
    14121427         ! 
    14131428#if ! defined key_xios2 
Note: See TracChangeset for help on using the changeset viewer.