anchor.diff:          if (zipAble.eq.0.and.T.le.75.0) then
             call pgsci(3)
             zipAble = 1
  -          write(*,*) "Milk at start makes coffee sipable after ",time
  +          !write(*,*) "Whiskey at start makes coffee sipable after ",time
  +          time1 = time
           end if
           !write(*,*) "Plot plot"
           call pgpt (1,time,T,-1)
           if (i.eq.2) then
             call pgmove(time, T)
             T = T - 5
  +          r = r_irish
             call pgsci(4)
             call pgdraw(time, T)
             call pgsci(2)
  
positive.code:       time = 0
        zipAble = 0
        n_steps = int(t_max/delta_t)
        T = T_ini
        call pgsci(6)
        do i=1,n_steps
  
negative.code:         call pgpt(23, time_stamps, temps_black, 2)
          return
        end
  