Skip to content

Feature/chgres ice check #35

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
GeorgeGayno-NOAA merged 2 commits into from
Dec 4, 2019
Merged

Feature/chgres ice check #35

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
a9081cc
feature/chgres_ice_check: This commit references #32.
GeorgeGayno-NOAA Nov 22, 2019
289d36a
feature/chgres_ice_check: This commit references #32.
GeorgeGayno-NOAA Nov 22, 2019
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
20 changes: 11 additions & 9 deletions sorc/chgres_cube.fd/model_grid.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -161,6 +161,8 @@ subroutine define_input_grid_gaussian(localpet, npets)
real(esmf_kind_r8), allocatable :: longitude(:,:)
real(esmf_kind_r8), pointer :: lat_src_ptr(:,:)
real(esmf_kind_r8), pointer :: lon_src_ptr(:,:)
real(esmf_kind_r8), pointer :: lat_corner_src_ptr(:,:)
real(esmf_kind_r8), pointer :: lon_corner_src_ptr(:,:)
real(esmf_kind_r8) :: deltalon
real(esmf_kind_r8), allocatable :: slat(:), wlat(:)

Expand Down Expand Up @@ -320,40 +322,40 @@ subroutine define_input_grid_gaussian(localpet, npets)
call error_handler("IN GridAddCoord", rc)

print*,"- CALL GridGetCoord FOR INPUT GRID X-COORD."
nullify(lon_src_ptr)
nullify(lon_corner_src_ptr)
call ESMF_GridGetCoord(input_grid, &
staggerLoc=ESMF_STAGGERLOC_CORNER, &
coordDim=1, &
farrayPtr=lon_src_ptr, rc=rc)
farrayPtr=lon_corner_src_ptr, rc=rc)
if(ESMF_logFoundError(rcToCheck=rc,msg=ESMF_LOGERR_PASSTHRU,line=__line__,file=__file__)) &
call error_handler("IN GridGetCoord", rc)

print*,"- CALL GridGetCoord FOR INPUT GRID Y-COORD."
nullify(lat_src_ptr)
nullify(lat_corner_src_ptr)
call ESMF_GridGetCoord(input_grid, &
staggerLoc=ESMF_STAGGERLOC_CORNER, &
coordDim=2, &
computationalLBound=clb, &
computationalUBound=cub, &
farrayPtr=lat_src_ptr, rc=rc)
farrayPtr=lat_corner_src_ptr, rc=rc)
if(ESMF_logFoundError(rcToCheck=rc,msg=ESMF_LOGERR_PASSTHRU,line=__line__,file=__file__)) &
call error_handler("IN GridGetCoord", rc)

print*,'bounds for corners ',localpet,clb(1),cub(1),clb(2),cub(2)

do j = clb(2), cub(2)
do i = clb(1), cub(1)
lon_src_ptr(i,j) = longitude(i,1) - (0.5_esmf_kind_r8*deltalon)
if (lon_src_ptr(i,j) > 360.0_esmf_kind_r8) lon_src_ptr(i,j) = lon_src_ptr(i,j) - 360.0_esmf_kind_r8
lon_corner_src_ptr(i,j) = longitude(i,1) - (0.5_esmf_kind_r8*deltalon)
if (lon_corner_src_ptr(i,j) > 360.0_esmf_kind_r8) lon_corner_src_ptr(i,j) = lon_corner_src_ptr(i,j) - 360.0_esmf_kind_r8
if (j == 1) then
lat_src_ptr(i,j) = 90.0_esmf_kind_r8
lat_corner_src_ptr(i,j) = 90.0_esmf_kind_r8
cycle
endif
if (j == jp1_input) then
lat_src_ptr(i,j) = -90.0_esmf_kind_r8
lat_corner_src_ptr(i,j) = -90.0_esmf_kind_r8
cycle
endif
lat_src_ptr(i,j) = 0.5_esmf_kind_r8 * (latitude(i,j-1)+ latitude(i,j))
lat_corner_src_ptr(i,j) = 0.5_esmf_kind_r8 * (latitude(i,j-1)+ latitude(i,j))
enddo
enddo

Expand Down
7 changes: 5 additions & 2 deletions sorc/chgres_cube.fd/surface.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -603,8 +603,11 @@ subroutine interp(localpet)
if (localpet == 0) then
do j = 1, j_target
do i = 1, i_target
if (data_one_tile(i,j) < 0.15) data_one_tile(i,j) = 0.0
if (data_one_tile(i,j) >= 0.15) mask_target_one_tile(i,j) = 2
if (data_one_tile(i,j) > 1.0_esmf_kind_r8) then
data_one_tile(i,j) = 1.0_esmf_kind_r8
endif
if (data_one_tile(i,j) < 0.15_esmf_kind_r8) data_one_tile(i,j) = 0.0_esmf_kind_r8
if (data_one_tile(i,j) >= 0.15_esmf_kind_r8) mask_target_one_tile(i,j) = 2
enddo
enddo
endif
Expand Down