Skip to content

fix bug in writing fesom.mesh.diag.nc #747

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.

Sign up for GitHub

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

Merged
JanStreffing merged 3 commits into from
Sep 2, 2025
Merged

fix bug in writing fesom.mesh.diag.nc #747

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
27a1060
fix bug in writing fesom.mesh.diag.nc, like the code was written unt…
patrickscholz Jul 24, 2025
bb80b71
Merge branch 'main' into workbench_bugfix_fesom.mesh.diag.nc
suvarchal Aug 19, 2025
b18f11c
Merge branch 'main' into workbench_bugfix_fesom.mesh.diag.nc
JanStreffing Sep 2, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 8 additions & 4 deletions src/io_mesh_info.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -530,7 +530,8 @@ subroutine write_mesh_info(partit, mesh)
allocate(rbuffer(elem2D))
do i=1, 3
call gather_elem(gradient_sca(i, 1:myDim_elem2D), rbuffer, partit)
call my_put_vara(ncid, gradient_sca_x_id, (/1, 4-i/), (/elem2D, 1/), rbuffer, partit) ! (4-i), NETCDF will permute otherwise
!!PS call my_put_vara(ncid, gradient_sca_x_id, (/1, 4-i/), (/elem2D, 1/), rbuffer, partit) ! (4-i), NETCDF will permute otherwise
call my_put_vara(ncid, gradient_sca_x_id, (/1, i/), (/elem2D, 1/), rbuffer, partit)
end do
deallocate(rbuffer)

Expand All @@ -539,7 +540,8 @@ subroutine write_mesh_info(partit, mesh)
allocate(rbuffer(elem2D))
do i=1, 3
call gather_elem(gradient_sca(i+3, 1:myDim_elem2D), rbuffer, partit)
call my_put_vara(ncid, gradient_sca_y_id, (/1, 4-i/), (/elem2D, 1/), rbuffer, partit)! (4-i), NETCDF will permute otherwise
!!PS call my_put_vara(ncid, gradient_sca_y_id, (/1, 4-i/), (/elem2D, 1/), rbuffer, partit)! (4-i), NETCDF will permute otherwise
call my_put_vara(ncid, gradient_sca_y_id, (/1, i/), (/elem2D, 1/), rbuffer, partit)
end do
deallocate(rbuffer)

Expand All @@ -548,7 +550,8 @@ subroutine write_mesh_info(partit, mesh)
allocate(rbuffer(elem2D))
do i=1, 3
call gather_elem(gradient_vec(i, 1:myDim_elem2D), rbuffer, partit)
call my_put_vara(ncid, gradient_vec_x_id, (/1, 4-i/), (/elem2D, 1/), rbuffer, partit) ! (4-i), NETCDF will permute otherwise
!!PS call my_put_vara(ncid, gradient_vec_x_id, (/1, 4-i/), (/elem2D, 1/), rbuffer, partit) ! (4-i), NETCDF will permute otherwise
call my_put_vara(ncid, gradient_vec_x_id, (/1, i/), (/elem2D, 1/), rbuffer, partit) ! (4-i), NETCDF will permute otherwise
end do
deallocate(rbuffer)

Expand All @@ -557,7 +560,8 @@ subroutine write_mesh_info(partit, mesh)
allocate(rbuffer(elem2D))
do i=1, 3
call gather_elem(gradient_vec(i+3, 1:myDim_elem2D), rbuffer, partit)
call my_put_vara(ncid, gradient_vec_y_id, (/1, 4-i/), (/elem2D, 1/), rbuffer, partit)! (4-i), NETCDF will permute otherwise
!!PS call my_put_vara(ncid, gradient_vec_y_id, (/1, 4-i/), (/elem2D, 1/), rbuffer, partit)! (4-i), NETCDF will permute otherwise
call my_put_vara(ncid, gradient_vec_y_id, (/1, i/), (/elem2D, 1/), rbuffer, partit)! (4-i), NETCDF will permute otherwise
end do
deallocate(rbuffer)

Expand Down