#
# power_vs_z
#
# get
#   field power vs z 
#

if test $# -ne 1 
then
  echo usage: `basename $0` ELFIN_outputfile
  exit
fi

dlen=`egrep 'Length of dispersion section' $1 |awk -F"=" '{print $2}'|awk '{printf "%e", $1}'`

grep '=' $1 |\
egrep 'z = |field    ' | egrep 'z = |power' | egrep 'z = |(W)' |\
awk -F"=" '
   BEGIN{ 
     und = 0 
     printf "# %16s %16s\n", "z", "Field_Power"
   }
   /z = / { z = $2 ; ok = 1 }
   /(W)/  {
            if( ok == 1 )
            { 
              if( z == 0 )
              {  
                und = und + 1 
                if( und == 2 ) { z1 = z0 + '${dlen}' }
              }  
              if( und == 2 ){ z = z1 + z }   
              power = $2  
              printf "%16.5e %16.5e\n" , z, power
              z0 = z 
            }
          } '