hskl.f

gsac程序包

	subroutine hskl(cossh,rsinsh,sinshr,TL,iwat,hl,atnb)
	implicit none
	complex*16 cossh, rsinsh, sinshr
	real TL
	integer iwat
	complex*16 hl(2,2)
	complex*16 atnb
	if(iwat.eq.0)then	
		hl(1,1) = cossh
		hl(2,1) = TL*rsinsh*atnb*atnb
		hl(1,2) = sinshr/(TL*atnb*atnb)
		hl(2,2) = cossh
	else
		hl(1,1) = dcmplx(1.0d+00,0.0d+00)
		hl(1,2) = dcmplx(0.0d+00,0.0d+00)
		hl(2,1) = dcmplx(0.0d+00,0.0d+00)
		hl(2,2) = dcmplx(1.0d+00,0.0d+00)
	endif
	return
	end