calc_ider_wye.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	calc_ider_wye (df, pt1, Tt1, xflow1, xflow2, pt2, Tt2, ichan_num, A1, A2, A_s, dh1, dh2, alpha, zeta_fac, kappa, R, ider, iflag)

Function/Subroutine Documentation

calc_ider_wye()

subroutine calc_ider_wye	(	real*8, dimension(6)	df,
		real*8	pt1,
		real*8	Tt1,
		real*8	xflow1,
		real*8	xflow2,
		real*8	pt2,
		real*8	Tt2,
		integer	ichan_num,
		real*8	A1,
		real*8	A2,
		real*8	A_s,
		real*8	dh1,
		real*8	dh2,
		real*8	alpha,
		real*8	zeta_fac,
		real*8	kappa,
		real*8	R,
		integer	ider,
		integer	iflag
	)

!
      implicit none
!
      integer ichan_num,ider,iflag
!
      real*8 
     &df(6),
     &pt1,
     &pt2,
     &tt1,
     &tt2,
     &xflow1,
     &xflow2,
     &a1,
     &a2,
     &a_s,
     &kappa,
     &r,
     &dh1,
     &dh2,
     &alpha,
     &calc_residual_wye,
     &eps,
     &h,
     &f0,
     &zeta_fac
!
      eps = 1.0e-4
!
      f0 = calc_residual_wye(pt1,tt1,xflow1,xflow2,pt2,
     &tt2,ichan_num,a1,a2,a_s,dh1,dh2,alpha,zeta_fac,kappa,
     &r,ider,iflag)
!
      h = eps*dabs(pt1)
      if(h.eq.0)then
         h = eps
      endif
      df(1) = (calc_residual_wye(pt1+h,tt1,xflow1,xflow2,pt2,
     &tt2,ichan_num,a1,a2,a_s,dh1,dh2,alpha,zeta_fac,kappa,r,ider,
     &iflag)-f0)/h
!
      h = eps*dabs(tt1)
      if(h.eq.0)then
         h = eps
      endif
      df(2) = (calc_residual_wye(pt1,tt1+h,xflow1,xflow2,pt2,
     &tt2,ichan_num,a1,a2,a_s,dh1,dh2,alpha,zeta_fac,kappa,r,ider,
     &iflag)-f0)/h
!
      h = eps*dabs(xflow1)
      if(h.eq.0)then
         h = eps
      endif
      df(3) = (calc_residual_wye(pt1,tt1,xflow1+h,xflow2,pt2,
     &tt2,ichan_num,a1,a2,a_s,dh1,dh2,alpha,zeta_fac,kappa,r,ider,
     &iflag)-f0)/h
!
      h = eps*dabs(xflow2)
      if(h.eq.0)then
         h = eps
      endif
      df(4) = (calc_residual_wye(pt1,tt1,xflow1,xflow2+h,pt2,
     &tt2,ichan_num,a1,a2,a_s,dh1,dh2,alpha,zeta_fac,kappa,r,ider,
     &iflag)-f0)/h
!
      h = eps*dabs(pt2)
      if(h.eq.0)then
         h = eps
      endif
      df(5) = (calc_residual_wye(pt1,tt1,xflow1,xflow2,pt2+h,
     &tt2,ichan_num,a1,a2,a_s,dh1,dh2,alpha,zeta_fac,kappa,r,ider,
     &iflag)-f0)/h
!
      h = eps*dabs(tt2)
      if(h.eq.0)then
         h = eps
      endif
      df(6) = (calc_residual_wye(pt1,tt1,xflow1,xflow2,pt2,
     &tt2+h,ichan_num,a1,a2,a_s,dh1,dh2,alpha,zeta_fac,kappa,r,ider,
     &iflag)-f0)/h
!
      return