Register to reply 
Fortran helpby rajpintoo
Tags: fortran 
Share this thread: 
#1
Mar1212, 06:42 PM

P: 9

MODULE space_data
REAL :: l !LENGTH OF THE BAR REAL :: temp !TEMPERATURE IN KELVIN REAL :: d_0 !CONSTANT DIFFUSION PARAMETER REAL :: q !ACTIVATION ENERGY REAL :: c_0 !INITIAL CONCENTRATION AT X=0 REAL :: r=8.31 !IDEAL GAS CONSTANT IN J/molK REAL :: c_r !DESIRED END CONCENTRATION INTEGER :: imax=30 !MAXIMUM ITERATIONS REAL :: eps=1E4 !MAXIUMUM ERROR REAL :: a, b !INTERVAL OF EVALUATION [A,B] REAL :: t REAL :: exp, df !, erf, c, c_avg END MODULE space_data PROGRAM carburization USE space_data IMPLICIT NONE CHARACTER :: choice REAL :: f, root INTEGER :: ercode DO WRITE(*,*) 'a) Enter new data' WRITE(*,*) 'b) Calculate time until desired concentration' WRITE(*,*) 'q) Quit' WRITE(*,'(A,\)') 'Enter choice: ' READ(*,*) choice SELECT CASE(choice) CASE('A','a') CALL call_data() CASE('B','b') CALL bisection(f, root, ercode) CALL results() CASE('Q','q') STOP END SELECT END DO CONTAINS SUBROUTINE call_data() USE space_data IMPLICIT NONE WRITE(*,'(A,\)') 'Enter initial end concentration (%): ' READ(*,*) c_0 WRITE(*,'(A,\)') 'Enter diffusion parameter (m^3/sec): ' READ(*,*) d_0 WRITE(*,'(A,\)') 'Enter activation energy (J): ' READ(*,*) q WRITE(*,'(A,\)') 'Enter temperature (K): ' READ(*,*) temp WRITE(*,'(A,\)') 'Enter length of bar (m): ' READ(*,*) l WRITE(*,'(A,\)') 'Enter desired average final concentration (%): ' READ(*,*) c_r WRITE(*,*) '' WRITE(*,*) ' The bisection method needs an initial search interval [A,B].' WRITE(*,*) '' DO WRITE(*,'(A,\)') 'Enter a A (Years):' READ(*,*) a WRITE(*,'(A,\)') 'Enter a B (Years):' READ(*,*) b IF(a < b) THEN EXIT ELSE WRITE(*,*) 'Error left end point needs to be smaller than right end point.' WRITE(*,*) 'Please enter again:' END IF END DO END SUBROUTINE call_data END PROGRAM carburization SUBROUTINE bisection(f, root, ercode) USE space_data IMPLICIT NONE INTEGER::I REAL,EXTERNAL::F REAL,INTENT(OUT)::ROOT !THE ROOT WE ARE TRYING TO FIND REAL::X1,X2,X3,F1,F2,F3,D,D01 !X VALUES, F(X) VALUES AND INTERVAL LENGTHS INTEGER,INTENT(OUT)::ERCODE !ERROR CODE WILL RETURN AN INT VALUE BASED ON OUTCOME !1 MEANS ITERATION CNT EXCEED !2 MEANS NO ROOT IN INTERVAL AB !3 UNKNOWN ERROR ! 0 ROOT FOUND, BISECTION METHOD SUCCESSFUL X1=A X3=B F1 = F(X1) F3 = F(X3) IF(F1*F3>0) THEN ERCODE = 2 !2 MEANS NO ROOT IN INTERVAL AB ELSE D=1 I=0 D01=AB DO X2 = (X1+X3)/2. F2 = F(X2) IF(D<EPS) THEN ROOT = X2 ERCODE = 0 !SUCCESS EXIT ELSE IF(I>IMAX) THEN ERCODE = 1 !1 MEANS ITERATION CNT EXCEED EXIT END IF IF(F1*F2<0) THEN !LEFT D = (X2X1)/D01 F3 = F2 X3 = X2 ELSE IF(F2*F3<0)THEN!RIGHT D = (X3X2)/D01 F1 = F2 X1 = X2 ELSE IF(F2 == 0) THEN ROOT = X2 ERCODE = 0 !SUCCESS EXIT ELSE ERCODE = 3 !3 UNKNOWN ERROR EXIT END IF I = I+1 END DO END IF END SUBROUTINE SUBROUTINE results() USE space_data IMPLICIT NONE WRITE(*, '(T5,A)') 'RESULTS:' WRITE(*, '(T5, 50(""))') WRITE(*,'(T5, "", 1X, A,T40, "", F12.3, 1X, "")') 'Initial end concentration (%)',c_0 WRITE(*,'(T5, "", 1X, A,T40, "", ES12.3, 1X, "")') 'Diffusion parameter (m^3/sec):',d_0 WRITE(*,'(T5, "", 1X, A,T40, "", ES12.3, 1X, "")') 'Activationenergy (J):',q WRITE(*,'(T5, "", 1X, A,T40, "", F12.2, 1X, "")') 'Temperature (K):',temp WRITE(*,'(T5, "", 1X, A,T40, "", F12.2, 1X, "")') 'Length of bar (m):',l WRITE(*,'(T5, "", 1X, A,T40, "", F12.3, 1X, "")') 'Final concentration (%):',c_r WRITE(*, '(T5, 50(""))') WRITE(*,'(T5, "", 1X, A,T40, "", F12.3, 1X, "")') 'Computed diffusion time (days):',t WRITE(*, '(T5, 50(""))') WRITE(*,*) END SUBROUTINE results REAL FUNCTION f() USE space_data IMPLICIT NONE f = c_avg  c_r END FUNCTION f REAL FUNCTION c () USE space_data IMPLICIT NONE REAL :: z, x df = d_0*(exp**(q/(r*temp))) x=l z = x/(2*df*t) c = (c_0/2.)*(1  ERF(z)) END FUNCTION c REAL FUNCTION c_avg () USE space_data IMPLICIT NONE REAL :: delta_x, idx INTEGER :: i delta_x = .2 DO i = 0, 5 idx = i*delta_x c_avg = c_avg + (1/6.)*( c(idx,t) ) END DO END FUNCTION c_avg 


#2
Mar1212, 06:44 PM

P: 9

Can someone just error check this and explain whats wrong the compiler just tells me random errors which I dont know



#3
Mar1312, 10:50 AM

Mentor
P: 21,397

Tell us what the errors are  they are probably not "random."



#4
Mar1312, 11:12 AM

P: 9

Fortran help
Hmm alriight it says that C_avg has no implicit type



#5
Mar1312, 12:31 PM

P: 20

Hey I think you're in my class. c_avg is a function but you defined it as just a REAL in the module. I didn't define any functions in the module because I kept getting errors. You also have to have input variables for your function's dummy variables. So c_avg would just be a function of time because we only have to calculate the concentration at the end of the bar (x = Length of bar). I don't have any compiler errors but I keep getting the wrong answer so I just have something wrong with the calculation.
REAL FUNCTION c_avg (t) USE space_data IMPLICIT NONE REAL, INTENT(IN):: t REAL :: delta_x, idx INTEGER :: i delta_x = .2 DO i = 0, 5 idx = i*delta_x c_avg = c_avg + (1/6.)*( c(idx,t) ) END DO END FUNCTION c_avg EDIT I just saw this in the module REAL :: exp, df !, erf, c, c_avg It might not be reading what's after the exclamation point because it counts it as a comment. 


#6
Mar1312, 12:57 PM

P: 9

K put in the t function stuff but I am still getting errors, how did you define your function within the main program? also what do you mean imput variables for dummy variables?



#7
Mar1312, 01:02 PM

P: 9

How exactly are you supposed to do the functions because thats what causing me the problems



#8
Mar1312, 01:06 PM

P: 9

oh for your computation problem was this the issue
The initial search interval [A,B] is given in years. The calculations require you to convert these two into seconds. Likewise, the answer you get will be in seconds. You have to convert it back to years before you print it out. 


#9
Mar1312, 01:51 PM

Mentor
P: 21,397

I'm not sure what c_r is. If it's another function, I didn't find it, but I didn't look that carefully at your code. If c_r is a function, you aren't calling it correctly, either.
I see problems in this code as well  t is undeclared. When you use IMPLICIT NONE, you HAVE TO declare all variables. The other problem is that you are calling your c function and passing two arguments to it, but the definition of this function shows that it has no parameters.



#10
Mar1312, 01:58 PM

P: 9

hmm so all I have to do in terms of putting a function within a function is to put the parenthesis afterwards but wont it thinks its an array or something?



#11
Mar1312, 02:06 PM

P: 9

MODULE space_data
REAL :: L !LENGTH OF THE BAR REAL :: temp !TEMPERATURE IN KELVIN REAL :: d_0 !CONSTANT DIFFUSION PARAMETER REAL :: q !ACTIVATION ENERGY REAL :: c_0 !INITIAL CONCENTRATION AT X=0 REAL :: r=8.31 !IDEAL GAS CONSTANT IN J/molK REAL :: c_r !DESIRED END CONCENTRATION INTEGER :: imax=30 !MAXIMUM ITERATIONS REAL :: eps=1E4 !MAXIUMUM ERROR REAL :: a, b !INTERVAL OF EVALUATION [A,B] END MODULE space_data PROGRAM carburization USE space_data IMPLICIT NONE CHARACTER :: choice REAL :: f, root INTEGER :: ercode DO WRITE(*,*) 'a) Enter new data' WRITE(*,*) 'b) Calculate time until desired concentration' WRITE(*,*) 'q) Quit' WRITE(*,'(A,\)') 'Enter choice: ' READ(*,*) choice SELECT CASE(choice) CASE('A','a') CALL call_data() CASE('B','b') CALL bisection(f, root, ercode) CALL results() CASE('Q','q') STOP END SELECT END DO CONTAINS SUBROUTINE call_data() USE space_data IMPLICIT NONE WRITE(*,'(A,\)') 'Enter initial end concentration (%): ' READ(*,*) c_0 WRITE(*,'(A,\)') 'Enter diffusion parameter (m^3/sec): ' READ(*,*) d_0 WRITE(*,'(A,\)') 'Enter activation energy (J): ' READ(*,*) q WRITE(*,'(A,\)') 'Enter temperature (K): ' READ(*,*) temp WRITE(*,'(A,\)') 'Enter length of bar (m): ' READ(*,*) l WRITE(*,'(A,\)') 'Enter desired average final concentration (%): ' READ(*,*) c_r WRITE(*,*) '' WRITE(*,*) ' The bisection method needs an initial search interval [A,B].' WRITE(*,*) '' DO WRITE(*,'(A,\)') 'Enter a A (Years):' READ(*,*) a WRITE(*,'(A,\)') 'Enter a B (Years):' READ(*,*) b IF(a < b) THEN EXIT ELSE WRITE(*,*) 'Error left end point needs to be smaller than right end point.' WRITE(*,*) 'Please enter again:' END IF END DO END SUBROUTINE call_data END PROGRAM carburization SUBROUTINE bisection(f, root, ercode) USE space_data IMPLICIT NONE INTEGER::I REAL,EXTERNAL::F REAL,INTENT(OUT)::ROOT !THE ROOT WE ARE TRYING TO FIND REAL::X1,X2,X3,F1,F2,F3,D,D01 !X VALUES, F(X) VALUES AND INTERVAL LENGTHS INTEGER,INTENT(OUT)::ERCODE !ERROR CODE WILL RETURN AN INT VALUE BASED ON OUTCOME !1 MEANS ITERATION CNT EXCEED !2 MEANS NO ROOT IN INTERVAL AB !3 UNKNOWN ERROR ! 0 ROOT FOUND, BISECTION METHOD SUCCESSFUL X1=A X3=B F1 = F(X1) F3 = F(X3) IF(F1*F3>0) THEN ERCODE = 2 !2 MEANS NO ROOT IN INTERVAL AB ELSE D=1 I=0 D01=AB DO X2 = (X1+X3)/2. F2 = F(X2) IF(D<EPS) THEN ROOT = X2 ERCODE = 0 !SUCCESS EXIT ELSE IF(I>IMAX) THEN ERCODE = 1 !1 MEANS ITERATION CNT EXCEED EXIT END IF IF(F1*F2<0) THEN !LEFT D = (X2X1)/D01 F3 = F2 X3 = X2 ELSE IF(F2*F3<0)THEN!RIGHT D = (X3X2)/D01 F1 = F2 X1 = X2 ELSE IF(F2 == 0) THEN ROOT = X2 ERCODE = 0 !SUCCESS EXIT ELSE ERCODE = 3 !3 UNKNOWN ERROR EXIT END IF I = I+1 END DO END IF END SUBROUTINE SUBROUTINE results() USE space_data IMPLICIT NONE WRITE(*, '(T5,A)') 'RESULTS:' WRITE(*, '(T5, 50(""))') WRITE(*,'(T5, "", 1X, A,T40, "", F12.3, 1X, "")') 'Initial end concentration (%)',c_0 WRITE(*,'(T5, "", 1X, A,T40, "", ES12.3, 1X, "")') 'Diffusion parameter (m^3/sec):',d_0 WRITE(*,'(T5, "", 1X, A,T40, "", ES12.3, 1X, "")') 'Activationenergy (J):',q WRITE(*,'(T5, "", 1X, A,T40, "", F12.2, 1X, "")') 'Temperature (K):',temp WRITE(*,'(T5, "", 1X, A,T40, "", F12.2, 1X, "")') 'Length of bar (m):',L WRITE(*,'(T5, "", 1X, A,T40, "", F12.3, 1X, "")') 'Final concentration (%):',c_r WRITE(*, '(T5, 50(""))') WRITE(*,'(T5, "", 1X, A,T40, "", F12.3, 1X, "")') 'Computed diffusion time (days):' WRITE(*, '(T5, 50(""))') WRITE(*,*) END SUBROUTINE results REAL FUNCTION Con() USE Space_data IMPLICIT NONE REAL :: t,z,x,DF,delta_x,idx INTEGER :: i delta_x=L/5 t=(BA)*31556926 DF=D_0*EXP(q/(R*temp)) DO i = 0,5 idx=i*delta_x z = idx/(2*DF*t) con = (c_0/2.)*(1  ERF(z)) END DO END FUNCTION REAL FUNCTION C_avg() USE Space_data IMPLICIT NONE REAL,EXTERNAL :: Con C_avg=(1/6)*Con() END FUNCTION REAL FUNCTION f() USE Space_data IMPLICIT NONE REAL,EXTERNAL :: C_avg f = C_avg()C_r END FUNCTION 


#12
Mar1312, 04:05 PM

P: 20

The t (time) is either A or B which you have to ask for as the input (10 and 50 in the pdf). When the program runs, it will evaluate the function(s) for these two times, but they have to be in seconds. We only care about the end of the bar so X or delta X will always just be the length of the bar, then these are only functions of time.
I don't want to post my code here but if you go to the lab on Thursday I can show you how I did the functions. Oh yeah, it says the time input is in years, but it's actually days, so when you convert A and B to seconds, make sure you convert from days to seconds, not years to seconds. 


#13
Mar1312, 04:05 PM

P: 9

So my code compiles but now i get an error when I put in data and enter choice b
it says exception access violation traceback not available 


#14
Mar1412, 10:37 AM

Mentor
P: 21,397

When you call a function there is more to it than just putting a pair of parentheses after the function name. The number of parameters in the function call has to agree with the number of parameters in the function definition. Also, the types of the parameters have to agree. 


#15
Mar1412, 10:42 AM

Mentor
P: 21,397

I'm guessing that the exception you're seeing results from attempting to access an array using an index that's out of range. 


#16
Mar1612, 01:50 PM

P: 1

Did anyone find a solution to this? my project has the same error.



#17
Mar1612, 06:51 PM

Mentor
P: 21,397

Did you try what I suggested a couple of posts ago?



Register to reply 
Related Discussions  
Fortran 77 help making an empty array (or blank list if they exist in fortran)  Programming & Computer Science  5  
FORTRAN 95: New to fortran, want to learn how to input a function  Programming & Computer Science  1  
Fortran vs C for CFD  Mechanical Engineering  2  
Accessing Fortran Modules within a Fortran library from Fortran  Programming & Computer Science  0  
Fortran v.s. Visual Fortran  Programming & Computer Science  0 