vname(1)='p'
        vname(2)='t'
        vname(3)='q'
        vname(4)='u'
        vname(5)='v'
        vname(6)='w'
        vname(7)='w2'
        vname(8)='wd'
        vname(9)='td'
        vname(10)='gust'
        vname(11)='vis'
        vname(12)='pblh'
        vname(13)='dist'
        vname(14)='wspd10m'
        vname(15)='td2m'
        vname(16)='mxtm'
        vname(17)='mitm'
        vname(18)='pmsl'
        vname(19)='howv'
        vname(20)='tcamt'
        vname(21)='lcbas'


        do n=1,nvars
           dfactor=1. ; if (n==11) scale=1./1609.344

           open (12,file=trim(vname)'_avgs_of_hdvar.txt',form='formatted')
           write(12,'(a11,a15,a1)') 'starttime="',trim(starttime),'"'
           write(12,'(a9,a1,a4,a1)') 'filetype=','"',trim(cpos1neg9),'"'

           write(12,'(a22,f6.2,a1)') trim(vname)'_bckgrmseavg="AVG =',vrmse1_hdvaravg(n)*scale, '"'
           write(12,'(a21,f6.2,a1)') trim(vname)'_anlrmseavg="AVG =',vrmse2_hdvaravg(n)*scale, '"'
           write(12,'(a23,f6.2,a1)') trim(vname)'_cressrmseavg="AVG =',vrmse3_hdvaravg(n)*scale, '"'

           write(12,'(a22,f6.2,a1)') trim(vname)'_bckgbiasavg="AVG =',vbias1_hdvaravg(n)*scale, '"'
           write(12,'(a21,f6.2,a1)') trim(vname)'_anlbiasavg="AVG =',vbias2_hdvaravg(n)*scale, '"'
           write(12,'(a23,f6.2,a1)') trim(vname)'_cressbiasavg="AVG =',vbias3_hdvaravg(n)*scale, '"'

           write(12,'(a22,f6.2,a1)') trim(vname)'_bckgabseavg="AVG =',vabse1_hdvaravg(n)*scale, '"'
           write(12,'(a21,f6.2,a1)') trim(vname)'_anlabseavg="AVG =',vabse2_hdvaravg(n)*scale, '"'
           write(12,'(a23,f6.2,a1)') trim(vname)'_cressabseavg="AVG =',vabse3_hdvaravg(n)*scale, '"'

           close(12)
        enddo