--- trunk/dyn3d/dynetat0.f	2014/03/05 14:57:53	82
+++ trunk/dyn3d/dynetat0.f	2014/03/06 17:35:22	85
@@ -146,6 +146,9 @@
 
     call NF95_INQ_VARID (ncid, "ps", varid)
     call NF95_GET_VAR(ncid, varid, ps)
+    ! Check that there is a single value at each pole:
+    call assert(ps(1, 1) == ps(2:, 1), "dynetat0 ps north pole")
+    call assert(ps(1, jjm + 1) == ps(2:, jjm + 1), "dynetat0 ps south pole")
 
     call NF95_CLOSE(ncid)