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

Separate removal of Antarctic river and ice runoff #103

Conversation

xylar
Copy link
Collaborator

@xylar xylar commented Jun 23, 2024

This is needed so so we can remove ice (solid) runoff but not river (liquid) runoff.

This merge also removes totalRemovedRiverRunoffFlux and totalRemovedIceRunoffFlux variables, as these were computed incorrectly (they were missing areaCell) and they are not currently used for anything within MPAS-Ocean or any diagnostics.

@@ -707,10 +707,11 @@ add_default($nl, 'config_flux_attenuation_coefficient_runoff');
# Namelist group: coupling #
############################

add_default($nl, 'config_remove_ais_river_runoff');
Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is a non-BFB change for simulations with DIB. @darincomeau and @cbegeman, do you agree that we no longer want to remove river runoff? Or should we make this change a separate PR?

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think we mostly likely want to make this change, but I'd like to see a 30-year B-case before deciding to make this the default.

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think that will be a lot to be able to afford for a SORRM run just to test this config option, but we can discuss it.

@xylar
Copy link
Collaborator Author

xylar commented Jun 23, 2024

@cbegeman and @darincomeau, I would appreciate your review of this work. It's a step needed before https://acme-climate.atlassian.net/wiki/spaces/PSC/pages/4210098268/Design+Document+Data+iceberg+and+ice-shelf+melt+flux+patterns+for+E3SM+spin-up+runs but it's not directly related to that work and seems like something we need for any simulations with DIB and either DISMF or PISMF.

@xylar
Copy link
Collaborator Author

xylar commented Jun 23, 2024

I haven't had a chance to test these changes yet, and probably won't have a chance to do that until Thursday so for now just a look at the code is what I'm after.

@xylar xylar force-pushed the ocn/separate-removial-ais-river-ice-runoff branch from 02a78c5 to c06108d Compare June 23, 2024 20:54
The totals are computed incorrectly (they are missing areaCell)
and are not currently used for anything within MPAS-Ocean or any
diagnostics.
This is needed so so we can remove ice (solid) runoff but not
river (liquid) runoff.
@xylar xylar force-pushed the ocn/separate-removial-ais-river-ice-runoff branch from c06108d to 5cb37fb Compare June 23, 2024 21:09
@xylar
Copy link
Collaborator Author

xylar commented Jun 24, 2024

I did test this after all with the e3sm_cryo_developer suite. I'm seeing segfaults in time_series_stats:

101: ==== backtrace (tid: 666453) ====
101:  0 0x0000000000012cf0 __funlockfile()  :0
101:  1 0x0000000005c62832 ocn_time_series_stats_mp_operate3r_avg_()  /lcrc/group/e3sm/ac.xylar/scratch/chrys/SMS_D_Ld1.TL319_IcoswISC30E3r5.GMPAS-JRA1p5-DIB-PISMF.chrysalis_intel.mpaso-jra_1958.20240623_155730_j02lqb/bld/cmake-bld/core_ocean/analysis_members/mpas_ocn_time_series_stats.f90:2542
101:  2 0x0000000005c50c2c ocn_time_series_stats_mp_typed_operate_()  /lcrc/group/e3sm/ac.xylar/scratch/chrys/SMS_D_Ld1.TL319_IcoswISC30E3r5.GMPAS-JRA1p5-DIB-PISMF.chrysalis_intel.mpaso-jra_1958.20240623_155730_j02lqb/bld/cmake-bld/core_ocean/analysis_members/mpas_ocn_time_series_stats.f90:2079
101:  3 0x0000000005c2654b ocn_time_series_stats_mp_ocn_compute_time_series_stats_()  /lcrc/group/e3sm/ac.xylar/scratch/chrys/SMS_D_Ld1.TL319_IcoswISC30E3r5.GMPAS-JRA1p5-DIB-PISMF.chrysalis_intel.mpaso-jra_1958.20240623_155730_j02lqb/bld/cmake-bld/core_ocean/analysis_members/mpas_ocn_time_series_stats.f90:383
101:  4 0x0000000005442aff ocn_analysis_driver_mp_ocn_compute_analysis_members_()  /lcrc/group/e3sm/ac.xylar/scratch/chrys/SMS_D_Ld1.TL319_IcoswISC30E3r5.GMPAS-JRA1p5-DIB-PISMF.chrysalis_intel.mpaso-jra_1958.20240623_155730_j02lqb/bld/cmake-bld/core_ocean/analysis_members/mpas_ocn_analysis_driver.f90:1394
101:  5 0x0000000005435f84 ocn_analysis_driver_mp_ocn_analysis_compute_()  /lcrc/group/e3sm/ac.xylar/scratch/chrys/SMS_D_Ld1.TL319_IcoswISC30E3r5.GMPAS-JRA1p5-DIB-PISMF.chrysalis_intel.mpaso-jra_1958.20240623_155730_j02lqb/bld/cmake-bld/core_ocean/analysis_members/mpas_ocn_analysis_driver.f90:831
101:  6 0x00000000031d4883 ocn_comp_mct_mp_ocn_run_mct_()  /lcrc/group/e3sm/ac.xylar/scratch/chrys/SMS_D_Ld1.TL319_IcoswISC30E3r5.GMPAS-JRA1p5-DIB-PISMF.chrysalis_intel.mpaso-jra_1958.20240623_155730_j02lqb/mpas-ocean/driver/ocn_comp_mct.f90:1053
101:  7 0x000000000049b235 component_mod_mp_component_run_()  /gpfs/fs1/home/ac.xylar/e3sm_work/E3SM/ocn/separate-removial-ais-river-ice-runoff/driver-mct/main/component_mod.F90:757
101:  8 0x0000000000453f36 cime_comp_mod_mp_cime_run_()  /gpfs/fs1/home/ac.xylar/e3sm_work/E3SM/ocn/separate-removial-ais-river-ice-runoff/driver-mct/main/cime_comp_mod.F90:3151
101:  9 0x0000000000482965 MAIN__()  /gpfs/fs1/home/ac.xylar/e3sm_work/E3SM/ocn/separate-removial-ais-river-ice-runoff/driver-mct/main/cime_driver.F90:153
101: 10 0x000000000041ab22 main()  ???:0
101: 11 0x000000000003ad85 __libc_start_main()  ???:0
101: 12 0x000000000041aa2e _start()  ???:0

So that obviously needs to be fixed.

@xylar
Copy link
Collaborator Author

xylar commented Jun 24, 2024

I reran SMS_D_Ld1.TL319_IcoswISC30E3r5.GMPAS-JRA1p5-DIB-PISMF.chrysalis_intel.mpaso-jra_1958 on its own and it was successful, so it doesn't seem like the failure above is reproducible.

I will run a test that actually includes removing ice runoff.

@cbegeman
Copy link
Collaborator

Should we run B-case tests both with and without @chloewhicker's snowcapping changes from E3SM-Project#6413? If I remember correctly, those significantly changed both liquid and solid runoff from AIS.

@cbegeman
Copy link
Collaborator

@xylar Let me know if you'd like me to run particular tests on different machines/compilers. The code looks good to me from visual inspection.

@xylar
Copy link
Collaborator Author

xylar commented Jun 24, 2024

Should we run B-case tests both with and without @chloewhicker's snowcapping changes from E3SM-Project#6413? If I remember correctly, those significantly changed both liquid and solid runoff from AIS.

I think that would be ideal but I don't think we will use that in any configuration other than SORRM. So I don't think it makes sense to do a long B-case run until we have a new version of the SORRM mesh and I'm not sure whether we think we can afford long B-case runs with a bunch of different configurations for comparison. We can certainly discuss this further. But, as an example, I've been sitting in the queue all day waiting to run a 1-hour test of SMS_D_Ld1.ne30pg2_SOwISC12to60E2r4.CRYO1850-DISMF.chrysalis_intel.

@cbegeman
Copy link
Collaborator

I think that would be ideal but I don't think we will use that in any configuration other than SORRM. So I don't think it makes sense to do a long B-case run until we have a new version of the SORRM mesh and I'm not sure whether we think we can afford long B-case runs with a bunch of different configurations for comparison. We can certainly discuss this further. But, as an example, I've been sitting in the queue all day waiting to run a 1-hour test of SMS_D_Ld1.ne30pg2_SOwISC12to60E2r4.CRYO1850-DISMF.chrysalis_intel.

@xylar Bummer. Thanks for the update.

@xylar
Copy link
Collaborator Author

xylar commented Jun 26, 2024

I ran a quick test that took 2 days in the queue to make sure river runoff was being removed and ice runoff was not. Results are in:

/lcrc/group/e3sm/ac.xasay-davis/scratch/chrys/SMS_D_Ld1.ne30pg2_SOwISC12to60E2r4.CRYO1850-DISMF.chrysalis_intel.20240624_034222_shla8m

As expected, I'm seeing:

$ ncdump SMS_D_Ld1.ne30pg2_SOwISC12to60E2r4.CRYO1850-DISMF.chrysalis_intel.20240624_034222_shla8m.mpaso.hist.am.conservationCheck.0001-01-01.nc

...
 accumulatedRemovedRiverRunoffFlux = 0, 0 ;

 accumulatedRemovedIceRunoffFlux = 0, 33623980.655436 ;

...

@xylar
Copy link
Collaborator Author

xylar commented Jun 26, 2024

I agree, we probably want this change only after E3SM-Project#6413 goes in.

@xylar xylar closed this Jul 1, 2024
@xylar xylar deleted the ocn/separate-removial-ais-river-ice-runoff branch July 1, 2024 14:10
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