Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fix various typos #17

Open
wants to merge 2 commits into
base: master
Choose a base branch
from
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 2 additions & 2 deletions CODE/bc_p.f90
Original file line number Diff line number Diff line change
Expand Up @@ -473,7 +473,7 @@ SUBROUTINE APPLY_BOUNDARY(N,XPER,YPER,ZPER,IPERIODICITY,XMPIELRANK)

!$OMP DO SCHEDULE (STATIC)
DO I=1,KMAXE ! For ALL ELEMENTS
if (ielem(n,i)%interior.eq.1)then ! THAT HAVE AT LEAST ONE UNKNWON NEIGHBOUR
if (ielem(n,i)%interior.eq.1)then ! THAT HAVE AT LEAST ONE UNKNOWN NEIGHBOUR
IF (ielem(n,i)%nofbc.GT.0)THEN ! THAT HAVE AT LEAST ESTABLISHED A BOUNDARY CONDITION CODE

DO J=1,ielem(n,i)%IFCA ! LOOP ALL THEIR FACES
Expand Down Expand Up @@ -574,7 +574,7 @@ SUBROUTINE APPLY_BOUNDARY(N,XPER,YPER,ZPER,IPERIODICITY,XMPIELRANK)
jj1=0
!$OMP DO SCHEDULE (STATIC)
DO I=1,KMAXE !> ALL ELEMENTS
if (ielem(n,i)%interior.eq.1)then ! THAT HAVE AT LEAST ONE UNKNWON NEIGHBOUR
if (ielem(n,i)%interior.eq.1)then ! THAT HAVE AT LEAST ONE UNKNOWN NEIGHBOUR
IF (ielem(n,i)%nofbc.GT.0)THEN ! THAT HAVE AT LEAST ESTABLISHED A BOUNDARY CONDITION CODE
DO J=1,ielem(n,i)%IFCA ! LOOP ALL THEIR BOUNDARY FACES
if (IELEM(N,I)%IBOUNDS(J).gt.0)then
Expand Down
20 changes: 10 additions & 10 deletions CODE/declarations.f90
Original file line number Diff line number Diff line change
Expand Up @@ -300,17 +300,17 @@ MODULE DECLARATION
!$OMP THREADPRIVATE(sb,IBFC,number_of_dog,WEIGHT_T2,lscqm1,B1_imp,DU1,DU2,DUMMY12,xxder,yyder,zzder,C1_imp,DUR,DUL,DURR,DULR,DUT1,B1T,DUMMY12T,SOURCE_T,UTMIN,UTMAX,ax,ay,az,nx,ny,nz,nnx,nny,nnz,VVA,DETA,VVA1,VVr1,VVr2,VVr3,VVR4,VVwg,VVnpox,VVnpoy,VVnpoz,VVwpox,VVwpoy,VVwpoz,VVnxi,VVneta,VVnzeta,VVxi,VVeta,VVzeta,VVnallx,VVnally,VVnallz,VVB,VVC,VVD,VVE,VVF,vvjacobsurf,vvjacobvolume)

REAL,ALLOCATABLE,DIMENSION(:)::gamma_IN,MP_IE,MP_A,MP_AR,MP_R,MP_R_IN,MP_A_IN,MP_PINF !MULTIPHASE COMPONENTS
REAL,ALLOCATABLE,DIMENSION(:)::CLEFT,lamc !VALUE OF THE RECONSTRUCTED SOLUTION REFERING TO MY CELL (I) CONSERVED VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::CLEFT_ROT !VALUE OF THE RECONSTRUCTED SOLUTION REFERING TO MY CELL (I) ROTATED WITH RESPECT TO ITS FACE CONSERVED VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::CRIGHT !VALUE OF THE RECONSTRUCTED SOLUTION REFERING TO THE NEIGHBOUR OF MY CELL (I) CONSERVED VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::CRIGHT_ROT !VALUE OF THE RECONSTRUCTED SOLUTION REFERING TO THE NEIGHBOUR OF MY CELL (I) ROTATED WITH RESPECT TO ITS FACE CONSERVED VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::CLEFT,lamc !VALUE OF THE RECONSTRUCTED SOLUTION REFERRING TO MY CELL (I) CONSERVED VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::CLEFT_ROT !VALUE OF THE RECONSTRUCTED SOLUTION REFERRING TO MY CELL (I) ROTATED WITH RESPECT TO ITS FACE CONSERVED VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::CRIGHT !VALUE OF THE RECONSTRUCTED SOLUTION REFERRING TO THE NEIGHBOUR OF MY CELL (I) CONSERVED VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::CRIGHT_ROT !VALUE OF THE RECONSTRUCTED SOLUTION REFERRING TO THE NEIGHBOUR OF MY CELL (I) ROTATED WITH RESPECT TO ITS FACE CONSERVED VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::SL !WAVE SPEED ESTIMATE OF LEFT STATE OF THE SOLUTION OF A FACE
REAL,ALLOCATABLE,DIMENSION(:)::SR !WAVE SPEED ESTIMATE OF RIGHT STATE OF THE SOLUTION OF A FACE
REAL,ALLOCATABLE,DIMENSION(:)::SM !WAVE SPEED ESTIMATE OF STAR REGION ACCORDING TO RIEMANN SOLVER FAN
REAL,ALLOCATABLE,DIMENSION(:)::CTURBL !VALUE OF THE RECONSTRUCTED SOLUTION OF ADDITIONAL EQUATIONS REFERING TO MY CELL (I) PRIMITIVE VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::CTURBR !VALUE OF THE RECONSTRUCTED SOLUTION OF ADDITIONAL EQUATIONS REFERING TO THE NEIGHBOUR OF MY CELL (I) PRIMITIVE VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::ROTVL !VALUE OF THE RECONSTRUCTED SOLUTION REFERING TO MY CELL (I) ROTATED WITH RESPECT TO ITS FACE PRIMITIVE VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::ROTVR !VALUE OF THE RECONSTRUCTED SOLUTION REFERING TO THE NEIGHBOUR OF MY CELL (I) ROTATED WITH RESPECT TO ITS FACE PRIMITIVE VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::CTURBL !VALUE OF THE RECONSTRUCTED SOLUTION OF ADDITIONAL EQUATIONS REFERRING TO MY CELL (I) PRIMITIVE VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::CTURBR !VALUE OF THE RECONSTRUCTED SOLUTION OF ADDITIONAL EQUATIONS REFERRING TO THE NEIGHBOUR OF MY CELL (I) PRIMITIVE VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::ROTVL !VALUE OF THE RECONSTRUCTED SOLUTION REFERRING TO MY CELL (I) ROTATED WITH RESPECT TO ITS FACE PRIMITIVE VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::ROTVR !VALUE OF THE RECONSTRUCTED SOLUTION REFERRING TO THE NEIGHBOUR OF MY CELL (I) ROTATED WITH RESPECT TO ITS FACE PRIMITIVE VARIABLE
REAL,ALLOCATABLE,DIMENSION(:)::HLLCFLUX !VALUE OF THE RIEMANN SOLUTION NORMAL TO FACE DIRECTION
REAL,ALLOCATABLE,DIMENSION(:)::RHLLCFLUX !VALUE OF THE RIEMANN SOLUTION AT CARTESIAN COORDINATES
REAL::L1NORM !L1 NORM OF SOLUTION FOR GRID CONVERGENCE STUDIES OF EULER AND LINEAR ADVECTION EQUATIONS
Expand Down Expand Up @@ -341,7 +341,7 @@ MODULE DECLARATION
REAL,ALLOCATABLE,DIMENSION(:,:)::IMPOFF_MF !IMPLICIT ONLY OFF DIAGONAL MATRIX D
REAL,ALLOCATABLE,DIMENSION(:,:,:)::IMPOFFT !IMPLICIT ONLY OFF DIAGONAL MATRIX D FOR TURBULENCE
REAL,ALLOCATABLE,DIMENSION(:,:)::IMPDIAGT !IMPLICIT ONLY DIAGONAL MATRIX D FOR TURBULENCE
REAL,ALLOCATABLE,DIMENSION(:,:)::SHT !IMPLICIT ONLY SORUCE TERM JACOBIAN
REAL,ALLOCATABLE,DIMENSION(:,:)::SHT !IMPLICIT ONLY SOURCE TERM JACOBIAN
REAL,ALLOCATABLE,DIMENSION(:,:)::INVTRI !INVERSE MATRIX FOR TRANSFORMATION
REAL,ALLOCATABLE,DIMENSION(:)::ROTVECT !ROTATED VECTOR OF CONSERVED VARIABLES
REAL,ALLOCATABLE,DIMENSION(:)::VECTCO ! VECTOR OF CONSERVED VARIABLES
Expand Down Expand Up @@ -689,7 +689,7 @@ MODULE DECLARATION
INTEGER::full !SPECIFIES IF SUFFICIENT NUMBER OF STENCILS ARE FOUND TO PROCEED WITH WENO FOR THIS CELL
INTEGER::INTERIOR !SPECIFIES IF THIS CELL HAS ANY SIDE BOUNDED (INTERIOR CELLS GET A VALUE OF 0, NON INTERIOR ONES GET A VALUE OF 1)
INTEGER::VDEC,HALO !NUMBER OF VOLUME DECOMPOSITIONS FOR EACH ELEMENT
INTEGER::MODE !SPECIFIES IF DECOMPOSITION CANBE AVOIDED FOR STRAIGHT SIDED ELEMENT (NOT IMPLEMENTED YET)
INTEGER::MODE !SPECIFIES IF DECOMPOSITION CAN BE AVOIDED FOR STRAIGHT SIDED ELEMENT (NOT IMPLEMENTED YET)
INTEGER::GGS ! SPECIFIES WITH WHAT ALGORITHM TO COMPUTE THE GRADIENTS (GREEN GAUSS, LEAST SQUARES OR BLEND OF THEM)
INTEGER::idegfree !DEGREES OF FREEDOM FOR POLYNOMIAL SELECTED
INTEGER::inumneighbours !NUMBER OF NEIGHBOURS IN EACH STENCIL
Expand Down
2 changes: 1 addition & 1 deletion CODE/gradients.f90
Original file line number Diff line number Diff line change
Expand Up @@ -1379,7 +1379,7 @@ END SUBROUTINE COMPUTE_GRADIENTS_MEAN_LSQ
!
SUBROUTINE COMPUTE_GRADIENTS_INNER_MEAN_LSQ_VISCOUS(N,ICONSIDERED,NUMBER_OF_DOG,NUMBER_OF_NEI)!check_all
!> @brief
!> This subroutine computes the gradients of the primitve variables of each interior cell using the least-squares
!> This subroutine computes the gradients of the primitive variables of each interior cell using the least-squares
IMPLICIT NONE
INTEGER,INTENT(IN)::N,ICONSIDERED,NUMBER_OF_DOG,NUMBER_OF_NEI
REAL,DIMENSION(nof_variables)::SOLS1,SOLS2
Expand Down
2 changes: 1 addition & 1 deletion CODE/grid_p.f90
Original file line number Diff line number Diff line change
Expand Up @@ -5985,7 +5985,7 @@ END SUBROUTINE STENCIILS
! !!!!!!!!!!!!!!!!!!!!!FOR VARIOUS ORDER OF ACCURACY!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
SUBROUTINE DETERMINE_SIZE(N,IORDER,ISELEM,ISELEMT,IOVERST,IOVERTO,ILX,NUMNEIGHBOURS,IDEGFREE,IMAXDEGFREE,IEXTEND)
!> @brief
!> This subroutine determines the degress of freedom, neighbours and polynomial order for each stencil of each cell
!> This subroutine determines the degrees of freedom, neighbours and polynomial order for each stencil of each cell
IMPLICIT NONE
INTEGER,INTENT(INOUT)::IORDER
INTEGER,ALLOCATABLE,DIMENSION(:),INTENT(INOUT)::ISELEMT
Expand Down
8 changes: 4 additions & 4 deletions CODE/io.f90
Original file line number Diff line number Diff line change
Expand Up @@ -12893,7 +12893,7 @@ SUBROUTINE CHECKPOINT(N)
CALL MPI_BARRIER(MPI_COMM_WORLD,IERROR)
! ioCPt1=MPI_WTIME()

ALLOCATE(DISPT(KMAXE),ARRAY2(KMAXE*(NOF_VARIABLES+turbulenceequations+passivescalar))) !I ALLOCATE IN MEMORY THE PATTERN OF ACCESS OF DATA IN TERMS OF DISPLACEMENT, AND IN TERMS OF BLOCKLENGTH, AND FINALY AN ARRAY WITH THIS PROCESSOR DATA
ALLOCATE(DISPT(KMAXE),ARRAY2(KMAXE*(NOF_VARIABLES+turbulenceequations+passivescalar))) !I ALLOCATE IN MEMORY THE PATTERN OF ACCESS OF DATA IN TERMS OF DISPLACEMENT, AND IN TERMS OF BLOCKLENGTH, AND FINALLY AN ARRAY WITH THIS PROCESSOR DATA

DO I=1,KMAXE
DISPT(I)=(XGO(I)-1)*(NOF_VARIABLES+turbulenceequations+passivescalar)
Expand Down Expand Up @@ -13263,7 +13263,7 @@ SUBROUTINE CHECKPOINT2D(N)
ICPUID=N
CALL MPI_BARRIER(MPI_COMM_WORLD,IERROR)

ALLOCATE(DISPT(KMAXE),ARRAY2(KMAXE*(NOF_VARIABLES+turbulenceequations+passivescalar))) !I ALLOCATE IN MEMORY THE PATTERN OF ACCESS OF DATA IN TERMS OF DISPLACEMENT, AND IN TERMS OF BLOCKLENGTH, AND FINALY AN ARRAY WITH THIS PROCESSOR DATA
ALLOCATE(DISPT(KMAXE),ARRAY2(KMAXE*(NOF_VARIABLES+turbulenceequations+passivescalar))) !I ALLOCATE IN MEMORY THE PATTERN OF ACCESS OF DATA IN TERMS OF DISPLACEMENT, AND IN TERMS OF BLOCKLENGTH, AND FINALLY AN ARRAY WITH THIS PROCESSOR DATA

DO I=1,KMAXE
DISPT(I)=(XGO(I)-1)*(NOF_VARIABLES+turbulenceequations+passivescalar)
Expand Down Expand Up @@ -13978,7 +13978,7 @@ END SUBROUTINE CHECKPOINTAV2d

SUBROUTINE PROBING
!> @brief
!> This subroutine writes the primitve variables at the probe positions
!> This subroutine writes the primitive variables at the probe positions
IMPLICIT NONE
INTEGER::INV
CHARACTER(LEN=120)::PROB,PROBFILE,PROC3
Expand Down Expand Up @@ -14026,7 +14026,7 @@ END SUBROUTINE PROBING

SUBROUTINE PROBING2D
!> @brief
!> This subroutine writes the primitve variables at the probe positions in 2D
!> This subroutine writes the primitive variables at the probe positions in 2D
IMPLICIT NONE
INTEGER::INV
CHARACTER(LEN=120)::PROB,PROBFILE,PROC3
Expand Down
4 changes: 2 additions & 2 deletions CODE/local_pt.f90
Original file line number Diff line number Diff line change
Expand Up @@ -1960,7 +1960,7 @@ END SUBROUTINE LOCALISE_STEN2d

subroutine direct_side(n)
!> @brief
!> This subroutine establishes the distance betwen cell centres for each face
!> This subroutine establishes the distance between cell centres for each face
implicit none
integer,intent(in)::n
integer::i,j,k,kmaxe,facexx,ixxfff
Expand Down Expand Up @@ -2063,7 +2063,7 @@ end subroutine direct_side

subroutine direct_side2d(n)
!> @brief
!> This subroutine establishes the distance betwen cell centres for each edge
!> This subroutine establishes the distance between cell centres for each edge
implicit none
integer,intent(in)::n
integer::i,j,k,kmaxe,facexx,ixxfff
Expand Down
4 changes: 2 additions & 2 deletions CODE/main.f90
Original file line number Diff line number Diff line change
Expand Up @@ -151,13 +151,13 @@ PROGRAM UCNS3D

CALL READ_INPUT(N,XMPIELRANK,XMPINRANK,XMPIE,XMPIN,IELEM,INODE,IMAXN,IMAXE,IBOUND,IMAXB,XMPINNUMBER,SCALER,INODER) !> Read the grid files and populate the allocated memory values for vertex coordinates and numbering

CALL DETERMINE_SIZE(N,IORDER,ISELEM,ISELEMT,IOVERST,IOVERTO,ILX,NUMNEIGHBOURS,IDEGFREE,IMAXDEGFREE,IEXTEND) !> Determing the stencil sizes, number of polynomial coefficients etc.
CALL DETERMINE_SIZE(N,IORDER,ISELEM,ISELEMT,IOVERST,IOVERTO,ILX,NUMNEIGHBOURS,IDEGFREE,IMAXDEGFREE,IEXTEND) !> Determining the stencil sizes, number of polynomial coefficients etc.

!$OMP PARALLEL DEFAULT(SHARED)
IF ((TURBULENCE.EQ.1).OR.(PASSIVESCALAR.GT.0))THEN
CALL ALLOCATETURB(N,EDDYFL,EDDYFR) !> Allocate memory for turbulence model or passive scalar variables
END IF
CALL GAUSSIANPOINTS(IGQRULES,NUMBEROFPOINTS,NUMBEROFPOINTS2) !> Establish the number of Gausian quadrature points for each element and face
CALL GAUSSIANPOINTS(IGQRULES,NUMBEROFPOINTS,NUMBEROFPOINTS2) !> Establish the number of Gaussian quadrature points for each element and face
CALL VERTALLOCATION(N,VEXT,LEFTV,RIGHTV,VISCL,LAML) !> Allocate the memory for the left and right state variables
CALL ALLOCATE1(N) !> Additional allocations
call QUADALLOC(QPOINTS,QPOINTS2D,WEQUA2D,WEQUA3D,NUMBEROFPOINTS,NUMBEROFPOINTS2) !>Allocate the memory required for the Gaussian integration rules
Expand Down
2 changes: 1 addition & 1 deletion CODE/parameters.f90
Original file line number Diff line number Diff line change
Expand Up @@ -159,7 +159,7 @@ SUBROUTINE READ_UCNS3D
CB2=0.622! 0.622 || Cb2
SIGMA=0.666666667! 0.666666667 || sigma
KAPPA=0.41 ! 0.41 || kappa
CW1=3.23886781677! 3.23886781677 || CW1 (0 for (CB1/kappa**2)+1/sigma(1+CB2) or predifined)
CW1=3.23886781677! 3.23886781677 || CW1 (0 for (CB1/kappa**2)+1/sigma(1+CB2) or predefined)
CW2=0.3! 0.3 || CW2
CW3=2.0 ! 2.0 || CW3
CV1=7.1! 7.1 || CV1
Expand Down
2 changes: 1 addition & 1 deletion CODE/reconstruct_p.f90
Original file line number Diff line number Diff line change
Expand Up @@ -3639,7 +3639,7 @@ END SUBROUTINE WENOWEIGHTS2d
! !---------------------------------------------------------------------------------------------!
! !---------------------------------------------------------------------------------------------!
! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! !!!!!!!!!!!!!!!!!!SUBROUTINE CALLED TO EMPLOY THE LEAST SQARES LINEAR INTERPOLATION!!!!!!!!!!!!
! !!!!!!!!!!!!!!!!!SUBROUTINE CALLED TO EMPLOY THE LEAST SQUARES LINEAR INTERPOLATION!!!!!!!!!!!!
! !!!!!!!!!!!!!!!!!!!!!!!FOR DETERMINING THE SLOPES OF EACH CELL IN EACH DIRECTION!!!!!!!!!!!!!!!
! !!!!!!!!!!!!!!!!IN A WEIGHTED AVERAGE WAY WITH THE INVERSE DISTANCE !!!!!!!!!!!!!!!!!!!!!!!!!!!
! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Expand Down
4 changes: 2 additions & 2 deletions CODE/translate.f90
Original file line number Diff line number Diff line change
Expand Up @@ -74,7 +74,7 @@ Subroutine Drive(interray)

Integer::ing2,jj,dimen,imaxe,imaxn,imaxb,dum,index10,zoneid,in1,dum1,bctypdum,iosx,ii,vrt1,vrt2,vrt3,vrt4,vct1,vct2,vct3,vct4,countb,eltype,ing,ibtr,nfin,icte
Integer::lexist,checkbrac,dum2,dum3,dum4,ierr,il,il1,il2,ispace1,ispace2
Integer :: countline,counto,countword,lengt,intsize,imaxnglobal,imaxeglobal,imin,imax,imaxfglobal,index1,ina,in2,ichen,iix,iiy,icountfc,corn,icorn
Integer :: countline,counto,countword,length,intsize,imaxnglobal,imaxeglobal,imin,imax,imaxfglobal,index1,ina,in2,ichen,iix,iiy,icountfc,corn,icorn
integer,dimension(4)::nod,nodx,cans,cand,cane,cang,canh,iiz,canf
integer,allocatable,dimension(:),intent(inout):: interray
integer,allocatable,dimension(:)::spaces,integerarray,ishape
Expand Down Expand Up @@ -1277,7 +1277,7 @@ Subroutine Drive(interray)
!
!
! read(82,"(A1,I2)",advance='NO')braco,dum
! ! if comment preceed index = 0
! ! if comment precede index = 0
! if (dum.eq.0) then
! read(82,"(1x)",advance='YES')
! endif
Expand Down
4 changes: 2 additions & 2 deletions PARAMETERS_README.txt
Original file line number Diff line number Diff line change
Expand Up @@ -78,7 +78,7 @@ PRANDTL= PRANDTL CONSTANT

Reynolds= REYNOLDS NUMBER
POSSIBLE VALUES= The value is defined as (Re=(RRES*UFREESTREAM*CHARLENGTH)/(VISC)),
and it is used to determing the freestream value of viscosity
and it is used to determining the freestream value of viscosity

CharLength= CHARACTERISTIC LENGTH
POSSIBLE VALUES= --> ANY VALUE
Expand Down Expand Up @@ -178,7 +178,7 @@ reSLIMIT= NORMALISED RESIDUAL CONVERGENCE CRITERION FOR STEADY STATE PROBL
A WIDE RANGE OF PROBLEMS
FOR OPTION (5, 10) A VALUE CLOSE TO 0.00001 MIGHT BE REQUIRED

iboundary= PRESENCE OF PERIDIC BOUNDARY IN THE DOMAIN
iboundary= PRESENCE OF PERIODIC BOUNDARY IN THE DOMAIN
POSSIBLE VALUES= 1 --> PERIODIC
0 --> NON PERIODIC

Expand Down
2 changes: 1 addition & 1 deletion README.txt
Original file line number Diff line number Diff line change
Expand Up @@ -69,7 +69,7 @@ The ucns3d code generates numerous files while running depending on the setup ch
a)history.txt: outlines the various stages of the initialisation and indicates the timestep size, iteration,
time, and message when writing solution or checkpoint files
b)STATISTICS.txt: outlines the cpu time taken per iteration for various procedures in the code while the code is running,
including cputime taken for reconstruction, fluxes, solution update, mpi boundary values exhange, mpi hallo exchange for stencils etc. Note that the collection of statistics has a noticeable computational footpring and should only be used for a limited number of iterations rather than production runs for the entire simulation.
including cputime taken for reconstruction, fluxes, solution update, mpi boundary values exchange, mpi hallo exchange for stencils etc. Note that the collection of statistics has a noticeable computational footpring and should only be used for a limited number of iterations rather than production runs for the entire simulation.
c) FORCE.dat : specifies in four columns the (iteration, time, Lift Coefficient, Drag Coefficient). You need to ensure that
the last two are normalised with respect to the reference area of your geometry.
d) RESTART.dat: A checkpoint file written when the simulation is stopped because one of the four criteria has been met
Expand Down
2 changes: 1 addition & 1 deletion README_TESTS.txt
Original file line number Diff line number Diff line change
Expand Up @@ -5,7 +5,7 @@ https://doi.org/10.5281/zenodo.3375432

These include 2D, and 3D test problems using Linear Advection, Euler, Navier-Stokes and RANS equations.

All the directores include the grid files (MESHES), the pointwise files for generating the meshes,
All the directories include the grid files (MESHES), the pointwise files for generating the meshes,
the parameter files UCNS3D.DAT (DATS), and sample run directories with some representative outputs.


Expand Down
2 changes: 1 addition & 1 deletion SCRIPTS/ARCHER_IB.pbs
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@
# Make sure any symbolic links are resolved to absolute path
#export PBS_O_WORKDIR=$(readlink -f $PBS_O_WORKDIR)
# export PBS_O_WORKDIR=/.../
# Change to the direcotry that the job was submitted from
# Change to the directory that the job was submitted from
cd $PBS_O_WORKDIR

# Set the number of threads to 1
Expand Down
2 changes: 1 addition & 1 deletion SCRIPTS/ARCHER_KNL.pbs
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@
# Make sure any symbolic links are resolved to absolute path
#export PBS_O_WORKDIR=$(readlink -f $PBS_O_WORKDIR)
# export PBS_O_WORKDIR=/.....././././.
# Change to the direcotry that the job was submitted from
# Change to the directory that the job was submitted from
cd $PBS_O_WORKDIR

# Set the number of threads to 1
Expand Down