diff --git a/src/gen_topo.f90 b/src/gen_topo.f90
index cae8e36..e1a208b 100644
--- a/src/gen_topo.f90
+++ b/src/gen_topo.f90
@@ -215,12 +215,12 @@ subroutine gen_topo(topo_file, out_file, grid_file, istripolar, offset)
     call handle_error(nf90_def_var(ncid_topo, 'depth_med', nf90_float, dids_topo, depth_med_id))
     call handle_error(nf90_def_var(ncid_topo, 'depth_all_med', nf90_float, dids_topo, depth_all_med_id))
     call handle_error(nf90_def_var(ncid_topo, 'sea_area_fraction', nf90_float, dids_topo, frac_id))
-    call handle_error(nf90_def_var(ncid_topo, 'geolon_t', nf90_float, dids, geolon_id, deflate_level=1, shuffle=.true.))
+    call handle_error(nf90_def_var(ncid_topo, 'geolon_t', nf90_float, dids_topo, geolon_id, deflate_level=1, shuffle=.true.))
     call handle_error(nf90_put_att(ncid_topo, geolon_id, 'long_name', 'tracer longitude'))
     call handle_error(nf90_put_att(ncid_topo, geolon_id, 'units', 'degrees_E'))
     call handle_error(nf90_put_var(ncid_topo, geolon_id, x_t))
 
-    call handle_error(nf90_def_var(ncid_topo, 'geolat_t', nf90_float, dids, geolat_id, deflate_level=1, shuffle=.true.))
+    call handle_error(nf90_def_var(ncid_topo, 'geolat_t', nf90_float, dids_topo, geolat_id, deflate_level=1, shuffle=.true.))
     call handle_error(nf90_put_att(ncid_topo, geolat_id, 'long_name', 'tracer latitude'))
     call handle_error(nf90_put_att(ncid_topo, geolat_id, 'units', 'degrees_N'))
     call handle_error(nf90_put_var(ncid_topo, geolat_id, y_t))