getabc.f

gsac程序包

	subroutine getabc(m,omg,wvn,a,b,c,d,e,f)
	implicit none
	integer m
	COMPLEX*16 omg,wvn
	COMPLEX*16 a, b, c, d, e, f
	integer NL
	parameter (NL=100)
	common/timodel/H(NL),TA(NL),TC(NL),TL(NL),TN(NL),TF(NL),
     1		TRho(NL),
     2		qa(NL),qb(NL),etap(NL),etas(NL), 
     3		frefp(NL), frefs(NL)
	real H, TA, TC, TF, TL, TN, TRho
	real qa, qb, etap, etas, frefp, frefs
	a = wvn * TF(m) / TC(m)
	b = 1.0/TC(m)
	c = - TRho(m)*omg*omg + wvn*wvn *(TA(m) -TF(m)*TF(m)/TC(m))
	d = - wvn
	e = 1.0/TL(m)
	f = - TRho(m)*omg*omg
	return
	end