Jump with parachute

globals [
  Fluft
  acceleration
  tyngde-acceleration
  NLoop
  A
  dt
  y
  vy
]

breed [men man]

to setup
  clear-all
  set tyngde-acceleration 9.82               ;Alle fysiske størrelse er i SI-enheder
  set acceleration tyngde-acceleration
  set NLoop 100
  set A 1.3 * 0.3
  set dt 1.0 / NLoop
  set y startshøjde
  set vy 0
  create-men 1
  [
    setxy 0 y ; start position
    set shape "person"
    set color red
    set size 8
  ]
  reset-ticks
  paint-sky
end 

; beregn ny position og fart

to go
  if (ticks > udløsningstid )
  [
    set A 10 * 5
    ask men
    [ set shape "parachute"
      set size 15
      set heading 0
    ]
  ]
  repeat Nloop [              ; beregningerne laves med et kortere tidsinterval end skærmen opdateres
       set Fluft 0.5 * 1.2 * A * vy * vy              ; Fluft = 0.5 · rho · A · cw · v^2
       set acceleration (tyngde-acceleration - Fluft / masse)
       set vy ( vy - acceleration * dt )
       set y ( y + vy * dt )
      ]
  if (y < 0 )     ; stop inden manden rammer jorden
    [
      ask patches with [ pycor < 0 ] [ set pcolor green ]
      print ticks
      stop
    ]
  ask men [setxy xcor ( y ) ]       ; hvis ikke jorden er ramt endnu flyttes manden
  paint-sky
  tick      ; ticks tæller antal sekunder
end 

to paint-sky
    ask patches [ set pcolor 100 - ( min list  ( y / 1000)  10 )]
end