diff --git a/core/quasigeostrophic.py b/core/quasigeostrophic.py
index 98598ae..fe90fea 100644
--- a/core/quasigeostrophic.py
+++ b/core/quasigeostrophic.py
@@ -114,8 +114,7 @@ def step(self, t, dt):
         if self.noslip:
             self.add_noslip(self.var.state)
 
-        # self.set_psi_from_pv()
-        
+        self.set_psi_from_pv()
         # 3/ diagnostic fields
         self.var.state[self.ipva] = self.var.state[self.ipv] - self.pvback
         self.var.state[self.ivor] = (self.var.state[self.ipva]