import numpy as np
import commands
import re

tunit=6.88027e+14/3.156e7/1e6

for n in [13,18,23,28,33,43,53,63,68,73]:

 fi=open('RD%04i/normal_galaxies_properties_luminosity_%04i'%(n,n),'r')
 lines1 = fi.readlines()
 fi.close()

 fo=open('RD%04i/normal_galaxies_properties_luminosity_noheader_%04i'%(n,n),'w')

 region=lines1[0].split()[1]
 result=[float(x) for x in re.findall("[0-9.]+",lines1[0])]

 for line1 in lines1[2:]:
    line1=line1.rstrip('\r\n') 
    fo.write(line1+" %s %4.1f %12.8f\n"%(region,result[0],result[1]) ) 
   
 fo.close()