anchor.diff:  !**********************************************************************************************************************************
   ! LICENSING
   ! Copyright (C) 2015-2016  National Renewable Energy Laboratory
  -! Copyright (C) 2016-2017  Envision Energy USA, LTD   
  -!   
  +! Copyright (C) 2016-2017  Envision Energy USA, LTD
  +!
   ! Licensed under the Apache License, Version 2.0 (the "License");
   ! you may not use this file except in compliance with the License.
   ! You may obtain a copy of the License at
  
positive.code:    CALL AllocAry(p%Mass0_QP,6,6,p%nqp*p%elem_total,'Mass0_QP',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
        if (ErrStat >= AbortErrLev) return
        
     IF(p%quadrature .EQ. GAUSS_QUADRATURE) THEN
  
         ! Compute sectional propertities ( 6 by 6 stiffness and mass matrices)
         ! at Gauss points
        
        !.... Compute the relative location of each quadrature node along the length of the curved blade ....
         CALL AllocAry(temp_ratio,p%nqp,p%elem_total,'temp_ratio',ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
negative.code: 
  
        ! Set the initial displacements
     CALL BD_QuadraturePointDataAt0(p)
        if (ErrStat >= AbortErrLev) then
           call cleanup()
           return
  