Changeset 10693 for NEMO/branches/2019

Timestamp:

2019-02-15T18:58:27+01:00 (6 years ago)

Author:

mathiot

Message:

add line to move across the north fold cut line intermediate positions/velocities/accelerations (ticket #2238 solution 2)

File:

• NEMO/branches/2019/fix_ticket2238_solution2/src/OCE/ICB/icblbc.F90 (modified) (1 diff)

NEMO/branches/2019/fix_ticket2238_solution2/src/OCE/ICB/icblbc.F90

@@ -708 +708 @@
                         !
                         ! moving across the cut line means both position and
-                        ! velocity must change
+                        ! velocity and acceleration must change for icb position but also for intermediaite RK 4th scheme value
                         ijglo = INT( ipts/nicbpack )
                         iiglo = ipts - nicbpack*ijglo
-                        pt%xi = iiglo - ( pt%xi - REAL(iine,wp) )
-                        pt%yj = ijglo - ( pt%yj - REAL(ijne,wp) )
-                        pt%uvel = -1._wp * pt%uvel
-                        pt%vvel = -1._wp * pt%vvel
+                        pt%xi = iiglo - ( pt%xi - REAL(iine,wp) ) ; pt%uvel = -1._wp * pt%uvel
+                        pt%yj = ijglo - ( pt%yj - REAL(ijne,wp) ) ; pt%vvel = -1._wp * pt%vvel
+                        pt%xRK1(1) = iiglo - ( pt%xRK1(1) - REAL(iine,wp) ) ; pt%xRK1(2:4) = -1._wp * pt%xRK1(2:4)
+                        pt%xRK2(1) = iiglo - ( pt%xRK2(1) - REAL(iine,wp) ) ; pt%xRK2(2:4) = -1._wp * pt%xRK2(2:4)
+                        pt%xRK3(1) = iiglo - ( pt%xRK3(1) - REAL(iine,wp) ) ; pt%xRK3(2:4) = -1._wp * pt%xRK3(2:4)
+                        pt%xRK4(1) = iiglo - ( pt%xRK4(1) - REAL(iine,wp) ) ; pt%xRK4(2:4) = -1._wp * pt%xRK4(2:4)
+                        pt%yRK1(1) = ijglo - ( pt%yRK1(1) - REAL(ijne,wp) ) ; pt%yRK1(2:4) = -1._wp * pt%yRK1(2:4)
+                        pt%yRK2(1) = ijglo - ( pt%yRK2(1) - REAL(ijne,wp) ) ; pt%yRK2(2:4) = -1._wp * pt%yRK2(2:4)
+                        pt%yRK3(1) = ijglo - ( pt%yRK3(1) - REAL(ijne,wp) ) ; pt%yRK3(2:4) = -1._wp * pt%yRK3(2:4)
+                        pt%yRK4(1) = ijglo - ( pt%yRK4(1) - REAL(ijne,wp) ) ; pt%yRK4(2:4) = -1._wp * pt%yRK4(2:4)
                         !
                         ! now remove berg from list and pack it into a buffer