Changeset 15081

Timestamp:

2021-07-05T16:26:10+02:00 (3 years ago)

Author:

clem

Message:

debug algo andreas2005 in the trunk

File:

• NEMO/trunk/src/OCE/SBC/sbcblk_algo_ice_an05.F90 (modified) (1 diff)

NEMO/trunk/src/OCE/SBC/sbcblk_algo_ice_an05.F90

@@ -267 +267 @@
 
             !! *** TABLE 1 of Andreas et al. 2005 ***
-            !! Smooth flow condition (R* <= 0.135):
-            zsmoot = 0.5_wp + SIGN( 0.5_wp, (0.135_wp   - zre) ) ! zre <= 0.135: zsmoot==1 ; otherwize: zsmoot==0
-            !! Transition (0.135 < R* < 2.5):
-            ztrans = 0.5_wp + SIGN( 0.5_wp, (2.49999_wp - zre) ) - zsmoot
-            !! Rough ( R* > 2.5):
-            zrough = 0.5_wp + SIGN( 0.5_wp, (zre - 2.5_wp) )
-
-            IF( (zsmoot+ztrans+zrough > 1.001_wp).OR.(zsmoot+ztrans+zrough < 0.999_wp) ) &
-               CALL ctl_stop( ' rough_leng_tq@mod_blk_ice_an05.f90 => something wrong with zsmoot, ztrans, zrough!' )
-
+            zsmoot = 0._wp ; ztrans = 0._wp ; zrough = 0._wp
+            IF    ( zre <= 0.135_wp ) THEN   ! Smooth flow condition (R* <= 0.135):
+               zsmoot = 1._wp
+            ELSEIF( zre < 2.5_wp )    THEN   ! Transition (0.135 < R* < 2.5)
+               ztrans = 1._wp
+            ELSE                             ! Rough ( R* > 2.5)
+               zrough = 1._wp
+            ENDIF
+               
             zlog  = LOG(zre)
             zlog2 = zlog*zlog