set ansupdn=n        # For spin-pol case plot both spins
    set ansneg=n        # for spin-pol case plot spin-dn negative
    set color=1        # BW/color setting
    set units=1        # units in Ry/eV
    set lsize=24         # pointsize for labels
    set count=4          # plotting N cases  (1-2/4)
    set lineint_batch=( 1 2 3 4 )  # set lines for plotting (must be N values)
                            # set N properties
    set label1=" total DOS"
    set linety1=1
    set linew1=1
    set label2="Am tot"
    set linety2=2
    set linew2=1
    set label3="N tot"
    set linety3=3
    set linew3=1
    set label4=" interstitial"
    set linety4=4
    set linew4=1
                            # repeat this N times for label2,3,4 ...
    set xmin=-5         
    set xmax=5          
    set ymin=
    set ymax=

    set out=/home/quanyin/.w2web/rupc04/tmp/365103-19115     # file to be used for plotting
    set noview=y             # set to skip preview and directly print