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

Changeset 14603


Ignore:
Timestamp:
2021-03-09T12:51:51+01:00 (3 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.4_mppmask/src/OCE/LBC/mppini.F90

    r14075 r14603  
    11701170      ! 
    11711171      INTEGER                           ::   inumsave                ! local logical unit 
    1172       REAL(wp), DIMENSION(jpiglo,kjcnt) ::   zbot, zbdy  
     1172      REAL(wp), DIMENSION(jpiglo,kjcnt) ::   zmsk, zbdy  
    11731173      !!---------------------------------------------------------------------- 
    11741174      ! 
     
    11761176      ! 
    11771177      IF( numbot /= -1 ) THEN 
    1178          CALL iom_get( numbot, jpdom_unknown, 'bottom_level', zbot, kstart = (/1,kjstr/), kcount = (/jpiglo, kjcnt/) ) 
     1178         CALL iom_get( numbot, jpdom_unknown, 'mppmask', zmsk, kstart = (/1,kjstr/), kcount = (/jpiglo, kjcnt/) ) 
    11791179      ELSE 
    1180          zbot(:,:) = 1.                         ! put a non-null value 
     1180         zmsk(:,:) = 1.                         ! put a non-null value 
    11811181      ENDIF 
    11821182 
    11831183       IF( numbdy /= -1 ) THEN                  ! Adjust with bdy_msk if it exists     
    11841184         CALL iom_get ( numbdy, jpdom_unknown, 'bdy_msk', zbdy, kstart = (/1,kjstr/), kcount = (/jpiglo, kjcnt/) ) 
    1185          zbot(:,:) = zbot(:,:) * zbdy(:,:) 
    1186       ENDIF 
    1187       ! 
    1188       ldoce(:,:) = zbot(:,:) > 0. 
     1185         zmsk(:,:) = zmsk(:,:) * zbdy(:,:) 
     1186      ENDIF 
     1187      ! 
     1188      ldoce(:,:) = zmsk(:,:) > 0. 
    11891189      numout = inumsave 
    11901190      ! 
Note: See TracChangeset for help on using the changeset viewer.