Changeset 2591 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90
- Timestamp:
- 2011-02-18T14:26:28+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90
r2528 r2591 32 32 !! NEMO/LIM3 3.3 , UCL - NEMO Consortium (2010) 33 33 !! $Id$ 34 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)34 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 35 35 !!---------------------------------------------------------------------- 36 36 CONTAINS … … 87 87 !! (11-2006) Vectorized by Xavier Fettweis (UCL-ASTR) 88 88 !! (04-2007) Energy conservation tested by M. Vancoppenolle 89 !!90 89 !!------------------------------------------------------------------ 91 !! * Arguments92 93 90 INTEGER , INTENT (in) :: & 94 91 kideb , & ! Start point on which the the computation is applied … … 173 170 zftrice, & ! solar radiation transmitted through the ice 174 171 zihic, zhsu 175 172 !!------------------------------------------------------------------ 176 173 ! 177 174 !------------------------------------------------------------------------------! … … 382 379 DO ji = kideb , kiut 383 380 ztcond_i(ji,layer) = rcdic + zbeta*( s_i_b(ji,layer) & 384 + s_i_b(ji,layer+1) ) / MIN(- zeps, &381 + s_i_b(ji,layer+1) ) / MIN(-2.0*zeps, & 385 382 t_i_b(ji,layer)+t_i_b(ji,layer+1)-2.0*rtt) 386 383 ztcond_i(ji,layer) = MAX(ztcond_i(ji,layer),zkimin) … … 389 386 ENDIF 390 387 391 IF ( thcon_i_swi .EQ. 1 ) THEN ! Pringle 388 IF ( thcon_i_swi .EQ. 1 ) THEN ! Pringle (0.011/2=0.0055) 392 389 DO layer = 1, nlay_i-1 393 390 DO ji = kideb , kiut 394 391 ztcond_i(ji,layer) = rcdic + 0.09*( s_i_b(ji,layer) & 395 + s_i_b(ji,layer+1) ) / MIN(- zeps, &392 + s_i_b(ji,layer+1) ) / MIN(-2.0*zeps, & 396 393 t_i_b(ji,layer)+t_i_b(ji,layer+1)-2.0*rtt) - & 397 0.0 11* ( t_i_b(ji,layer) + t_i_b(ji,layer+1) - 2.0*rtt )394 0.0055* ( t_i_b(ji,layer) + t_i_b(ji,layer+1) - 2.0*rtt ) 398 395 ztcond_i(ji,layer) = MAX(ztcond_i(ji,layer),zkimin) 399 396 END DO
Note: See TracChangeset
for help on using the changeset viewer.