Lines Matching refs:dz
98 PetscReal dx, dz; member
314 z = (j - grid->jlid - 0.5) * grid->dz; in HorizVelocity()
331 z = (j - grid->jlid) * grid->dz; in VertVelocity()
347 z = (j - grid->jlid - 0.5) * grid->dz; in Pressure()
413 eps22 = (wN - wS) / grid->dz; in CalcSecInv()
414 eps12 = 0.5 * ((uN - uS) / grid->dz + (wE - wW) / grid->dx); in CalcSecInv()
466 PetscScalar dx = grid->dx, dz = grid->dz; in XMomentumResidual() local
475 TS = param->potentialT * TInterp(x, i, j - 1) * PetscExpScalar((j - 1.0) * dz * z_scale); in XMomentumResidual()
477 else TN = param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * z_scale); in XMomentumResidual()
478 TW = param->potentialT * x[j][i].T * PetscExpScalar((j - 0.5) * dz * z_scale); in XMomentumResidual()
479 TE = param->potentialT * x[j][i + 1].T * PetscExpScalar((j - 0.5) * dz * z_scale); in XMomentumResidual()
487 etaN = Viscosity(TN, epsN, dz * (j + 0.5), param); in XMomentumResidual()
488 etaS = Viscosity(TS, epsS, dz * (j - 0.5), param); in XMomentumResidual()
489 etaW = Viscosity(TW, epsW, dz * j, param); in XMomentumResidual()
490 etaE = Viscosity(TE, epsE, dz * j, param); in XMomentumResidual()
494 else dudzN = etaN * (x[j + 1][i].u - x[j][i].u) / dz; in XMomentumResidual()
495 dudzS = etaS * (x[j][i].u - x[j - 1][i].u) / dz; in XMomentumResidual()
500 + (dudxE - dudxW) / dx + (dudzN - dudzS) / dz; in XMomentumResidual()
506 residual += (dudxE - dudxW) / dx + (dwdxN - dwdxS) / dz; in XMomentumResidual()
517 PetscScalar dx = grid->dx, dz = grid->dz; in ZMomentumResidual() local
528 TN = param->potentialT * x[j + 1][i].T * PetscExpScalar((j + 0.5) * dz * z_scale); in ZMomentumResidual()
529 TS = param->potentialT * x[j][i].T * PetscExpScalar((j - 0.5) * dz * z_scale); in ZMomentumResidual()
530 TW = param->potentialT * TInterp(x, i - 1, j) * PetscExpScalar(j * dz * z_scale); in ZMomentumResidual()
532 else TE = param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * z_scale); in ZMomentumResidual()
540 etaN = Viscosity(TN, epsN, dz * (j + 1.0), param); in ZMomentumResidual()
541 etaS = Viscosity(TS, epsS, dz * (j + 0.0), param); in ZMomentumResidual()
542 etaW = Viscosity(TW, epsW, dz * (j + 0.5), param); in ZMomentumResidual()
543 etaE = Viscosity(TE, epsE, dz * (j + 0.5), param); in ZMomentumResidual()
545 dPdz = (x[j + 1][i].p - x[j][i].p) / dz; in ZMomentumResidual()
546 dwdzN = etaN * (x[j + 1][i].w - x[j][i].w) / dz; in ZMomentumResidual()
547 dwdzS = etaS * (x[j][i].w - x[j - 1][i].w) / dz; in ZMomentumResidual()
548 if (i == ilim) dwdxE = etaE * (x[j][i].u - x[j + 1][i].u) / dz; in ZMomentumResidual()
554 + (dwdzN - dwdzS) / dz + (dwdxE - dwdxW) / dx; in ZMomentumResidual()
557 dudzE = etaE * (x[j + 1][i].u - x[j][i].u) / dz; in ZMomentumResidual()
558 dudzW = etaW * (x[j + 1][i - 1].u - x[j][i - 1].u) / dz; in ZMomentumResidual()
560 residual += (dwdzN - dwdzS) / dz + (dudzE - dudzW) / dx; in ZMomentumResidual()
577 dwdz = (wN - wS) / grid->dz; in ContinuityResidual()
587 PetscScalar dx = grid->dx, dz = grid->dz; in EnergyResidual() local
593 dTdzN = (x[j + 1][i].T - x[j][i].T) / dz; in EnergyResidual()
594 dTdzS = (x[j][i].T - x[j - 1][i].T) / dz; in EnergyResidual()
598 residual = ((dTdzN - dTdzS) / dz + /* diffusion term */ in EnergyResidual()
600 dx * dz / param->peclet; in EnergyResidual()
625 fE = uE * TE * dz; in EnergyResidual()
626 fW = uW * TW * dz; in EnergyResidual()
630 …scAbsScalar(uE) * (-x[j][i + 2].T + 3.0 * (x[j][i + 1].T - x[j][i].T) + x[j][i - 1].T) / 8.0) * dz; in EnergyResidual()
631 …scAbsScalar(uW) * (-x[j][i + 1].T + 3.0 * (x[j][i].T - x[j][i - 1].T) + x[j][i - 2].T) / 8.0) * dz; in EnergyResidual()
672 return (uN - uS) / grid->dz + (wE - wW) / grid->dx; in ShearStress()
680 PetscScalar dx = grid->dx, dz = grid->dz; in XNormalStress() local
690 TC = param->potentialT * x[j][i].T * PetscExpScalar((j - 0.5) * dz * z_scale); in XNormalStress()
692 etaC = Viscosity(TC, epsC, dz * j, param); in XNormalStress()
701 TC = param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * z_scale); in XNormalStress()
703 etaC = Viscosity(TC, epsC, dz * (j + 0.5), param); in XNormalStress()
719 PetscScalar dz = grid->dz; in ZNormalStress() local
730 TC = param->potentialT * x[j][i].T * PetscExpScalar((j - 0.5) * dz * z_scale); in ZNormalStress()
732 etaC = Viscosity(TC, epsC, dz * j, param); in ZNormalStress()
740 TC = param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * z_scale); in ZNormalStress()
742 etaC = Viscosity(TC, epsC, dz * (j + 0.5), param); in ZNormalStress()
749 return 2.0 * etaC * (wN - wS) / dz - pC; in ZNormalStress()
785 grid->dz = grid->dx * PetscTanReal(param->slab_dip); /* km */ in SetParams()
786 grid->nj = (PetscInt)(param->depth / grid->dz + 3.0); /* gridpoints*/ in SetParams()
787 param->depth = grid->dz * (grid->nj - 2); /* km */ in SetParams()
876 grid->dz = grid->dz / param->L; /* dim'less */ in SetParams()
877 grid->jlid = (PetscInt)(param->lid_depth / grid->dz); /* gridcells */ in SetParams()
878 grid->jfault = (PetscInt)(param->fault_depth / grid->dz); /* gridcells */ in SetParams()
879 param->lid_depth = grid->jlid * grid->dz; /* dim'less */ in SetParams()
880 param->fault_depth = grid->jfault * grid->dz; /* dim'less */ in SetParams()
907 …= %g, %g km\n", grid->ni, grid->nj, (double)(grid->dx * param->L), (double)(grid->dz * param->L))); in ReportParams()
992 x[j][i].T = PetscMin(grid->dz * (j - 0.5), 1.0); in Initialize()
1039 PetscCall(VecSetValue(pars, 3, (PetscScalar)grid->dz, INSERT_VALUES)); in DoOutput()
1099 PetscReal eps, /* dx,*/ dz, T, epsC, TC; in ViscosityField() local
1116 /* dx = grid->dx; */ dz = grid->dz; in ViscosityField()
1125 …T = PetscRealPart(param->potentialT * x[j][i].T * PetscExpScalar((j - 0.5) * dz * param->z_scale)); in ViscosityField()
1127 …TC = PetscRealPart(param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * param->z_scale)); in ViscosityField()
1136 v[j][i].p = Viscosity(T, eps, dz * (j - 0.5), param); in ViscosityField()
1137 v[j][i].T = Viscosity(TC, epsC, dz * j, param); in ViscosityField()
1214 if (plate == PLATE_LID) z = (j - 0.5) * user->grid->dz; in PlateModel()
1215 else z = (j - 0.5) * user->grid->dz * param->cb; /* PLATE_SLAB */ in PlateModel()