BagaBaga
Full Member | Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Сбоственно, вот (только вместо 200 и 100 метки 2 и 1)
Есть два DO с меткой 2. Я немного не пойму, сколько раз выполняется внутренний цикл.
Код:
SUBROUTINE CONDUC(L,NL2,IDT,NTR,RADS,NH,PAR,STT,SFT,SFF,KPAR
* ,TETMT)
DIMENSION RADS(NH),PAR(KPAR,NTR,NL2),STT(IDT,NTR,NL2)
* ,SFT(IDT,NTR,NL2),SFF(IDT,NTR,NL2),TETMT(NL2)
DATA E/1.60219E-20/,OE/1.76E7/,OI/3.11E2/,G10/-.30356/,
*RE/6371.02E5/,CI1/4.23E-10/,CI2/4.28E-10/,CI3/2.58E-10/,
*CE1/1.82E-10/,CE11/3.6E-2/,CE2/2.33E-11/,CE21/1.21E-4/,
*CE3/2.8E-10/
IE=(NL2+1)/2
DO 2 I=1,NL2
TET=RAD(TETMT(I))
CT=COS(TET)
ST=SIN(TET)
SK=SQRT(1.+3.*CT*CT)
SI=(CT+CT)/SK
CI=ST/SK
SIS=SI*SI
CIS=CI*CI
DO 2 J=1,NTR
ROQ=(RE/(RE+RADS(J)))**3
B=-G10*ROQ*SK
EE=E/B
OME=OE*B
OMI=OI*B
CM=PAR(4,J,I)
CO2=PAR(1,J,I)
CN2=PAR(2,J,I)
CO=PAR(3,J,I)
TN=PAR(5,J,I)
B=SQRT(TN)
FI=CI1*CO2+CI2*CN2+CI3*CO
FE=CE1*(1.+CE11*B)*B*CO2+CE2*(1.-CE21*TN)*TN*CN2+CE3*B*CO
BU=OME*OME
CU=OMI*OMI
CO2=FI*FI
CN2=FE*FE
B=1./(BU+CN2)
CO=1./(CU+CO2)
SP=CM*EE*(OMI*FI*CO+OME*FE*B)
SH=CM*EE*(CU*CO-BU*B)
S0=CM*EE*(OMI/FI+OME/FE)
IF(I.NE.IE)GOTO1
STT(L,J,I)=S0
SFT(L,J,I)=0.
SFF(L,J,I)=SP+SH*SH/SP
GOTO2
1 CONTINUE
B=1./(S0*SIS+SP*CIS)
STT(L,J,I)=S0*SP*B
SFT(L,J,I)=S0*SH*SI*B
SFF(L,J,I)=(S0*SP*SIS+(SP*SP+SH*SH)*CIS)*B
2 CONTINUE
RETURN
END
|
PS
увидел GOTO1, но все равно не понял, зачем общее continue у внешнего и вложенного циклов. |
