1C Copyright (c) 2017, Lawrence Livermore National Security, LLC. Produced at 2C the Lawrence Livermore National Laboratory. LLNL-CODE-734707. All Rights 3C reserved. See files LICENSE and NOTICE for details. 4C 5C This file is part of CEED, a collection of benchmarks, miniapps, software 6C libraries and APIs for efficient high-order finite element and spectral 7C element discretizations for exascale applications. For more information and 8C source code availability see http://github.com/ceed 9C 10C The CEED research is supported by the Exascale Computing Project (17-SC-20-SC) 11C a collaborative effort of two U.S. Department of Energy organizations (Office 12C of Science and the National Nuclear Security Administration) responsible for 13C the planning and preparation of a capable exascale ecosystem, including 14C software, applications, hardware, advanced system engineering and early 15C testbed platforms, in support of the nation's exascale computing imperative. 16 17C Dimension file to be included 18C 19C HCUBE array dimensions 20C 21 parameter (lx1=8) 22 parameter (lelg=10000, lelt=1024) 23 24 parameter (ldim=3) 25 parameter (lp=32) 26 27 parameter (ly1=lx1,lz1=lx1,lelv=lelt) 28 parameter (lxd=lx1+2,lyd=lxd,lzd=lxd) 29 parameter (lelx=1,lely=1,lelz=1) 30c 31 parameter (lzl=3 + 2*(ldim-3)) 32c 33 parameter (lx2=lx1) 34 parameter (ly2=ly1) 35 parameter (lz2=lz1) 36 parameter (lx3=lx2) 37 parameter (ly3=ly2) 38 parameter (lz3=lz2) 39c 40c parameter (lpelv=lelv,lpelt=lelt,lpert=3) ! perturbation 41c parameter (lpx1=lx1,lpy1=ly1,lpz1=lz1) ! array sizes 42c parameter (lpx2=lx2,lpy2=ly2,lpz2=lz2) 43c 44 parameter (lpelv=1,lpelt=1,lpert=1) ! perturbation 45 parameter (lpx1=1,lpy1=1,lpz1=1) ! array sizes 46 parameter (lpx2=1,lpy2=1,lpz2=1) 47c 48c 49c parameter (lbelv=lelv,lbelt=lelt) ! MHD 50c parameter (lbx1=lx1,lby1=ly1,lbz1=lz1) ! array sizes 51c parameter (lbx2=lx2,lby2=ly2,lbz2=lz2) 52c 53 parameter (lbelv=1,lbelt=1) ! MHD 54 parameter (lbx1=1,lby1=1,lbz1=1) ! array sizes 55 parameter (lbx2=1,lby2=1,lbz2=1) 56c 57C LX1M=LX1 when there are moving meshes; =1 otherwise 58 parameter (lx1m=1,ly1m=1,lz1m=1) 59 parameter (ldimt= 2) ! 3 passive scalars + T 60 parameter (ldimt1=ldimt+1) 61 parameter (ldimt3=ldimt+3) 62c 63c Note: In the new code, LELGEC should be about sqrt(LELG) 64c 65 PARAMETER (LELGEC = 1) 66 PARAMETER (LXYZ2 = 1) 67 PARAMETER (LXZ21 = 1) 68c 69 PARAMETER (LMAXV=LX1*LY1*LZ1*LELV) 70 PARAMETER (LMAXT=LX1*LY1*LZ1*LELT) 71 PARAMETER (LMAXP=LX2*LY2*LZ2*LELV) 72 PARAMETER (LXZ=LX1*LZ1) 73 PARAMETER (LORDER=4) 74 PARAMETER (MAXOBJ=4,MAXMBR=LELT*6,lhis=100) 75C 76C Common Block Dimensions 77C 78 PARAMETER (LCTMP0 =2*LX1*LY1*LZ1*LELT) 79 PARAMETER (LCTMP1 =4*LX1*LY1*LZ1*LELT) 80C 81C The parameter LVEC controls whether an additional 42 field arrays 82C are required for Steady State Solutions. If you are not using 83C Steady State, it is recommended that LVEC=1. 84C 85 PARAMETER (LVEC=1) 86C 87C Uzawa projection array dimensions 88C 89 parameter (mxprev = 02) 90 parameter (lgmres = 02) 91C 92C Split projection array dimensions 93C 94 parameter(lmvec = 1) 95 parameter(lsvec = 1) 96 parameter(lstore=lmvec*lsvec) 97c 98c NONCONFORMING STUFF 99c 100 parameter (maxmor = lelt) 101C 102C Array dimensions 103C 104 COMMON/DIMN/NELV,NELT,NX1,NY1,NZ1,NX2,NY2,NZ2 105 $,NX3,NY3,NZ3,NDIM,NFIELD,NPERT,NID 106 $,NXD,NYD,NZD 107