Commits

Anonymous committed 87cd717

Changed initial perm and porosity values to -999 to catch uninitialized values.

Comments (0)

Files changed (2)

src/pflotran/material.F90

   allocate(material_property)
   material_property%id = 0
   material_property%name = ''
-  material_property%permeability = 0.d0
+  ! initialize to -999.d0 to catch bugs
+  material_property%permeability = -999.d0
   material_property%isotropic_permeability = PETSC_TRUE
   material_property%vertical_anisotropy_ratio = 0.d0
   material_property%permeability_scaling_factor = 0.d0
   material_property%permeability_min_scale_fac = 1.d0
   material_property%permeability_dataset_name = ''
   nullify(material_property%permeability_dataset)
-  material_property%porosity = 0.d0
+  ! initialize to -999.d0 to catch bugs
+  material_property%porosity = -999.d0
   material_property%porosity_dataset_name = ''
   nullify(material_property%porosity_dataset)
   material_property%tortuosity = 1.d0

src/pflotran/timestepper.F90

   do
       
     call PetscTime(log_start_time, ierr)
-
+    print *, 'FLOW STEP SIZE: ', option%flow_dt/3600.d0/24.d0/365.d0, option%time/3600.d0/24.d0/365.d0
     select case(option%iflowmode)
       case(MPH_MODE,TH_MODE,THC_MODE,THMC_MODE,IMS_MODE,MIS_MODE,FLASH2_MODE,G_MODE)
         call SNESSolve(solver%snes, PETSC_NULL_OBJECT, field%flow_xx, ierr)
 !       The below is set within RTUpdateAuxVarsPatch() when PETSC_TRUE,PETSC_TRUE,* are passed
 !       patch%aux%RT%aux_vars_up_to_date = PETSC_TRUE 
     endif
+        print *, 'TRAN STEP SIZE: ', option%tran_dt/3600.d0/24.d0/365.d0
+
     if (realization%reaction%use_log_formulation) then
       call VecCopy(field%tran_xx,field%tran_log_xx,ierr)
       call VecLog(field%tran_log_xx,ierr)