Skip to content
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.

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

Add TUV-x height grid updates #132

Closed

Conversation

mattldawson
Copy link
Collaborator

@mattldawson mattldawson commented Oct 9, 2024

Updates TUV-x height grid with CAM-SIMA state data before calculating rate constants.

Uses geopontential height values to calculate actual heights above sea level, as currently done in the CAMChem gas-phase chemistry module. If there is a better way to get these height values, we can update this PR.

closes #94

@mattldawson mattldawson changed the title add TUV-x height grid updates and run function add TUV-x height grid updates Oct 11, 2024
@mattldawson mattldawson marked this pull request as ready for review October 11, 2024 19:16
@boulderdaze boulderdaze changed the base branch from main to development October 11, 2024 22:37
@boulderdaze boulderdaze changed the title add TUV-x height grid updates Add TUV-x height grid updates Oct 15, 2024
Comment on lines +46 to +50
height_grid => create_height_grid( vertical_layer_dimension, &
vertical_interface_dimension, errmsg, errcode )
if (errcode /= 0) return
call grids%add( height_grid, error )
if (has_error_occurred( error, errmsg, errcode )) return
Copy link
Collaborator

Choose a reason for hiding this comment

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

I think the error checking function call can be consistent unless there is reason to not to?

Suggested change
height_grid => create_height_grid( vertical_layer_dimension, &
vertical_interface_dimension, errmsg, errcode )
if (errcode /= 0) return
call grids%add( height_grid, error )
if (has_error_occurred( error, errmsg, errcode )) return
height_grid => create_height_grid( vertical_layer_dimension, &
vertical_interface_dimension, errmsg, errcode )
if (has_error_occurred( error, errmsg, errcode )) return
call grids%add( height_grid, error )
if (has_error_occurred( error, errmsg, errcode )) return


type(grid_t), intent(inout) :: height_grid
real(kind_phys), intent(in) :: host_midpoints(:) ! km
real(kind_phys), intent(in) :: host_edges(:) ! km
Copy link
Collaborator

Choose a reason for hiding this comment

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

Is the term edge equivalent to interface? If so, can we use interface across the code base to enhance the consistency?

Suggested change
real(kind_phys), intent(in) :: host_edges(:) ! km
real(kind_phys), intent(in) :: host_interface(:) ! km

@boulderdaze
Copy link
Collaborator

Looks like I can't commit the code to Matt's fork, which is expected. I am going to create a new branch in this repository for this PR so we can push the commits

@boulderdaze
Copy link
Collaborator

Opened the same PR. Closes

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Update height grid in TUV-x before calculating rate constants
2 participants