Skip to content

DELWAQ-1242 Retrieving Latitude from D-Flow FM when using online DELWAQ#805

Open
michelle-deltares wants to merge 3 commits intomainfrom
all/feature/DELWAQ-1242-retrieving-Latitude-from-D-Flow-FM-when-using-online-DELWAQ
Open

DELWAQ-1242 Retrieving Latitude from D-Flow FM when using online DELWAQ#805
michelle-deltares wants to merge 3 commits intomainfrom
all/feature/DELWAQ-1242-retrieving-Latitude-from-D-Flow-FM-when-using-online-DELWAQ

Conversation

@michelle-deltares
Copy link
Copy Markdown
Contributor

What was done

  • For spherical models, the latitude of faces/cells is made available as a segment function in integrated FM-Delwaq calculations when Latitude is mentioned in the sub-file

Evidence of the work done

  • Clear from the issue description

Tests

  • New test added: c116_latitude_from_fm_to_waq

Documentation

  • Documentation updated \

Issue link

https://issuetracker.deltares.nl/browse/DELWAQ-1242

@michelle-deltares michelle-deltares added enhancement New feature or request kernel::fm Changes related to the fm kernel type::feature Adds a feature labels Apr 14, 2026
real(hp), allocatable, dimension(:, :) :: waqoutputs !< waq outputs, dim(noout,Ndkx)

integer :: isfsurf !< pointer to surface segment function
integer :: isflat !< pointer to latitude segment function
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

I suggest using the name "isflatitude" - now I read the name as "is flat", whereas it is not zero, if the model is NOT flat.

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

Are you telling me that earth isn't flat ;-).
I was in doubt, and I did notice it had become 'flat', but didn't notice the consequence of the combination with 'is'. :-D. I'll update this.

else
num_spatial_time_fuctions = num_spatial_time_fuctions + 1
isflat = num_spatial_time_fuctions
call realloc(sfunname, num_spatial_time_fuctions, keepExisting=.true., fill='latitude')
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Why not reuse the paramter clatitude?

Copy link
Copy Markdown
Contributor Author

@michelle-deltares michelle-deltares Apr 16, 2026

Choose a reason for hiding this comment

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

Oh, that is a nice suggestion. That could be don on other locations as well, but I think not always... I think I had some issues with string comparison, which I resolved by using a parameter.
I think this has to do with the function call... Since I don't want to use the parameter only for this one, and also do not want to change it everywhere else, I leave it like it is. I'll focus on the weird fact that two unrelated testcases are now failing. Restart runs, while the 'normal' version have no issue....

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

enhancement New feature or request kernel::fm Changes related to the fm kernel type::feature Adds a feature

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants