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 8562 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd_ent.F90 – NEMO

Ignore:
Timestamp:
2017-09-25T21:11:19+02:00 (7 years ago)
Author:
clem
Message:

cosmetics only

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd_ent.F90

    r8534 r8562  
    7878      !  1) Cumulative integral of old enthalpy * thickness and layers interfaces 
    7979      !-------------------------------------------------------------------------- 
    80       zeh_cum0(:,0:nlay_i+2) = 0._wp  
    81       zh_cum0 (:,0:nlay_i+2) = 0._wp 
     80      zeh_cum0(1:nidx,0) = 0._wp  
     81      zh_cum0 (1:nidx,0) = 0._wp 
    8282      DO jk0 = 1, nlay_i+2 
    8383         DO ji = 1, nidx 
     
    9696 
    9797      ! new layers interfaces 
    98       zh_cum1(:,0:nlay_i) = 0._wp 
     98      zh_cum1(1:nidx,0) = 0._wp 
    9999      DO jk1 = 1, nlay_i 
    100100         DO ji = 1, nidx 
     
    103103      END DO 
    104104 
    105       zeh_cum1(:,0:nlay_i) = 0._wp  
     105      zeh_cum1(1:nidx,0) = 0._wp  
    106106      ! new cumulative q*h => linear interpolation 
    107107      DO jk0 = 1, nlay_i+1 
     
    117117      END DO 
    118118      ! to ensure that total heat content is strictly conserved, set: 
    119       zeh_cum1(:,nlay_i) = zeh_cum0(:,nlay_i+2)  
     119      zeh_cum1(1:nidx,nlay_i) = zeh_cum0(1:nidx,nlay_i+2)  
    120120 
    121121      ! new enthalpies 
Note: See TracChangeset for help on using the changeset viewer.