C
C Подпрограмма вычисления интерпол. параметров (трубка-трубка)
C для узла в экватор. плоскости
C
SUBROUTINE DRQ0(UX,U,KPART,PAR,NTSL,NL,MSUM,P1,NTPL)
INCLUDE 'PARAMETR.INC'
DIMENSION NTSL(NL),P1(KPART),U(NL),MSUM(NL)
*,PAR(KPART,NTPL),VN1(KPART0)
COMMON /OUTBRD/ COUT,VOUT,TIOUT,TEOUT,RMAXT
IF(UX.LT.U(NL) ) THEN
C Линия полярнее самой экваториальной:
DO 1 K=1,NL
IF(NTSL(K).NE.NTSL(1)) THEN
NI=K
GO TO 2
1 ENDIF
C NI - номер первой замкнутой линии:
2 IF(UX.LT.U(NI)) THEN
C линия полярнее первой замкнутой, т.е. в области разомкнутых -
C - интерполяция между внеш.гран.усл. и 1-й замкн.линией:
C VN1 - внешние граничные условия:
UN1=1./RMAXT
DO 3 M=1,3
VN1(M)=COUT
VN1(M+3)=VOUT
3 CONTINUE
VN1(7)=TIOUT
VN1(8)=TEOUT
C VN2 - вершина линии NI:
L=MSUM(NI)+NTSL(NI)/2+1
CALL INTERP(UN1,U(NI),UX,VN1,PAR(1,L),P1,KPART)
ELSE
C IF(UX.GE.U(NI)) THEN...
C линия между замкнутыми линиями сетки -
C - интерполяция между соседними линиями:
CALL FIND( NL,UX,U,NI)
L1=MSUM(NI)+NTSL(NI)/2+1
NI1=NI+1
L2=MSUM(NI1)+NTSL(NI1)/2+1
CALL INTERP(U(NI),U(NI1),UX,PAR(1,L1),PAR(1,L2),P1,KPART)
ENDIF
ELSE
C IF(UX.GE.U(NL) ) THEN...
C линия экваториальнее самой экваториальной - ей присваиваются
C значения узла в экватор.плоскости с самой нижней линии сетки:
L=MSUM(NL)+NTSL(NL)/2+1
DO 4 M=1,KPART
P1(M)=PAR(M,L)
4 CONTINUE
ENDIF
RETURN
END
C
|
