Lines Matching refs:grid
104 GridInfo *grid; member
128 GridInfo grid; in main() local
147 PetscCall(SetParams(¶m, &grid)); in main()
148 PetscCall(ReportParams(¶m, &grid)); in main()
153 …PetscCall(DMDACreate2d(comm, grid.bx, grid.by, grid.stencil, grid.ni, grid.nj, PETSC_DECIDE, PETSC… in main()
167 user->grid = &grid; in main()
309 GridInfo *grid = user->grid; in HorizVelocity() local
313 x = (i - grid->jlid) * grid->dx; in HorizVelocity()
314 z = (j - grid->jlid - 0.5) * grid->dz; in HorizVelocity()
326 GridInfo *grid = user->grid; in VertVelocity() local
330 x = (i - grid->jlid - 0.5) * grid->dx; in VertVelocity()
331 z = (j - grid->jlid) * grid->dz; in VertVelocity()
343 GridInfo *grid = user->grid; in Pressure() local
346 x = (i - grid->jlid - 0.5) * grid->dx; in Pressure()
347 z = (j - grid->jlid - 0.5) * grid->dz; in Pressure()
358 GridInfo *grid = user->grid; in CalcSecInv() local
359 PetscInt ilim = grid->ni - 1, jlim = grid->nj - 1; in CalcSecInv()
368 if (j <= grid->jlid) return EPS_ZERO; in CalcSecInv()
383 if (j == grid->jlid + 1) uS = 0.0; in CalcSecInv()
387 if (j < grid->jlid) return EPS_ZERO; in CalcSecInv()
401 if (j == grid->jlid) { in CalcSecInv()
412 eps11 = (uE - uW) / grid->dx; in CalcSecInv()
413 eps22 = (wN - wS) / grid->dz; in CalcSecInv()
414 eps12 = 0.5 * ((uN - uS) / grid->dz + (wE - wW) / grid->dx); in CalcSecInv()
465 GridInfo *grid = user->grid; in XMomentumResidual() local
466 PetscScalar dx = grid->dx, dz = grid->dz; in XMomentumResidual()
470 PetscInt jlim = grid->nj - 1; in XMomentumResidual()
516 GridInfo *grid = user->grid; in ZMomentumResidual() local
517 PetscScalar dx = grid->dx, dz = grid->dz; in ZMomentumResidual()
521 PetscInt ilim = grid->ni - 1; in ZMomentumResidual()
569 GridInfo *grid = user->grid; in ContinuityResidual() local
574 dudx = (uE - uW) / grid->dx; in ContinuityResidual()
577 dwdz = (wN - wS) / grid->dz; in ContinuityResidual()
586 GridInfo *grid = user->grid; in EnergyResidual() local
587 PetscScalar dx = grid->dx, dz = grid->dz; in EnergyResidual()
588 PetscInt ilim = grid->ni - 1, jlim = grid->nj - 1, jlid = grid->jlid; in EnergyResidual()
645 GridInfo *grid = user->grid; in ShearStress() local
646 PetscInt ilim = grid->ni - 1, jlim = grid->nj - 1; in ShearStress()
649 if (j <= grid->jlid || i < j || i == ilim || j == jlim) return EPS_ZERO; in ShearStress()
661 if (j == grid->jlid + 1) uS = 0.0; in ShearStress()
672 return (uN - uS) / grid->dz + (wE - wW) / grid->dx; in ShearStress()
679 GridInfo *grid = user->grid; in XNormalStress() local
680 PetscScalar dx = grid->dx, dz = grid->dz; in XNormalStress()
681 PetscInt ilim = grid->ni - 1, jlim = grid->nj - 1, ivisc; in XNormalStress()
683 if (i < j || j <= grid->jlid) return EPS_ZERO; in XNormalStress()
718 GridInfo *grid = user->grid; in ZNormalStress() local
719 PetscScalar dz = grid->dz; in ZNormalStress()
720 PetscInt ilim = grid->ni - 1, jlim = grid->nj - 1, ivisc; in ZNormalStress()
723 if (i < j || j <= grid->jlid) return EPS_ZERO; in ZNormalStress()
758 PetscErrorCode SetParams(Parameter *param, GridInfo *grid) in SetParams() argument
780 PetscCall(PetscOptionsGetInt(NULL, NULL, "-jfault", &grid->jfault, NULL)); in SetParams()
781 grid->ni = 82; in SetParams()
782 PetscCall(PetscOptionsGetInt(NULL, NULL, "-ni", &grid->ni, NULL)); in SetParams()
784 grid->dx = param->width / ((PetscReal)(grid->ni - 2)); /* km */ in SetParams()
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()
788 grid->inose = 0; /* gridpoints*/ in SetParams()
789 PetscCall(PetscOptionsGetInt(NULL, NULL, "-inose", &grid->inose, NULL)); in SetParams()
790 grid->bx = DM_BOUNDARY_NONE; in SetParams()
791 grid->by = DM_BOUNDARY_NONE; in SetParams()
792 grid->stencil = DMDA_STENCIL_BOX; in SetParams()
793 grid->dof = 4; in SetParams()
794 grid->stencil_width = 2; in SetParams()
795 grid->mglevels = 1; in SetParams()
875 grid->dx = grid->dx / param->L; /* dim'less */ 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()
881 grid->corner = grid->jlid + 1; /* gridcells */ in SetParams()
892 PetscErrorCode ReportParams(Parameter *param, GridInfo *grid) in ReportParams() argument
907 …etscInt_FMT " [dx,dz] = %g, %g km\n", grid->ni, grid->nj, (double)(grid->dx * param->L), (do… in ReportParams()
908 …jlid = %3" PetscInt_FMT " jfault = %3" PetscInt_FMT " \n", grid->jlid, grid->jfault)); in ReportParams()
964 GridInfo *grid; in Initialize() local
974 grid = user->grid; in Initialize()
982 else if (j <= grid->jlid) x[j][i].u = 0.0; in Initialize()
986 else if (j <= grid->jlid) x[j][i].w = 0.0; in Initialize()
989 if (i < j || j <= grid->jlid) x[j][i].p = 0.0; in Initialize()
992 x[j][i].T = PetscMin(grid->dz * (j - 0.5), 1.0); in Initialize()
1006 GridInfo *grid; in DoOutput() local
1018 grid = user->grid; in DoOutput()
1036 PetscCall(VecSetValue(pars, 0, (PetscScalar)grid->ni, INSERT_VALUES)); in DoOutput()
1037 PetscCall(VecSetValue(pars, 1, (PetscScalar)grid->nj, INSERT_VALUES)); in DoOutput()
1038 PetscCall(VecSetValue(pars, 2, (PetscScalar)grid->dx, INSERT_VALUES)); in DoOutput()
1039 PetscCall(VecSetValue(pars, 3, (PetscScalar)grid->dz, INSERT_VALUES)); in DoOutput()
1044 PetscCall(VecSetValue(pars, 8, (PetscScalar)grid->jlid, INSERT_VALUES)); in DoOutput()
1046 PetscCall(VecSetValue(pars, 10, (PetscScalar)grid->jfault, INSERT_VALUES)); in DoOutput()
1096 GridInfo *grid; in ViscosityField() local
1105 grid = user->grid; in ViscosityField()
1116 /* dx = grid->dx; */ dz = grid->dz; in ViscosityField()
1118 ilim = grid->ni - 1; in ViscosityField()
1119 jlim = grid->nj - 1; in ViscosityField()
1195 GridInfo *grid = user->grid; in SlabVel() local
1199 else if (j <= grid->jfault) return 0.0; in SlabVel()
1204 else if (j <= grid->jfault) return 0.0; in SlabVel()
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()
1292 GridInfo *grid = user->grid; in FormFunctionLocal() local
1315 … else if (j <= grid->jlid || (j < grid->corner + grid->inose && i < grid->corner + grid->inose)) { in FormFunctionLocal()
1345 …} else if (j <= grid->jlid || (j < grid->corner + grid->inose && i < grid->corner + grid->inose)) { in FormFunctionLocal()
1373 …if (i < j || j <= grid->jlid || (j < grid->corner + grid->inose && i < grid->corner + grid->inose)… in FormFunctionLocal()