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 2591 – NEMO

Changeset 2591


Ignore:
Timestamp:
2011-02-18T14:26:28+01:00 (13 years ago)
Author:
gm
Message:

v3.3: #799 , correct a LIM3 small bug (limthd_dif)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

    r2528 r2591  
    3232   !! NEMO/LIM3 3.3 , UCL - NEMO Consortium (2010) 
    3333   !! $Id$ 
    34    !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
     34   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3535   !!---------------------------------------------------------------------- 
    3636CONTAINS 
     
    8787      !!           (11-2006) Vectorized by Xavier Fettweis (UCL-ASTR) 
    8888      !!           (04-2007) Energy conservation tested by M. Vancoppenolle 
    89       !! 
    9089      !!------------------------------------------------------------------ 
    91       !! * Arguments 
    92  
    9390      INTEGER , INTENT (in) ::  & 
    9491         kideb ,  &  ! Start point on which the  the computation is applied 
     
    173170         zftrice,            &  ! solar radiation transmitted through the ice 
    174171         zihic, zhsu 
    175  
     172      !!------------------------------------------------------------------ 
    176173      ! 
    177174      !------------------------------------------------------------------------------! 
     
    382379               DO ji = kideb , kiut 
    383380                  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,     & 
    385382                     t_i_b(ji,layer)+t_i_b(ji,layer+1)-2.0*rtt) 
    386383                  ztcond_i(ji,layer)   = MAX(ztcond_i(ji,layer),zkimin) 
     
    389386         ENDIF 
    390387 
    391          IF ( thcon_i_swi .EQ. 1 ) THEN ! Pringle 
     388         IF ( thcon_i_swi .EQ. 1 ) THEN ! Pringle (0.011/2=0.0055) 
    392389            DO layer = 1, nlay_i-1 
    393390               DO ji = kideb , kiut 
    394391                  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,      & 
    396393                     t_i_b(ji,layer)+t_i_b(ji,layer+1)-2.0*rtt) - & 
    397                      0.011* ( 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 )   
    398395                  ztcond_i(ji,layer) = MAX(ztcond_i(ji,layer),zkimin) 
    399396               END DO 
Note: See TracChangeset for help on using the changeset viewer.