CTSM icon indicating copy to clipboard operation
CTSM copied to clipboard

Receive ozone from atmosphere

Open billsacks opened this issue 6 years ago • 27 comments

(Moved from https://trello.com/c/fHXntN8l/150-ozone-receive-from-atmosphere)

  • [x] add ozone stream to datm (Danica had diffs in datm_comp_mod in SourceMods, as well as some things she did for adding a new stream for her cases)

  • [x] add coupling hooks; to do this, incorporate diffs from clm_atmlnd.F90, clm_cpl_indices.F90, lnd_comp_mct.F90 in Danica's original branch (o3_amax_pft); also add coupling field seq_flds_mod (Danica had diffs in SourceMods)

  • [x] See bug report 2294 (now #134)

  • [x] determine how coupling to CAM will happen (see wjs comments below - 10-11-2016 and 12-27-16)

    • (2022-03-25) This is done: see https://github.com/ESCOMP/CAM/pull/404 and https://github.com/ESCOMP/CMEPS/pull/212 and https://github.com/ESCOMP/CESM_share/pull/4. Briefly: the atm component will always send some ozone concentration. It will set a flag saying whether it's providing subdaily ozone values or a multi-day average. If the latter, CTSM will apply anomalies on top of the provided values, which CTSM will read from an anomaly file.
  • [x] add a forc_ozone history field (O3CONC_IN): see histFldsMod in Danica's original branch (o3_amax_pft)

  • [x] make ozone a required field from atm.

    • (2022-03-25) This is done: see https://github.com/ESCOMP/CMEPS/pull/212 and https://github.com/ESCOMP/CAM/pull/404
  • [x] Work with Francis for getting from CAM

    • (2022-03-25) This is done: see https://github.com/ESCOMP/CAM/pull/404
  • [ ] add ozone as a cpl hist field (written by coupler, and then read by datm when running in cplhist mode)

  • [x] Read and apply diurnal ozone anomalies if the atmosphere forcing is smoothed

billsacks avatar Feb 12 '18 23:02 billsacks

Original comment from @billsacks 2016-10-11

For coupling to CAM: It seems like ozone will only be available when running with chemistry. I can see a few options here:

Have CAM (and datm) set a flag like atm_provides_ozone in the coupler in initialization (in seq_infodata?). Then, if ozone is on in CLM, then CLM checks whether this flag is set; if not, aborts. (This would depend on atm init happening before lnd init, I think.) This would be simple if it works, but probably isn't the most robust solution.

Can we use the drv_flds_in functionality? CAM can add ozone to its list of fields if running with chemistry; datm can always add ozone to its list of fields. Then CLM can abort if use_ozone is true but the ozone field index isn't set.

Have ozone controlled via an xml variable. Set USE_OZONE=TRUE to turn on ozone; this (a) turns on the CLM use_ozone namelist variable; and (b) tells CAM/DATM to send ozone. CAM should check at build-namelist time to confirm that USE_OZONE is only true if it is running with chemistry.

@ekluzek thinks that it makes sense to use drv_flds_in (i.e., option 2). CAM will add the ozone field if it's running with a chemistry configuration that supports it. When running CLM offline, we'll either (a) have CLM add ozone to drv_flds_in when use_ozone is .true., or (b) have datm add ozone to drv_flds_in always. Then CLM should check at runtime: if use_ozone is true, but the ozone field has index = 0, then abort.

But upon further reflection, I think it makes sense to require CAM to always send some ozone concentration, even when it's running without chemistry. If it's running without chemistry, it could just send a global, constant value for now (as has been the case in CLM up until now). This seems important so that you can at least run some flavor of ozone in CLM when running without CAM-chem. That would also simplify the coupling logic - because we wouldn't need to anything with drv_flds_in, etc.

(12-5-16) In talking with Danica: She'd prefer if you got an error message if you're running CAM in a configuration without ozone. But she'd be okay with CAM sending a constant - she suggests it just sends 0 (an alternative would be for it to send some garbage value that's obviously bad).

billsacks avatar Feb 12 '18 23:02 billsacks

Original comment from @billsacks 2016-10-11

Currently Danica is using an hourly "mask" (8760 values on the file - one for each hour of the year) that is added to the monthly values. She does the addition in CLM. If we wanted to use that functionality, then we'd want to:

  1. Move the addition to datm, rather than doing the addition in CLM

  2. Probably we'd want to turn off the ozone read unless it's needed (e.g., see the idea of a USE_OZONE xml variable in my recent comment), because reading those hourly data could have significant performance overhead.

For the first pass, we'll probably skip the hourly "mask", just supporting monthly-average ozone values.

billsacks avatar Feb 12 '18 23:02 billsacks

Original comment from @billsacks 2016-11-30

Mariana has done the necessary xml and buildnml mods for datm on this branch:

https://github.com/ESMCI/cime/tree/mvertens/datm_ozone

billsacks avatar Feb 12 '18 23:02 billsacks

Original comment from @billsacks 2016-12-27

Update on the hourly "mask", from Danica: It turns out that CAM, when run without CHEM, can supply monthly average ozone values. So we'll probably want to apply the "mask" in CLM after all, rather than in datm - since this "mask" will be applied both when running with datm and when running with CAM without CHEM. But in order for this to work right, CLM will need to know when it's being given monthly values, and when it's being given hourly / half-hourly data already. Maybe have datm and CAM set a field in seq_infodata to communicate this to CLM?

billsacks avatar Feb 12 '18 23:02 billsacks

cc @danicalombardozzi

billsacks avatar Feb 12 '18 23:02 billsacks

That sounds good to me if it's a viable option. I should probably recalculate the mask after the new CAM-CHEM simulations are complete so that we are using the most up-to-date chemistry.

danicalombardozzi avatar Feb 12 '18 23:02 danicalombardozzi

@danicalombardozzi to be honest: I just moved this over blindly from trello without even re-reading some of this. Hopefully once the never-ending release crunch is done, we can finally get this in!

billsacks avatar Feb 13 '18 00:02 billsacks

~See also #134 , which may be a duplicate of this.~ Edit: #134 is a piece of this, not a duplicate.

billsacks avatar Jan 14 '20 16:01 billsacks

Notes from @danicalombardozzi's email on updating this for CESM3.

CAM folks have put the pieces in place to send ozone from CAM to the coupler. We need to connect the pieces in CLM. The outstanding tasks include:

  1. Connect O3 from the coupler into CLM (for coupled simulations)
  2. Add a new forcing stream to datm (for land-only simulations)
  3. Generate files for hourly ozone concentrations (needed for both datm and coupled simulations without chemistry) (Louisa)

Please note that @lkemmons has agreed to help with generating the appropriate files for 3, and @mvertens still help with 1 and 2. There are two university groups who are actively trying to do this for their own needs, and incorporating this capability into CLM will save time and energy for future users.

wwieder avatar Apr 04 '22 15:04 wwieder

@billsacks, @ekluzek and @lkemmons, is there more from CTSM that needs to happen for integration with CAM-chem? or is this issue done effectively (aside from documentation)?

wwieder avatar Nov 15 '22 13:11 wwieder

Most of this is done thanks to work by @adrifoster a few months ago. Everything should work correctly now when running with CAM-chem, where we have prognostic ozone values every time step. The main remaining piece – which @adrifoster started but hasn't completed – is applying diurnal anomalies to smoothed ozone values to get better ozone forcings when running with either DATM or CAM without chemistry.

billsacks avatar Nov 15 '22 13:11 billsacks

Seems like it may be time to revisit this. How much of a priority is this for CESM3?

wwieder avatar Jan 24 '24 20:01 wwieder

It would be nice to complete this - I know of users wanting it.

lkemmons avatar Feb 20 '24 02:02 lkemmons

Thanks, @lkemmons. I'm not really clear what it would take to get this wrapped up for CESM3, we can discuss at an upcoming CTSM software meeting.

wwieder avatar Feb 20 '24 15:02 wwieder

From the CAM-Chem group, how does prioritization stack up compared to the soil NOx work, #2290?

wwieder avatar Feb 20 '24 16:02 wwieder

Soil NO emissions are more important for CAM-chem, in my opinion.

On Tue, Feb 20, 2024 at 9:33 AM will wieder @.***> wrote:

From the CAM-Chem group, how does prioritization stack up compared to the soil NOx work, #2290 https://github.com/ESCOMP/CTSM/pull/2290?

— Reply to this email directly, view it on GitHub https://github.com/ESCOMP/CTSM/issues/270#issuecomment-1954599416, or unsubscribe https://github.com/notifications/unsubscribe-auth/AH5BH7LEYFAJCZBPBIHHTU3YUTF4ZAVCNFSM4EQL5IEKU5DIOJSWCZC7NNSXTN2JONZXKZKDN5WW2ZLOOQ5TCOJVGQ2TSOJUGE3A . You are receiving this because you were mentioned.Message ID: @.***>

lkemmons avatar Feb 20 '24 17:02 lkemmons

We think this only matters for offline I cases at this point so we don't think this matters for CAM-Chem.

But, @wwieder will setup a meeting to discuss with @lkemmons and @adrifoster to discuss.

ekluzek avatar Feb 22 '24 17:02 ekluzek

We think this only matters for offline I cases at this point so we don't think this matters for CAM-Chem.

That matches my recollection, too.

billsacks avatar Feb 22 '24 17:02 billsacks

That is correct - the connection to CAM is complete. The last thing that needs to be completed is temporal downscaling of the ozone streams for offline cases.

On Thu, Feb 22, 2024 at 10:46 AM Bill Sacks @.***> wrote:

We think this only matters for offline I cases at this point so we don't think this matters for CAM-Chem.

That matches my recollection, too.

— Reply to this email directly, view it on GitHub https://github.com/ESCOMP/CTSM/issues/270#issuecomment-1959954572, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADE42IRK74IYASMVRALCBADYU6AARAVCNFSM4EQL5IEKU5DIOJSWCZC7NNSXTN2JONZXKZKDN5WW2ZLOOQ5TCOJVHE4TKNBVG4ZA . You are receiving this because you were mentioned.Message ID: @.***>

adrifoster avatar Feb 22 '24 18:02 adrifoster

Is the CAM connection on by default for B-cases?

The notes above state that the temporal downscaling of ozone streams for I-cases is something that @lkemmons will help with, although that was 2 years ago! Is this something you're still interested / willing to do? I'd assume this would be something we'd want to generate from cplhist files of CESM3 simulations. Maybe we need to make sure this output is generated in cplhist output next time we generate these data?

wwieder avatar Feb 22 '24 18:02 wwieder

Presumably this is not on by default in B or F-cases, because that would mean that we are doing ozone damage by default. Or, maybe my question is whether or not turning on ozone damage on land would 'automatically' trigger passing the ozone fields from the atmosphere in coupled cases or whether you need to turn on both ozone damage and the passing of ozone from atmosphere to land.

On Thu, Feb 22, 2024 at 11:15 AM will wieder @.***> wrote:

Is the CAM connection on by default for B-cases?

The notes above state that the temporal downscaling of ozone streams for I-cases is something that @lkemmons https://github.com/lkemmons will help with, although that was 2 years ago! Is this something you're still interested / willing to do? I'd assume this would be something we'd want to generate from cplhist files of CESM3 simulations. Maybe we need to make sure this output is generated in cplhist output next time we generate these data?

— Reply to this email directly, view it on GitHub https://github.com/ESCOMP/CTSM/issues/270#issuecomment-1960003688, or unsubscribe https://github.com/notifications/unsubscribe-auth/AFABYVFRLIHVI3PJO55MOZTYU6DKXAVCNFSM4EQL5IEKU5DIOJSWCZC7NNSXTN2JONZXKZKDN5WW2ZLOOQ5TCOJWGAYDAMZWHA4A . You are receiving this because you are subscribed to this thread.Message ID: @.***>

dlawrenncar avatar Feb 22 '24 18:02 dlawrenncar

I think it is only on when you use the ozone-specific namelist options for any case, but I can check on this.

@wwieder I think @lkemmons already helped by providing the downscaling data, so we/I just need to put it in CTSM. I started this a long time ago so I think it could be fairly straightforward to finish.

adrifoster avatar Feb 22 '24 18:02 adrifoster

Oh, that would be great if this is relatively straightforward to wrap up, @adrifoster.

wwieder avatar Feb 22 '24 18:02 wwieder

I even wrote a design doc about it!

adrifoster avatar Feb 22 '24 18:02 adrifoster

I didn't look carefully at what I was responding to. I am happy to help provide updated ozone fields if necessary, or provide instructions if this is something that should be produced routinely.

On Thu, Feb 22, 2024 at 11:45 AM Adrianna Foster @.***> wrote:

I even wrote a design doc https://github.com/ESCOMP/CTSM/blob/06a62712a5f3784921f0981181e3b5037cfbdb8e/doc/design/design_doc_diurnal_ozone.md about it!

— Reply to this email directly, view it on GitHub https://github.com/ESCOMP/CTSM/issues/270#issuecomment-1960047605, or unsubscribe https://github.com/notifications/unsubscribe-auth/AH5BH7JF547NXUI4YUXZRSDYU6G25AVCNFSM4EQL5IEKU5DIOJSWCZC7NNSXTN2JONZXKZKDN5WW2ZLOOQ5TCOJWGAYDINZWGA2Q . You are receiving this because you were mentioned.Message ID: @.***>

lkemmons avatar Feb 22 '24 20:02 lkemmons

We think this only matters for offline I cases at this point so we don't think this matters for CAM-Chem.

That matches my recollection, too.

Replying to my own comment as well as some of the other recent comments: Actually, I now think / remember that it's more subtle than this: I think that, in F/B cases with CAM-Chem – so with prognostic ozone – the work @adrifoster did before she went on maternity leave put everything needed in place. I think the coupling is already set up so that ozone is always being passed from CAM to CTSM, so I think all that's needed in those cases is to turn on the ozone damage namelist option. (Though I'm not sure if we ever got to the bottom of https://github.com/ESCOMP/CTSM/issues/962.) But the application of diurnal anomalies is needed both for offline (I) compsets and for F/B compsets without chemistry, in which case CAM passes smoothed ozone concentrations without diurnal variability.

billsacks avatar Feb 22 '24 22:02 billsacks

@billsacks Do we still want to do this?

  • [ ] add ozone as a cpl hist field (written by coupler, and then read by datm when running in cplhist mode)

I believe it's the last thing to do in the list at the top of this issue. I'm not sure what is involved in this.

also are we still worried about #962 above?

adrifoster avatar Apr 26 '24 21:04 adrifoster