- Timestamp:
- 2015-06-22T16:40:58+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/restart_datestamp/NEMOGCM/NEMO/LIM_SRC_3/limthd_ent.F90
r5420 r5462 25 25 USE sbc_oce ! Surface boundary condition: ocean fields 26 26 USE ice ! LIM variables 27 USE par_ice ! LIM parameters28 27 USE thd_ice ! LIM thermodynamics 29 28 USE limvar ! LIM variables … … 87 86 88 87 !-------------------------------------------------------------------------- 89 ! 1) Cumulative integral of old enthalpy * thic nkess and layers interfaces88 ! 1) Cumulative integral of old enthalpy * thickness and layers interfaces 90 89 !-------------------------------------------------------------------------- 91 90 zqh_cum0(:,0:nlay_i+2) = 0._wp … … 103 102 ! new layer thickesses 104 103 DO ji = kideb, kiut 105 zhnew(ji) = SUM( h_i_old(ji,0:nlay_i+1) ) / REAL( nlay_i )104 zhnew(ji) = SUM( h_i_old(ji,0:nlay_i+1) ) * r1_nlay_i 106 105 ENDDO 107 106 … … 133 132 DO jk1 = 1, nlay_i 134 133 DO ji = kideb, kiut 135 rswitch = 1._wp - MAX( 0._wp , SIGN( 1._wp , - zhnew(ji) + epsi10 ) )136 qnew(ji,jk1) = rswitch * ( zqh_cum1(ji,jk1) - zqh_cum1(ji,jk1-1) ) / MAX( zhnew(ji), epsi 10 )134 rswitch = MAX( 0._wp , SIGN( 1._wp , zhnew(ji) - epsi20 ) ) 135 qnew(ji,jk1) = rswitch * ( zqh_cum1(ji,jk1) - zqh_cum1(ji,jk1-1) ) / MAX( zhnew(ji), epsi20 ) 137 136 ENDDO 138 137 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.