made changes in calcflows: add subroutine for calculating flows by OX and OY in the same time, change name of variables 'hyj' -> 'hjy'
add subroutine for calculating R of flows
replace calculation of flows by new subroutine instead of calling two subroutines
replace CalculateFRX and CalculateFRY by calling a new subroutine CalculateFR, refactor CalculateRXY