diff --git a/src/bmi_main_lgar.cxx b/src/bmi_main_lgar.cxx
index 37b0ac9..26c0ed8 100644
--- a/src/bmi_main_lgar.cxx
+++ b/src/bmi_main_lgar.cxx
@@ -163,7 +163,7 @@ int main(int argc, char *argv[])
 
   // do final mass balance
   model_state.global_mass_balance();
-
+  model_state.Finalize();
   if (outdata_fptr) {
     fclose(outdata_fptr);
     fclose(outlayer_fptr);