1 c---------------------------------------------------------------------- 2 c 3 c This file contains the common blocks and the data declaration needed 4 c for the routines. 5 c 6 c Input variables that have been previously declared in common_c.h have to be 7 c re-declared here, in a consistant block. 8 c 9 c Zdenek Johan, Winter 1991. (Fortran 90) 10 c---------------------------------------------------------------------- 11 12 IMPLICIT REAL*8 (a-h,o-z) 13 c 14 c.... parameters IF YOU CHANGE THES YOU HAVE TO CHANGE THEM IN 15 c common_c.h ALSO 16 c 17 parameter ( MAXBLK = 50000, MAXTS = 100) 18 parameter ( MAXSH = 32, NSD = 3 , NSDSQ = 9) 19 c 20 c The five types of region topology are 1= Tet, 2=Hex, 3= Wedge (tri-start), 21 c 4= Wedge (quad-first) 5=pyramid 22 c 23 c The two types of face topology are 1= tri, 2=quad 24 c 25 parameter ( MAXTOP = 6, MAXSURF=1000 ) 26 27 c the common block nomodule holds all the things which have been removed 28 c from different modules 29 30 integer seqsize, stepseq 31 integer consrv_sclr_conv_vel 32 integer spongecontinuity, spongemomentum1, spongemomentum2 33 integer spongeenergy, spongemomentum3 34 integer*8 nshgt,minowned, maxowned 35 common /amgvarr/strong_eps,ramg_eps,ramg_relax,ramg_trunc, 36 & ramg_chebyratio 37 common /amgvari/irun_amg,irun_amg_prec, 38 & iamg_verb, 39 & iamg_neg_sten,iamg_nlevel, 40 & iamg_c_solver, 41 & iamg_init, 42 & iamg_setup_frez, 43 & iamg_interp,maxnev,maxncv,iamg_smoother,mlsdeg, 44 & iamg_reduce 45 46 common /nomodule/ bcttimescale,ValueListResist(0:MAXSURF), 47 & rhovw,thicknessvw, evw, rnuvw, rshearconstantvw, betai, 48 & icardio, itvn, ipvsq, numResistSrfs, nsrflistResist(0:MAXSURF), 49 & numImpSrfs, nsrflistImp(0:MAXSURF),impfile, 50 & numRCRSrfs, nsrflistRCR(0:MAXSURF),ircrfile, 51 & ideformwall, iwallmassfactor, iwallstiffactor, iviscflux 52 common /sequence/ seqsize, stepseq(100) 53 common /workfc/ master, numpe, myrank 54 common /fronts/ maxfront, nlwork 55 common /newdim/ nshgt, minowned,maxowned, numper, nshg0 56 common /timer4/ birth, death, comtim 57 common /extrat/ ttim(100) 58 common /spongevar/ zoutSponge, radSponge, zinSponge, 59 & grthOSponge,grthISponge,betamax, 60 & spongecontinuity, spongemomentum1, spongemomentum2, 61 & spongeenergy, spongemomentum3 62 common /turbvar/ eles,ylimit(3,9), rampmdot(2,3), 63 & rmutarget, pzero, wtavei, 64 & dtavei, dke, fwr1, flump, DES_SA_hmin, 65 & ierrcalc, ihessian, itwmod, ngaussf,idim, 66 & nlist, nintf(MAXTOP) 67 common /turbvari/iRANS, iLES, idistcalc, isubmod, ifproj, 68 & i2filt, modlstats, idis, nohomog, 69 & ierrsmooth, iramp 70 common /mpistats/iISend, iISendScal, iIRecv, iIRecvScal, 71 & iWaitAll,iWaitAllScal, iAllR, iAllRScal, 72 & impistat, impistat2, rmpitmr, 73 & rISend, rISendScal, rIRecv, rIRecvScal, 74 & rWaitAll, rWaitAllScal, rAllR, rAllRScal, 75 & rCommu, rCommuScal 76 77 common /memstats/rheap,rheapavail,rstack,rstackavail,rshared, 78 & rpersist,rguard,rmmap,rmemstats 79 80 common /spebcvr/ irscale, intpres, plandist, 81 & thetag, ds, tolerence, radcyl, rbltin, rvscal 82 83 common /sclrs/ scdiff(5),tdecay,nsclr,isclr,nsolt,nosource, 84 & consrv_sclr_conv_vel 85 c 86 c.... common blocks 87 c 88 parameter (MAXQPT = 125) 89 c 90 c.... common blocks for hierarchic basis functions 91 c 92 common /intpt/ Qpt (MAXTOP ,4,MAXQPT), Qwt (MAXTOP ,MAXQPT), 93 & Qptb(MAXTOP,4,MAXQPT), Qwtb(MAXTOP,MAXQPT), 94 & nint(MAXTOP), nintb(MAXTOP), 95 & ngauss, ngaussb, intp, 96 & maxnint 97 98 c nsrflist is a binary switch that tells us if a given srfID should be 99 c included in the consistent flux calculation. It starts from zero 100 c since we need to be able to handle/ignore surfaces with no srfID attached 101 c 102 c flxID(numfluxes,nIDs+1) 103 c numfluxes = area, mass, fx, fy, fz, heat, scalar_flux_{1,2,3,4} 104 c nIDs currently set to MAXSURF, each surface has its own 105 c 106 common /aerfrc/ flxID(10,0:MAXSURF), Force(3),HFlux, 107 & nsrflist(0:MAXSURF), isrfIM, 108 & flxIDsclr(4,MAXSURF), 109 & irankfilesforce(0:MAXSURF) 110 c 111 common /astore/ a(100000) 112 c 113 common /blkdat/ lcblk(10,MAXBLK+1), lcblkb(10,MAXBLK+1) 114 c 115 common /mbndnod/ mnodeb(9,8,3) 116 c 117 integer, target :: numnp, numel, numelb, numpbc, nen, nfaces, 118 & numflx, ndof, nelblk, nelblb,ntopsh, nlwork, 119 & nedof, 120 & nshg, nnz, nflow, 121 & nfath, ncorpsize, iownnodes, usingpetsc 122 123 common /conpar/ numnp, numel, numelb, numpbc, nen, nfaces, 124 & numflx, ndof, iALE, icoord, navier, 125 & irs, iexec, necho, ichem, iRK, nedof, 126 & nshg, nnz, istop, nflow, nnz_tot, idtn, 127 & ncorpsize, iownnodes, usingpetsc 128 c......................NASA Inlet Woolwine................................. 129 common /nasavari/ BCdtKW, tsBase 130 common /nasavar/ PresFreq, PresAmp, AlphaFreq, AlphaAmp 131 c.......................................................................... 132 c........................................................................... 133 common /ctrlvari/ iI2Binlet, isetOutPres, isetInitial 134 135 common /Ductvari/ BlowingVelDuct, 136 & BlowingIniMdotDuct, 137 & BlowingFnlMdotDuct, 138 & suctionVbottom, 139 & suctionVside_lower, 140 & suctionVside_upper, 141 & suctionVtop, 142 & blowerVelocity, 143 & blowerTemperature, 144 & blowerEV, 145 & isetOutletID, 146 & isetInitial_Duct, 147 & isetInlet_Duct, 148 & isetSuctionID_Duct, 149 & isetBlowerID_Duct, 150 & iDuctgeometryType, 151 & iStraightPrint, 152 & isetEV_IC_BC, 153 & isetEVramp, 154 & isetBlowing_Duct, 155 & ifixBlowingVel_Duct, 156 & nBlowingStepsDuct 157 real*8 inletVelX 158 common /ctrlvar/ inletVelX, outPres1, 159 & xvel_ini, yvel_ini, zvel_ini, 160 & temp_ini, pres_ini, evis_ini 161 162 common /Ductvar/ evis_IC_BC, 163 & EVrampXmin, 164 & EVrampXmax, 165 & EVrampMin, 166 & EVrampMax 167 c........................................................................... 168 169 c 170 common /levlset/ epsilon_ls, epsilon_lsd, dtlset, iLSet, 171 & ivconstraint, iExpLSSclr1, iExpLSSclr2 172 173 c 174 common /shpdat/ nshape, nshapeb, maxshb, 175 & nshl, nshlb,nfath, ntopsh, nsonmax 176 c 177 common /datpnt/ mshp, mshgl, mwght, mshpb, mshglb, mwghtb, 178 & mmut, mrhot, mxst 179 c 180 common /melmcat/ mcsyst, melCat, nenCat(8,3), nfaCat(8,3) 181 c 182 common /elmpar/ lelCat, lcsyst, iorder, nenb, 183 & nelblk, nelblb, ndofl, nsymdl, nenl, nfacel, 184 & nenbl, intind, mattyp 185 c 186 187 integer EntropyPressure 188 189 common /genpar/ E3nsd, I3nsd, nsymdf, ndofBC, ndiBCB, ndBCB, 190 & Jactyp, jump, ires, iprec, iprev, ibound, 191 & idiff, lhs, itau, ipord, ipred, lstres, 192 & iepstm, dtsfct, taucfct, ibksiz, iabc, isurf, 193 & idflx, Bo, EntropyPressure, irampViscOutlet, 194 & istretchOutlet, iremoveStabTimeTerm, iLHScond 195 196 c 197 integer :: svLSType, svLSFlag 198 common /inpdat/ epstol(6), Delt(MAXTS), CFLfl(MAXTS), 199 & CFLsl(MAXTS), nstep(MAXTS), niter(MAXTS), 200 & impl(MAXTS), rhoinf(MAXTS), 201 & LHSupd(6), loctim(MAXTS), deltol(MAXTS,2), 202 & leslib, svLSFlag, svLSType 203 c 204 common /intdat/ intg(2,MAXTS), intpt(3), intptb(3) 205 c 206 common /mintpar/ indQpt(3,3,4), numQpt(3,3,4), 207 & intmax 208 c 209 common /mio / iin, igeom, ipar, ibndc, imat, iecho, 210 & iout, ichmou, irstin, irstou, ihist, iflux, 211 & ierror, itable, iforce, igraph, itime 212 c 213 c /* common /andres/ fwr1,ngaussf,idim,nlist */ 214 215 character*80 fin, fgeom, fpar, fbndc, fmat, fecho, 216 & frstin, frstou, fhist, ferror, ftable, fforce, 217 & fgraph, ftime, iotype 218 common /mioname/ fin, fgeom, fpar, fbndc, fmat, fecho, 219 & frstin, frstou, fhist, ferror, ftable, fforce, 220 & fgraph, ftime 221 c 222 common /itrpar/ eGMRES, lGMRES, lGMRESs, iKs, iKss, ntotGM, ntotGMs 223 c 224 REAL*8 Nh, Msh 225 common /mmatpar/ pr, Planck, Stefan, Nh, Rh, Rgas, 226 & gamma, gamma1, s0, const, xN2, xO2, 227 & yN2, yO2, Msh(5), cpsh(5),s0sh(5),h0sh(5), 228 & Rs(5), cps(5), cvs(5), h0s(5), Trot(5),sigs(5), 229 & Tvib(5),g0s(5), dofs(5),ithm 230 c 231 logical mexist 232 common /matdat/ datmat(3,7,MAXTS), matflg(6,MAXTS), 233 & nummat, mexist 234 235 integer input_mode, output_mode 236 common /outpar/ ro, vel, temper, press, entrop, ntout, 237 & ioform, iowflux, iofieldv, iotype, ioybar, 238 & nstepsincycle, nphasesincycle, 239 & ncycles_startphaseavg, ivort, icomputevort, 240 & nsynciofiles, nsynciofieldswriterestart, 241 & iv_rankpercore, iv_corepernode, 242 & input_mode, output_mode 243 244 common /point / mbeg, mend, mprec 245 c 246 common /precis/ epsM, iabres 247 c 248 common /propar/ npro 249 c 250 common /resdat/ resfrt, resfrts 251 c 252 common /solpar/ imap, ivart, iDC, iPcond, Kspace, nGMRES, 253 & iconvflow, iconvsclr, idcsclr(2) 254 c 255 common /msympar/ indsym(5,5) 256 c 257 common /timdat/ time, CFLfld, CFLsld, Dtgl, Dtmax, alpha, 258 & etol, lstep, ifunc, itseq, istep, iter, 259 & nitr, almi, alfi, gami, flmpl, flmpr, 260 & dtol(2), iCFLworst, lskeep 261 c 262 common /timpar/ LCtime, ntseq 263 c 264 common /incomp/ numeqns(100), minIters, maxIters, 265 & iprjFlag, nPrjs, ipresPrjFlag, nPresPrjs, 266 & prestol, statsflow(6), statssclr(6), 267 & iverbose 268 c 269 character(8) :: ccode(13) 270 common /mtimer1/ ccode 271 c 272 integer flops, gbytes, sbytes 273 common /mtimer2/ flops, gbytes, sbytes, iclock, icd, icode, 274 & icode2, icode3 275 c 276 common /timer3/ cpu(11), cpu0(11), nacess(11) 277 c 278 character*80 title, ititle 279 common /title / title, ititle 280 c 281 character*8 machin 282 parameter ( machin = 'RS/6000 ' ) 283 parameter ( machfl = 4 ) 284 285 parameter 286 & ( zero = 0.0000000000000000000000000000000d0, 287 & pt125 = 0.1250000000000000000000000000000d0, 288 & pt25 = 0.2500000000000000000000000000000d0, 289 & pt33 = 0.3333333333333333333333333333333d0, 290 & pt39 = 0.3968502629920498686879264098181d0, 291 & pt5 = 0.5000000000000000000000000000000d0, 292 & pt57 = 0.5773502691896257645091487805020d0, 293 & pt66 = 0.6666666666666666666666666666667d0, 294 & pt75 = 0.7500000000000000000000000000000d0, 295 & one = 1.0000000000000000000000000000000d0, 296 & sqt2 = 1.4142135623730950488016887242097d0, 297 & onept5 = 1.5000000000000000000000000000000d0, 298 & two = 2.0000000000000000000000000000000d0, 299 & three = 3.0000000000000000000000000000000d0, 300 & four = 4.0000000000000000000000000000000d0, 301 & five = 5.0000000000000000000000000000000d0, 302 & pi = 3.1415926535897932384626433832795d0, 303 & inv1024sq = 9.5367431640625e-7) 304 305 c 306 c---------------------------------------------------------------------- 307 c 308 c.... element pointers 309 c 310 c mmat (MAXBLK) : pointer to interior element material number 311 c mmatb (MAXBLK) : pointer to boundary element material number 312 c mien (MAXBLK) : pointer to ien array 313 c mienb (MAXBLK) : pointer to ienb array 314 c miBCB (MAXBLK) : pointer to iBCB array 315 c mDt (MAXBLK) : pointer to Dt array 316 c mDC (MAXBLK) : pointer to DC array 317 c mBCB (MAXBLK) : pointer to BCB array 318 c mstiff (MAXBLK) : pointer to stiff array 319 c 320 c---------------------------------------------------------------------- 321 c 322 c.... common /aerfrc/ : aerodynamic forces 323 c 324 c Force(3) : components of the aerodynamic forces 325 c HFlux : total heat flux 326 c 327 c---------------------------------------------------------------------- 328 c 329 c.... common /astore/ : the dynamic memory allocation area 330 c 331 c a(...) : the blank array used for front-end data storage 332 c 333 c---------------------------------------------------------------------- 334 c 335 c.... common /blkdat/ : blocking data 336 c 337 c lcblk (10,MAXBLK+1) : blocking data for the interior elements 338 c lcblkb (10,MAXBLK+1) : blocking data for the boundary elements 339 c 340 c---------------------------------------------------------------------- 341 c 342 c.... common /bndnod/ : boundary nodes of boundary elements 343 c 344 c mnodeb (9,8,3) : boundary nodes of each element category and dimension 345 c 346 c---------------------------------------------------------------------- 347 c 348 c.... common /conpar/ : input constants 349 c 350 c numnp : number of nodal points 351 c numel : number of elements 352 c numelb : number of boundary elements 353 c numpbc : number of nodes having a boundary condition 354 c nen : maximum number of element nodes 355 c nfaces : maximum number of element faces 356 c nsd : number of space dimensions 357 c numflx : number of flux boundary nodes 358 c ndof : number of degrees of freedom per node 359 c iALE : ALE formulation flag 360 c icoord : coordinate system flag 361 c navier : Navier-Stokes calculation flag 362 c irs : restart option 363 c iexec : execute flag 364 c necho : input echo parameter 365 c ichem : equilibrium chemistry flag (for outchem.step dump) 366 c iRK : Runge-Kutta flag 367 c nshg : global number of shape functions (degrees of freedom, 368 c or equations). Computed from the specified p-order, 369 c the number of edges, and the number of faces (in the 370 c entire mesh) 371 c 372 c---------------------------------------------------------------------- 373 c 374 c.... common /datpnt/ : front-end data pointers 375 c 376 c mshp : pointer to shape-functions 377 c mshgl : pointer to local-grad-shape-functions 378 c mwght : pointer to quadrature weights 379 c mshpb : pointer to shape-functions of boundary elements 380 c mshglb : pointer to local-grad-shape-functions of bound. elem. 381 c mwghtb : pointer to quadrature weights of bound. elements 382 c mmut : pointer to table mu = mu (p,T) 383 c mrhot : pointer to table rho = rho (p,T) 384 c mxst : pointer to table xs = xs (p,T) 385 c 386 c---------------------------------------------------------------------- 387 c 388 c.... common /elmcat/ : element category information 389 c 390 c mcsyst : maximum number of element coordinate system 391 c melCat : maximum number of element categories 392 c nenCat (8,3) : number of nodes for each category and dimension 393 c nfaCat (8,3) : number of faces for each category and dimension 394 c 395 c---------------------------------------------------------------------- 396 c 397 c.... common /elmpar/ : element parameters 398 c 399 c lelCat : element category (P1, Q1, P2, Q2, etc.) 400 c lcsyst : element coordinate system 401 c iorder : element order (=k for Pk and Qk) 402 c nenb : number of element nodes per boundary sides 403 c maxsh : total number integration points 404 c maxshb : total number integration points of boundary elements 405 c nelblk : number of element blocks 406 c nelblb : number of boundary element blocks 407 c ndofl : number of degrees of freedom (for current block) 408 c nsymdl : number of d.o.f for symm. storage (for current block) 409 c nenl : number of element nodes (for current block) 410 c nfacel : number of element faces (for current block) 411 c nenbl : number of boundary element nodes 412 c intind : integration data index 413 c nintg : number of integration points 414 c mattyp : material type ( = 0 for fluid; = 1 for solid ) 415 c 416 c---------------------------------------------------------------------- 417 c 418 c.... common /genpar/ : control parameters 419 c 420 c E3nsd : NSD .eq. 3 flag; 0. for 2D, 1. for 3D 421 c I3nsd : NSD .eq. 3 flag; 0 for 2D, 1 for 3D 422 c nsymdf : number of d.o.f.'s in symm. storage (= ndof*(ndof+1)/2) 423 c ndofBC : dimension size of the boundary condition array BC 424 c ndiBCB : dimension size of the boundary condition array iBCB 425 c ndBCB : dimension size of the boundary condition array BCB 426 c Jactyp : Jacobian type flag 427 c jump : jump term computation flag 428 c ires : residual type computation flag 429 c iprec : block-diagonal preconditioner flag 430 c iprev : ypl array allocation flag 431 c ibound : boundary element flag 432 c idiff : diffusive flux vector flag 433 c ( = 0 not used; = 1 global reconstruction ) 434 c itau : type of tau to be used 435 c iLHScond : add contributiosn from the heat flux BC to the LHS 436 c tangency matrix. 437 c 438 c---------------------------------------------------------------------- 439 c 440 c.... common /inpdat/ : time sequence input data 441 c 442 c epstol (MAXTS) : tolerance for GMRES solvers 443 c Delt (MAXTS) : global time step 444 c CFLfl (MAXTS) : CFL number for fluid flow 445 c CFLsl (MAXTS) : CFL number for structural heating 446 c nstep (MAXTS) : number of time steps 447 c niter (MAXTS) : number of iterations per time step 448 c impl (MAXTS) : solver flag 449 c iturb (MAXTS) : turbulence model flag 450 c rhoinf (MAXTS) : time integration spectral radius paramter 451 c (0=Gears 1= trapezoidal rule) 452 c LHSupd (MAXTS) : LHS/preconditioner update 453 c loctim (MAXTS) : local time stepping flag 454 c 455 c---------------------------------------------------------------------- 456 c 457 c.... common /intdat/ : integration data 458 c 459 c intg (2,MAXTS) : integration parameters 460 c intpt (3) : integration pointers 461 c intptb(3) : integration pointers of boundary elements 462 c 463 c---------------------------------------------------------------------- 464 c 465 c.... common /shpdat/ : hierarchic shape function quadrature data 466 c 467 c Qpt (3,MAXQPT) : interior element quadrature points (xi,eta,zeta) 468 c Qwt (MAXQPT) : interior element quad. weights 469 c Qptb (2,MAXQPT) : boundary element quad. pnts. 470 c Qwtb (MAXQPT) : boundary element quad. weights 471 c nshape : number of interior element shape functions 472 c nshapeb : " " boundary " " " 473 c ngauss : number of interior element integration points 474 c ngaussb : " " boundary " " " 475 c---------------------------------------------------------------------- 476 c 477 c.... common /intpar/ : integration parameters 478 c 479 c Qpt (4,*) : xi, eta, zeta, weight of quadrature points 480 c indQpt(3,3,4) : index to quadrature points for a given rule 481 c numQpt(3,3,4) : number of quadrature points for a given rule 482 c intmax : number of allowable spatial integ. points per nsd 483 c 484 c---------------------------------------------------------------------- 485 c 486 c.... common /io / : io channels 487 c 488 c iin : input (main parameters) [INPUT.DAT] 489 c igeom : input (problem geometry) [GEOM.DAT] 490 c ipar : in/out (spectral mapping) [PARTITION.DAT] 491 c ibndc : input (problem boundary cond.) [BC.DAT] 492 c imat : input (element material types) [MATERIAL.DAT] 493 c iecho : output (echo of input) [ECHO.DAT] 494 c iout : output (result output) [OUTPUT.lstep] 495 c ichmou : output (chemistry output) [OUTCHM.lstep] 496 c irstin : input (input restart) [RESTAR.INP] 497 c irstou : output (output restart) [RESTAR.OUT] 498 c ihist : output (history output) [HISTOR.DAT] 499 c iflux : output (boundary flux) [FLUX.lstep] 500 c ierror : output (error messages) [ERROR.DAT] 501 c itable : input (equilibrium chemistry) [TABLE.DAT] 502 c iforce : output (aerodynamic forces) [FORCES.DAT] 503 c 504 c---------------------------------------------------------------------- 505 c 506 c.... common /ioname/ : io file names 507 c 508 c fin : input.dat 509 c fgeom : geom.dat 510 c fpar : partition.dat 511 c fbndc : bc.dat 512 c fmat : material.dat 513 c fecho : echo.dat 514 c frstin : restar.inp 515 c frstou : restar.out 516 c fhist : histor.dat 517 c ferror : error.dat 518 c ftable : table.dat 519 c fforce : forces.dat 520 c 521 c---------------------------------------------------------------------- 522 c 523 c.... common /itrpar/ : Preconditioned GMRES parameters 524 c 525 c eGMRES : finite difference interval 526 c lGMRES : number of GMRES cycles 527 c iKs : current Krylov vector 528 c ntotGM : total number of GMRES iterations 529 c 530 c---------------------------------------------------------------------- 531 c 532 c.... common /itrpnt/ : Preconditioned GMRES array pointers 533 c 534 c mHBrg : pointer to Hessenberg matrix 535 c meBrg : pointer to Hessenberg's RHS matrix 536 c myBrg : pointer to minimize solution matrix 537 c mRcos : pointer to Rotation Cosine of QR algorithm 538 c mRsin : pointer to Rotation Sine of QR algorithm 539 c 540 c---------------------------------------------------------------------- 541 c 542 c.... common /matpar/ : material constants 543 c 544 c pr : Prandtl number 545 c Planck : Planck's constant 546 c Stefan : Stefan's constant (for radiation) 547 c Nh : Avogadro's number 548 c Rh : universal gas constant 549 c Rgas : specific gas constant 550 c gamma : specific heat ratio 551 c gamma1 : gamma - 1 552 c s0 : reference specific entropy 553 c const : special constant 554 c xN2 : mole fraction of diatomic nitrogen 555 c xO2 : mole fraction of diatomic oxygen 556 c yN2 : mole fraction of diatomic nitrogen 557 c yO2 : mole fraction of diatomic oxygen 558 c Msh (5) : molar mass of species 559 c cpsh (5) : molar heat at constant pressure of species 560 c s0sh (5) : molar reference entropy of species 561 c h0sh (5) : molar heat of formation of species 562 c Rs (5) : specific gas constant of species 563 c cps (5) : specific heat at constant pressure of species 564 c cvs (5) : specific heat at constant volume of species 565 c h0s (5) : specific heat of formation of species 566 c Trot (5) : characteristic rotational temperature of species 567 c sigs (5) : symmetry factor of species 568 c Tvib (5) : characteristic vibrational temperature of species 569 c g0s (5) : ground degeneracy of electronic energy 570 c dofs (5) : degrees of freedom of species 571 c ithm : thermodynamic property flag 572 c 573 c---------------------------------------------------------------------- 574 c 575 c.... common /matdat/ : material data 576 c 577 c datmat (3,5,2) : material data 578 c matflg (5,100) : material type flag 579 c nummat : number of materials 580 c mexist : flag indicating the presence of MATERIAL.DAT 581 c 582 c---------------------------------------------------------------------- 583 c 584 c.... common /outpar/ : output parameters 585 c 586 c ro : density rescaling factor for output 587 c vel : velocity rescaling factor for output 588 c temper : temperature rescaling factor for output 589 c press : pressure rescaling factor for output 590 c entrop : entropy rescaling factor for output 591 c ntout : number of steps between consecutive printouts 592 c ioform : output I/O format 593 c 594 c---------------------------------------------------------------------- 595 c 596 c.... common /point / : dynamic storage pointer management data 597 c 598 c mbeg : pointer to the beginning of the free storage 599 c mend : pointer to the end of the storage 600 c mprec : precision of the floating point data 601 c 602 c---------------------------------------------------------------------- 603 c 604 c.... common /precis/ : finite difference interval data 605 c 606 c epsM : square root of machine precision 607 c iabres : absolute value residual flag 608 c 609 c---------------------------------------------------------------------- 610 c 611 c....common /propar/ : processor related information 612 c 613 c npro : number of virtual processors for the current block 614 c 615 c---------------------------------------------------------------------- 616 c 617 c....common /resdat/ : residual statistics data 618 c 619 c resfrt : first residual of convergence 620 c 621 c---------------------------------------------------------------------- 622 c 623 c.... common /solpar/ : solution parameters 624 c 625 c imap : permutation mapping flag 626 c ivart : variational formulation type 627 c iDC : DC type 628 c iPcond : type of preconditioner 629 c Kspace : dimension of Krylov space 630 c nGMRES : maximum number of GMRES iterations 631 c 632 c---------------------------------------------------------------------- 633 c 634 c.... common /sympar/ : symmetric storage parameters 635 c 636 c indsym (5,5) : mapping from 2D storage to symmetric one 637 c 638 c---------------------------------------------------------------------- 639 c 640 c.... common /timdat/ : time data 641 c 642 c time : current run time 643 c CFLfld : CFL number for fluid flow 644 c CFLsld : CFL number for structural heating 645 c Dtgl : inverse of global time step 646 c Dtmax : maximum delta-time 647 c alpha : trapezoidal rule parameter 648 c etol : epsilon tolerance for GMRES 649 c lstep : current time step 650 c ifunc : func. eval. counter (=niter*(lstep-lstep0) + iter) 651 c itseq : sequence number 652 c istep : step number (reseted at the beginning of the run) 653 c iter : iteration number 654 c nitr : number of multi-corrector iterations for this sequence 655 c 656 c---------------------------------------------------------------------- 657 c 658 c.... common /timpar/ : time integration parameters 659 c 660 c LCtime : local time stepping flag 661 c ntseq : number of time sequences 662 c 663 c---------------------------------------------------------------------- 664 c 665 c.... common /timer1/ : timer parameters 666 c.... common /timer2/ : timer parameters 667 c.... common /timer3/ : timer parameters 668 c 669 c ccode(13) : timing entities codes 670 c flops : flop counter 671 c gbytes : byte counter for gather operation 672 c sbytes : byte counter for scatter operation 673 c iclock : wall-clock time (in milliseconds) 674 c icd : number of timing entities 675 c icode : current timer code 676 c icode2 : last timer code 677 c icode3 : next-to-last timer code 678 c cpu(11) : cpu time of each entity 679 c cpu0(11) : initial cpu time of each entity 680 c nacess(11) : number of times each entity is accessed 681 c 682 c---------------------------------------------------------------------- 683 c 684 c.... common /title / : problem title 685 c 686 c title : problem title 687 c ititle : problem title (with form feed) 688 c 689 c---------------------------------------------------------------------- 690 c 691 c.... common /avging / : nfath 692 c 693 c nfath : total number of global fathers over which certain 694 c quantities will be averaged 695 c 696 c---------------------------------------------------------------------- 697 c 698 c.... parameters : machine data 699 c 700 c machin : machine type 701 c (set parameter) 702 c machfl : single precision floating point lenght in bytes 703 c (set parameter) 704 c 705 c---------------------------------------------------------------------- 706 c 707 c.... parameters : useful constants 708 c 709 c zero : 0.0 710 c pt125 : 0.125 711 c pt25 : 0.25 712 c pt33 : 0.33 (1/3) 713 c pt39 : 2^(-4/3) 714 c pt5 : 0.5 715 c pt57 : 1/sqrt(3) 716 c pt66 : 0.66 (2/3) 717 c pt75 : 0.75 718 c one : 1.0 719 c sqt2 : sqrt(2) 720 c onept5 : 1.5 721 c two : 2.0 722 c three : 3.0 723 c four : 4.0 724 c five : 5.0 725 c pi : the magical number :-) 726 c 727 c---------------------------------------------------------------------- 728 c 729 c Zdenek Johan, Winter 1991. 730 c 731 c---------------------------------------------------------------------- 732