Fortran Program for Regional-Residual Separation by Finite Element Method
Appendix I 259 Appendix I Fortran Program for Regional-Residual Separation by Finite Element Method C AG-BOUGUER,GN-REGIONAL,D-RESIDUAL C PROGRAM FOR SEPARATION OF REGIONAL & RESIDUAL USING FEM METHOD C PROGRAM NEEDS THE GRAVITY VALUES ALSO IN ADDITION TO THE A(I)& B(I)VALUES DIMENSION G(12),X(12),Y(12) DIMENSION A(200000),B(200000),GN(200000),AG(200000),D(200000) OPEN(UNIT=3,FILE=‘F1-8.DAT’, STATUS=‘OLD’) OPEN(UNIT=4,FILE=‘A1.DAT’, STATUS=‘NEW’) OPEN(UNIT=11,FILE=‘A2.DAT’, STATUS=‘NEW’) OPEN(UNIT=12,FILE=‘A3.DAT’, STATUS=‘NEW’) OPEN(UNIT=13,FILE=‘A4.DAT’, STATUS=‘NEW’) READ (3,*) N,NN WRITE (*,*)‘GIVE OUTPUT DATA FILE NAME AS UNIT 4’ WRITE (4,*)‘NUMBER OF NODES{A(I)& B(I)} AND DATA(8 or 12) POINTS’ WRITE (4,*)N,NN WRITE (4,*)‘INPUT VALUES X,Y,AND GRAVITY DATA’ DO 10 I=1,NN READ (3,*) X(I),Y(I),G(I) 10 WRITE (4,*) X(I),Y(I),G(I) WRITE (4,*)‘INPUT VALUES OF A(I) & B(I)’ DO 20 I=1,N READ (3,*) A(I),B(I),AG(I) 20 WRITE (4,*)A(I),B(I),AG(I) WRITE (4,*)‘OUTPUT X , Y , REGIONAL GRAVITY & RESIDUAL VALUES’ DO 100 I=1,N A1=1+A(I) A2=1-A(I) AA=1-A(I)*A(I) K. Mallick et al., Bouguer Gravity Regional and Residual Separation: Application to Geology 259 and Environment, DOI 10.1007/978-94-007-0406-0, © Capital Publishing Company 2012 260 Bouguer Gravity Regional and Residual Separation B1=1+B(I) B2=1-B(I) BB=1-B(I)*B(I) C1=(9*AA)/32 C2=(9*BB)/32 C3=(-10+9*(2-AA-BB))/32 AN1=0.25*A2*B2*(A2+B2-3) AN3=0.25*A1*B2*(B2-A2-1) AN5=0.25*A1*B1*(1-A2-B2) AN7=0.25*A2*B1*(A2-B2-1) AN2=0.5*AA*B2 AN6=0.5*AA*B1 AN4=0.5*BB*A1 AN8=0.5*BB*A2 T1=AN1*G(1)+AN2*G(2)+AN3*G(3)+AN4*G(4)
[Show full text]