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