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.

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

In filter_topo added extrapolation to halo points #571

Merged
merged 11 commits into from
Sep 17, 2021
Merged

In filter_topo added extrapolation to halo points #571

Changes from 1 commit
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
b602411
Add extrapolation for geolat and geolon to halo points. This commit r…
RatkoVasic-NOAA Aug 20, 2021
dc72965
Add comments for doxygen.
RatkoVasic-NOAA Aug 23, 2021
06cc530
Fix bugs in FV3_zs_filter and two_delta_filter.
RatkoVasic-NOAA Aug 25, 2021
1ec6689
Fix oro and mask fields on halo regions in regional application.
RatkoVasic-NOAA Aug 25, 2021
f87f1f5
Move new routine 'ftst_fill_regional_halo' to utils.F90
GeorgeGayno-NOAA Aug 26, 2021
88fa6d9
Comment out call to routine in new unit test.
GeorgeGayno-NOAA Aug 26, 2021
852405d
Fixed coefficient a1 in subroutine two_delta_filter. Corners were use…
RatkoVasic-NOAA Sep 1, 2021
9b73442
Add test for routine fill_regional_halo.
RatkoVasic-NOAA Sep 10, 2021
8313442
Fix 'module use' in test file.
RatkoVasic-NOAA Sep 10, 2021
bab5f9b
Add deallocate.
RatkoVasic-NOAA Sep 10, 2021
fe4072a
Merge branch 'develop' into filter_topo
RatkoVasic-NOAA Sep 16, 2021
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
39 changes: 39 additions & 0 deletions sorc/grid_tools.fd/filter_topo.fd/filter_topo.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -748,6 +748,9 @@ subroutine read_grid_file(regional)
call fill_cubic_grid_halo(geolat_c, geolat_c, ng, 1, 1, 1, 1)
if(.not. nested) call fill_bgrid_scalar_corners(geolon_c, ng, npx, npy, isd, jsd, XDir)
if(.not. nested) call fill_bgrid_scalar_corners(geolat_c, ng, npx, npy, isd, jsd, YDir)
else
call fill_regional_halo(geolon_c, ng)
call fill_regional_halo(geolat_c, ng)
endif

!--- compute grid cell center
Expand Down Expand Up @@ -1144,6 +1147,42 @@ subroutine fill_cubic_grid_halo(data, data2, halo, ioff, joff, sign1, sign2)

end subroutine fill_cubic_grid_halo

!> This routine extrapolate geolat_c and geolon_c halo points for the
!! regional standalone grid. Halo points are needed for dxc and dyc
!! calculation.
!!
RatkoVasic-NOAA marked this conversation as resolved.
Show resolved Hide resolved
!! @param[in] data ???
!! @param[in] halo ???
!! @author Ratko
subroutine fill_regional_halo(data, halo)
RatkoVasic-NOAA marked this conversation as resolved.
Show resolved Hide resolved
integer, intent(in) :: halo
real, dimension(1-halo:,1-halo:,:), intent(inout) :: data
integer :: h, i_st, i_ed, j_st, j_ed

i_st=1
i_ed=npx
j_st=1
j_ed=npy

do h = 1, halo
data(i_st:i_ed, j_st-1 , :) = 2* data(i_st:i_ed, j_st , :) - data(i_st:i_ed, j_st+1 , :)! north
data(i_st:i_ed, j_ed+1 , :) = 2* data(i_st:i_ed, j_ed , :) - data(i_st:i_ed, j_ed-1 , :)! south
data(i_st-1 , j_st:j_ed, :) = 2* data(i_st , j_st:j_ed, :) - data(i_st+1 , j_st:j_ed, :)! east
data(i_ed+1 , j_st:j_ed, :) = 2* data(i_ed , j_st:j_ed, :) - data(i_ed-1 , j_st:j_ed, :)! west

data(i_st-1, j_st-1, :) = (data(i_st-1, j_st, :) + data(i_st, j_st-1, :))*0.5 !NW Corner
data(i_ed+1, j_st-1, :) = (data(i_ed+1, j_st, :) + data(i_ed, j_st-1, :))*0.5 !NE Corner
data(i_st-1, j_ed+1, :) = (data(i_st-1, j_ed, :) + data(i_st, j_ed+1, :))*0.5 !SW Corner
data(i_ed+1, j_ed+1, :) = (data(i_ed+1, j_ed, :) + data(i_ed, j_ed+1, :))*0.5 !SE Corner

i_st=i_st-1
i_ed=i_ed+1
j_st=j_st-1
j_ed=j_ed+1
enddo

end subroutine fill_regional_halo

!> ???
!!
!! @param[in] is ???
Expand Down