Subroutine QUICK(Con_f,Fi_U,Fi_C,Fi_D,Fi_DD,Delta_f)
If(Con_f.GT.0.)then
Delta_f= 0.375 * Fi_D - 0.25 * Fi_C - 0.125 * Fi_U
else
Delta_f= -0.375 * Fi_DD + 0.25 * Fi_D + 0.125 * Fi_C
End If
Return End