invalid character, no implicit type, unclassifiable statemen

------------------------------------------------------------------------------ * C Driver program * IMPLICIT NONE * C Declaration of variables: * INTEGER ITEST, CELLS, N, NFREQ, NTMAXI * REAL CFLCOE, DOMLEN, DT, TIME, TIMEOUT, TIMETO * COMMON / DATAIN/ CFLCOE, DOMLEN, ITEST, CELLS, & NFREQ, NTMAXI, TIMEOUT COMMON /DELTAT/ DT * DATA TIMETO /1.0E-07/ * C Parameters of problem are read in form C file 'b1god.ini' * CALL READER * C initial conditions are set up * CALL INITIA(DOMLEN, ITEST, CELLS) * WRITE(6,)'---------------------------------------------------------------' WRITE(6,)' Time step N TIME TIMEOUT' WRITE(6,)'---------------------------------------------------------------' * C Time marching procedure * TIME = 0.0 * DO 10 N = 1, NTMAXI * C Boundary conditions are set * CALL BCONDI(CELLS) * C Courant-Friedrichs-Lewy (CFL) condition imposed * CALL CFLCON(CFLCOE, CELLS, TIME, TIMEOUT) * TIME = TIME + DT * C Intercell numerical fluxes are computed * CALL FLUXES(CELLS) * C Solution is updated according to c conservative formula * CALL UPDATE(CELLS) * IF(MOD(N,NFREQ).EQ.0)WRITE(6,20)N, TIME * C Check output time * IF(ABS(TIME - TIMEOUT).LE.TIMETO)THEN * C Solution is written to 'numer.out? at C specified time TIMEOUT * CALL OUTPUT(CELLS) * WRITE(6,)'---------------------------------------------------------' WRITE(6,)' Number of time steps = ',N * STOP ENDIF * 10 CONTINUE * 20 FORMAT(I12,6X, F12.7) * END * ------------------------------------------------------------------------------ * SUBROUTINE READER * C Purpose: to read initial parameters of the problem * IMPLICIT NONE * C Declaration of variables * INTEGER ITEST, CELLS, NFREQ, NTMAXI * REAL CFLCOE, DOMLEN, TIMEOUT * COMMON /DATAIN/ CFLCOE, DOMLEN, ITEST, CELLS, NFREQ, & NTMAXI, TIMEOUT * OPEN(UNIT = 1,FILE = 'b1god.ini',STATUS = 'UNKNOWN') * READ(1,)CFLCOE ! Courant number coefficient READ(1,)DOMLEN ! Domain length READ(1,)ITEST ! Test problem READ(1,)CELLS ! Number of cells in domain READ(1,)NFREQ ! Output frequency to screen READ(1,)NTMAXI ! Maximum number of time steps READ(1,)TIMEOUT ! Output time * CLOSE(1) * WRITE(6,)'--------------------------------' WRITE(6,)'Data read in is echoed to screen' WRITE(6,)'--------------------------------' WRITE(6,)'CFLCOE = ',CFLCOE WRITE(6,)'DOMLEN = ',DOMLEN WRITE(6,)'ITEST = ',ITEST WRITE(6,)'CELLS = ',CELLS WRITE(6,)'NFREQ = ',NFREQ WRITE(6,)'NTMAXI = ',NTMAXI WRITE(6,)'TIMEOUT = ',TIMEOUT WRITE(6,*)'--------------------------------' * RETURN END * ------------------------------------------------------------------------------ * SUBROUTINE INITI(DOMLEN, ITEST, CELLS) * C Purpose: to set initial conditions for solution U C and initialise other variables. There are C two choices of initial conditions, C determined by ITEST * C Local variables: * C Name Description

• ====          ===========
  

C DX Spatial mesh size C I Variable in do loop C ITEST Defines test problem C FLUX Array for intercell fluxes C U Array for numerical solution C XPOS Position along x-axis C XLEFT Left diaphragm C XMIDDL Middle diaphragm C XRIGHT Right diaphragm * IMPLICIT NONE * C Declaration of variables * INTEGER I, ITEST, CELLS, IDIM * REAL DOMLEN, DX, FLUX, U, XLEFT, XPOS, XMIDDL, & XRIGHT * PARAMETER (IDIM = 1000) * DIMENSION FLUX(0:IDIM + 1), U(0:IDIM + 1) * COMMON /DELTAX/ DX COMMON /FLUXFS/ FLUX COMMON /SOLUTI/ U * C Calculate mesh size DX * DX = DOMLEN/REAL(CELLS) * C Initialise arrays * DO 10 I = 0, IDIM + 1 FLUX(I) = 0.0 U(I) = 0.0 10 CONTINUE * IF(ITEST.EQ.1)THEN * C Test 1: smooth profile * XPOS = -1.0 * DO 20 I = 1, CELLS XPOS = XPOS + 2.0/REAL(CELLS) U(I) = EXP(-8.0XPOSXPOS) 20 CONTINUE * ELSE * C Test 2: square waves * XLEFT = 0.1DOMLEN XMIDDL = 0.5DOMLEN XRIGHT = 0.9DOMLEN * DO 30 I = 1, CELLS * XPOS = (REAL(I)-1.0)DX * IF(XPOS.LT.XLEFT)THEN U(I) = -1.0 ENDIF * IF(XPOS.GE.XLEFT.AND.XPOS.LE.XMIDDL)THEN U(I) = 1.0 ENDIF * IF(XPOS.GT.XMIDDL.AND.XPOS.LE.XRIGHT)THEN U(I) = 0.0 ENDIF * IF(XPOS.GT.XRIGHT)THEN U(I) = -1.0 ENDIF * 30 CONTINUE * ENDIF * RETURN END * ------------------------------------------------------------------------------ * SUBROUTINE BCONDI(CELLS) * C Purpose: to apply boundary conditions * IMPLICIT NONE * C Declaration of variables * INTEGER CELLS, IDIM * REAL U * PARAMETER (IDIM = 1000) * DIMENSION U(0:IDIM + 1) * COMMON /SOLUTI/ U * C Left boundary, periodic boundary condition * U(0) = U(CELLS) * C Right boundary, periodic boundary condition * U(CELLS + 1) = U(1) * RETURN END * ------------------------------------------------------------------------------ * SUBROUTINE CFLCON(CFLCOE, CELLS, TIME, TIMEOUT) * C Purpose: to apply the CFL condition to compute a C stable time step DT * IMPLICIT NONE * C Declaration of variables * INTEGER I, CELLS, IDIM * REAL CFLCOE, DT, DX, SMAX, TIME, TIMEOUT, U * PARAMETER (IDIM = 1000) * DIMENSION U(0:IDIM + 1) * COMMON /SOLUTI/ U COMMON /DELTAT/ DT COMMON /DELTAX/ DX * SMAX = -1.0E+06 * C Find maximum characteristic speed * DO 10 I = 0, CELLS + 1 IF(ABS(U(I)).GT.SMAX)SMAX = ABS(U(I)) 10 CONTINUE * DT = CFLCOEDX/SMAX * C Check size of DT to avoid exceeding output time * IF((TIME + DT).GT.TIMEOUT)THEN * C Recomplete DT * DT = TIMEOUT - TIME ENDIF * RETURN END * ------------------------------------------------------------------------------ * SUBROUTINE UPDATE(CELLS) * C Purpose: to update the solution to a new time level C using the explicit conservative formula * IMPLICIT NONE * C@@Declaration of variables * INTEGER I, CELLS, IDIM * REAL DT, DX, DTODX, FLUX, U * PARAMETER (IDIM = 1000) * DIMENSION U(0:IDIM + 1), FLUX(0:IDIM + 1) * COMMON /DELTAT/ DT COMMON /DELTAX/ DX COMMON /FLUXFS/ FLUX COMMON /SOLUTI/ U * DTODX = DT/DX * DO 10 I = 1, CELLS U(I) = U(I) + DTODX(FLUX(I-1) - FLUX(I)) 10 CONTINUE * RETURN END * ------------------------------------------------------------------------------ * SUBROUTINE OUTPUT(CELLS) * C Purpose: to output the solution at specified time C TIMEOUT * IMPLICIT NONE * C Declaration of variables * INTEGER I, CELLS, IDIM * REAL DX, U, XPOS * PARAMETER (IDIM = 1000) * DIMENSION U(0:IDIM + 1) * COMMON /DELTAX/ DX COMMON /SOLUTI/ U * OPEN(UNIT = 1,FILE = 'numer.out',STATUS = 'UNKNOWN') * DO 10 I = 1, CELLS XPOS = REAL(I)DX WRITE(1,20)XPOS, U(I) 10 CONTINUE * CLOSE(1) * 20 FORMAT(2(4X, F10.5)) * RETURN END * ------------------------------------------------------------------------------ * SUBROUTINE FLUXES(CELLS) * C Purpose: to compute intercell fluxes according to C the Godunov first-order upwind method, C in conjunction with the exact Riemann C solver * IMPLICIT NONE * C Declaration of variables * INTEGER I, CELLS, IDIM * REAL FLUX, U, UL, UR, USTAR * PARAMETER (IDIM = 1000) * DIMENSION FLUX(0:IDIM + 1), U(0:IDIM + 1) * COMMON /FLUXFS/ FLUX COMMON /SOLUTI/ U * C Compute intercell flux FLUX(I), I = 0, CELLS C Solution of Riemann problem RP(I, I+1) is stored C in FLUX(I) * DO 10 I =@0, CELLS * C Define states UL (Left) and UR (Right) for local C Riemann problem RP(UL, UR) * UL = U(I) UR = U(I+1) * C Solve the Riemann problem RP(UL, UR) exactly * CALL RIEMANN(UL, UR, USTAR) * C Compute Godunov intercell flux * FLUX(I) = 0.5USTAR*USTAR * 10 CONTINUE * RETURN END * ------------------------------------------------------------------------------ * SUBROUTINE RIEMANN(UL, UR, USTAR) * C Purpose: to solve the Riemann problem for the inviscid C BurgerS equation exactly. * C Local variables: * C Name Description

• ==== ===========

C UL Left data state C UR Right data state C S Shock speed C USTAR Sampled state * IMPLICIT NONE * REAL S, UL, UR, USTAR * IF(UL.GT.UR)THEN * C Solution is a shock wave C Compute shock speed S * S = 0.5*(UL + UR) * C Sample the state along the t-axis * IF(S.GE.0.0)THEN USTAR = UL ELSE USTAR = UR ENDIF * ELSE * C Solution is a rarefaction wave. C There are 3 cases: * IF(UL.GE.0.0)THEN * C Right supersonic rarefaction * USTAR = UL ENDIF * IF(UR.LE.0.0)THEN * C Left supersonic rarefaction * USTAR = UR ENDIF * IF(UL.LE.0.0.AND.UR.GE.0.0)THEN * C Transonic rarefaction * USTAR = 0.0 ENDIF * ENDIF * RETURN END * ------------------------------------------------------------------------------ *

Posting long source codes in-line in this forum does not work.

There is a limit on the size of a post, which your code exceeded. As you can see, there is a missing portion of your code at the end.

If you post code without [code]...[/code] markers, leading blanks in the code get stripped, leaving the code in a state unfit to be compiled.

It is best to upload the code in a Zip file to a cloud storage site such as Dropbox, make the file accessible to the public, and post a link to the file here.

I suspect that the problem with your code is that it has some lines that exceed 72 in length, but let's hold off the diagnosis until we see an uncorrupted file.

I'd recommend that you don't use unit 1 to open a file as this can conflict with the screen input/output. I'd suggest you use unit 11. Also unit 6 is not the file unit for the screen with FTN95 (it is with some other compilers)

With FTN95, either use the following to a file opened on unit 6 WRITE(6,)'---------------------------------------------------------' WRITE(6,)' Number of time steps = ',N

or use the following if you want to write to the screen WRITE(,)'---------------------------------------------------------' WRITE(,)' Number of time steps = ',N

As a general rule, I only open a file with a unit number > 10

It is difficult to tell if you are using free or fixed format code layout, as you appear to be mixing both. Free format can be preferable for new code, so either change: * C Time marching procedure *

to

!* !C Time marching procedure !*

or use the fixed format code layout restrictions for

• comments starting with C or * in column 1 only
• statement numbers are in columns 1 to 5
• continuation in column 6 and
• code in columns 7 to 72

free format is much easier !!

I attach link for this complete file

https://www.dropbox.com/s/5v8m0qo29kcxzts/Dropbox.zip?dl=0

I hope someone can fix it, i have tried so much time but always error

There are many lines in the source file (which is contained in the Zip file) that do not obey the rules of fixed format Fortran source. Specifically, many lines begin before column 7. A few lines contains Format strings that continue beyond column 72. Either shorten these or use the /WI compiler option.

Line 395 in Subroutine FLUXES appears to be corrupted:

DO 10 I =ü@0, CELLS

Replacing this line by

DO 10 I=1,CELLS

and reformatting the source code to make statements conform to fixed format rules allowed the file to be compiled and linked.

Actually I have follow this program from reference book but also so many error.

If any one can fix it I will appreciate

This the link for book reference and program file https://www.dropbox.com/s/v0xtavw3wdde493/Dropbox%20%281%29.zip?dl=0

If it is fixed format (old) Fortran code then you should save it in a file with the extension .for or use /FIXED on the FTN95 command line.

Can you give an example?

I typed this program online in this link https://www.tutorialspoint.com/compile_fortran_online.php

sorry,

there is no way I will look at a .php file That I don't know about

Riogurky:

I think that you are attempting to do something that is beyond the reach of your capabilities. CFD is an application area that requires a high level of education, mathematical ability and experience. I recognize that you are attempting to learn by running some codes from a textbook, but you are unlikely to succeed at this point in time.

The second Zip file that you linked to contains (i) a copy of the buggy source code and (ii) a PDF document on how to use Dropbox!

The link to Tutorialspoint takes one to the entry page of a Cloud site for compiling and running short programs using Gfortran, and is irrelevant to what you are trying to do if you have a Windows PC and have FTN95 installed on it.

Please spend some time educating yourself regarding Fortran and Fortran programming on Windows. It is far wiser to do that instead of lurching about in the dark and asking others to fix the buggy code for you.

I am sorry mecej4, actually the reference link is: https://www.dropbox.com/s/vheugh95762qp3g/Fortran%20program%20for%20Gudunov%20Method.pdf?dl=0

I learn and curious about that, thank you so much.