diff --git a/parallel.py b/parallel.py
index 567bab8..b4da8ba 100644
--- a/parallel.py
+++ b/parallel.py
@@ -246,7 +246,10 @@ class PPP:
         dim = (ie-ib+1+2*self.nghx,
                je-jb+1+2*self.nghy,
                ke-kb+1+2*self.nghz)
-        return Field3d.allocate(dim,origin,self.spacing,fill=fill,pseudo=pseudo,dtype=dtype)
+        if pseudo:
+            return Field3d.pseudo_field(dim,origin,self.spacing)
+        else:
+            return Field3d.allocate(dim,origin,self.spacing,fill=fill,dtype=dtype)
 
     def qcriterion(self,key_out=None,from_pressure=False,keep_derivatives=False):
         '''Computes Q-criterion'''