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

halo update bugs in ACCESS driver? #70

Open
aekiss opened this issue Mar 31, 2023 · 7 comments · May be fixed by #71
Open

halo update bugs in ACCESS driver? #70

aekiss opened this issue Mar 31, 2023 · 7 comments · May be fixed by #71

Comments

@aekiss
Copy link
Contributor

aekiss commented Mar 31, 2023

Hi @russfiedler, following on from #68, these also look wrong to me (should be field_type_scalar, otherwise their signs will be flipped on the tripole in halo updates)

call ice_HaloUpdate(um_tmlt, halo_info, field_loc_center,field_type_vector)
call ice_HaloUpdate(um_bmlt, halo_info, field_loc_center,field_type_vector)

and

call ice_HaloUpdate(um_swflx, halo_info,field_loc_center,field_type_vector)
call ice_HaloUpdate(um_lwflx, halo_info,field_loc_center,field_type_vector)
call ice_HaloUpdate(um_shflx, halo_info,field_loc_center,field_type_vector)
call ice_HaloUpdate(um_press, halo_info,field_loc_center,field_type_vector)
call ice_HaloUpdate(um_co2, halo_info, field_loc_center, field_type_vector)
call ice_HaloUpdate(um_wnd, halo_info, field_loc_center, field_type_vector)

@aekiss aekiss linked a pull request Apr 3, 2023 that will close this issue
@aekiss
Copy link
Contributor Author

aekiss commented Apr 3, 2023

I've made a PR to fix this. @aidanheerdegen does ACCESS-CM2 use this same code?

@aidanheerdegen
Copy link
Contributor

I've made a PR to fix this. @aidanheerdegen does ACCESS-CM2 use this same code?

ACCESS-OM2 uses the auscom driver I believe, and yes I think the access driver is used for ACCESS-CM2. @MartinDix would have a better idea, and would probably be very interested.

@MartinDix
Copy link

The CM2 driver has these problems. I hope it's not as catastrophic as it sounds.

I've started a new run with a fix.

@aekiss
Copy link
Contributor Author

aekiss commented Apr 5, 2023

The bug was introduced in Dec 2014 - does that mean it doesn't affect previous versions of ACCESS-CM and ACCESS-ESM?
49b80b6#diff-fc19414c40fff80d4d66454635b302d42cec694887445d8c275eaf2a219f9285

@MartinDix
Copy link

A CM2 run with a fixed driver gives identical results !! ??

Does that mean these halos aren't actually used?

@aekiss
Copy link
Contributor Author

aekiss commented Apr 5, 2023

Hm, maybe we got lucky then? Would we ever need to take horizontal finite differences of these scalar fields?

@MartinDix
Copy link

I think the fluxes would only ever be used in the column.

Could test by setting the halos of these variables to NaNs.

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 a pull request may close this issue.

3 participants