8000 RNITPD division by zero in PODS.for · Issue #435 · DSSAT/dssat-csm-os · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

RNITPD division by zero in PODS.for #435

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
fabiooliveira72 opened this issue Aug 22, 2024 · 0 comments
Open

RNITPD division by zero in PODS.for #435

fabiooliveira72 opened this issue Aug 22, 2024 · 0 comments

Comments

@fabiooliveira72
Copy link
Contributor

Error of division by zero in PODS, equation:
RNITPD = (RNITP*0.01-PROLFF*0.16)/(FNINL-PROLFF*0.16)
Reason:
The division by zero happens because FNINL is assigned in DEMAND.for with:
FNINL = PROLFI * 0.16
This value is not changed during the simulation and will be passed to PODS subroutine. PROLFI and PROLFF comes from the Species file. The values of this variables were set equal as showed below:

!*PLANT COMPOSITION VALUES
 .200  .200  .200  .028  .028  .028   PROLFI,PROLFG,PROLFF,PROSTI,PROSTG,PROSTF

The division by zero happens because the FNINL - PROLFF*0.16 will result in 0.0. To bypass this problem the value of PROLFI was increased in 0.001, because it is usually higher than PROLFF based on other CROPGRO models. However, this issue is also present in all CROPGRO-based models and should be added a protection in PODS.for.
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
@fabiooliveira72
0