diff --git a/src/outcdf.f90 b/src/outcdf.f90
index 98c0ba6..f1a8418 100644
--- a/src/outcdf.f90
+++ b/src/outcdf.f90
@@ -327,8 +327,10 @@ subroutine write_netcdf(header)
 
     if (debug) print *, 'Defined TGLF variables'
 
-    call neasyf_write(ncid, "rpts", rpts, dim_ids=dimids2d, units='m')
-    call neasyf_write(ncid, "zpts", zpts, dim_ids=dimids2d, units='m')
+    call neasyf_write(ncid, "rpts", rpts, dim_ids=dimids2d, units='m', &
+         long_name="Flux surface R points")
+    call neasyf_write(ncid, "zpts", zpts, dim_ids=dimids2d, units='m', &
+         long_name="Flux surface Z points")
     call neasyf_write(ncid, "bppts", bppts, dim_ids=dimids2d, units='T')
 
     call write_input_parameters(ncid)