anchor.diff:  
      AbortErrLev     = ErrID_Fatal                                        ! Until we read otherwise from the FAST input file, we abort only on FATAL errors
      t_global        = t_initial - 20.                                    ! initialize this to a number < t_initial for error message in ProgAbort
  -   NextJacCalcTime = t_global  - 1.                                     ! We want to calculate the Jacobian on the first step (I just made it negative to ensure it is <= t_global
  +   calcJacobian    = .TRUE.                                             ! we need to calculate the Jacobian
  +   NextJacCalcTime = t_global                                           ! We want to calculate the Jacobian on the first step
  +   
      
         ! ... Initialize NWTC Library (open console, set pi constants) ...
      CALL NWTC_Init( ProgNameIN=FAST_ver%Name, EchoLibVer=.FALSE. )       ! sets the pi constants, open console for output, etc...
  
positive.code:   
        t_global_next = t_initial + (n_t_global+1)*p_FAST%DT  ! = t_global + p_FAST%dt
                         
        !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
        ! Step 1.a: Extrapolate Inputs (gives predicted values at t+dt
        ! 
  
negative.code:       END DO
        CALL SD_CopyInput (SD_Input(1),  u_SD,  MESH_NEWCOPY, Errstat, ErrMsg) ! do this to initialize meshes/allocatable arrays for output of ExtrapInterp routine
           CALL CheckError( ErrStat, 'Message from SD_CopyInput (u_SD): '//NewLine//ErrMsg )      
                                 
           
           ! Initialize predicted states for j_pc loop:
        CALL SD_CopyContState   ( x_SD,  x_SD_pred, MESH_NEWCOPY, Errstat, ErrMsg)
           CALL CheckError( ErrStat, 'Message from SD_CopyContState (init): '//NewLine//ErrMsg )
  