diff --git a/src/gotm/gotm.F90 b/src/gotm/gotm.F90
index 2e0ee875..161ecee0 100644
--- a/src/gotm/gotm.F90
+++ b/src/gotm/gotm.F90
@@ -731,6 +731,7 @@ subroutine integrate_gotm()
 !BOC
    if (i == 0 .and. .not. restart) then
       LEVEL1 'saving initial conditions'
+      call do_diagnostics(nlev)
       call output_manager_save(julianday,int(fsecondsofday),int(mod(fsecondsofday,_ONE_)*1000000),0)
       STDERR LINE
       LEVEL1 'integrate_gotm'