anchor.diff:           
            DO I=1,6
               
  -            IF (I < 4 ) THEN           
  +            IF (I < 4 ) THEN
                  
  -               outDistribMesh%Force(I,J   ) = D_F_BF(I,J) 
  +               outDistribMesh%Force(I,J   ) = D_F_BF(I,J)
                  
               ELSE
                  
  
positive.code:       
        
           ! Attach the external distributed buoyancy loads to the distributed mesh so they can be transferred to the WRP
           
        DO J = 1, outDistribMesh%Nnodes
           
           DO I=1,6
              
              IF (I < 4 ) THEN           
                 
                 outDistribMesh%Force(I   ,J) = D_F_B(I,J) 
              
              ELSE
                 
                 outDistribMesh%Moment(I-3,J) = D_F_B(I,J)
                 
              END IF
              
           END DO  ! DO I
           
        END DO ! DO J
        
   
  
negative.code:       
        
           ! 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
           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%Moment(I-3,J) = L_F_B(I,J)
                 
                 END IF
              
              END DO  ! DO I
           end if              ! <= MSL2SWL check
        END DO ! DO J
        
           ! Remap for the lumped to WRP mesh transfer       
  