How can I calculate time dilation between our clock, and clock in satellite? We know satellite is moving , so there should be time dilation due to motion, also we know due to Earth's Gravitational field there should be time dilation, what should I do, just calculate both of them and add? 
 A: Yes, just add. From GR the satellite clock goes faster, but then SR slows it down a bit, though it's still faster than on the ground.
Taylor and Wheeler's book "Exploring Black Holes" has a short chapter on GPS doing just this. The first edition (which I read) is out of print, but there is an online second edition, which I didn't read, but which I see does have a chapter on GPS.
A: You should not really separate special relativistic time dilation and general relativistic time dilation into two separate effects and you should read John Rennie's excellent answer to  Is gravitational time dilation different from other forms of time dilation? 
As an approximation you can use the expression for time dilation of an object in a circular orbit in the Schwarzschild metric.
$$ d\tau_S = \left(1 - \frac{3GM}{c^2 r_S}\right)^{1/2}\ dt,$$
where $d\tau$ represents a tick of the clock in its reference frame (the proper time), $dt$ is the tick length as observed by a distant observer in asymptotically flat space-time and $r$ is the orbital radius. This formula assumes a spherically symmetric gravitational potential and that the satellite orbits outside the Earth!
For an observer on the surface of the Earth, then they are not in orbit, so the above equation does not apply. In fact to a good approximation they are equivalent to a stationary observer in the Schwarzschild metric at distance $r$ from the centre. This approximation becomes very good for an observer at the north or south poles and is roughly true at the equator because the speed of an observer at the equator is about 15 times less than that of a satellite in orbit just above the Earth's surface.
$$d \tau_{\rm O} \simeq \left(1 - \frac{2GM}{c^2r_O}\right)^{1/2} \ dt,$$
where $d\tau_O$ would be the ticks measured at the surface of the Earth and $r_O$ is the Earth radius.
If we divide one equation by the other, you can get the expression for the relative time-dilation. 
$$ \frac{d\tau_O}{d\tau_S} = \left( \frac{1 - 2GM/c^2r_O}{1- 3GM/c^2r_S}\right)^{1/2}$$
