I wanted to add a new term in this parameterization.
subroutine hsfluxpv(nz,n,bl,bw,albr,rs,swddif,emr,rl,tr,tair,sigma,hspv,eppv,pv_frac_roof,uout,vout,tpv,dt)
!
implicit none
!
! Input variables
!
integer,intent(in) :: nz !Maximum number of vertical levels in the urban...