diff --git a/geom-exp.cpp b/geom-exp.cpp
index 1c5ed51e..44769f69 100644
--- a/geom-exp.cpp
+++ b/geom-exp.cpp
@@ -1252,6 +1252,12 @@ int read_geom_args() {
     shift(); gp::param.second = argi();
     set_variation(eVariation::goldberg);
     }
+  else if(argis("-gph")) {
+    PHASEFROM(2);
+    shift(); int x = argi();
+    shift(); int y = argi();
+    gp::whirl_set({x, y});
+    }
   else if(argis("-unrectified")) {
     PHASEFROM(2);
     set_variation(eVariation::unrectified);
diff --git a/goldberg.cpp b/goldberg.cpp
index 9e96939d..a365b0bb 100644
--- a/goldberg.cpp
+++ b/goldberg.cpp
@@ -761,7 +761,7 @@ EX namespace gp {
     return eVariation::goldberg;
     }
   
-  void whirl_set(loc xy) {
+  EX void whirl_set(loc xy) {
     xy = internal_representation(xy);
     if(xy.second && xy.second != xy.first && nonorientable) {
       addMessage(XLAT("This does not work in non-orientable geometries"));