Skip to content

Commit

Permalink
Part of the fix for issue NOAA-EMC#54.
Browse files Browse the repository at this point in the history
  • Loading branch information
kshedstrom committed Apr 6, 2022
1 parent 6f48b9d commit 24f412d
Showing 2 changed files with 4 additions and 4 deletions.
2 changes: 1 addition & 1 deletion src/framework/MOM_horizontal_regridding.F90
Original file line number Diff line number Diff line change
@@ -860,7 +860,7 @@ subroutine horiz_interp_and_extrap_tracer_fms_id(fms_id, Time, conversion, G, t
do k=1,kd
do j=js,je
do i=is,ie
tr_z(i,j,k)=data_in(i,j,k)
tr_z(i,j,k)=data_in(i,j,k) * conversion
if (.not. ans_2018) mask_z(i,j,k) = 1.
if (abs(tr_z(i,j,k)-missing_value) < abs(roundoff*missing_value)) mask_z(i,j,k) = 0.
enddo
6 changes: 3 additions & 3 deletions src/parameterizations/vertical/MOM_ALE_sponge.F90
Original file line number Diff line number Diff line change
@@ -143,7 +143,7 @@ module MOM_ALE_sponge

!> This subroutine determines the number of points which are within sponges in this computational
!! domain. Only points that have positive values of Iresttime and which mask2dT indicates are ocean
!! points are included in the sponges. It also stores the target interface heights. This
!! points are included in the sponges. It also stores the target interface heights.
subroutine initialize_ALE_sponge_fixed(Iresttime, G, GV, param_file, CS, data_h, nz_data, &
Iresttime_u_in, Iresttime_v_in)

@@ -965,7 +965,7 @@ subroutine apply_ALE_sponge(h, dt, G, GV, US, CS, Time)
if (CS%time_varying_sponges) then
nz_data = CS%Ref_val_u%nz_data
! Interpolate from the external horizontal grid and in time
call horiz_interp_and_extrap_tracer(CS%Ref_val_u%id, Time, 1.0, G, sp_val, mask_z, z_in, &
call horiz_interp_and_extrap_tracer(CS%Ref_val_u%id, Time, US%m_s_to_L_T, G, sp_val, mask_z, z_in, &
z_edges_in, missing_value, CS%reentrant_x, CS%tripolar_N, .false., &
spongeOnGrid=CS%SpongeDataOngrid, m_to_Z=US%m_to_Z,&
answers_2018=CS%hor_regrid_answers_2018)
@@ -1014,7 +1014,7 @@ subroutine apply_ALE_sponge(h, dt, G, GV, US, CS, Time)
deallocate(sp_val, mask_u, mask_z, hsrc)
nz_data = CS%Ref_val_v%nz_data
! Interpolate from the external horizontal grid and in time
call horiz_interp_and_extrap_tracer(CS%Ref_val_v%id, Time, 1.0, G, sp_val, mask_z, z_in, &
call horiz_interp_and_extrap_tracer(CS%Ref_val_v%id, Time, US%m_s_to_L_T, G, sp_val, mask_z, z_in, &
z_edges_in, missing_value, CS%reentrant_x, CS%tripolar_N, .false., &
spongeOnGrid=CS%SpongeDataOngrid, m_to_Z=US%m_to_Z,&
answers_2018=CS%hor_regrid_answers_2018)

0 comments on commit 24f412d

Please sign in to comment.