ERROR: Variable ID 134 not supported for the process model: RICHARDS

[Eugenio Pescimoro]

Dear PFLOTRAN developers,

I am using the following PLFOTRAN input file to run a diffusion simulation:

#================================================

SIMULATION
  SIMULATION_TYPE SUBSURFACE
  PROCESS_MODELS
    SUBSURFACE_FLOW flow
      MODE RICHARDS
    /
    SUBSURFACE_TRANSPORT transport
      MODE GIRT
      OPTIONS
      /
    /
  /
END

SUBSURFACE


CHEMISTRY
  PRIMARY_SPECIES
    TRACER
  /
  OUTPUT
    TOTAL
    TRACER
  /
END
#=========================== discretization ===================================
GRID
  TYPE unstructured_explicit full_mesh_vol_area.uge
  GRAVITY 0.d0 0.d0 0.d0
END


#=========================== fluid properties =================================
FLUID_PROPERTY
  DIFFUSION_COEFFICIENT 4.5d-6
END

DATASET Permeability
  FILENAME dfn_properties.h5
END

#=========================== material properties ==============================
MATERIAL_PROPERTY soil1
  ID 1
  POROSITY 0.25
  TORTUOSITY 0.2d0
  CHARACTERISTIC_CURVES default
  PERMEABILITY
    DATASET Permeability
  /
END

#=========================== characteristic curves ============================
CHARACTERISTIC_CURVES default
  SATURATION_FUNCTION VAN_GENUCHTEN
    M 0.5d0
    ALPHA  1.d-3
    LIQUID_RESIDUAL_SATURATION 0.01d0
    MAX_CAPILLARY_PRESSURE 1.d8
  /
  PERMEABILITY_FUNCTION MUALEM_VG_LIQ
    M 0.5d0
    LIQUID_RESIDUAL_SATURATION 0.1d0
  /
END

#=========================== output options ===================================
OUTPUT
  PERIODIC TIME 100 y ## Comment this out if you don't want to visualize the results
#  FORMAT TECPLOT BLOCK
  PRINT_PRIMAL_GRID
  FORMAT VTK
  ACKNOWLEDGE_VTK_FLAW
  MASS_BALANCE
  VARIABLES
    LIQUID_PRESSURE
    PERMEABILITY_X
    PERMEABILITY_Y
    PERMEABILITY_Z
    POROSITY
    # SOLUTE_CONCENTRATION
  /
END

#=========================== times ============================================
TIME
  INITIAL_TIMESTEP_SIZE  1.d-1 y
  FINAL_TIME 1.d3 y
  MAXIMUM_TIMESTEP_SIZE 1.d2 y
  #STEADY_STATE
END

#=========================== regions ==========================================
REGION All
  COORDINATES
    -5.d20 -5.d20 -5.d20
    5.d20 5.d20 5.d20
  /
END 

REGION inflow
  FILE boundary_bottom.ex
END

REGION outflow
  FILE boundary_top.ex
END

TRANSPORT_CONDITION initial
  TYPE DIRICHLET_ZERO_GRADIENT 
    CONSTRAINT_LIST
    0.d0 initial
  /
END

TRANSPORT_CONDITION inject
  TYPE DIRICHLET_ZERO_GRADIENT
    CONSTRAINT_LIST
    0.d0 inlet
  /
END

CONSTRAINT initial
  CONCENTRATIONS
    TRACER 1.d-20 T
  /
END

CONSTRAINT inlet
  CONCENTRATIONS
    TRACER 1.d0 T
  /
END

#=========================== flow conditions ==================================
FLOW_CONDITION initial
  TYPE
     LIQUID_PRESSURE dirichlet 
  /
  LIQUID_PRESSURE 1.01325d6
END


FLOW_CONDITION outflow 
  TYPE 
     LIQUID_PRESSURE dirichlet
  /
  LIQUID_PRESSURE 1.01325d6
END

FLOW_CONDITION inflow
  TYPE
    LIQUID_PRESSURE dirichlet
  /
  LIQUID_PRESSURE 1.01325d6
END

#=========================== condition couplers ===============================
# initial condition
INITIAL_CONDITION
  FLOW_CONDITION initial
  TRANSPORT_CONDITION initial
  REGION All
END

BOUNDARY_CONDITION INFLOW 
  FLOW_CONDITION inflow
  TRANSPORT_CONDITION inject
  REGION inflow
END

BOUNDARY_CONDITION OUTFLOW
  FLOW_CONDITION outflow
  TRANSPORT_CONDITION initial 
  REGION outflow
END

#=========================== stratigraphy couplers ============================
STRATA
  REGION All 
  MATERIAL soil1
END

END_SUBSURFACE

I am interested in the solute concentration at the outlet of the domain [m^-3] however when I add the SOLUTE_CONCENTRATION keyword at the "output options" I do get the following error message:

"ERROR: Variable ID 134 not supported for the process model: RICHARDS. Please look at variables.F90 to match this integer ID with a variable parameter constant.  If you feel this is in error, please email this message to pflotr...@googlegroups.com."

Do you know how to troubleshoot it?

Best wishes,

Eugenio Pescimoro Postdoc Instituto de Diagnóstico Ambiental y Estudios del Agua (IDAEA-CSIC) CONSEJO SUPERIOR DE INVESTIGACIONES CIENTÍFICAS ✆ +34 93 400 61 00 (extension 437728) @ eugenio....@idaea.csic.es ✉ C/ Jordi Girona, 18-26, 08034, Barcelona

[Hammond, Glenn E]

Eugenio,

 

SOLUTE_CONCENTRATION under the outer OUTPUT block is for flow modes (e.g., the tightly coupled tracer in ZFLOW). You should not include that card for RICHARDS mode as the OUTPUT block under CHEMISTRY will output solution concentrations from the transport side of the code.

 

Glenn

 

From: 'Eugenio Pescimoro' via pflotran-dev <pflotr...@googlegroups.com>
Date: Tuesday, June 3, 2025 at 3:23 AM
To: pflotr...@googlegroups.com <pflotr...@googlegroups.com>
Subject: [pflotran-dev: 6337] ERROR: Variable ID 134 not supported for the process model: RICHARDS

Check twice before you click! This email originated from outside PNNL.

 

Eugenio Pescimoro Postdoc Instituto de Diagnóstico Ambiental y Estudios del Agua (IDAEA-CSIC) CONSEJO SUPERIOR DE INVESTIGACIONES CIENTÍFICAS ✆ +34 93 400 61 00 (extension 437728) @ eugenio....@idaea.csic.es ✉ C/ Jordi Girona, 18-26, 08034, Barcelona

 

--
You received this message because you are subscribed to the Google Groups "pflotran-dev" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pflotran-dev...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/pflotran-dev/1954297262.2172623.1748946225197.JavaMail.zimbra%40cid.csic.es.