Changeset 10693


Ignore:
Timestamp:
2019-02-15T18:58:27+01:00 (18 months ago)
Author:
mathiot
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/fix_ticket2238_solution2/src/OCE/ICB/icblbc.F90

    r10679 r10693  
    708708                        ! 
    709709                        ! moving across the cut line means both position and 
    710                         ! velocity must change 
     710                        ! velocity and acceleration must change for icb position but also for intermediaite RK 4th scheme value 
    711711                        ijglo = INT( ipts/nicbpack ) 
    712712                        iiglo = ipts - nicbpack*ijglo 
    713                         pt%xi = iiglo - ( pt%xi - REAL(iine,wp) ) 
    714                         pt%yj = ijglo - ( pt%yj - REAL(ijne,wp) ) 
    715                         pt%uvel = -1._wp * pt%uvel 
    716                         pt%vvel = -1._wp * pt%vvel 
     713                        pt%xi = iiglo - ( pt%xi - REAL(iine,wp) ) ; pt%uvel = -1._wp * pt%uvel 
     714                        pt%yj = ijglo - ( pt%yj - REAL(ijne,wp) ) ; pt%vvel = -1._wp * pt%vvel 
     715                        pt%xRK1(1) = iiglo - ( pt%xRK1(1) - REAL(iine,wp) ) ; pt%xRK1(2:4) = -1._wp * pt%xRK1(2:4) 
     716                        pt%xRK2(1) = iiglo - ( pt%xRK2(1) - REAL(iine,wp) ) ; pt%xRK2(2:4) = -1._wp * pt%xRK2(2:4) 
     717                        pt%xRK3(1) = iiglo - ( pt%xRK3(1) - REAL(iine,wp) ) ; pt%xRK3(2:4) = -1._wp * pt%xRK3(2:4) 
     718                        pt%xRK4(1) = iiglo - ( pt%xRK4(1) - REAL(iine,wp) ) ; pt%xRK4(2:4) = -1._wp * pt%xRK4(2:4) 
     719                        pt%yRK1(1) = ijglo - ( pt%yRK1(1) - REAL(ijne,wp) ) ; pt%yRK1(2:4) = -1._wp * pt%yRK1(2:4) 
     720                        pt%yRK2(1) = ijglo - ( pt%yRK2(1) - REAL(ijne,wp) ) ; pt%yRK2(2:4) = -1._wp * pt%yRK2(2:4) 
     721                        pt%yRK3(1) = ijglo - ( pt%yRK3(1) - REAL(ijne,wp) ) ; pt%yRK3(2:4) = -1._wp * pt%yRK3(2:4) 
     722                        pt%yRK4(1) = ijglo - ( pt%yRK4(1) - REAL(ijne,wp) ) ; pt%yRK4(2:4) = -1._wp * pt%yRK4(2:4) 
    717723                        ! 
    718724                        ! now remove berg from list and pack it into a buffer 
Note: See TracChangeset for help on using the changeset viewer.