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

fv_diagnostics update #332

Open
wants to merge 6 commits into
base: main
Choose a base branch
from
Open

fv_diagnostics update #332

wants to merge 6 commits into from

Conversation

tmerlis
Copy link

@tmerlis tmerlis commented Apr 2, 2024

Description

fv_diagnostics.F90 and fv_diagnostics.h are updated to
i) include new diagnostics related to height (vertically integrated) and its horizontal fluxes (both 3d and vertically integrate), as these are components of the moist static energy
ii) correct the logic of existing vertically integrated horizontal T, q fluxes (previously, having just those in a diag table without the 3d version wouldn't work)

Fixes # (issue)

How Has This Been Tested?

The revision compiles and runs on on Princeton's Stellar platform with Intel Fortran Compiler.

There is now output if the diag table only has vertically integrated q and T fluxes.

Checklist:

Please check all whether they apply or not

  • My code follows the style guidelines of this project
  • I have performed a self-review of my own code
  • I have commented my code, particularly in hard-to-understand areas
  • I have made corresponding changes to the documentation
  • My changes generate no new warnings
  • Any dependent changes have been merged and published in downstream modules

Copy link
Contributor

@lharris4 lharris4 left a comment

Choose a reason for hiding this comment

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

I think these are fine. The product variables are very useful for computing eddy fluxes and are a neat addition. I have a question and one requested fix.

@abrooks1085 @bensonr are there any steps that need to be done to get the CI working properly here?

Thanks,
Lucas

@@ -1140,7 +1150,10 @@ subroutine fv_diag_init(Atm, axes, Time, npx, npy, npz, p_ref)
'Vertically Integrated Snow', 'kg/m**2', missing_value=missing_value )
id_intqg = register_diag_field ( trim(field), 'intqg', axes(1:2), Time, &
'Vertically Integrated Graupel', 'kg/m**2', missing_value=missing_value )

id_inthght = register_diag_field ( trim(field), 'inthght', axes(1:2), Time, &
Copy link
Contributor

Choose a reason for hiding this comment

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

Are inthght and inttemp supposed to be proportional to the vertically-integrated potential and internal energies?

Copy link
Author

Choose a reason for hiding this comment

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

Yes, that's correct. This could be overt by multiplying by the relevant constants 'grav' and 'cp_air' that are available in the module and changing the unit accordingly.

tools/fv_diagnostics.F90 Outdated Show resolved Hide resolved
Copy link
Contributor

@lharris4 lharris4 left a comment

Choose a reason for hiding this comment

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

@tmerlis thanks for taking care of these. Looks good to me.

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.

3 participants