Addition of Water Vapor Age Distribution Moments Tracers

[PhilippeBoulangerPhysics]

We have developed passive tracers for calculating the water vapor age spectrum in an Eulerian model. Our method resolves the temporal and spatial variability of the water vapor age distribution, and allows us to calculate the age spectrum of the precipitating water as well. The theory behind our method and the results of our climate experiments will be submitted to EGU's Weather and Climate Dynamics in the near future.

Our method was implemented by adding passive tracers to track the moments of the water vapor age distribution.

The age_moments module updates these tracers, and the tracers can be turned on/off at will. (Default is off)

Passed trip tests (axisymmetric, bucket_model, frierson, giant_planet, realistic_continents_fixed_sst, realistic_continents_variable_qflux, variable_co2_grey, variable_co2_rrtm).

Philippe Boulanger and Robert Fajber
Department of Atmospheric and Oceanic Sciences, McGill University, Montreal, Quebec Canada
Brace Water Center, McGill University, Ste Anne de Bellevue, Quebec Canada

[sit23]

Dear Philippe - thanks very much for submitting this P/R. There's a lot of really interesting extensions to Isca here, which is great. I don't quite have time to do a full review for now, but just a few thoughts after a quick scan:

• Please could you remove any references to your username or your local directory from the python scripts you've provided? We try and keep these as agnostic to the individual user as possible.
• Please could you remove the .sh runscripts from the experiment folders? These might be necessary at a local level, but this could be confusing for users to find a python runscript and a .sh runscript to run the python.
• I'm very grateful that you've put in lots of land masks, but Adding a script to create smoothed ERA5 topography files #288 has only just merged a script to do this. So I'd prefer if those weren't present.
• I see you've modified some of the test cases to include your new tracer functionality. What I'd prefer is for the test cases themselves to stay as they are, and for you to include a seperate new test that specifically shows off your functionality. This could be where you have e.g. a MiMA test case and a MiMA + tracer test case. That way we can tell very specifically what's different / added.

I hope this is all OK. Once you've made those changes I'll look into the rest of the code in more detail. Thanks again, Stephen

[PhilippeBoulangerPhysics]

Hi,

I have made the changes you mentioned. I cleaned up any references to my username or my local directory and I added the test case you mentioned.

Furthermore, in my original pull request, I have failed to mention that I added a 6th surface albedo option to modify the albedo over a specific rectangular region centered at point (albedo_cntr_lat,albedo_cntr_lon) extending to +/- albedo_wdth_lat and +/- albedo_wdth_lon.

I hope the code is now OK.

Best,

Philippe