diff --git a/epfd_script.py b/epfd_script.py
index dadf0ab..4d157ae 100644
--- a/epfd_script.py
+++ b/epfd_script.py
@@ -102,11 +102,11 @@ def interp(lev, data, lev_int):
 	V_res_int = np.copy(V_res) 
 	O_res_int = np.copy(O_res)
 
-
+Fdiv_int /= np.reshape(np.cos(phi),(1,1,nlat))   
 Fp_int = Fp_int*np.reshape(np.cos(phi),(1,1,nlat))   
 Fphi_int = Fphi_int/a
 
-Fp_int = Fp_int/1.0e4
+Fp_int = Fp_int/1.0e5
 Fphi_int = Fphi_int/PI
 
 if scale_by_sqrt_p: