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 13126 – NEMO

Changeset 13126


Ignore:
Timestamp:
2020-06-17T18:05:30+02:00 (4 years ago)
Author:
cguiavarch
Message:

Read mppmask from domain_cfg instead of bottom_level

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.2_mppmask/src/OCE/LBC/mppini.F90

    r12658 r13126  
    11491149      ! 
    11501150      INTEGER                           ::   inumsave                ! local logical unit 
    1151       REAL(wp), DIMENSION(jpiglo,kjcnt) ::   zbot, zbdy  
     1151      REAL(wp), DIMENSION(jpiglo,kjcnt) ::   zmsk, zbdy  
    11521152      !!---------------------------------------------------------------------- 
    11531153      ! 
     
    11551155      ! 
    11561156      IF( numbot /= -1 ) THEN 
    1157          CALL iom_get( numbot, jpdom_unknown, 'bottom_level', zbot, kstart = (/1,kjstr/), kcount = (/jpiglo, kjcnt/) ) 
     1157         CALL iom_get( numbot, jpdom_unknown, 'mppmask', zmsk, kstart = (/1,kjstr/), kcount = (/jpiglo, kjcnt/) ) 
    11581158      ELSE 
    1159          zbot(:,:) = 1.                         ! put a non-null value 
     1159         zmsk(:,:) = 1.                         ! put a non-null value 
    11601160      ENDIF 
    11611161 
    11621162       IF( numbdy /= -1 ) THEN                  ! Adjust with bdy_msk if it exists     
    11631163         CALL iom_get ( numbdy, jpdom_unknown, 'bdy_msk', zbdy, kstart = (/1,kjstr/), kcount = (/jpiglo, kjcnt/) ) 
    1164          zbot(:,:) = zbot(:,:) * zbdy(:,:) 
    1165       ENDIF 
    1166       ! 
    1167       ldoce(:,:) = zbot(:,:) > 0. 
     1164         zmsk(:,:) = zmsk(:,:) * zbdy(:,:) 
     1165      ENDIF 
     1166      ! 
     1167      ldoce(:,:) = zmsk(:,:) > 0. 
    11681168      numout = inumsave 
    11691169      ! 
Note: See TracChangeset for help on using the changeset viewer.