function sofie_at_zmax, z,data,zmax,zcut ;---------------------------------------------------------------- ; Routine extracts data at PMC zmax, where zmax > zcut ; ; Input: ; z........altitude (km), arr(z) ; data.....the data array, arr(z,event) ; zmax.....altitude of peak PMC mass density, arr(event) ; zcut.....consider only PMCs with zmax >= zcut (km), scalar ; ; Output: ; data_zm..the data at Zmax, arr(event) ; ; Example call for water vapor (h2o): h2o_zm = sofie_at_zmax(z,h2o,zmax,80.) ; ; Source: Mark Hervig, GATS ; ;---------------------------------------------------------------- nevent = n_elements(data(0,*)) data_zm = fltarr(nevent) for i = 0,nevent-1 do begin if (zmax(i) ge zcut) then begin m = min( abs(z - zmax(i)), im) data_zm(i) = data(im,i) endif endfor ;- done return, data_zm end