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 у внешнего и вложенного циклов. |