- Timestamp:
- 2017-09-25T21:11:19+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd_ent.F90
r8534 r8562 78 78 ! 1) Cumulative integral of old enthalpy * thickness and layers interfaces 79 79 !-------------------------------------------------------------------------- 80 zeh_cum0( :,0:nlay_i+2) = 0._wp81 zh_cum0 ( :,0:nlay_i+2) = 0._wp80 zeh_cum0(1:nidx,0) = 0._wp 81 zh_cum0 (1:nidx,0) = 0._wp 82 82 DO jk0 = 1, nlay_i+2 83 83 DO ji = 1, nidx … … 96 96 97 97 ! new layers interfaces 98 zh_cum1( :,0:nlay_i) = 0._wp98 zh_cum1(1:nidx,0) = 0._wp 99 99 DO jk1 = 1, nlay_i 100 100 DO ji = 1, nidx … … 103 103 END DO 104 104 105 zeh_cum1( :,0:nlay_i) = 0._wp105 zeh_cum1(1:nidx,0) = 0._wp 106 106 ! new cumulative q*h => linear interpolation 107 107 DO jk0 = 1, nlay_i+1 … … 117 117 END DO 118 118 ! 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) 120 120 121 121 ! new enthalpies
Note: See TracChangeset
for help on using the changeset viewer.