#!/usr/bin/env ruby

require "narray"
#include NumRu

CpDry = 820.0
Grav  = 8.87
GammaD = Grav/CpDry

temp   = NArray.float(26)
height = NArray.float(26)

temp[0]=178.76712328767124 ; height[0]=94.54545454545455
temp[1]=180.82191780821918 ; height[1]=89.35064935064935
temp[2]=184.93150684931507 ; height[2]=84.5021645021645
temp[3]=193.15068493150685 ; height[3]=80.0
temp[4]=203.4246575342466  ; height[4]=75.84415584415585
temp[5]=217.8082191780822  ; height[5]=71.34199134199135
temp[6]=236.3013698630137  ; height[6]=67.18614718614718
temp[7]=252.73972602739724 ; height[7]=64.06926406926407
temp[8]=273.28767123287673 ; height[8]=60.60606060606061
temp[9]=295.8904109589041  ; height[9]=56.45021645021645
temp[10]=320.54794520547944; height[10]=52.98701298701299
temp[11]=343.1506849315068 ; height[11]=49.523809523809526
temp[12]=365.75342465753425 ; height[12]=46.40692640692641
temp[13]=392.4657534246575  ; height[13]=42.5974025974026
temp[14]=419.17808219178085 ; height[14]=38.78787878787879
temp[15]=445.8904109589041 ; height[15]=35.324675324675326
temp[16]=472.6027397260274 ; height[16]=31.861471861471863
temp[17]=499.3150684931507 ; height[17]=28.3982683982684
temp[18]=528.0821917808219 ; height[18]=24.935064935064936
temp[19]=556.8493150684931 ; height[19]=21.125541125541126
temp[20]=587.6712328767123 ; height[20]=17.316017316017316
temp[21]=614.3835616438356 ; height[21]=13.852813852813853
temp[22]=645.2054794520548 ; height[22]=10.043290043290042
temp[23]=669.8630136986301 ; height[23]=6.926406926406926
temp[24]=698.6301369863014 ; height[24]=3.463203463203463
temp[25]=725.3424657534247 ; height[25]=0.3463203463203463

for i in 0..24
  tempbar = (temp[i]+temp[i+1])/2
  gamma   = (temp[i]-temp[i+1])/(height[i]-height[i+1])*1e-3
  n2      = Grav/tempbar * (GammaD + gamma)
  printf( "%d %f %e %e\n", i, (height[i]+height[i+1])/2, gamma, n2 )
end