Fortran programe.docx
Content
Fortran programe
c surface area and volume of box
open(unit=1,file='c:\msdev\box.dat')
open(unit=2,file='c:\msdev\result.dat')
read(1,*)a,b,c
sur=2.0*(a*b+a*c+b*c)
vol=a*b*c
write(2,10)a,b,c,sur,vol
10
format(5x,'dimension of box'/7x,3(f6.2,15x)//
+6x,'surface area =',f10.1//6x,'volume=',f16.2)
stop
end
c programe fr calculating the area and premeter f a triangle
open(unit=1,file='c:\msdev\tri.dat')
open(unit=2,file='c:\msdev\result of tri.dat')
read (1,*)a,b,c
pre=a+b+c
s=pre/2
area=sqrt(s*(s-a)*(s-b)*(s-c))
write(2,10)a,b,c,pre,area
10 format(15x,'trianle specificaion'//15x,3(f5.2,15x)//
+15x,'premeter of triangle=
',15x,f7.3,10x,'meter'/
+15x,'surface area of triangle= ',15x,f7.3,10x,'meter^2')
stop
end
c
c **** triangle area for more triangles ******
dimension a(20),b(20),c(20),s(20),p(20),area(20)
n=5
open(unit=1,file='c:\msdev\mretri.dat')
open(unit=2,file='c:\msdev\result.dat')
read(1,*)(a(20),b(20),c(20),i=1,n)
do 10 i=1,n
p(i)=a(i)+b(i)+c(i)
s(i)=p(i)/2.0
area(i)=sqrt(s(i)*(s(i)-a(i))*(s(i)-b(i))*(s(i)-c(i)))
write(*,*)i,a(i),b(i),c(i),p(i),s(i),area(i)
write(2,15)i,a(i),b(i),c(i),p(i),s(i),area(i)
15
format(2x,i2,2x,3(f10.2,2x),2x,f10.2,2x,f10.2,2x,f10.2)
10
continue
stop
end
C THE PROGRAM USED TO CALCULATE THE AMOUNT OF HEAT
CONDUCTED
C
ACROSS THE RECTANGULAR FIN
C
CALCULATE THE TEMPERATURE DISTRIBUTION
DIMENSION T(100),TT(100),X(100)
OPEN(UNIT=1,FILE='C:/MSDEV/AS1.DAT')
OPEN(UNIT=2,FILE='C:/MSDEV/AS2.DAT')
OPEN(UNIT=3,FILE='C:/MSDEV/AS3.DAT')
READ(1,*)TZERO,H,AK,THICK,TAIR,AL
WRITE(2,102)TZERO,H,AK,THICK,TAIR,AL
102
FORMAT(1h1,/,
+'Temperature at the base of the fin =',F7.2,' deg.c',//,
+'Convective heat transfer coeff
+'Theramal conductivity
+'Thickness of the fin
+'Temp of the ambient air
+'Length of the fin
=',F7.2,' watts/sq.m.,k',//,
=',F7.2,' watts/m.k.',//,
=',F5.2,' mm',//,
=',F6.2,' deg.c',//,
=',F6.2,' cms',//)
AL=AL/100.0
THICK=THICK/1000.0
C
READ (1,*)M
C
IF(M.EQ.0) GOTO 500
DELX=AL/(M-1.)
A1=H*DELX**2/(AK*THICK)
READ(1,*)TOL
MM=M-1
C------INITIAL TRIAL VALUE
C------------------------DO 203 I=1,M
X(I)=DELX*(I-1)
203
T(I)=TZERO+(TAIR-TZERO)*((I-1.)*DELX/AL)
KOUNT=0
C------ITERATED VALUE
C------------------300
TT(1)=TZERO
C --------EQUATION 5-2a
C-----------------------
DO 204 I=2,MM
204
TT(I)=(2.*A1*TAIR +T(I-1)+T(I+1))/(2.*(A1+1.))
C
C --------EQUATION 5-5a
C---------------------TT(M)=((A1+H*DELX/AK)*TAIR+T(MM))/(1.+A1+H*DELX/AK)
KOUNT=KOUNT+1
DO 205 I=2,M
205
IF(ABS((TT(I)-T(I))/T(I)).GT.TOL) GOTO 12
GOTO 302
12
DO 333 I=1,M
333
T(I)=TT(I)
GO TO 300
302
CONTINUE
C
C
call polyft(x,tt,z,p,q,c,m,mm,mn,std,1.)
C
qfin1=-c(2)*Ak*thick*100.
QFIN2=AK*THICK*(T(1)-T(2))/DELX+(.75*T(1)+.25*T(2)TAIR)*H*DELX
WRITE(3,604)
604
FORMAT(3X,'I',8X,'X(I)',4X,'T(I)',7X,'TT(I)')
WRITE(*,*)(I,X(I),T(I),TT(I),I=1,M)
WRITE(3,605)(I,X(I),T(I),TT(I),I=1,M)
605
FORMAT(I4,3x,f8.2,3x,f8.2,3x,f8.2)
Tol=Tol*100.0
WRITE(2,555) TOL
555
FORMAT(/,'the last two iterated values are within ',f7.2,'%',//)
WRITE(2,607) KOUNT
607
FORMAT(//,'solution was obtained after',I7,2X,'iterations',//)
WRITE(2,606)QFIN2
606
FORMAT('RATE OF HT.TR.(Finite Diff.))',f12.2,'watts')
STOP
END
C
THE PROGRAM USED TO CALCULATE THE AMOUNT OF HEAT
CONDUCTED
C
ACROSS THE PLAT
C
CALCULATE THE TEMPERATURE DISTRIBUTION
DIMENSION T(1000,1000),TT(1000,1000),X(1000,1000),Y(1000,1000)
OPEN(UNIT=1,FILE='C:/MSDEV/PROJ1.DAT')
OPEN(UNIT=2,FILE='C:/MSDEV/PROJ2.DAT')
OPEN(UNIT=3,FILE='C:/MSDEV/PROJ3.DAT')
READ(1,*)((T(M,N),I=1,M),J=1,N )
READ(1,*)TZERO,H,AK,THICK,TAIR,XA,YB
WRITE(2,102)TZERO,H,AK,THICK,TAIR,XA,YB
102
FORMAT(1h1,/,
+'Temperature at the lower plat =',F7.2,' deg.c',//,
+'Convective heat transfer coeff
+'Theramal conductivity
+'Thickness of the fin
+'Temp of the ambient air
=',F7.2,' watts/sq.m.,k',//,
=',F7.2,' watts/m.k.',//,
=',F5.2,' mm',//,
=',F6.2,' deg.c',//,
+'Length of the plat
=',F6.2,' cms',//,
+'width of the plat
=',F6.2,' cms',//)
XA=XA/100.0
YB=YB/100.0
THICK=THICK/1000.0
C
READ (1,*)M,N
C
IF(M.EQ.0) GOTO 500
DELX=XA/(M-1.)
DELY=YB/(N-1.)
B1=(H*DELX)/AK
READ(1,*)TOL
MM=M-1
NN=N-1
C------INITIAL TRIAL VALUE
C------------------------DO 203 I=1,M
DO 203 J=1,N
X(I)=DELX*(I-1)
Y(J)=DELY*(J-1)
203
T(I)=TZERO+(TAIR-TZERO)*((I-1.)*DELX/XA)
T(J)=TZERO+(TAIR-TZERO)*((J-1.)*DELX/YB)
KOUNT=0
C------ITERATED VALUE
C------------------300
TT(1,1)=TZERO
C --------EQUATION 5-2a
C----------------------DO 204 I=2,MM
DO 204 J=2,NN
204
TT(I,J)=(T(I-1,J)+T(I+1,J)+T(I,J-1))/4
C
C --------EQUATION 5-5a
C---------------------C---------FOR INSULATED NODE
C---------LEFT NODES
DO 2000 I=1
DO 2000 J=2,NN
2000 TT(I,J)=(T(I,J+1)+T(I,J-1)+(2*T(I+1,J)))/4
C---------RIGHT NODES
DO 3000 I=M
DO 3000 J=2,NN
3000 TT(I,J)=(T(I,J+1)+T(I,J-1)+(2*T(I-1,J)))/4
C---------BOTTOM NODES
DO 4000
I=2,3
DO 4000
J=1
4000 TT(I,J)=(T(I+1,J)+T(I-1,J)+(2*T(I,J+1)))/4
DO 5000 I=4,7
DO 5000 J=1
5000 TT(I,J)=TZERO
DO 6000 I=8,9
DO 6000 J=1
6000 TT(I,J)=(T(I+1,J)+T(I-1,J)+(2*T(I,J+1)))/4
C---------FOR CONVECTION NODES
C---------UPPER NODES
DO 7000 I=2,MM
DO 7000 J=N
7000 TT(I,J)=(T(I,J-1)+((T(I+1,J)+T(I-1,J))/2)+(B1*TAIR))/(2+B1)
C---------FOR INSULATED CORNER NODE
C---------BOTTOM RIGHT NODES
TT(1,M)=(T(I,J+1)+T(I-1,J))/2
C---------BOTTOM LEFT NODES
TT(1,1)=(T(I,J+1)+T(I+1,J))/2
C---------FOR JNSULATED & CONVECTION CORNER NODE
C---------RIGHT UPPER NODE
TT(M,N)=(T(I-1,J)+T(I,J-1)+(B1*TAIR))/(2+B1)
C---------LEFT UPPER NODE
TT(1,N)=(T(I+1,J)+T(I,J-1)+(B1*TAIR))/(2+B1)
KOUNT=KOUNT+1
DO 205 I=2,M
DO 205 J=2,N
205
IF(ABS((TT(I,J)-T(I,J))/T(I,J)).GT.TOL) GOTO 12
GOTO 302
12
DO 333 I=1,M
DO 333 J=1,N
333
T(I,J)=TT(I,J)
GO TO 300
302
CONTINUE
C
C
call polyft(x,tt,z,p,q,c,m,mm,mn,std,1.)
C
qfin1=-c(2)*Ak*thick*100.
C
QFIN2=AK*THICK*(T(1)-T(2))/DELX+(.75*T(1)+.25*T(2)TAIR)*H*DELX
C
WRITE(3,604)
C604 FORMAT(3X,'I',3X, 'J' ,8X,'X(I)',8X,'Y(J),4X,'T(I,J)',7X,'TT(I,J)
C +')
C
WRITE(*,*)(I,J,X(I),Y(J),T(I,J),TT(I,J),I=1,M,J=1,N)
WRITE(3,605)(I,J,X(I),Y(J),T(I,J),TT(I,J),I=1,M,J=1,N)
605
FORMAT(I4,3x,f8.2,3x,f8.2,3x,f8.2)
Tol=Tol*100.0
WRITE(2,555) TOL
555
FORMAT(/,'the last two iterated values are within ',f7.2,'%',//)
WRITE(2,607) KOUNT
607
FORMAT(//,'solution was obtained after',I7,2X,'iterations',//)
WRITE(2,606)QFIN2
606
FORMAT('RATE OF HT.TR.(Finite Diff.))',f12.2,'watts')
STOP
END
c surface area and volume of box
open(unit=1,file='c:\msdev\box.dat')
open(unit=2,file='c:\msdev\result.dat')
read(1,*)a,b,c
sur=2.0*(a*b+a*c+b*c)
vol=a*b*c
write(2,10)a,b,c,sur,vol
10
format(5x,'dimension of box'/7x,3(f6.2,15x)//
+6x,'surface area =',f10.1//6x,'volume=',f16.2)
stop
end
c programe fr calculating the area and premeter f a triangle
open(unit=1,file='c:\msdev\tri.dat')
open(unit=2,file='c:\msdev\result of tri.dat')
read (1,*)a,b,c
pre=a+b+c
s=pre/2
area=sqrt(s*(s-a)*(s-b)*(s-c))
write(2,10)a,b,c,pre,area
10 format(15x,'trianle specificaion'//15x,3(f5.2,15x)//
+15x,'premeter of triangle=
',15x,f7.3,10x,'meter'/
+15x,'surface area of triangle= ',15x,f7.3,10x,'meter^2')
stop
end
c
c **** triangle area for more triangles ******
dimension a(20),b(20),c(20),s(20),p(20),area(20)
n=5
open(unit=1,file='c:\msdev\mretri.dat')
open(unit=2,file='c:\msdev\result.dat')
read(1,*)(a(20),b(20),c(20),i=1,n)
do 10 i=1,n
p(i)=a(i)+b(i)+c(i)
s(i)=p(i)/2.0
area(i)=sqrt(s(i)*(s(i)-a(i))*(s(i)-b(i))*(s(i)-c(i)))
write(*,*)i,a(i),b(i),c(i),p(i),s(i),area(i)
write(2,15)i,a(i),b(i),c(i),p(i),s(i),area(i)
15
format(2x,i2,2x,3(f10.2,2x),2x,f10.2,2x,f10.2,2x,f10.2)
10
continue
stop
end
C THE PROGRAM USED TO CALCULATE THE AMOUNT OF HEAT
CONDUCTED
C
ACROSS THE RECTANGULAR FIN
C
CALCULATE THE TEMPERATURE DISTRIBUTION
DIMENSION T(100),TT(100),X(100)
OPEN(UNIT=1,FILE='C:/MSDEV/AS1.DAT')
OPEN(UNIT=2,FILE='C:/MSDEV/AS2.DAT')
OPEN(UNIT=3,FILE='C:/MSDEV/AS3.DAT')
READ(1,*)TZERO,H,AK,THICK,TAIR,AL
WRITE(2,102)TZERO,H,AK,THICK,TAIR,AL
102
FORMAT(1h1,/,
+'Temperature at the base of the fin =',F7.2,' deg.c',//,
+'Convective heat transfer coeff
+'Theramal conductivity
+'Thickness of the fin
+'Temp of the ambient air
+'Length of the fin
=',F7.2,' watts/sq.m.,k',//,
=',F7.2,' watts/m.k.',//,
=',F5.2,' mm',//,
=',F6.2,' deg.c',//,
=',F6.2,' cms',//)
AL=AL/100.0
THICK=THICK/1000.0
C
READ (1,*)M
C
IF(M.EQ.0) GOTO 500
DELX=AL/(M-1.)
A1=H*DELX**2/(AK*THICK)
READ(1,*)TOL
MM=M-1
C------INITIAL TRIAL VALUE
C------------------------DO 203 I=1,M
X(I)=DELX*(I-1)
203
T(I)=TZERO+(TAIR-TZERO)*((I-1.)*DELX/AL)
KOUNT=0
C------ITERATED VALUE
C------------------300
TT(1)=TZERO
C --------EQUATION 5-2a
C-----------------------
DO 204 I=2,MM
204
TT(I)=(2.*A1*TAIR +T(I-1)+T(I+1))/(2.*(A1+1.))
C
C --------EQUATION 5-5a
C---------------------TT(M)=((A1+H*DELX/AK)*TAIR+T(MM))/(1.+A1+H*DELX/AK)
KOUNT=KOUNT+1
DO 205 I=2,M
205
IF(ABS((TT(I)-T(I))/T(I)).GT.TOL) GOTO 12
GOTO 302
12
DO 333 I=1,M
333
T(I)=TT(I)
GO TO 300
302
CONTINUE
C
C
call polyft(x,tt,z,p,q,c,m,mm,mn,std,1.)
C
qfin1=-c(2)*Ak*thick*100.
QFIN2=AK*THICK*(T(1)-T(2))/DELX+(.75*T(1)+.25*T(2)TAIR)*H*DELX
WRITE(3,604)
604
FORMAT(3X,'I',8X,'X(I)',4X,'T(I)',7X,'TT(I)')
WRITE(*,*)(I,X(I),T(I),TT(I),I=1,M)
WRITE(3,605)(I,X(I),T(I),TT(I),I=1,M)
605
FORMAT(I4,3x,f8.2,3x,f8.2,3x,f8.2)
Tol=Tol*100.0
WRITE(2,555) TOL
555
FORMAT(/,'the last two iterated values are within ',f7.2,'%',//)
WRITE(2,607) KOUNT
607
FORMAT(//,'solution was obtained after',I7,2X,'iterations',//)
WRITE(2,606)QFIN2
606
FORMAT('RATE OF HT.TR.(Finite Diff.))',f12.2,'watts')
STOP
END
C
THE PROGRAM USED TO CALCULATE THE AMOUNT OF HEAT
CONDUCTED
C
ACROSS THE PLAT
C
CALCULATE THE TEMPERATURE DISTRIBUTION
DIMENSION T(1000,1000),TT(1000,1000),X(1000,1000),Y(1000,1000)
OPEN(UNIT=1,FILE='C:/MSDEV/PROJ1.DAT')
OPEN(UNIT=2,FILE='C:/MSDEV/PROJ2.DAT')
OPEN(UNIT=3,FILE='C:/MSDEV/PROJ3.DAT')
READ(1,*)((T(M,N),I=1,M),J=1,N )
READ(1,*)TZERO,H,AK,THICK,TAIR,XA,YB
WRITE(2,102)TZERO,H,AK,THICK,TAIR,XA,YB
102
FORMAT(1h1,/,
+'Temperature at the lower plat =',F7.2,' deg.c',//,
+'Convective heat transfer coeff
+'Theramal conductivity
+'Thickness of the fin
+'Temp of the ambient air
=',F7.2,' watts/sq.m.,k',//,
=',F7.2,' watts/m.k.',//,
=',F5.2,' mm',//,
=',F6.2,' deg.c',//,
+'Length of the plat
=',F6.2,' cms',//,
+'width of the plat
=',F6.2,' cms',//)
XA=XA/100.0
YB=YB/100.0
THICK=THICK/1000.0
C
READ (1,*)M,N
C
IF(M.EQ.0) GOTO 500
DELX=XA/(M-1.)
DELY=YB/(N-1.)
B1=(H*DELX)/AK
READ(1,*)TOL
MM=M-1
NN=N-1
C------INITIAL TRIAL VALUE
C------------------------DO 203 I=1,M
DO 203 J=1,N
X(I)=DELX*(I-1)
Y(J)=DELY*(J-1)
203
T(I)=TZERO+(TAIR-TZERO)*((I-1.)*DELX/XA)
T(J)=TZERO+(TAIR-TZERO)*((J-1.)*DELX/YB)
KOUNT=0
C------ITERATED VALUE
C------------------300
TT(1,1)=TZERO
C --------EQUATION 5-2a
C----------------------DO 204 I=2,MM
DO 204 J=2,NN
204
TT(I,J)=(T(I-1,J)+T(I+1,J)+T(I,J-1))/4
C
C --------EQUATION 5-5a
C---------------------C---------FOR INSULATED NODE
C---------LEFT NODES
DO 2000 I=1
DO 2000 J=2,NN
2000 TT(I,J)=(T(I,J+1)+T(I,J-1)+(2*T(I+1,J)))/4
C---------RIGHT NODES
DO 3000 I=M
DO 3000 J=2,NN
3000 TT(I,J)=(T(I,J+1)+T(I,J-1)+(2*T(I-1,J)))/4
C---------BOTTOM NODES
DO 4000
I=2,3
DO 4000
J=1
4000 TT(I,J)=(T(I+1,J)+T(I-1,J)+(2*T(I,J+1)))/4
DO 5000 I=4,7
DO 5000 J=1
5000 TT(I,J)=TZERO
DO 6000 I=8,9
DO 6000 J=1
6000 TT(I,J)=(T(I+1,J)+T(I-1,J)+(2*T(I,J+1)))/4
C---------FOR CONVECTION NODES
C---------UPPER NODES
DO 7000 I=2,MM
DO 7000 J=N
7000 TT(I,J)=(T(I,J-1)+((T(I+1,J)+T(I-1,J))/2)+(B1*TAIR))/(2+B1)
C---------FOR INSULATED CORNER NODE
C---------BOTTOM RIGHT NODES
TT(1,M)=(T(I,J+1)+T(I-1,J))/2
C---------BOTTOM LEFT NODES
TT(1,1)=(T(I,J+1)+T(I+1,J))/2
C---------FOR JNSULATED & CONVECTION CORNER NODE
C---------RIGHT UPPER NODE
TT(M,N)=(T(I-1,J)+T(I,J-1)+(B1*TAIR))/(2+B1)
C---------LEFT UPPER NODE
TT(1,N)=(T(I+1,J)+T(I,J-1)+(B1*TAIR))/(2+B1)
KOUNT=KOUNT+1
DO 205 I=2,M
DO 205 J=2,N
205
IF(ABS((TT(I,J)-T(I,J))/T(I,J)).GT.TOL) GOTO 12
GOTO 302
12
DO 333 I=1,M
DO 333 J=1,N
333
T(I,J)=TT(I,J)
GO TO 300
302
CONTINUE
C
C
call polyft(x,tt,z,p,q,c,m,mm,mn,std,1.)
C
qfin1=-c(2)*Ak*thick*100.
C
QFIN2=AK*THICK*(T(1)-T(2))/DELX+(.75*T(1)+.25*T(2)TAIR)*H*DELX
C
WRITE(3,604)
C604 FORMAT(3X,'I',3X, 'J' ,8X,'X(I)',8X,'Y(J),4X,'T(I,J)',7X,'TT(I,J)
C +')
C
WRITE(*,*)(I,J,X(I),Y(J),T(I,J),TT(I,J),I=1,M,J=1,N)
WRITE(3,605)(I,J,X(I),Y(J),T(I,J),TT(I,J),I=1,M,J=1,N)
605
FORMAT(I4,3x,f8.2,3x,f8.2,3x,f8.2)
Tol=Tol*100.0
WRITE(2,555) TOL
555
FORMAT(/,'the last two iterated values are within ',f7.2,'%',//)
WRITE(2,607) KOUNT
607
FORMAT(//,'solution was obtained after',I7,2X,'iterations',//)
WRITE(2,606)QFIN2
606
FORMAT('RATE OF HT.TR.(Finite Diff.))',f12.2,'watts')
STOP
END
