### read every100m spectrum power to output data Dspectrum bin #### you have to calc every 100m spectrum(about each particle) before exec this program. 1000m高度は別途計算 ===> modified to include. 12/27'8 include 1000m-height INPUT: read file name fixed "/home/koshida/jobs/2003/mrr/prg/melting/E100/every100.csv" OUTPUT: out_eta.csv --> sent to 'every100NN6.xls' FORMAT nyu1(0.10m/s) nyu2(0.29m/s) ............ nyu64(12.10m/s) ------------------------------------------------ 2/24'9 Two file must be fit for analysis date. nddd.csv ept.txt :: copy from Excel file. ------------------------------------------------ 11/24'8 out_eta_50m.csv :: every 50m pitch 9/15'8- fort.85 -- fort.83 fort.83 col-AL;ratioM/R,,wrong...col-Q repair mie2r ------------------------------------------------ 6/26'8-- add reverse order condition 6/28'8 --- now calclate 31-bin --- had better to use64bin. (corresponding to 2004/12/29case) 2008年07月04日 (金) 05時37分56秒 JST - thunder.. output ::out_eta.csv ;;intersected data input ::E100/every100.csv;; calculated previously nddd.csv ;; N(D) -- respect to nyu. N(D)(m-3/(0.19m/s)) convert N(D) [m-3/mm] following steps N(D)*0.19[m/s] -> devided D_rain.step every100NN2.xls effective from 12 to 42 @2004.12.29 case. S_ze ok! but S_Ze_nyu NG. error in dD/dV 8/25'8 mie/rayleigh calb 9/13'8 input file.name check ~/mrr/prg/melting/E100/every100.csv "every100.csv" file.format slightly changed -> source re-write "nddd.csv" re-stored 9/14'8 calc.result=(2.19, 3.1549) @100m request (2.19, 146670) |--> col-O;etae: 8.706536 should be 5.892963823 in fort.83 |--> col-I;ndd 483.31 should be 327.125206 in fort.84 REPAIR :: dD exter porate should be favorable. but Notime.