Commit 5407349d authored by Stefan Kerkemeier's avatar Stefan Kerkemeier
Browse files

remove unnecessary files

parent f1a933c4
#!/bin/sh
# Copyright 1998-2019 Lawrence Livermore National Security, LLC and other
# HYPRE Project Developers. See the top-level COPYRIGHT file for details.
#
# SPDX-License-Identifier: (Apache-2.0 OR MIT)
#=============================================================================
# This script prints the hypre version number, date, and time.
# It currently inspects the 'configure' file for this info.
#=============================================================================
case $1 in
-h|-help)
echo
echo "$0 [options]"
echo " -h|-help - prints usage information"
echo " -number - prints the version number"
echo " -date - prints the version day"
echo " -time - prints the version day and time"
echo
exit;;
esac
# NOTE: In order to call this script from other directories,
# we need to get the path info from the command line
VPATH=`dirname $0`
VFILE="${VPATH}/../configure"
NUMBER=`grep "HYPRE_VERSION=" $VFILE | cut -d= -f 2 | sed 's/"//g'`
DATE=`grep "HYPRE_DATE=" $VFILE | cut -d= -f 2 | sed 's/"//g'`
TIME=`grep "HYPRE_TIME=" $VFILE | cut -d= -f 2 | sed 's/"//g'`
# this is the no-option print line
VPRINT=`echo hypre Version $NUMBER Date: $DATE`
# this defines the print lines for the various options
case $1 in
-number)
VPRINT=$NUMBER;;
-date)
VPRINT=$DATE;;
-time)
VPRINT=$TIME;;
esac
# print the version information
echo $VPRINT
examples @ 35bac752
Subproject commit 35bac75238bf3e7abb6f621615be1f5b3b2bed04
# Ignore everything in this directory?
*
# Except this file
!.gitignore
Steady isentropic vortex solution to Euler equations. Calorically perfect gas.
Widely used for simple validation in the 1990s. See
Shu (1998) MRUGESH WHAT REFERENCE IS THIS
Yee, Vinokur, Djomehri (1999) "Low-dissipative high-order shock-capturing methods
using characteristic-based filters"
J. Comp. Phys. 150, pp. 199-238
CMT-Nek_2d_isentropic_vortex.pdf is an old document summarizing our test efforts with the
OLD, INCORRECTLY-DEALIASED version of the code. The case represented here in pvort.usr
and pvort.rea corresponds to the lower-left-hand corner of Fig. 6, p. 12. It is a periodic
2D square of 4-by-4 elements with curved edges between them.
The variable "bvort" is beta in CMT-Nek_2d_isentropic_vortex.pdf. It represents vortex strength.
Maximum Mach number, referenced somewhat clumsily to the freestream, is
exp(0.5)*beta/(2*sqrt(2)*pi)/c_sound.
It is set to 5 in useric.
There is a base flow that advects the vortex through the periodic domain. Userchk dumps
an error measurment as a function of time, with the following line appearing in the logfile
'L2norms_sqrt ',time,err(rho),err(rho*u),err(rho*v),err(rho*w),err(rho*e)
TO DO
-stabilize the code correctly. filtering does not ever force the error to drop or plateau
(neither did our old dealiasing effort; see Fig. 10 of CMT-Nek_2d_isentropic_vortex.pdf).
-correct the error measurement. It is not the quotient of two L2 norms, and it is not meaningful for
rho*v.
L2err_norms 0.10000000E+01 0.16282454E-06 0.13085757E-06 0.13524302E-06 0.00000000E+00 0.52143596E-06
L2err_norms 0.20000000E+01 0.23907328E-06 0.18972729E-06 0.20414780E-06 0.00000000E+00 0.46727824E-06
L2err_norms 0.30000000E+01 0.23173725E-06 0.29259952E-06 0.20754394E-06 0.00000000E+00 0.44411979E-06
L2err_norms 0.40000000E+01 0.29042242E-06 0.38543722E-06 0.28891266E-06 0.00000000E+00 0.71206093E-06
L2err_norms 0.50000000E+01 0.28360368E-06 0.36137926E-06 0.32522802E-06 0.00000000E+00 0.57468093E-06
L2err_norms 0.60000000E+01 0.30810915E-06 0.44263109E-06 0.40753937E-06 0.00000000E+00 0.85348250E-06
L2err_norms 0.70000000E+01 0.32961038E-06 0.51922627E-06 0.42701081E-06 0.00000000E+00 0.87381191E-06
L2err_norms 0.80000000E+01 0.30444863E-06 0.49995765E-06 0.46721256E-06 0.00000000E+00 0.77739276E-06
L2err_norms 0.90000000E+01 0.35646304E-06 0.55411221E-06 0.49129136E-06 0.00000000E+00 0.10338507E-05
L2err_norms 0.10000000E+02 0.33708978E-06 0.55441730E-06 0.50826328E-06 0.00000000E+00 0.87627426E-06
****** PARAMETERS *****
2.60999990 NEKTON VERSION
2 DIMENSIONAL RUN
118 PARAMETERS FOLLOW
1.00000 p001 DENSITY
-100.000 p002 VISCOS
0.00000 p003 : BETAG
0.00000 p004 : GTHETA
0.00000 p005 : PGRADX
0.00000 p006
1.00000 p007 RHOCP
1.00000 p008 CONDUCT
0.00000 p009
0.00000 p010 FINTIME
10000.0 p011 NSTEPS
-1.000000E-03 p012 DT
0.00000 p013 IOCOMM
0.00000 p014 IOTIME
1000.00 p015 IOSTEP
0.00000 p016 PSSOLVER: 0=default
1.00000 p017
0.500000E-01 p018 GRID < 0 --> # cells on screen
-1.00000 p019 INTYPE
7.00000 p020 NORDER
0.100000E-06 p021 DIVERGENCE
0.100000E-08 p022 HELMHOLTZ
0 p23 NPSCAL
0.100000E-09 p024 TOLREL
0.100000E-04 p025 TOLABS
0.50000 p026 COURANT/NTAU
3.00000 p027 TORDER
0.00000 p028 TORDER: mesh velocity (0: p28=p27)
0.00000 p029 = magnetic visc if > 0, = -1/Rm if < 0
0.00000 p030 > 0 ==> properties set in uservp()
0.00000 p031 NPERT: #perturbation modes
0.00000 p032 #BCs in re2 file, if > 0
0.00000 p033 :
0.00000 p034 :
0.00000 p035 :
0.00000 p036 : XMAGNET
0.00000 p037 : NGRIDS
0.00000 p038 : NORDER2
0.00000 p039 : NORDER3
0.00000 p040
0.00000 p041 1-->multiplicative SEMG
0.00000 p042 0=gmres/1=pcg
0.00000 p043 0=semg/1=schwarz
0.00000 p044 0=E-based/1=A-based prec.
0.00000 p045 Relaxation factor for DTFS
0.00000 p046 reserved
0.00000 p047 vnu: mesh matieral prop.
0.00000 p048 :
0.00000 p049 :
0.00000 p050 :
0.00000 p051
0.00000 p052 IOHIS
0.00000 p053
0.00000 p054 fixed flow rate dir: |p54|=1,2,3=x,y,z
0.00000 p055 vol.flow rate (p54>0) or Ubar (p54<0)
0.00000 p056 :
0.00000 p057 :
0.00000 p058
0.00000 p059 !=0 --> full Jac. eval. for each el.
0.00000 p060 !=0 --> init. velocity to small nonzero
0.00000 p061
0.00000 p062 >0 --> force byte_swap for output
0.00000 p063 =8 --> force 8-byte output !change to 8 for full_restart
0.00000 p064 =1 --> perturbation restart
1.00000 p065 #iofiles (eg, 0 or 64); <0 --> sep. dirs
6.00000 p066 output : <0=ascii, else binary
6.00000 p067 restart: <0=ascii, else binary
0.00000 p068 iastep: freq for avg_all (0=iostep)
0.00000 p069 : frequency of srf dump
0.00000 p070 :
0.00000 p071 :
0.00000 p072 :
0.00000 p073
0.00000 p074 verbose Helmholtz
0.00000 p075 : eps for Laurette Tucke
0.00000 p076 : critical Rayleigh Numb
0.00000 p077 : sqrt Taylor number
0.00000 p078 :
0.00000 p079 :
0.00000 p080 :
0.00000 p081 :
0.00000 p082 : p82 nxc
0.00000 p083
0.00000 p084 !=0 --> sets initial timestep if p12>0
0.00000 p085 dt ratio if p84 !=0, for timesteps>0
0.00000 p086 reserved
0.00000 p087 :
0.00000 p088 : p88 1 = no prec for E
0.00000 p089 :
0.00000 p090 :
0.00000 p091 :
0.00000 p092
0.00000 p093 Number of previous pressure solns saved
0.00000 p094 start projecting velocity after p94 step
0.00000 p095 start projecting pressure after p95 step
0.00000 p096 :
0.00000 p097 :
0.00000 p098
3.00000 p099 dealiasing: <0--> off/3--> old/4--> new
0.00000 p100
0.00000 p101 Number of additional modes to filter
1.00000 p102 Dump out divergence at each time step
0.00000 p103 weight of stabilizing filter (.01)
0.00000 p104
0.00000 p105
0.00000 p106
0.00000 p107 !=0 --> add to h2 array in hlmhotz eqn
0.00000 p108
0.00000 p109
0.00000 p110
0.00000 p111
0.00000 p112
0.00000 p113
0.00000 p114
0.00000 p115
0.00000 p116 !=0: x elements for fast tensor product
0.00000 p117 !=0: y elements for fast tensor product
0.00000 p118 !=0: z elements for fast tensor product
4 Lines of passive scalar data follows2 CONDUCT; 2RHOCP
1.00000 1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000
13 LOGICAL SWITCHES FOLLOW
T IFFLOW
F IFHEAT
T IFTRAN
T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)
F F T T T T T T T T T T IFTMSH (IF mesh for this field is T mesh)
F IFAXIS
F IFSTRS
T IFSPLIT
F IFMGRID
F IFMODEL
F IFKEPS
F IFMVBD
F IFCHAR
20.0000 20.0000 -3.00000 -10.0000 XFAC,YFAC,XZERO,YZERO
**MESH DATA** 1st line is X of corner 1,2,3,4. 2nd line is Y.
16 2 16 NEL,NDIM,NELV
ELEMENT 1 [ 1a] GROUP 0
0.00000 2.50000 2.50000 0.00000
-5.00000 -5.00000 -2.50000 -2.50000
ELEMENT 2 [ 1b] GROUP 0
2.50000 5.00000 5.00000 2.50000
-5.00000 -5.00000 -2.50000 -2.50000
ELEMENT 3 [ 1c] GROUP 0
5.00000 7.50000 7.50000 5.00000
-5.00000 -5.00000 -2.50000 -2.50000
ELEMENT 4 [ 1d] GROUP 0
7.50000 10.0000 10.0000 7.50000
-5.00000 -5.00000 -2.50000 -2.50000
ELEMENT 5 [ 1e] GROUP 0
0.00000 2.50000 2.50000 0.00000
-2.50000 -2.50000 0.00000 0.00000
ELEMENT 6 [ 1f] GROUP 0
2.50000 5.00000 5.00000 2.50000
-2.50000 -2.50000 0.00000 0.00000
ELEMENT 7 [ 1g] GROUP 0
5.00000 7.50000 7.50000 5.00000
-2.50000 -2.50000 0.00000 0.00000
ELEMENT 8 [ 1h] GROUP 0
7.50000 10.0000 10.0000 7.50000
-2.50000 -2.50000 0.00000 0.00000
ELEMENT 9 [ 1i] GROUP 0
0.00000 2.50000 2.50000 0.00000
0.00000 0.00000 2.50000 2.50000
ELEMENT 10 [ 1j] GROUP 0
2.50000 5.00000 5.00000 2.50000
0.00000 0.00000 2.50000 2.50000
ELEMENT 11 [ 1k] GROUP 0
5.00000 7.50000 7.50000 5.00000
0.00000 0.00000 2.50000 2.50000
ELEMENT 12 [ 1l] GROUP 0
7.50000 10.0000 10.0000 7.50000
0.00000 0.00000 2.50000 2.50000
ELEMENT 13 [ 1m] GROUP 0
0.00000 2.50000 2.50000 0.00000
2.50000 2.50000 5.00000 5.00000
ELEMENT 14 [ 1n] GROUP 0
2.50000 5.00000 5.00000 2.50000
2.50000 2.50000 5.00000 5.00000
ELEMENT 15 [ 1o] GROUP 0
5.00000 7.50000 7.50000 5.00000
2.50000 2.50000 5.00000 5.00000
ELEMENT 16 [ 1p] GROUP 0
7.50000 10.0000 10.0000 7.50000
2.50000 2.50000 5.00000 5.00000
***** CURVED SIDE DATA *****
0 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE
***** BOUNDARY CONDITIONS *****
***** FLUID BOUNDARY CONDITIONS *****
P 1 1 13.0000 3.00000 0.00000 0.00000 0.00000
E 1 2 2.00000 4.00000 0.00000 0.00000 0.00000
E 1 3 5.00000 1.00000 0.00000 0.00000 0.00000
P 1 4 4.00000 2.00000 0.00000 0.00000 0.00000
P 2 1 14.0000 3.00000 0.00000 0.00000 0.00000
E 2 2 3.00000 4.00000 0.00000 0.00000 0.00000
E 2 3 6.00000 1.00000 0.00000 0.00000 0.00000
E 2 4 1.00000 2.00000 0.00000 0.00000 0.00000
P 3 1 15.0000 3.00000 0.00000 0.00000 0.00000
E 3 2 4.00000 4.00000 0.00000 0.00000 0.00000
E 3 3 7.00000 1.00000 0.00000 0.00000 0.00000
E 3 4 2.00000 2.00000 0.00000 0.00000 0.00000
P 4 1 16.0000 3.00000 0.00000 0.00000 0.00000
P 4 2 1.00000 4.00000 0.00000 0.00000 0.00000
E 4 3 8.00000 1.00000 0.00000 0.00000 0.00000
E 4 4 3.00000 2.00000 0.00000 0.00000 0.00000
E 5 1 1.00000 3.00000 0.00000 0.00000 0.00000
E 5 2 6.00000 4.00000 0.00000 0.00000 0.00000
E 5 3 9.00000 1.00000 0.00000 0.00000 0.00000
P 5 4 8.00000 2.00000 0.00000 0.00000 0.00000
E 6 1 2.00000 3.00000 0.00000 0.00000 0.00000
E 6 2 7.00000 4.00000 0.00000 0.00000 0.00000
E 6 3 10.0000 1.00000 0.00000 0.00000 0.00000
E 6 4 5.00000 2.00000 0.00000 0.00000 0.00000
E 7 1 3.00000 3.00000 0.00000 0.00000 0.00000
E 7 2 8.00000 4.00000 0.00000 0.00000 0.00000
E 7 3 11.0000 1.00000 0.00000 0.00000 0.00000
E 7 4 6.00000 2.00000 0.00000 0.00000 0.00000
E 8 1 4.00000 3.00000 0.00000 0.00000 0.00000
P 8 2 5.00000 4.00000 0.00000 0.00000 0.00000
E 8 3 12.0000 1.00000 0.00000 0.00000 0.00000
E 8 4 7.00000 2.00000 0.00000 0.00000 0.00000
E 9 1 5.00000 3.00000 0.00000 0.00000 0.00000
E 9 2 10.0000 4.00000 0.00000 0.00000 0.00000
E 9 3 13.0000 1.00000 0.00000 0.00000 0.00000
P 9 4 12.0000 2.00000 0.00000 0.00000 0.00000
E 10 1 6.00000 3.00000 0.00000 0.00000 0.00000
E 10 2 11.0000 4.00000 0.00000 0.00000 0.00000
E 10 3 14.0000 1.00000 0.00000 0.00000 0.00000
E 10 4 9.00000 2.00000 0.00000 0.00000 0.00000
E 11 1 7.00000 3.00000 0.00000 0.00000 0.00000
E 11 2 12.0000 4.00000 0.00000 0.00000 0.00000
E 11 3 15.0000 1.00000 0.00000 0.00000 0.00000
E 11 4 10.0000 2.00000 0.00000 0.00000 0.00000
E 12 1 8.00000 3.00000 0.00000 0.00000 0.00000
P 12 2 9.00000 4.00000 0.00000 0.00000 0.00000
E 12 3 16.0000 1.00000 0.00000 0.00000 0.00000
E 12 4 11.0000 2.00000 0.00000 0.00000 0.00000
E 13 1 9.00000 3.00000 0.00000 0.00000 0.00000
E 13 2 14.0000 4.00000 0.00000 0.00000 0.00000
P 13 3 1.00000 1.00000 0.00000 0.00000 0.00000
P 13 4 16.0000 2.00000 0.00000 0.00000 0.00000
E 14 1 10.0000 3.00000 0.00000 0.00000 0.00000
E 14 2 15.0000 4.00000 0.00000 0.00000 0.00000
P 14 3 2.00000 1.00000 0.00000 0.00000 0.00000
E 14 4 13.0000 2.00000 0.00000 0.00000 0.00000
E 15 1 11.0000 3.00000 0.00000 0.00000 0.00000
E 15 2 16.0000 4.00000 0.00000 0.00000 0.00000
P 15 3 3.00000 1.00000 0.00000 0.00000 0.00000
E 15 4 14.0000 2.00000 0.00000 0.00000 0.00000
E 16 1 12.0000 3.00000 0.00000 0.00000 0.00000
P 16 2 13.0000 4.00000 0.00000 0.00000 0.00000
P 16 3 4.00000 1.00000 0.00000 0.00000 0.00000
E 16 4 15.0000 2.00000 0.00000 0.00000 0.00000
***** NO THERMAL BOUNDARY CONDITIONS *****
0 PRESOLVE/RESTART OPTIONS *****
7 INITIAL CONDITIONS *****
C Default
C Default
C Default
C Default
C Default
C Default
C Default
***** DRIVE FORCE DATA ***** BODY FORCE, FLOW, Q
4 Lines of Drive force data follow
C
C
C
C
***** Variable Property Data ***** Overrrides Parameter data.
1 Lines follow.
0 PACKETS OF DATA FOLLOW
***** HISTORY AND INTEGRAL DATA *****
0 POINTS. Hcode, I,J,H,IEL
***** OUTPUT FIELD SPECIFICATION *****
6 SPECIFICATIONS FOLLOW
F COORDINATES
T VELOCITY
T PRESSURE
T TEMPERATURE
F TEMPERATURE GRADIENT
0 PASSIVE SCALARS
***** OBJECT SPECIFICATION *****
0 Surface Objects
0 Volume Objects
0 Edge Objects
0 Point Objects
subroutine userchk
include 'SIZE'
include 'TOTAL'
ifxyo = .true.
if (istep.gt.1) ifxyo=.false.
return
end
!-----------------------------------------------------------------------
subroutine cmtchk
include 'SIZE'
include 'TOTAL'
include 'CMTDATA'
scale = 1.
nxyz=nx1*ny1*nz1
if (istep.gt.0)then
if (mod(istep,1000).eq.0)
> call l2norm_2d_vortex
endif
! sigh. need settime to bend to my will
time=time_cmt
return
end
!-----------------------------------------------------------------------
subroutine cmt_userEOS(ix,iy,iz,eg)
include 'SIZE'
include 'NEKUSE'
include 'PARALLEL'
include 'CMTDATA'
include 'PERFECTGAS'
integer e,eg
e=gllel(eg)
cp=cpgref
cv=cvgref
temp=e_internal/cv
asnd=MixtPerf_C_GRT(gmaref,rgasref,temp)
pres=MixtPerf_P_DRT(rho,rgasref,temp)
return
end
c-----------------------------------------------------------------------
subroutine uservp (ix,iy,iz,eg)
include 'SIZE'
include 'TOTAL'
include 'NEKUSE'
include 'CMTDATA'
integer e,f,eg
c e = gllel(eg)
mu=0.0
lambda=-2.0/3.0*mu
udiff =0.
utrans=0.
nu_s=0.0
return
end
c-----------------------------------------------------------------------
subroutine userf (ix,iy,iz,eg)
include 'SIZE'
include 'TOTAL'
include 'NEKUSE'
integer e,f,eg
c e = gllel(eg)
c Note: this is an acceleration term, NOT a force!
c Thus, ffx will subsequently be multiplied by rho(x,t).
ffx = 0.0
ffy = 0.0
ffz = 0.0
return
end
c-----------------------------------------------------------------------
subroutine userq (ix,iy,iz,eg)
include 'SIZE'
include 'TOTAL'
include 'NEKUSE'
integer e,f,eg
c e = gllel(eg)
qvol = 0.0
return
end
c-----------------------------------------------------------------------
subroutine userbc (ix,iy,iz,iside,ieg)
include 'SIZE'
include 'NEKUSE'
include 'CMTDATA'
integer f,e,eg
flux=0.0
molarmass=molmass
return
end
c-----------------------------------------------------------------------
subroutine useric (ix,iy,iz,ieg)
include 'SIZE'
include 'TOTAL'
include 'NEKUSE'
include 'CMTDATA'
include 'PERFECTGAS'
integer e,eg, eqnum
real r2
real bvort,uxinf,uyinf,tempinf
common /VORTXSTRN/ bvort,uxinf,uyinf,tempinf
xctr = 5.0!22.50
yctr = 0.0!2.5
x=x-xctr
y=y-yctr
r2=(x**2+y**2)
bvort = 5.0
pi = 4.0*atan(1.0)
phi= 1.0
uxinf = 1.0
uyinf = 0.0
Tempinf = 1.0
rho = (Tempinf-(gmaref-1.0)*bvort**2/(16.0*gmaref
> *pi**2)* exp(2.0*(1.0-r2)))**(1.0/(gmaref-1.0))
ux = uxinf-bvort*exp((1.0-r2))*y/(2.0*pi)
uy = uyinf+bvort*exp((1.0-r2))*x/(2.0*pi)
uz = 0.0
temp=Tempinf-(gmaref-1.0)*bvort**2/(16.0*gmaref*pi**2)*
> exp(2.0*(1.0-r2))
pres=MixtPerf_P_DRT(rho,rgasref,temp)
cp=cpgref
cv=cvgref
e_internal=cvgref*temp
return
end
c-----------------------------------------------------------------------
subroutine usrdat
include 'SIZE'
include 'TOTAL'
include 'CMTBCDATA'
include 'CMTDATA'
include 'CMTTIMERS'
include 'PERFECTGAS'
molmass = 8314.3
muref = 0.0
coeflambda = -2.0/3.0
suthcoef = 1.0
reftemp = 1.0
prlam = 0.72
pinfty = 1.0
gmaref = 1.4
rgasref = MixtPerf_R_M(molmass,dum)
cvgref = rgasref/(gmaref-1.0)
cpgref = MixtPerf_Cp_CvR(cvgref,rgasref)
gmaref = MixtPerf_G_CpR(cpgref,rgasref)
res_freq = 5000000
flio_freq = 10000
nbc=0
do i=1,2 ! get density output
call add_temp(f2tbc,nbc,1)
enddo
igeom=2
call setup_topo
return
end
c-----------------------------------------------------------------------
subroutine usrdat2
include 'SIZE'
include 'TOTAL'
include 'CMTBCDATA'
include 'CMTDATA'
call fix_geom
outflsub=.true.
IFRESTART=.false.
IFCNTFILT=.false.
! Avoid overstepping bx bounds in nekuic with ldimt>3, lbelt=1
ifldmhd=9999
return
end
c-----------------------------------------------------------------------
subroutine usrdat3
include 'SIZE'
include 'TOTAL'
c
return
end
!-----------------------------------------------------------------------
subroutine cmt_usrflt(rmult)
include 'SIZE'
real rmult(nx1)
real alpfilt
integer sfilt, kut
real eta, etac
call rone(rmult,nx1)
alpfilt=36.0 ! H&W 5.3
kut=max((nx1/2),nx1-2)
kut=nx1/2
sfilt=8
etac=real(kut)/real(nx1)
do i=kut,nx1
eta=real(i)/real(nx1)
rmult(i)=exp(-alpfilt*((eta-etac)/(1.0-etac))**sfilt)
enddo
return