anchor.diff:        InitOut%CouplingScheme = ExplicitLoose
      END IF
      
  -   if (InitInp%Linearize) then
  -            
  -      ! Otherwise, if the module does allow linearization, return the appropriate Jacobian row/column names here:   
  -      ! Allocate and set these variables: InitOut%LinNames_y, InitOut%LinNames_x, InitOut%LinNames_xd, InitOut%LinNames_z, InitOut%LinNames_u 
  -      
  -      CALL AllocAry( InitOut%LinNames_y, 6+p%NumOuts, 'LinNames_y', ErrStat2, ErrMsg2 )
  -         CALL CheckError( ErrStat2, ErrMsg2 )
  -         IF (ErrStat >= AbortErrLev) RETURN
  -         
  -      do i=1,size(Indx_Y_BlPitchCom)
  -         InitOut%LinNames_y(Indx_Y_BlPitchCom(i)) = 'BlPitchCom('//trim(num2lstr(i))//')'
  -      end do
  -      InitOut%LinNames_y(Indx_Y_YawMom)  = 'YawMom'
  -      InitOut%LinNames_y(Indx_Y_GenTrq)  = 'GenTrq'
  -      InitOut%LinNames_y(Indx_Y_ElecPwr) = 'ElecPwr'
  -      do i=1,p%NumOuts
  -         InitOut%LinNames_y(i+Indx_Y_WrOutput) = p%OutParam(i)%Name
  -      end do
  -      
  -      
  -      CALL AllocAry( InitOut%LinNames_u, 3, 'LinNames_u', ErrStat2, ErrMsg2 )
  -         CALL CheckError( ErrStat2, ErrMsg2 )
  -         IF (ErrStat >= AbortErrLev) RETURN
  -
  -      InitOut%LinNames_u(Indx_u_Yaw    ) = 'Yaw'
  -      InitOut%LinNames_u(Indx_u_YawRate) = 'YawRate'
  -      InitOut%LinNames_u(Indx_u_HSS_Spd) = 'HSS_Spd'
  -                          
  -   end if
  -   
         !............................................................................................
         ! tip brakes - this may be added back, later, so we'll keep these here for now
         !............................................................................................      
  
positive.code:    PUBLIC :: SrvD_JacobianPConstrState           ! Routine to compute the Jacobians of the output (Y), continuous- (X), discrete-
                                                   !   (Xd), and constraint-state (Z) equations all with respect to the constraint
                                                   !   states (z)
     
     
  CONTAINS
  
negative.code:       END IF      
     
     END IF
     
     
     
        !............................................................................................
        ! Set Init outputs for linearization (after TMD, in case we ever add the TMD to the linearization features):
        !............................................................................................
     
     if (InitInp%Linearize) then
              
        ! If the module does allow linearization, return the appropriate Jacobian row/column names here:   
        ! Allocate and set these variables: InitOut%LinNames_y, InitOut%LinNames_x, InitOut%LinNames_xd, InitOut%LinNames_z, InitOut%LinNames_u 
        
        CALL AllocAry( InitOut%LinNames_y, 6+p%NumOuts, 'LinNames_y', ErrStat2, ErrMsg2 )
           CALL CheckError( ErrStat2, ErrMsg2 )
           IF (ErrStat >= AbortErrLev) RETURN
           
        do i=1,size(Indx_Y_BlPitchCom)
           InitOut%LinNames_y(Indx_Y_BlPitchCom(i)) = 'BlPitchCom('//trim(num2lstr(i))//'), rad'
        end do
        InitOut%LinNames_y(Indx_Y_YawMom)  = 'YawMom, Nm'
        InitOut%LinNames_y(Indx_Y_GenTrq)  = 'GenTrq, Nm'
        InitOut%LinNames_y(Indx_Y_ElecPwr) = 'ElecPwr, W'
        do i=1,p%NumOuts
           InitOut%LinNames_y(i+Indx_Y_WrOutput) = trim(p%OutParam(i)%Name)//', '//p%OutParam(i)%Units
        end do
        
        
        CALL AllocAry( InitOut%LinNames_u, 3, 'LinNames_u', ErrStat2, ErrMsg2 )
           CALL CheckError( ErrStat2, ErrMsg2 )
           IF (ErrStat >= AbortErrLev) RETURN
  
        InitOut%LinNames_u(Indx_u_Yaw    ) = 'Yaw, rad'
        InitOut%LinNames_u(Indx_u_YawRate) = 'YawRate, rad/s'
        InitOut%LinNames_u(Indx_u_HSS_Spd) = 'HSS_Spd, rad/s'
                            
     end if
     
     
        !............................................................................................
        ! Clean up the local variables:
  