Do_Upwind

input: int qt, int *nmaploc, int *nmapadj,Coord *nc, Coord *tc, Coord *pc, double **dat, double **adj_dat, double **dat_fl, double **adj_dat_fl, double **f, double **f_adj

output:

Description: This function does the upwind,the initial data for the face fa is in edges [0 ...Nfields][ 0 ...qt-1], adj_edges,the output is put to edges_fl int. flux, adj_edges_fl - ext. flux, f - upwinded flux.