What size is that array?

earth-system-radiation · Jan 4, 2024 · c6524e8 · c6524e8
1 parent 080fad2
commit c6524e8
Showing 1 changed file with 8 additions and 5 deletions.
diff --git a/tests/rte_unit_tests.F90 b/tests/rte_unit_tests.F90
@@ -60,11 +60,11 @@ program rte_unit_tests
   !
   real(wp), parameter :: sigma = 5.670374419e-8_wp, & ! Stefan-Boltzmann constant 
                          D     = 1.66_wp              ! Diffusivity angle, from single-angle RRTMGP solver
-  real(wp), dimension(    ncol), parameter :: sfc_t     = [(285._wp, icol = 1,ncol/2), & 
-                                                           (310._wp, icol = 1,ncol/2)]
-  real(wp), dimension(    ncol), parameter :: lw_total_tau = [0.1_wp, 1._wp, 10._wp, 50._wp, &
-                                                              0.1_wp, 1._wp, 10._wp, 50._wp] ! Would be nice to parameterize 
-  real(wp), dimension(1:1,ncol), parameter :: sfc_emis   = 1._wp
+  real(wp), dimension(  ncol), parameter :: sfc_t     = [(285._wp, icol = 1,ncol/2), & 
+                                                         (310._wp, icol = 1,ncol/2)]
+  real(wp), dimension(  ncol), parameter :: lw_total_tau = [0.1_wp, 1._wp, 10._wp, 50._wp, &
+                                                            0.1_wp, 1._wp, 10._wp, 50._wp] ! Would be nice to parameterize 
+  real(wp), dimension(1,ncol), parameter :: sfc_emis   = 1._wp
 
   type(ty_optical_props_1scl) :: lw_atmos 
   type(ty_source_func_lw)     :: lw_sources
@@ -79,6 +79,7 @@ program rte_unit_tests
 
   ! ------------------------------------------------------------------------------------------------------
   top_at_1 = .true. 
+  print *, "sizes of sfc_emis", size(sfc_emis,1), size(sfc_emis,2)
   ! ------------------------------------------------------------------------------------------------------
   ! 
   ! Longwave tests
@@ -323,6 +324,8 @@ subroutine lw_clear_sky_subset
     call stop_on_err(atmos_subset%init(lw_atmos))
     fluxes%flux_up => up
     fluxes%flux_dn => dn
+    print *, "sizes of sfc_emis in the subroutine", size(sfc_emis,1), size(sfc_emis,2)
+
     do i = 1, 2
       colS = ((i-1) * ncol/2) + 1
       colE = i * ncol/2