forums.silverfrost.com

[yuvadon]

This is my final year project...I have to run the program using silverfrost ftn95 in order to complete my project.But as I compiled it shows to much of errors.I have no idea about it.So, could you guys please help me...Your assistance is highly appreciated.Thank you.( I need to submit the result by this Thursday)

so this is the code.

c#######################################################################
c
c SPH: Smoothed Particle Hydrodynamics
c
c#######################################################################
c
program SPH
c
implicit real*8(a-h,o-z)
c
include 'parameter.f'
c
parameter(
& TOTAL_TIME=5.0D-6, ! �V�~�����[�V�����̍ő厞��
& MAXLOOP=100000) ! �V�~�����[�V�����̍ő�J��Ԃ���
c
common/model/ DISA,XMAX,XMIN,XSIZE,YMAX,YMIN,YSIZE,
& NUMX1,NUMY1,NUMX2,NUMY2,NUMP,NP1,NP2
c
common/particle/ X(NDIM,MAXP), X0(NDIM,MAXP)
c
common/property/ RHO1,RHO2,CS1,CS2,G1,G2,
& RHO(MAXP),PMASS(MAXP),CS(MAXP)
c
common/composite/ CMAT(3,3,2),MATER(MAXP)
c
common/boundary/
& V_BC_GIVEN(NDIM,MAXBC),V_LOAD_GIVEN(NDIM,MAXBC),
& I_BC_GIVEN(NDIM,MAXBC),I_LOAD_GIVEN(NDIM,MAXBC),
& N_BC_GIVEN(NDIM),N_LOAD_GIVEN(NDIM),SPEED
c
common/calc/ V(NDIM,MAXP),VISC(NEIGHBOR,MAXP),P(MAXP),E(MAXP)
c
common/param/ H0,H,DT,PDT
c
common/neighbor_table/
& NEIGH(NEIGHBOR,MAXP),NNEIGH(MAXP),
& NEIGH0(NEIGHBOR,MAXP),NNEIGH0(MAXP),JCONNECT0(NEIGHBOR,MAXP),
& RIJ(NDIM,NEIGHBOR,MAXP),RIJ1(NEIGHBOR,MAXP),RIJ2(NEIGHBOR,MAXP),
& RIJ0(NEIGHBOR,MAXP)
c
common/kernel/ GRADW(NDIM,NEIGHBOR,MAXP)
c
common/stress/
& SIGM(3,MAXP),HENSAS(3,MAXP),HENSASOLD(3,MAXP),STRNDOT(3,MAXP)
c
common/vardot/ DRHODT(MAXP),DVDT(NDIM,MAXP),DEDT(MAXP)
c
common/rkgill/ DX(NDIM,MAXP),DV(NDIM,MAXP),DR(MAXP),DE(MAXP)
c
common/yield/ SGMY(MAXP),EPSP(MAXP),T(MAXP),JYIELD(MAXP)
c
common/damage/ DMG(MAXP)
c
c==== output files
open(2,file='rst.txt')
open(3,file='out_damage.plt')
open(4,file='rst_damage.csv')
open(7,file='out_stress.plt')
open(8,file='out_pstress.plt')
open(100,file='check_coordinates.csv')
c
c==== input
call inp(angle_inc)
c
call set_param
c
call set_constitutive_matrix
c
c==== initialize
do ip=1,NUMP
X(1:NDIM,ip) = X0(1:NDIM,ip) ! ���W
E(ip) = 0.0D0 ! �G�l���M�[
do k=1,3
SIGM(k,ip) = 0.0D0 ! ����
HENSAS(k,ip) = 0.0D0 ! �΍�����
HENSASOLD(k,ip) = 0.0D0
STRNDOT(k,ip) = 0.0D0 ! �Ђ��ݑ��x
enddo
c
DX(1:NDIM,ip) = 0.0D0 ! �����Q�E�N�b�^�@
DV(1:NDIM,ip) = 0.0D0 ! �ł̎��Ԑϕ��Ŏg��
DR(ip) = 0.0D0
DE(ip) = 0.0D0
c
SGMY(ip) = YIELD0 ! �~������
EPSP(ip) = 0.0D0 ! �����Y���Ђ���
T(ip) = 300.0D0 ! ���x
c
DMG(ip) = 0.0D0 ! �����p�����[�^
enddo
c
c---- define boundary conditions
call definebc
c
V(1:NDIM,1:NUMP) = 0.0D0 ! ���x�̏�����
c
do ip=NP1+1,NUMP ! ���đ̂̏���
V(1,ip) = SPEED * dcos(angle_inc)
V(2,ip) = SPEED * dsin(angle_inc)
enddo

[IanLambley]

Tidied up a bit and unicode characters removed.

[IanLambley]