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

Correcting units on accumulated variables in ocean conservation check #6257

Merged
merged 1 commit into from
Mar 13, 2024

Conversation

darincomeau
Copy link
Member

Corrects units on accumulated variables in ocean conservation check, which have been multiplied by area.

[BFB]

@darincomeau
Copy link
Member Author

Copying comment over from #6229

here's my comment from the Ocean Discussion PR for more context:

@mark-petersen can you check the proposed unit changes in the conservation check analysis member's Registry here E3SM-Ocean-Discussion@285be91?

@cbegeman and I looked into this and see the flux variables being multiplied by m^2 here:
E3SM/components/mpas-ocean/src/analysis_members/mpas_ocn_conservation_check.F

Line 954 in 50485d0

sumArray( 1) = sumArray( 1) + areaCell(iCell) * rainFlux(iCell)

and not divided out again. Also when we look at the values that are written out to the conservation check analysis member, we can back out values in kg m^-2 s^-1 that appear in the ocn.log file from this analysis member by scaling by 1/area again.
The reason we're looking here is we want to use these accumulated*Flux quantities to produce a net AIS runoff time series plot for polar simulations to compare to global ssh anomaly, and to us everything is consistent if these quantities are actually in units kg s^-1.

Also, many thanks @mark-petersen for your work on this analysis member in the past. It's made verification of this coupler budget accounting much easier.

Copy link
Contributor

@cbegeman cbegeman left a comment

Choose a reason for hiding this comment

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

I agree that this change is needed.

Copy link
Contributor

@mark-petersen mark-petersen left a comment

Choose a reason for hiding this comment

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

@darincomeau Thanks, this change is correct. Looking at the output:

  MASS CONSERVATION CHECK

 MASS FLUXES
 MPAS-Ocean name          kg/s (F)       coupler name     short name     kg/m^2/s*1e6 (F/A)
 frazilFreshwaterFlux    -9.60609187E+07 o2x_Fioo_frazil  wfreeze             -0.18831740
 seaIceFreshwaterFlux     3.80124701E+09 x2o_Fioi_meltw   wmelt                7.45194768
 rainFlux                 1.35659405E+10 x2o_Faxa_rain    wrain               26.59460915
 snowFlux                 3.95339636E+08 x2o_Faxa_snow    wsnow                0.77502206
 icebergFreshwaterFlux    9.56158749E+07 x2o_Fioi_bergw   wberg                0.18744493
 evaporationFlux         -1.59515036E+10 x2o_Foxx_evap    wevap              -31.27125659
 riverRunoffFlux          6.97379066E+07 x2o_Foxx_rofl    wrunoff              0.13671388

those variables are clearly in kg/s, based on both the labels and magnitude of the numbers.

jonbob added a commit that referenced this pull request Mar 12, 2024
…6257)

Correcting units on accumulated variables in ocean conservation check

Corrects units on accumulated variables in ocean conservation check,
which have been multiplied by area.

[BFB]
@jonbob
Copy link
Contributor

jonbob commented Mar 12, 2024

Passes:

  • ERS_Ld5.T62_oQU120.CMPASO-NYF.chrysalis_intel
  • ERP_Ld3.ne30pg2_r05_IcoswISC30E3r5.WCYCL1850.chrysalis_intel.allactive-pioroot1

merged to next

@jonbob jonbob merged commit 93e511d into E3SM-Project:master Mar 13, 2024
10 checks passed
@jonbob
Copy link
Contributor

jonbob commented Mar 13, 2024

merged to master

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
BFB PR leaves answers BFB mpas-ocean
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants