Changeset 230 for trunk/NEMO

Timestamp:

2005-03-18T20:45:28+01:00 (19 years ago)

Author:

opalod

Message:

CT : UPDATE168 : add a test to avoid under-float running error due to a too small e-fold

File:

• trunk/NEMO/OPA_SRC/TRA/traqsr.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/TRA/traqsr.F90

@@ -320 +320 @@
       IF( .NOT. lk_sco ) THEN
          ! z-coordinate with or without partial step : same before last ocean w-level everywhere
+         gdsr(:) = 0.e0
          DO jk = 1, jpk
             zdp1 = -fsdepw(1,1,jk)
             gdsr(jk) = ro0cpr * (  rabs  * EXP( zdp1/xsi1 ) + (1.-rabs) * EXP( zdp1/xsi2 )  )
+            IF ( gdsr(jk) <= 1.e-10 ) EXIT
          END DO
          indic = 0