Per passare una sottomatrice, in f77 si può fare:
call horizontal_diffusion(conc(1,1,1,spec),1,h_diff_coeff)
in f90 c'è la possibilità di scrivere
call horizontal_diffusion(conc(:,:,:,spec),1,h_diff_coeff)