added the NTV calculation to the multiple species fortran code /HS