87   PetscInt  ivisc, adv_scheme, ibound, output_ivisc;  member
230   if (param->ivisc == VISC_CONST && !param->stop_solve) {  in UpdateSolution()
231     param->ivisc = VISC_CONST;  in UpdateSolution()
242   if (param->ivisc >= VISC_DIFN && !param->stop_solve) {  in UpdateSolution()
424   PetscInt    iVisc     = param->ivisc;  in Viscosity()
474   if (param->ivisc == VISC_DIFN || param->ivisc >= VISC_DISL) { /* viscosity is T-dependent */  in XMomentumResidual()
480     if (param->ivisc >= VISC_DISL) { /* olivine dislocation creep */  in XMomentumResidual()
502   if (param->ivisc != VISC_CONST) {  in XMomentumResidual()
527   if (param->ivisc == VISC_DIFN || param->ivisc >= VISC_DISL) { /* viscosity is T-dependent */  in ZMomentumResidual()
533     if (param->ivisc >= VISC_DISL) { /* olivine dislocation creep */  in ZMomentumResidual()
556   if (param->ivisc != VISC_CONST) {  in ZMomentumResidual()
681   PetscInt    ilim = grid->ni - 1, jlim = grid->nj - 1, ivisc;  in XNormalStress()  local
685   ivisc   = param->ivisc;  in XNormalStress()
691     if (ivisc >= VISC_DISL) epsC = CalcSecInv(x, i, j, CELL_CENTER, user);  in XNormalStress()
702     if (ivisc >= VISC_DISL) epsC = CalcSecInv(x, i, j, CELL_CORNER, user);  in XNormalStress()
720   PetscInt    ilim = grid->ni - 1, jlim = grid->nj - 1, ivisc;  in ZNormalStress()  local
725   ivisc   = param->ivisc;  in ZNormalStress()
731     if (ivisc >= VISC_DISL) epsC = CalcSecInv(x, i, j, CELL_CENTER, user);  in ZNormalStress()
741     if (ivisc >= VISC_DISL) epsC = CalcSecInv(x, i, j, CELL_CORNER, user);  in ZNormalStress()
816   param->ivisc        = 3;    /* 0=isovisc, 1=difn creep, 2=disl creep, 3=full */  in SetParams()
833   PetscCall(PetscOptionsGetInt(NULL, NULL, "-ivisc", &param->ivisc, NULL));  in SetParams()
836   param->output_ivisc = param->ivisc;  in SetParams()
912     if (param->ivisc == VISC_CONST) {  in ReportParams()
915     } else if (param->ivisc == VISC_DIFN) {  in ReportParams()
918     } else if (param->ivisc == VISC_DISL) {  in ReportParams()
921     } else if (param->ivisc == VISC_FULL) {  in ReportParams()
948 …if (param->output_ivisc != param->ivisc) PetscCall(PetscPrintf(PETSC_COMM_WORLD, "                …  in ReportParams()
1019   ivt   = param->ivisc;  in DoOutput()
1021   param->ivisc = param->output_ivisc;  in DoOutput()
1049       PetscCall(VecSetValue(pars, 13, (PetscScalar)param->ivisc, INSERT_VALUES));  in DoOutput()
1085   param->ivisc = ivt;  in DoOutput()
1106   ivt          = param->ivisc;  in ViscosityField()
1107   param->ivisc = param->output_ivisc;  in ViscosityField()
1144   param->ivisc = ivt;  in ViscosityField()