- Timestamp:
- 2015-02-11T16:32:03+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90
r5078 r5079 138 138 DO jj = 1, jpj 139 139 DO ji = 1, jpi 140 z2d(ji,jj) = ( tm_i(ji,jj) - rt t) * zswi(ji,jj)140 z2d(ji,jj) = ( tm_i(ji,jj) - rt0 ) * zswi(ji,jj) 141 141 END DO 142 142 END DO … … 149 149 DO jj = 1, jpj 150 150 DO ji = 1, jpi 151 z2d(ji,jj) = z2d(ji,jj) + zswi(ji,jj) * ( t_su(ji,jj,jl) - rt t) * a_i(ji,jj,jl) / MAX( at_i(ji,jj) , epsi06 )151 z2d(ji,jj) = z2d(ji,jj) + zswi(ji,jj) * ( t_su(ji,jj,jl) - rt0 ) * a_i(ji,jj,jl) / MAX( at_i(ji,jj) , epsi06 ) 152 152 END DO 153 153 END DO … … 265 265 rswitch = MAX( 0._wp , SIGN( 1._wp , a_i(ji,jj,jl) - epsi06 ) ) 266 266 zei(ji,jj,jl) = zei(ji,jj,jl) + 100.0* & 267 ( - tmut * s_i(ji,jj,jk,jl) / MIN( ( t_i(ji,jj,jk,jl) - rt t), - epsi06 ) ) * &267 ( - tmut * s_i(ji,jj,jk,jl) / MIN( ( t_i(ji,jj,jk,jl) - rt0 ), - epsi06 ) ) * & 268 268 rswitch * r1_nlay_i 269 269 END DO … … 352 352 CALL histwrite( kid, "iicethic", kt, icethi , jpi*jpj, (/1/) ) 353 353 CALL histwrite( kid, "iiceconc", kt, at_i , jpi*jpj, (/1/) ) 354 CALL histwrite( kid, "iicetemp", kt, tm_i - rt t, jpi*jpj, (/1/) )354 CALL histwrite( kid, "iicetemp", kt, tm_i - rt0 , jpi*jpj, (/1/) ) 355 355 CALL histwrite( kid, "iicevelu", kt, u_ice , jpi*jpj, (/1/) ) 356 356 CALL histwrite( kid, "iicevelv", kt, v_ice , jpi*jpj, (/1/) )
Note: See TracChangeset
for help on using the changeset viewer.