Skip to content

Add SAVE_UNMASKED_GEOM_FILE option #28

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
yichengt900 merged 1 commit into from
May 21, 2025
Merged

Add SAVE_UNMASKED_GEOM_FILE option #28

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
14 changes: 10 additions & 4 deletions src/framework/MOM_domains.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -128,6 +128,7 @@ subroutine MOM_domains_init(MOM_dom, param_file, symmetric, static_memory, &
character(len=200) :: mask_table ! The file name and later the full path to the diag table
character(len=64) :: inc_nm ! The name of the memory include file
character(len=200) :: mesg ! A string to use for error messages
logical :: do_unmasked_ocn_geo ! If true, create an unmasked domain for the ocean geometry file.

integer :: nip_parsed, njp_parsed
character(len=8) :: char_xsiz, char_ysiz, char_niglobal, char_njglobal
Expand Down Expand Up @@ -444,10 +445,15 @@ subroutine MOM_domains_init(MOM_dom, param_file, symmetric, static_memory, &

! Create an unmasked domain if requested. This is used for writing out unmasked ocean geometry.
if (present(MOM_dom_unmasked) .and. mask_table_exists) then
call MOM_define_layout(n_global, PEs_used, layout_unmasked)
call create_MOM_domain(MOM_dom_unmasked, n_global, n_halo, reentrant, tripolar_N, layout_unmasked, &
domain_name=domain_name, symmetric=symmetric, thin_halos=thin_halos, &
nonblocking=nonblocking)
call get_param(param_file, mdl, "SAVE_UNMASKED_GEOM_FILE", do_unmasked_ocn_geo, &
"If true and a mask table is used, create an unmasked MOM domain that is "//&
"saved in the ocean geometry file.", default=.false.)
if (do_unmasked_ocn_geo) then
call MOM_define_layout(n_global, PEs_used, layout_unmasked)
call create_MOM_domain(MOM_dom_unmasked, n_global, n_halo, reentrant, tripolar_N, layout_unmasked, &
domain_name=domain_name, symmetric=symmetric, thin_halos=thin_halos, &
nonblocking=nonblocking)
endif
endif

call create_MOM_domain(MOM_dom, n_global, n_halo, reentrant, tripolar_N, layout, &
Expand Down