Changeset 1046


Ignore:
Timestamp:
2008-06-02T12:24:22+02:00 (12 years ago)
Author:
cetlod
Message:

local computation of freezing point in trcadv_cen2.F90 routine, see ticket:181

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/TRP/trcadv_cen2.F90

    r941 r1046  
    137137         zfui1, zfvj1                          !    "         " 
    138138#endif 
     139#if defined key_lim3 || defined key_lim2 
     140      REAL(wp) ::                           & 
     141         ztfreez                               ! freezing point 
     142#endif 
    139143      CHARACTER (len=22) :: charout 
    140144      !!---------------------------------------------------------------------- 
     
    166170         DO jj = 1, jpj 
    167171            DO ji = 1, jpi 
     172#if defined key_lim3 || defined key_lim2 
     173               ztfreez = ( - 0.0575 + 1.710523e-3 * SQRT( sn(ji,jj,1) )   & 
     174                 &                  - 2.154996e-4 *       sn(ji,jj,1)   ) * sn(ji,jj,1) 
     175 
    168176               zind(ji,jj,jk) = MAX (   & 
    169177                  rnfmsk(ji,jj) * rnfmsk_z(jk),      &  ! near runoff mouths (& closed sea outflows) 
    170178                  upsmsk(ji,jj)                      &  ! some of some straits 
    171 #if defined key_lim3 || defined key_lim2 
    172179                  !                                     ! below ice covered area (if tn < "freezing"+0.1 ) 
    173                   , MAX(  0., SIGN( 1., fzptn(ji,jj) + 0.1 - tn(ji,jj,jk) )  ) * tmask(ji,jj,jk)   & 
    174 #endif 
     180                  , MAX(  0., SIGN( 1., ztfreez + 0.1 - tn(ji,jj,jk) )  ) * tmask(ji,jj,jk)   & 
    175181                  &                  ) 
     182 
     183#else 
     184               zind(ji,jj,jk) = MAX (   & 
     185                  rnfmsk(ji,jj) * rnfmsk_z(jk),      &  ! near runoff mouths (& closed sea outflows) 
     186                  upsmsk(ji,jj)                      &  ! some of some straits 
     187                  &                  ) 
     188#endif 
    176189            END DO 
    177190         END DO 
Note: See TracChangeset for help on using the changeset viewer.