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

Changeset 15046


Ignore:
Timestamp:
2021-06-23T12:46:01+02:00 (3 years ago)
Author:
clem
Message:

trunk: bug found by ecmwf: do not do ice remapping if, by some mystery and in very rare occasion, hR=hL=hice.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/trunk/src/ICE/iceitd.F90

    r14997 r15046  
    378378            ! 
    379379            ! Compute coefficients of g(eta) = g0 + g1*eta 
    380             zdhr = 1._wp / (phR(ji) - phL(ji)) 
     380            IF( phR(ji) > phL(ji) ) THEN   ;   zdhr = 1._wp / (phR(ji) - phL(ji)) 
     381            ELSE                           ;   zdhr = 0._wp ! if hR=hL=hice => no remapping 
     382            ENDIF 
     383            !!zdhr = 1._wp / (phR(ji) - phL(ji)) 
    381384            zwk1 = 6._wp * paice(ji) * zdhr 
    382385            zwk2 = ( phice(ji) - phL(ji) ) * zdhr 
Note: See TracChangeset for help on using the changeset viewer.