anchor.diff:    ! - The model then solves for the tendencies at the current step
     !   before solving for the fields at the next time step.
   
  -  print "(A)", "Initialized"
  +  cur_time = time()
  +  print "(A, I0, A)", "Initialized in ", cur_time - start_time, " seconds."
  +  last_report_time = cur_time
   
     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     !!! MAIN LOOP OF THE MODEL STARTS HERE                                  !!!
  
positive.code:         wind_x, wind_y, nx, ny, layers, &
          n, nwrite, avwrite, RedGrav, DumpWind)
  
    end do
  
    print "(A, i10.10)", "Run finished at time step ", n
    return
  end subroutine model_run
  
  
negative.code:         "Running an n-layer configuration of size ", &
          nx, "x", ny, "x", layers, " by ", nTimeSteps, " time steps."
    end if
    last_report_time = start_time
  
    nwrite = int(dumpFreq/dt)
    avwrite = int(avFreq/dt)
  