anchor.diff:        time = 0
         zipAble = 0
         n_steps = int(t_max/delta_t)
  +      r = r_black
         T = T_ini
         call pgsci(6)
         do i=1,n_steps
  
positive.code:         call pgpt(23, time_stamps, temps_black, 2)
          return
        end
  
        function getR(t_stamps, temps_coffee)
        real temps_coffee(23)
        getR = 1
        end
  
negative.code:         if (zipAble.eq.0.and.T.le.75.0) then
            call pgsci(3)
            zipAble = 1
            !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)
  