O1.5.4a Prescribe aerosol concentrations

[trontrytel]

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: @szy21
• Collaborators: @Sbozzolo @sriharshakandala @trontrytel
• Reviewers:

Results

Run AMIP simulations with aerosol concentrations read from files.

Tasks

Preview Give feedback

1. Add aerosol data to artifacts
2. Aerosol prescribing enabled in Atmos
3. add aerosol optics in rrtmgp
4. Add aerosol radiative effects #3140
   4 of 4
   radiation ☀️ +1
   
5. Add support for time interpolation with jump ClimaUtilities.jl#74
   +0
   
6. Use multiple aerosols for radiation #3274
   +0
   
7. use aerosols in the AMIP simulation
8. Update aerosol SW lookup table RRTMGP.jl#539
   bug +1
   

SDI Revision Log

CC

@tapios @simonbyrne @cmbengue

[trontrytel]

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:

Seasalt:

Dust:

[trontrytel]

The data is now on CliMA box under Aerosol_RCP45 folder

[trontrytel]

Remember to take a look at

ClimaAtmos.jl/src/parameterized_tendencies/radiation/radiation.jl

Lines 85 to 97 in b26b78b

	# interpolate the ozone concentrations to our initial pressures
	pressure2ozone = Spline1D(
	input_center_pressure,
	input_center_volume_mixing_ratio_o3,
	)
	if device isa ClimaComms.CUDADevice
	fv = Fields.field_values(ᶜp)
	fld_array = DA(pressure2ozone.(Array(parent(fv))))
	data = DataLayouts.rebuild(fv, fld_array)
	ᶜvolume_mixing_ratio_o3_field = Fields.Field(data, axes(ᶜp))
	else
	ᶜvolume_mixing_ratio_o3_field = @. FT(pressure2ozone(ᶜp))
	end

[Sbozzolo]

Is there anything else needed to close this issue?

[szy21]

Is there anything else needed to close this issue?

Sea salt is not added, but let's close this. I will add sea salt after updating rrtmgp.