subroutine xljs(Qt,Ht,n11,D1,Yt)
real Qt,Ht,n11,D1
c
real Zz1(5000),Pz1(5000),vz1(5000)
real vz1x(5000),vz1y(5000),vz1z(5000)
real Zz2(5000),Pz2(5000),vz2(5000)
real vz2x(5000),vz2y(5000),vz2z(5000)
c
real Zy1(5000),Py1(5000),vy1(5000)
real vy1x(5000),vy1y(5000),vy1z(5000)
real Zy2(5000),Py2(5000),vy2(5000)
real vy2x(5000),vy2y(5000),vy2z(5000)
c
real Zr1(5000),Pr1(5000),vr1(5000)
real vr1x(5000),vr1y(5000),vr1z(5000)
real Zr2(5000),Pr2(5000),vr2(5000)
real vr2x(5000),vr2y(5000),vr2z(5000)
c
real Zw1(5000),Pw1(5000),vw1(5000)
real vw1x(5000),vw1y(5000),vw1z(5000)
real Zw2(5000),Pw2(5000),Vw2(5000)
real vw2x(5000),vw2y(5000),vw2z(5000)
c
real Hsl,Hy,Hr,Hw
real Ys,Yv,Ymd,Ymu,Yt
real Dp,miu,sa,cfo,wo,Rdm,Ru,Lu
real s1,s2,s3,s4,s5,s6,s7,s8,sp1,sp2
integer ncnum,numr,numz1,numz2
integer numy1,numy2,numr1,numr2,numw1,numw2
character xz,name,fluency |