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

MLBalance returns NaN's when not using the model's resolution for B #1124

Open
guillaumevernieres opened this issue Jan 21, 2025 · 4 comments
Assignees
Labels
bug Something is not working SOCA Sea-ice, Ocean, and Coupled Assimilation

Comments

@guillaumevernieres
Copy link
Contributor

Description

There is something wrong in the MLBalance saber block. It returns NaN's in the first inner iteration. I suspect something related to the trajectory.

To Reproduce:

Run a variational application using the MLBalance block at a different resolution than the deterministic.

@guillaumevernieres guillaumevernieres added bug Something is not working SOCA Sea-ice, Ocean, and Coupled Assimilation labels Jan 21, 2025
@guillaumevernieres guillaumevernieres self-assigned this Jan 21, 2025
@travissluka
Copy link
Collaborator

you are using this parameter, right? (

change background resolution: true
)

@guillaumevernieres
Copy link
Contributor Author

you are using this parameter, right? (

soca/test/testinput/3dvar_lowres.yml

Line 56 in 3dd464a

change background resolution: true
)

yes

@guillaumevernieres
Copy link
Contributor Author

Skipping the mask/ghost cells like it is done in the other saber blocks does the trick. I'll try to squeeze in a code cleanup and PR in between v17 work duties.

@travissluka
Copy link
Collaborator

Skipping the mask/ghost cells like it is done in the other saber blocks does the trick. I'll try to squeeze in a code cleanup and PR in between v17 work duties.

great. You probably had junk in your halos. For future reference, if you are using information in the ghost cells, make sure you always call halo exchange before accessing them. atlas will only do an actual halo exchange if an individual field has been previously marked dirty. Related: if you change this to skip ghost cells, be sure to call set_dirty() on any field you write to

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something is not working SOCA Sea-ice, Ocean, and Coupled Assimilation
Projects
None yet
Development

No branches or pull requests

2 participants