:::::::::::::: tor0.f :::::::::::::: PROGRAM TOR0 REAL*4 IN,K,ETA,ETA1,ETA2,LETA,DLETA,NY,HRO,R,R1,R2,DR,T REAL*4 T1,T2,DT,RB,PI,GRAD,FI,FI1,F12,FI0,OM,VR,F1,F2,F3 REAL*4 RV(20),C,C1,C2,C3,RV1,RV2,RV3,LK,DR0,D,DF,R12,DTUT REAL*4 DE(20),QE(20),EKIN,EKINV,EPOT,EPOTV,TAP,TUT,TP,OM0 REAL*4 DD,RK,EPOT0 INTEGER I,IE,IR,IT,J,JR,JT,N,NR,NT,M,MR,MT,NETA C C GAGNLEGIR STUDLAR SEM KOMA FYRIR C PI=3.1415269 DTUT=1. GRAD=180./PI LK=.0254 C LENGD KEFLIS ^MF C RB=.04 DR0=.005 C TEGRUNARBIL ^MF C C LESA INN STUDLA SEM REIKNA A MED C WRITE(6,*) ' RADIUS BIKARS ^MF=' READ(5,*) RB WRITE(6,2) 2 FORMAT(' LESA INN HRO OLIU OG SNUNINGSSTUDUL: ',$) READ(5,*) HRO,K WRITE(6,5) 5 FORMAT(' LESA INN ETA1,ETA2- BILID OG FJOLDA GILDA:') READ(5,*) ETA1,ETA2,NETA WRITE(6,10) 10 FORMAT(' LESA INN RADIUS KEFLIS, HVERFITREGDU OG SNHORN: ',$) READ(5,*) RK,IN,FI0 FI0=FI0/GRAD C BREYTA I RADIANA TP=2.*PI*SQRT(IN/K) EPOT0=.5*K*FI0**2 R1=RK R12=R1**2 WRITE(1,16) HRO,K,RK 16 FORMAT(' HRO OLIU=',F4.0,' KG/M**3 SNUNINGSSTUDULL=',E10.3 1 ,'N*M/RAD RADIUS KEFLIS=',E10.3,'M') WRITE(1,17) IN,TP,FI0*GRAD 17 FORMAT(' HVERFITREGDA=',E10.3,'KG*M**2 SVEIFLUTIMI=',F6.3 1 ,'S BYRJUNARHORN=',F5.0,'GRADUR') WRITE(1,18) EPOT0 18 FORMAT(' BYRJUNARORKA (EPOT)= ',E11.4,' JOULE'///) C C UNDIRBUA LYKKJU I ETA C DLETA=LOG(ETA2/ETA1)/(NETA-1) LETA=LOG(ETA1)-DLETA DO 100 IE=1,NETA TUT=DTUT NR=INT((RB-R1)/DR0)+1 DR=(RB-R1)/(NR-1) R2=R1+DR LETA=LETA+DLETA ETA=EXP(LETA) C C UNDIRBUA TEGRUNINA DT OG DR C NY=ETA/HRO C KINEMATISK SEIGJA DT=0.4*DR**2/NY C TIMASKREF '0.4'<0.5 WRITE(6,*)' DT: ',DT C1=2.-K*DT**2/IN C2=2.*PI*LK*ETA*RK*DT**2/IN DD=NY*DT/DR**2 D=NY*DT*.5/DR DF=NY*DT C C HER KEMUR T TEGRUNIN C F2=FI0 F3=FI0 DO 105 I=1,NR C UPPHAFSHRADI = 0 RV(I)=0. 105 CONTINUE T=0. 200 T=T+DT F1=F2 F2=F3 F3=-F1+C1*F2+C2*((4.*RV(2)-3.*RV(1)-RV(3))/(2.*DR)-2.*RV(1)/R1) OM=(3.*F3-4.*F2+F1)/(2.*DT) RV2=R12*OM RV(1)=RV2 RV3=RV(2) R=R1 C C SVO KEMUR TEGRUNIN YFIR R C DO 110 IR=2,NR-1 RV1=RV2 RV2=RV3 RV3=RV(IR+1) R=R+DR RV(IR)=RV2+DD*(RV1-(RV2+RV2)+RV3)-D*(RV3-RV1)/R 110 CONTINUE C ENDIR R LYKKJU C C UTSKRIFT FYRIR HVERJAR DTUT SEKUNDUR c IF(T-TUT) 210,220,220 220 TUT=TUT+DTUT DO 115 I=1,NR DE(I)=PI*HRO*LK*RV(I)**2/(RK+(I-1)*DR) 115 CONTINUE C CALL QSF(DR,DE,QE,NR) EKINV=QE(NR) EKIN=.5*IN*OM**2 EPOT=.5*K*F3**2 TAP=EPOT0-(EPOT+EKIN) WRITE(1,21) 21 FORMAT('0',' TIMI^SF HORN^GF OMEGA EPOT EKIN 1 TAP EKIN-VOKVI ^JF') WRITE(1,22) T,F3*GRAD,OM,EPOT,EKIN,TAP,EKINV 22 FORMAT(7E10.3/) WRITE(1,23) (RK+(I-1)*DR,I=1,NR) 23 FORMAT(' R= ',10E10.3) WRITE(1,24) (RV(I),I=1,NR) 24 FORMAT(' R*V= ',10E10.3) WRITE(1,26) (RV(I)/(RK+(I-1)*DR),I=1,NR) 26 FORMAT(' V= ',10E10.3) C C ATHUGA HVORT KEFLI SE AD SNUA VID C 210 IF(F2.GE.F1 .AND. F3.GT.F2 .OR. F2.LE.F1 .AND. F3.LT.F2) 1 GO TO 200 WRITE(6,15) C PRENTA HAUS WRITE(1,15) 15 FORMAT(5X,'ETA ^KG/M*SF',5X,'HORN ^GRADF ',5X,' TIMI ^SEKF ') WRITE(6,20) ETA,F3*GRAD,T WRITE(1,20) ETA,F3*GRAD,T WRITE(6,25) (RV(I),I=1,NR) 20 FORMAT(3(7X,F7.2)///) 25 FORMAT(10E10.3) 100 CONTINUE C ENDIR ETA LYKKJU STOP END