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

Code additions for alternate g-point cloud optics #300

Conversation

mjiacono
Copy link
Collaborator

This PR includes code revisions to input and use alternate g-point cloud optics input data. The all-sky example test case was revised to demonstrate the new capability.

Longwave and shortwave flux profile comparisons (between band and g-point clouds optics) were completed for the tropical and mid-latitude winter atmospheric profiles using two cloud cases, one an all liquid low cloud and one an all ice mid-level cloud. The graphics below illustrate impacts of up to one percent in the shortwave liquid cloud case, with smaller impacts in the longwave and in the ice cloud case.

flxupdn_lwsw_diff_cld_bnd_gpt_mlw
flxupdn_lwsw_diff_cld_bnd_gpt_trp

A comparison was also completed between shortwave flux calculations using the RRTMGP band and g-point cloud optics and the LBLRTM shortwave cloudy flux calculations of Lu et al. (2011; https://doi.org/10.1175/JAS-D-10-05001.1). The LBLRTM fluxes from Table 3 of that paper and the comparable results from RRTMGP are shown in the table below. This comparison is not optimized, since small differences in TSI, cloud physical properties, gaseous absorption coefficients, or other parameters may be present.

LBL_comp_band_g-point_cloud_optics

Iacono and others added 5 commits July 26, 2023 18:23
… API (earth-system-radiation#284)

Accumulated changes and bug fixes. Changes LW answers (data repo is alsoupdated).
- New LW weights and secants (Change LW quadrature angles earth-system-radiation#282)
- A single source function on levels (Simplify LW source functions earth-system-radiation#250)
- Kernel API header files in Fortran and C (Add kernel API earth-system-radiation#272)
- Refactored two-stream, fixing RTE shortwave kernel not vectorizing earth-system-radiation#215 (Re-vectorize SW two-stream earth-system-radiation#275)
- Buxfixes, change to internal build systems and continuous integration (thanks to @skosukhin for ongoing help)
@RobertPincus
Copy link
Member

@mjiacono Two requests.

  1. Please open a matching PR that adds the new cloud absorption files to the rrtmgp-data repo.
  2. I propose that, rather than changing any existing data files, you add a new file with the cloud properties by g-point. It should have the same format as the original file (the same dataset names) and differ only in having tables defined on g-points rather than on bands. This will keep code changes to a minimum.

@RobertPincus
Copy link
Member

See also earth-system-radiation/rrtmgp-data#7

@RobertPincus RobertPincus changed the base branch from main to hotfix-relabel-ice-sizes December 23, 2024 21:17
@RobertPincus RobertPincus deleted the branch earth-system-radiation:hotfix-relabel-ice-sizes January 14, 2025 23:13
@mjiacono mjiacono deleted the feature-cloud-optics-gpoint branch February 3, 2025 22:52
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.

2 participants