Changeset 14259 for NEMO/branches/UKMO/NEMO_4.0.3_icb_speed_limit/src/OCE/ICB/icbdyn.F90

Timestamp:

2021-01-05T12:11:16+01:00 (3 years ago)

Author:

davestorkey

Message:

UKMO/NEMO_4.0.3_icb_speed_limit : Bug fix and extra write-outs for error trap in icbutl_bilin_e.

File:

• NEMO/branches/UKMO/NEMO_4.0.3_icb_speed_limit/src/OCE/ICB/icbdyn.F90 (modified) (5 diffs)

NEMO/branches/UKMO/NEMO_4.0.3_icb_speed_limit/src/OCE/ICB/icbdyn.F90

@@ -97 +97 @@
             !                                         !**   A1 = A(X1,V1)
             CALL icb_accel( berg , zxi1, ze1, zuvel1, zuvel1, zax1,     &
-                 &                   zyj1, ze2, zvvel1, zvvel1, zay1, zdt_2 )
+                 &                 zyj1, ze2, zvvel1, zvvel1, zay1, zdt_2 )
             !
             zu1 = zuvel1 / ze1                           !**   V1 in d(i,j)/dt
@@ -119 +119 @@
             !                                         !**   A2 = A(X2,V2)
             CALL icb_accel( berg , zxi2, ze1, zuvel2, zuvel1, zax2,    &
-                 &                   zyj2, ze2, zvvel2, zvvel1, zay2, zdt_2 )
+                 &                 zyj2, ze2, zvvel2, zvvel1, zay2, zdt_2 )
             !
             zu2 = zuvel2 / ze1                           !**   V2 in d(i,j)/dt
@@ -140 +140 @@
             !                                         !**   A3 = A(X3,V3)
             CALL icb_accel( berg , zxi3, ze1, zuvel3, zuvel1, zax3,    &
-                 &                   zyj3, ze2, zvvel3, zvvel1, zay3, zdt )
+                 &                 zyj3, ze2, zvvel3, zvvel1, zay3, zdt )
             !
             zu3 = zuvel3 / ze1                           !**   V3 in d(i,j)/dt
@@ -153 +153 @@
             IF( rn_speed_limit > 0._wp ) THEN  ! Limit speed of bergs based on a CFL criteria (if asked)
                zuvel_rk4 = zuvel1 + 2.*(zuvel2 + zuvel3) + zuvel4
-               zvvel_rk4 = zvvel1 + 2.*(zvvel2 + zvvel3) + zuvel4
+               zvvel_rk4 = zvvel1 + 2.*(zvvel2 + zvvel3) + zvvel4
                CALL icb_speed_limit(ze1, ze2, zdt_6, zuvel_rk4, zvvel_rk4, rn_speed_limit, l_loop, zvel_factor, icount, 3)
                IF(l_loop) CYCLE
@@ -163 +163 @@
             !                                         !**   A4 = A(X4,V4)
             CALL icb_accel( berg , zxi4, ze1, zuvel4, zuvel1, zax4,    &
-                 &                   zyj4, ze2, zvvel4, zvvel1, zay4, zdt )
+                 &                 zyj4, ze2, zvvel4, zvvel1, zay4, zdt )
 
             zu4 = zuvel4 / ze1                           !**   V4 in d(i,j)/dt