anchor.diff:      type(json_value),intent(in),pointer :: json_command
       character(len=:),allocatable :: cval
       real :: CL_target,alpha_temp
  -    real :: delta,relaxation,alpha,maxres,residual,ans1,ans2,gradient
  +    real :: delta,relaxation,maxres,residual,ans1,ans2,gradient
       integer :: iter,maxiter
       type(json_value),pointer    :: p_root
       character(100) :: filename
  
positive.code:         t%alpha = t%alpha - relaxation*residual/gradient
  
          call sf_target_CL_residual(t,CL_target,residual)
          write(*,*) alpha*180.0/pi, t%GL(t%nrealwings+1,3), residual
  #else
          t%alpha = t%alpha + (CL_target - t%GL(t%nrealwings+1, 3)) / t%GL(t%nrealwings+1, 3)%dx(1)
          t%alpha%dx(1) = t%GL(t%nrealwings+1, 3)%dx(1) / t%GL(t%nrealwings+1, 3)%dx(1)
  
negative.code:     write(*,*)
      write(*,*) '   alpha[deg]                CL                        CL_Residual'
      call sf_target_CL_residual(t,CL_target,residual)
      write(*,*) t%alpha*180.0/pi,t%GL(t%nrealwings+1,3), residual
  #else
      write(*,*) '                   convergence = ',maxres
  
  