anchor.diff:        
         
            ! Transfer the loads from the lumped mesh to the (0,0,0) point mesh
  -         
  +
  +         ! Because of wave stretching and user-supplied waves, we may have loads above the still water line (SWL) which will be used
  +         ! in the hydrodynamics for conditions where the wave height is > SWL.  So we now need to check that the vertical position
  +         ! is <= SWL for this summary file calculation.
  +
         DO J = 1, outLumpedMesh%Nnodes
  -          
  -         DO I=1,6
  -            
  -            IF (I < 4 ) THEN           
  +         if ( outLumpedMesh%Position(3,J) <= MSL2SWL ) then 
  +            DO I=1,6
  +            
  +               IF (I < 4 ) THEN           
                  
  -               outLumpedMesh%Force(I   ,J) = L_F_B(I,J) 
  -            
  -            ELSE
  +                  outLumpedMesh%Force(I   ,J) = L_F_B(I,J)
  +            
  +               ELSE
                  
  -               outLumpedMesh%Moment(I-3,J) = L_F_B(I,J)
  +                  outLumpedMesh%Moment(I-3,J) = L_F_B(I,J)
                  
  -            END IF
  -            
  -         END DO  ! DO I
  -         
  +               END IF
  +            
  +            END DO  ! DO I
  +         end if              ! <= MSL2SWL check
         END DO ! DO J
         
            ! Remap for the lumped to WRP mesh transfer       
  
positive.code:          
           DO I=1,6
              
              IF (I < 4 ) THEN           
                 
                 outDistribMesh%Force(I,J   ) = D_F_BF(I,J) 
                 
              ELSE
                 
  
negative.code:          
           DO I=1,6
              
              IF (I < 4 ) THEN
                 
                 outLumpedMesh%Force(I,J) = L_F_BF(I,J)
              
              ELSE
                 
  