swrefindx_col1(isize,k) = swrefindx(iclm,k,jclm,isize,3) ! at 600 nm lwrefindx_col(isize,k,:) = lwrefindx(iclm,k,jclm,isize,:) radius_core_col(isize,k) = radius_core ...