anchor.diff:  	cont_rms_old = cont_rms
   enddo
   call system_clock(time_clicks)
  -write(*,*)real(time_clicks-start_time)/real(time_scale),'-->Continuity converged to ',cont_rms
  +write(*,"(F5.2,A,3X,A,ES10.3)")real(time_clicks-start_time)/real(time_scale),'s','-->Continuity converged to ',cont_rms
   
   !	PREPARE DATA FOR GRAPHICS
   call system_clock(time_clicks)
  -write(*,*)real(time_clicks-start_time)/real(time_scale),'Preparing data for graphics.'
  +write(*,"(F5.2,A,3X,A)")real(time_clicks-start_time)/real(time_scale),'s','Preparing data for graphics.'
   allocate(u_g(x_steps,y_steps))
   allocate(v_g(x_steps,y_steps))
   allocate(x_g(x_steps))
  
positive.code: 
  !	SETUP MESH
  call system_clock(time_clicks)
  write(*,*)real(time_clicks-start_time)/real(time_scale),'Setting up mesh.'
  dx = x_size/real(x_steps,wp)
  dy = y_size/real(y_steps,wp)
  
  
negative.code: 
  !	BEGIN SOLUTION
  call system_clock(time_clicks)
  write(*,"(F5.2,A,3X,A)")real(time_clicks-start_time)/real(time_scale),'s','Solving navier stokes equations.'
  cont_rms_old = 0.0_wp
  do k=1,3000
  	!	--	MASS FLOW RATES
  