Attempt to modify solar radiation downscaling to apply to topounits

lanl · Jun 10, 2024 · 1942e06 · 1942e06
1 parent ac694a3
commit 1942e06
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/components/elm/src/main/surfrdMod.F90 b/components/elm/src/main/surfrdMod.F90
@@ -1631,6 +1631,7 @@ subroutine surfrd_get_topo_for_solar_rad(domain,filename)
 ! !USES:
     use domainMod , only : domain_type
     use fileutils , only : getfil
+    use topounit_varcon ,  only : max_topounits
 
 ! !ARGUMENTS:
     implicit none
@@ -1642,6 +1643,7 @@ subroutine surfrd_get_topo_for_solar_rad(domain,filename)
 !
 ! !REVISION HISTORY:
 ! Created by Dalei Hao
+! Modified for multiple topounits by Rich Fiorella
 !
 ! !LOCAL VARIABLES:
 !EOP
@@ -1681,7 +1683,7 @@ subroutine surfrd_get_topo_for_solar_rad(domain,filename)
     beg = domain%nbeg
     end = domain%nend
 
-    allocate(latc(beg:end),lonc(beg:end))
+    allocate(latc(beg:end),lonc(beg:end),max_topounits)
 
     call ncd_io(ncid=ncid, varname='LONGXY', flag='read', data=lonc, &
          dim1name=grlnd, readvar=readvar)
@@ -1716,7 +1718,7 @@ subroutine surfrd_get_topo_for_solar_rad(domain,filename)
          dim1name=grlnd, readvar=readvar)
     If (.not. readvar) call endrun( trim(subname)//' ERROR: SINSL_SINAS  NOT on fsurdat file' )
 
-    deallocate(latc,lonc)
+    deallocate(latc,lonc,max_topounits)
 
     call ncd_pio_closefile(ncid)