call fortran with matrix like LBNL