anchor.diff:          "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)
  
positive.code:   ! - The model then solves for the tendencies at the current step
    !   before solving for the fields at the next time step.
  
    print "(A)", "Initialized"
  
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!! MAIN LOOP OF THE MODEL STARTS HERE                                  !!!
  
negative.code:   double precision, dimension(:,:),   allocatable :: wind_x
    double precision, dimension(:,:),   allocatable :: wind_y
  
    ! Time
    integer*8 :: start_time, last_report_time, cur_time
  
    allocate(dhdt(0:nx+1, 0:ny+1, layers))
    allocate(dhdtold(0:nx+1, 0:ny+1, layers))
    allocate(dhdtveryold(0:nx+1, 0:ny+1, layers))
  