1. PFLOTRAN
  2. Untitled project
  3. pflotran-dev

Commits

Peter Lichtner  committed c003dce

Surface armoring continued...

  • Participants
  • Parent commits cc489c7
  • Branches default

Comments (0)

Files changed (4)

File src/pflotran/option.F90

View file
     
     PetscReal :: minimum_hydrostatic_pressure
     
-!   PetscBool :: update_mnrl_surf_with_porosity
-    
     PetscBool :: jumpstart_kinetic_sorption
     PetscBool :: no_checkpoint_kinetic_sorption
     PetscBool :: no_restart_kinetic_sorption

File src/pflotran/reaction.F90

View file
         reaction%update_mineral_surface_area = PETSC_TRUE
       case('UPDATE_MNRL_SURF_AREA_WITH_POR')
         reaction%update_mnrl_surf_with_porosity = PETSC_TRUE
+      case('UPDATE_ARMOR_MINERAL_SURFACE')
+        reaction%update_armor_mineral_surface = PETSC_TRUE
       case('MOLAL','MOLALITY')
         reaction%initialize_with_molality = PETSC_TRUE
       case('ACTIVITY_H2O','ACTIVITY_WATER')

File src/pflotran/reaction_aux.F90

View file
     PetscReal :: minimum_porosity
     PetscBool :: update_mineral_surface_area
     PetscBool :: update_mnrl_surf_with_porosity
-    
+    PetscBool :: update_armor_mineral_surface
+
+
     PetscBool :: use_sandbox
     
   end type reaction_type
   reaction%minimum_porosity = 0.d0
   reaction%update_mineral_surface_area = PETSC_FALSE
   reaction%update_mnrl_surf_with_porosity = PETSC_FALSE
+  reaction%update_armor_mineral_surface = PETSC_FALSE
   reaction%use_sandbox = PETSC_FALSE
 
   ReactionCreate => reaction

File src/pflotran/realization.F90

View file
 
   if (reaction%update_mineral_surface_area) then
     porosity_scale = 1.d0
-!   if (option%update_mnrl_surf_with_porosity) then
     if (reaction%update_mnrl_surf_with_porosity) then
       ! placing the get/restore array calls within the condition will
       ! avoid improper access.