jsv5D.,
'STD_ATM
'Written by HEW 4JAN07
'Computes middle latitude standard atmosphere P & T to 40 km
'Every 1/2 km


P0 = 1013.25 'Standard atm sfc pressure
T0 = 288.15  'Standard atm sfc
LR0  =  6.5  'Lapse rate from sfc to 11 km
LR20 = -1.0  'Lapse rate from 20 to 32 km
LR32 = -2.8  'Lapse rate from 32 to 47 km
GMR = 34.1632

'Loop from surface to 11 km. Lapsing 6.5 C/km
for i = 1 to 23 do
cell("Z",i) = 0.5*(i-1)
cell("T",i) = T0 - LR0*cell("Z",i)
cell("P",i) = P0*(cell("T",i)/T0)^(GMR/LR0)
end for 


'Loop from 11.5 km to 20 km. Isothermal
for j = 24 to 41 do
cell("Z",j) = 0.5*(j-1)
cell("T",j) = cell("T",23)
cell("P",j) = cell("P",23)*exp(-GMR*(cell("Z",j)-11)/cell("T",23))
end for


'Loop from 20.5 km to 32 km. Lapsing -1 C/ km
for k = 42 to 65 do
cell("Z",k) = 0.5*(k-1)
cell("T",k) = cell("T",41) - LR20*(cell("Z",k)-20)
cell("P",k) = cell("P",41)*(cell("T",k)/cell("T",41))^(GMR/LR20)
end for 


'Loop from 32.5 km to 47 km. Lapsing -2.8 C/ km
for m = 66 to 95 do
cell("Z",m) = 0.5*(m-1)
cell("T",m) = cell("T",65) - LR32*(cell("Z",m)-32)
cell("P",m) = cell("P",65)*(cell("T",m)/cell("T",65))^(GMR/LR32)
end for 


'Loop from 47.5 km to 51 km. Isothermal
for n = 96 to 103 do
cell("Z",n) = 0.5*(n-1)
cell("T",n) = cell("T",95)
cell("P",n) = cell("P",95)*exp(-GMR*(cell("Z",n)-47)/cell("T",95))
end for