Skip to content

Update int/sdirk.f90 to the same version as was in int/beuler.f90#130

Merged
yantosca merged 3 commits intodevfrom
feature/update-sdirk
Feb 19, 2025
Merged

Update int/sdirk.f90 to the same version as was in int/beuler.f90#130
yantosca merged 3 commits intodevfrom
feature/update-sdirk

Conversation

@yantosca
Copy link
Copy Markdown
Contributor

@yantosca yantosca commented Feb 18, 2025

This is the companion PR to #129. In this PR we do the following:

  • Moved int/beuler.f90 to int/sdirk.f90, since the code in int/beuler.f90 was a newer version of the SDIRK integrator module
  • Removed int/beuler.def
  • Updated documentation to instruct users wishing to select the Backwards Euler method to use #INTEGRATOR sdirk and set ICNTRL(3) = 6
  • Only use 101*Roundoff in the computation for Hstart if the selected integration method is Backward Euler
  • Trimmed trailing whitespace, indented code

int/sdirk.f90
- Original version has been deleted
- int/beuler.f90 is now renamed to int/sdirk.f90

int/beuler.f90
int/beuler.def
- Removed

docs/source/tech_info/07_numerical_methods.rst
- Added documentation instructing user to select #INTEGRATOR sdirk
  with ICNTRL(3) = 6 for Backwards Euler integration

CHANGELOG.md
- Updated accordingly

Signed-off-by: Bob Yantosca <[email protected]>
int/sdirk.f90
- Indented & lined up code for readability
- Added () to subroutines taking no arguments
- Test if the method is Backward Euler (BEL) before multiplying
  the RoundOff variable by 101 when computing Hstart
- Trimmed trailing whitespace

Signed-off-by: Bob Yantosca <[email protected]>
@yantosca yantosca added feature New feature or request integrators Related to numerical integrators labels Feb 18, 2025
@yantosca yantosca added this to the 3.2.0 milestone Feb 18, 2025
@yantosca yantosca requested a review from RolfSander February 18, 2025 23:33
@yantosca yantosca self-assigned this Feb 18, 2025
Copy link
Copy Markdown
Contributor

@RolfSander RolfSander left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good.

@yantosca yantosca merged commit 3bfd682 into dev Feb 19, 2025
1 check passed
@yantosca yantosca deleted the feature/update-sdirk branch February 19, 2025 15:25
@yantosca yantosca linked an issue Feb 19, 2025 that may be closed by this pull request
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

feature New feature or request integrators Related to numerical integrators

Projects

None yet

Development

Successfully merging this pull request may close these issues.

[BUG/ISSUE] Potential bug in beuler.f90

2 participants