calc_ider_tee.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	calc_ider_tee (df, pt1, Tt1, xflow1, xflow2, pt2, Tt2, A1, A2, zeta_fac, kappa, R, ider, iflag)

Function/Subroutine Documentation

calc_ider_tee()

subroutine calc_ider_tee	(	real*8, dimension(6)	df,
		real*8	pt1,
		real*8	Tt1,
		real*8	xflow1,
		real*8	xflow2,
		real*8	pt2,
		real*8	Tt2,
		real*8	A1,
		real*8	A2,
		real*8	zeta_fac,
		real*8	kappa,
		real*8	R,
		integer	ider,
		integer	iflag
	)

!
      implicit none
!
      integer ider,iflag
!
      real*8 
     &df(6),
     &pt1,
     &pt2,
     &tt1,
     &tt2,
     &xflow1,
     &xflow2,
     &a1,
     &a2,
     &kappa,
     &r,
     &calc_residual_tee,
!     Accuracy of the numerical derivatives
     &eps,
!     Step for the derivative
     &h,
     &f0,
     &zeta_fac
!
      eps = 1.0e-8
!
      f0 = calc_residual_tee(pt1,tt1,xflow1,xflow2,pt2,
     &tt2,a1,a2,zeta_fac,kappa,r,ider,iflag)
!
      h = eps*dabs(pt1)
      if(h.eq.0)then
         h = eps
      endif
      df(1) = (calc_residual_tee(pt1+h,tt1,xflow1,xflow2,pt2,
     &tt2,a1,a2,zeta_fac,kappa,r,ider,iflag)-f0)/h
!
      h = eps*dabs(tt1)
      if(h.eq.0)then
         h = eps
      endif
      df(2) = (calc_residual_tee(pt1,tt1+h,xflow1,xflow2,pt2,
     &tt2,a1,a2,zeta_fac,kappa,r,ider,iflag)-f0)/h
!
      h = eps*dabs(xflow1)
      if(h.eq.0)then
         h = eps
      endif
      df(3) = (calc_residual_tee(pt1,tt1,xflow1+h,xflow2,pt2,
     &tt2,a1,a2,zeta_fac,kappa,r,ider,iflag)-f0)/h
!
      h = eps*dabs(xflow2)
      if(h.eq.0)then
         h = eps
      endif
      df(4) = (calc_residual_tee(pt1,tt1,xflow1,xflow2+h,pt2,
     &tt2,a1,a2,zeta_fac,kappa,r,ider,iflag)-f0)/h
!
      h = eps*dabs(pt2)
      if(h.eq.0)then
         h = eps
      endif
      df(5) = (calc_residual_tee(pt1,tt1,xflow1,xflow2,pt2+h,
     &tt2,a1,a2,zeta_fac,kappa,r,ider,iflag)-f0)/h
!
      h = eps*dabs(tt2)
      if(h.eq.0)then
         h = eps
      endif
      df(6) = (calc_residual_tee(pt1,tt1,xflow1,xflow2,pt2,
     &tt2+h,a1,a2,zeta_fac,kappa,r,ider,iflag)-f0)/h
!
      return