Ignore:
Timestamp:
2018-05-18T11:52:45+02:00 (2 years ago)
Author:
frrh
Message:

Introduce and apply a 3D interior mask for T, U, V and W grid diagnostics.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_diag_bitcomp/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90

    r9321 r9608  
    14011401         ! mask land points, keep values on coast line -> specific mask for U, V and W points 
    14021402         SELECT CASE ( cdgrd ) 
    1403          CASE('T')   ;   zmask(:,:,:)       = tmask(:,:,:) 
    1404          CASE('U')   ;   zmask(2:jpim1,:,:) = tmask(2:jpim1,:,:) + tmask(3:jpi,:,:)   ;   CALL lbc_lnk( zmask, 'U', 1. ) 
    1405          CASE('V')   ;   zmask(:,2:jpjm1,:) = tmask(:,2:jpjm1,:) + tmask(:,3:jpj,:)   ;   CALL lbc_lnk( zmask, 'V', 1. ) 
    1406          CASE('W')   ;   zmask(:,:,2:jpk  ) = tmask(:,:,1:jpkm1) + tmask(:,:,2:jpk)   ;   zmask(:,:,1) = tmask(:,:,1) 
     1403         CASE('T')    
     1404            zmask(:,:,:)       = tmask_i_3d(:,:,:) 
     1405         CASE('U') 
     1406            zmask(2:jpim1,:,:) = tmask_i_3d(2:jpim1,:,:) + tmask_i_3d(3:jpi,:,:)    
     1407            CALL lbc_lnk( zmask, 'U', 1. ) 
     1408         CASE('V')  
     1409            zmask(:,2:jpjm1,:) = tmask_i_3d(:,2:jpjm1,:) + tmask_i_3d(:,3:jpj,:)    
     1410            CALL lbc_lnk( zmask, 'V', 1. ) 
     1411         CASE('W')    
     1412            zmask(:,:,2:jpk  ) = tmask_i_3d(:,:,1:jpkm1) + tmask_i_3d(:,:,2:jpk)    
     1413            zmask(:,:,1) = tmask_i_3d(:,:,1) 
    14071414         END SELECT 
    14081415         ! 
Note: See TracChangeset for help on using the changeset viewer.