8000 Recharge array multiplier input file issue · Issue #577 · pypest/pyemu · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

Recharge array multiplier input file issue #577

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom
Open
iaped opened this issue Feb 11, 2025 · 0 comments
Open

Recharge array multiplier input file issue #577

iaped opened this issue Feb 11, 2025 · 0 comments

Comments

@iaped
Copy link
iaped commented Feb 11, 2025

I am trying to use PEST++ GLM to calibrate a recharge multiplier array in my MODFLOW 6 input file. However, when I run GLM it says:

 WARNING: jacobian matrix has no non-zeros - the parameter pertubations
         have no effect on the control file observations.
         This usually means something is not setup correctly.
save_binary() WARNING: matrix has no non-zeros...

I think this is because PEST is not opening the template file when it gets to the OPEN/CLOSE part of the recharge file. I am using zone calibration as I have 5 different zones of recharge. The file I am using as my input/template is rchMult.dat. My recharge input file has the following format:

BEGIN options
  READASARRAYS
  auxiliary  rchMult
  AUXMULTNAME  rchmult
END options

BEGIN period  1
  recharge
    CONSTANT  4.20779661E-04
  rchMult
    OPEN/CLOSE  'rchMult.dat'  FACTOR  1.0
END period  1

BEGIN period  2
  recharge
    CONSTANT  8.29733333E-04
  rchMult
    OPEN/CLOSE  'rchMult.dat'  FACTOR  1.0
END period  2

Any advice on why this is not working? I can't seem to figure out why PEST wont read this file and I'm a little confused by the naming conventions written by pf.add_parameters. This it what my add parameters looks like:

f= os.path.join('rchMult.dat')
pf.add_parameters(
    f,
    par_type="zone",
    par_name_base="rchMult.dat",
    pargp= 'rchMult.dat',
    initial_value= 0.1,
    zone_array=zones,
    upper_bound=0.2,
    lower_bound=0.001,
    transform= "log"
)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
1 participant
@iaped
0