anchor.diff:  
   !	BEGIN SOLUTION
   call system_clock(time_clicks)
  -write(*,*)real(time_clicks-start_time)/real(time_scale),'Solving navier stokes equations.'
  +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
  
positive.code: 	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
  
  !	PREPARE DATA FOR GRAPHICS
  call system_clock(time_clicks)
  write(*,*)real(time_clicks-start_time)/real(time_scale),'Preparing data for graphics.'
  allocate(u_g(x_steps,y_steps))
  allocate(v_g(x_steps,y_steps))
  allocate(x_g(x_steps))
  
negative.code: 
  !	-- 	Read in namelists
  call system_clock(time_clicks)
  write(*,"(F5.2,A,3X,A)")real(time_clicks-start_time)/real(time_scale),'s','Reading in namelist data.'
  namelist_file = file_open(adjustl(trim(namelist_file_name)))
  read(namelist_file,FLUID_PROPERTIES,IOSTAT=io_stat)
  read(namelist_file,BOUNDARY_CONDITIONS,IOSTAT=io_stat)
  