Skip to content

Add read-in capability for urban 2d fraction #131

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
cenlinhe merged 1 commit into from
Jan 5, 2024
+64 −1
Merged

Add read-in capability for urban 2d fraction #131

Show file tree
Hide file tree
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
11 changes: 10 additions & 1 deletion hrldas/IO_code/module_NoahMP_hrldas_driver.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -292,6 +292,15 @@ subroutine land_driver_ini(NTIME_out,wrfits,wrfite,wrfjts,wrfjte)
call READ_TILE_DRAIN_MAP(NoahmpIO%TDINPUT_FLNM, NoahmpIO%XSTART, NoahmpIO%XEND, NoahmpIO%YSTART, NoahmpIO%YEND, NoahmpIO%TD_FRACTION)
endif

!------------------------------------------------------------------------
! For SF_URBAN_PHYSICS > 0 read 2D map of urban parameters
!------------------------------------------------------------------------

NoahmpIO%FRC_URB2D = 0.0
if (NoahmpIO%SF_URBAN_PHYSICS > 0) then
call READ_URBAN_MAP(NoahmpIO%HRLDAS_SETUP_FILE, NoahmpIO%XSTART, NoahmpIO%XEND, NoahmpIO%YSTART, NoahmpIO%YEND, NoahmpIO%FRC_URB2D)
endif

!----------------------------------------------------------------------
! Initialize Model State
!----------------------------------------------------------------------
Expand Down Expand Up @@ -1014,7 +1023,7 @@ subroutine land_driver_exe(itime)
NoahmpIO%drelr_urb2d, NoahmpIO%drelb_urb2d, NoahmpIO%drelg_urb2d, & !H urban
NoahmpIO%flxhumr_urb2d, NoahmpIO%flxhumb_urb2d, NoahmpIO%flxhumg_urb2d, & !H urban
NoahmpIO%julday, NoahmpIO%yr, & !H urban
NoahmpIO%frc_urb2d, NoahmpIO%utype_urb2d, & !I urban
NoahmpIO%FRC_URB2D, NoahmpIO%utype_urb2d, & !I urban
NoahmpIO%chs, NoahmpIO%chs2, NoahmpIO%cqs2, & !H
NoahmpIO%num_urban_ndm, NoahmpIO%urban_map_zrd, NoahmpIO%urban_map_zwd, & !I multi-layer urban
NoahmpIO%urban_map_gd, & !I multi-layer urban
Expand Down
54 changes: 54 additions & 0 deletions hrldas/IO_code/module_hrldas_netcdf_io.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -806,6 +806,60 @@ subroutine read_tile_drain_map(tdinput_flnm, &

end subroutine read_tile_drain_map

!---------------------------------------------------------------------------------------------------------

subroutine read_urban_map(wrfinput_flnm, &
xstart, xend, ystart, yend, &
urban_fraction)

implicit none
character(len=*), intent(in) :: wrfinput_flnm
integer, intent(in) :: xstart, xend, ystart, yend
real, dimension(xstart:xend,ystart:yend), intent(out) :: urban_fraction
character(len=24) :: name
integer :: ierr,iret
integer :: ncid, varid
integer :: rank

#ifdef _PARALLEL_
call MPI_COMM_RANK(MPI_COMM_WORLD, rank, ierr)
if (ierr /= MPI_SUCCESS) stop "MPI_COMM_RANK"
#else
rank = 0
#endif


! Open the NetCDF file.
if (rank == 0) write(*,'("wrfinput_flnm: ''", A, "''")') trim(wrfinput_flnm)
#ifdef _PARALLEL_
ierr = nf90_open_par(wrfinput_flnm, NF90_NOWRITE, MPI_COMM_WORLD, MPI_INFO_NULL, ncid)
#else
ierr = nf90_open(wrfinput_flnm, NF90_NOWRITE, ncid)
#endif
if (ierr /= 0) then
write(*,'("read_urban_map: Problem opening wrfinput file: ''", A, "''")') trim(wrfinput_flnm)
#ifdef _PARALLEL_
call mpi_finalize(ierr)
if (ierr /= 0) write(*, '("Problem with MPI_finalize.")')
#endif
stop
endif

! Get Urban Fraction
name = "FRC_URB2D"
iret = nf90_inq_varid(ncid, name, varid)
if (iret == 0) then
ierr = nf90_get_var(ncid, varid, urban_fraction, start=(/xstart,ystart/), count=(/xend-xstart+1,yend-ystart+1/))
else
write(*,*) "MODULE_HRLDAS_NETCDF_IO: Problem finding variable '"//trim(name)//"' in NetCDF file. Using default table values."
endif

! Close the NetCDF file
ierr = nf90_close(ncid)
if (ierr /= 0) stop "MODULE_NOAHLSM_HRLDAS_INPUT: read_urban_map: NF90_CLOSE"

end subroutine read_urban_map

!---------------------------------------------------------------------------------------------------------

subroutine read_3d_soil(spatial_filename,xstart, xend,ystart, yend, &
Expand Down