anchor.diff:           ErrStat = ErrID_Fatal
            RETURN
         END IF
  -      
  +      m%L_F_B = 0.0_ReKi
         ALLOCATE ( m%L_F_D(3,y%LumpedMesh%Nnodes), STAT = ErrStat )
         IF ( ErrStat /= ErrID_None ) THEN
            ErrMsg  = ' Error allocating space for L_F_D array.'
            ErrStat = ErrID_Fatal
            RETURN
         END IF
  +      m%L_F_D = 0.0_ReKi
         ALLOCATE ( m%L_F_I(6,y%LumpedMesh%Nnodes), STAT = ErrStat )
         IF ( ErrStat /= ErrID_None ) THEN
            ErrMsg  = ' Error allocating space for L_F_I array.'
            ErrStat = ErrID_Fatal
            RETURN
         END IF
  +      m%L_F_I = 0.0_ReKi
         !ALLOCATE ( m%L_F_DP(6,y%LumpedMesh%Nnodes), STAT = ErrStat )
         !IF ( ErrStat /= ErrID_None ) THEN
         !   ErrMsg  = ' Error allocating space for L_F_DP array.'
  
positive.code:                y%DistribMesh%Force(I,J) = m%D_F_AM(I,J) + m%D_F_D(I,J)  + m%D_F_I(I,J) + m%D_F_B(I,J) +  p%D_F_MG(I,J) + p%D_F_BF(I,J)
              ELSE
                 m%D_F_B(I,J) = elementWaterState * p%D_F_B(I,J)
                 y%DistribMesh%Moment(I-3,J) =   m%D_F_AM(I,J) + m%D_F_B(I,J) + p%D_F_BF(I,J)     
              END IF
           END DO  ! DO I
           
  
negative.code:          ErrStat = ErrID_Fatal
           RETURN
        END IF
        m%L_FV = 0.0_ReKi
        ALLOCATE ( m%L_FA(3,y%LumpedMesh%Nnodes), STAT = ErrStat )
        IF ( ErrStat /= ErrID_None ) THEN
           ErrMsg  = ' Error allocating space for L_FA array.'
           ErrStat = ErrID_Fatal
           RETURN
        END IF
        m%L_FA = 0.0_ReKi
        ALLOCATE ( m%L_FDynP(y%LumpedMesh%Nnodes), STAT = ErrStat )
        IF ( ErrStat /= ErrID_None ) THEN
           ErrMsg  = ' Error allocating space for L_FDynP array.'
           ErrStat = ErrID_Fatal
           RETURN
        END IF
        m%L_FDynP = 0.0_ReKi
        ALLOCATE ( m%L_F_AM(6,y%LumpedMesh%Nnodes), STAT = ErrStat )
        IF ( ErrStat /= ErrID_None ) THEN
           ErrMsg  = ' Error allocating space for L_F_AM array.'
           ErrStat = ErrID_Fatal
           RETURN
        END IF
        m%L_F_AM = 0.0_ReKi
        
           ! Define initial guess for the system inputs here:
  
  