trontrytel commented 7 months ago

The Climate Modeling Alliance

Software Design Issue 📜

Purpose

Prescribe the aerosol concentrations based on RCP 4.5 from CMIP 5. The files are here: https://tntcat.iiasa.ac.at/RcpDb/dsd?Action=htmlpage&page=download

Benefits

We need aerosol optical properties for radiation. (We also need to code a way to compute optical properties from aerosol concentrations, see SDI)
We need aerosol precursors for cloud droplet activation and ice nucleation in 2-moment microphysics scheme. (The 2-moment scheme is under development. See SDI and SDI)

People

Lead: @trontrytel
Collaborators: @LenkaNovak @szy21
Reviewers:

Results

Run AMIP simulations with aerosol concentrations read from files.

Tasks for prescribing data from ClimaCoupler (this was abandoned for the solution below)

[x] Visualize the RCP 4.5 data
[x] Store the data on Caltech box
https://github.com/CliMA/ClimaCoupler.jl/pull/506
~Convert (lat, lon, p, t) fields to (lat, lon, z, t) fields~
~Run spherical simulations with prescribed aerosols (no coupling between aerosol and anything) and test that our output diagnostics reproduce the prescribed files~

Tasks for prescribing data from ClimaAtmos

- [x] [Aerosol prescribing enabled in Atmos](https://github.com/CliMA/ClimaAtmos.jl/pull/2815)
- [ ] [add aerosol optics in rrtmgp](https://github.com/CliMA/RRTMGP.jl/issues/460)
- [ ] use aerosols in from the coupler

SDI Revision Log

CC

@tapios @simonbyrne @cmbengue

trontrytel commented 7 months ago

The RCP 4.5 data is in lon(0, 357.5), lat(-90, 90), sigma pressure (3.5, 992.6), time(56605, 56939) format.

The dataset does not include surface height. So in order to convert it to height I would have to look for a corresponding topography file. I think it would be easier to just try interpolating along pressure inside the Coupler?

Here is a couple of plots from the lowest level (fist column Jan-Jun, second column Jul-Dec)

Temperature: