Dmitry Alexeev suggests adding (:,:,:) in vr()

earth-system-radiation · Jan 9, 2024 · a79eec4 · a79eec4
1 parent 4652ce7
commit a79eec4
Showing 1 changed file with 7 additions and 7 deletions.
diff --git a/tests/mo_testing_utils.F90 b/tests/mo_testing_utils.F90
@@ -237,7 +237,7 @@ subroutine increment_with_nstr(atmos)
   end subroutine increment_with_nstr 
   ! ----------------------------------------------------------------------------
   !
-  ! Vertically reverse optical
+  ! Vertically reverse optical properties
   ! 
   subroutine vr(atmos, sources)
     class(ty_optical_props_arry), intent(inout) :: atmos
@@ -252,16 +252,16 @@ subroutine vr(atmos, sources)
 
     select type (atmos)
       type is (ty_optical_props_2str)
-        atmos%ssa = atmos%ssa(:,nlay:1:-1,:)
-        atmos%g   = atmos%g  (:,nlay:1:-1,:)
+        atmos%ssa(:,:,:) = atmos%ssa(:,nlay:1:-1,:)
+        atmos%g  (:,:,:)= atmos%g  (:,nlay:1:-1,:)
       type is (ty_optical_props_nstr)
-        atmos%ssa = atmos%ssa(:,nlay:1:-1,:)
-        atmos%p   = atmos%p(:,:,nlay:1:-1,:)
+        atmos%ssa(:,:,:) = atmos%ssa(:,nlay:1:-1,:)
+        atmos%p(:,:,:,:) = atmos%p(:,:,nlay:1:-1,:)
     end select    
 
     if(present(sources)) then 
-      sources%lev_source = sources%lev_source(:,nlay+1:1:-1,:)
-      sources%lay_source = sources%lay_source(:,nlay  :1:-1,:)
+      sources%lev_source(:,:,:) = sources%lev_source(:,nlay+1:1:-1,:)
+      sources%lay_source(:,:,:) = sources%lay_source(:,nlay  :1:-1,:)
     end if                           
   end subroutine vr
   ! ----------------------------------------------------------------------------