Pression atmosphérique

rem ============================================================================
rem  Evolution de la pression atmosphérique moyenne avec l altitude
rem                Par Papydall
rem ============================================================================
' Calcul de la pression atmosphérique P en fonction de l'altitude Z
' avec la formule internationnale du nivellement barométrique.
' On prend comme altitude de référence Z0, le niveau moyen de la mer.
' On prend pour l'atmosphère un état moyen défini par l'atmosohère normalisée
' type OACI (Organisation de l'Aviation Civile Internationale), c-à-d :
' Température 15°C = 288.15 K;
' Pression 1013.25 hPa;
' Gradient vertical de température = 0.65 K pour 100 m.
rem  La formule : P(Z) = 1013.25*power(((1-(0.0065*z)/288.15)),5.255)
' donne la pression exprimée en hectopascals (hPa) à l'altitude Z exprimée en mètres.
rem ============================================================================
Init()
Tracer_Repere()
Tracer_Courbe()
end
rem ============================================================================
SUB Init()
    dim p : ' Pression atmosphérique
    width 0,700 : height 0,600 : color 0,0,0,0
    picture 10 : full_space 10 : print_target_is 10 : 2d_target_is 10 : color 10,0,0,0 : font_bold 10
    alpha 20 : top 20, 500 : left 20,10 : font_bold 20 : font_size 20,14 : font_color 20,255,0,255
    caption 20,"Evolution de la pression atmosphérique moyenne avec l'altitude" + chr$(13)+"entre 0 et 20 000 m"
END_SUB
rem ============================================================================
SUB P(Z)
    p = 1013.25*power((1-(0.0065*z)/288.15),5.255)
END_SUB
rem ============================================================================
SUB Tracer_Repere()
    dim_local i
    2d_pen_color 255,255,0  : 2d_pen_width 2
    print_locate 10,05 : print " Altitude en Km "
    2d_line 50,20,50,450 : 2d_poly_to 650,450 : print_locate 35,450 : print "0"
    for i = 1 to 20
        2d_line 50,450-20*i,600,450-20*i : print_locate 33,442-20*i : print str$(i)
    next i
    print_locate 585,460 : print " Pression en hPa "
    for i = 100 to 1000 step 100
        2d_line 50+i/2,450,50+i/2,50 : print_locate 40+i/2,460 : print str$(i)
    next i
    print_locate 220,442-180 : print " Nécessité de respirer de l'oxygène "
    print_locate 140,442-260 : print " Nécessité d'une préssurisation "
    print_locate 080,442-400 : print " Ebulution du sang "
END_SUB
rem ============================================================================
SUB Tracer_Courbe()
    dim_local z,x,y
    2d_pen_color 0,255,255 : p(0) : 2d_poly_from 50+p/2, 450
    for z = 0 to 21000 step 10
        p(z) : x = 50 + p/2 : y = 450-z/50 : 2d_poly_to x,y
    next z
END_SUB
rem ============================================================================