import os
import mf
import w2
import ADtest as AD
import cPickle as pickle
 
dtg='2009090412'
adps=['/w21/dat/tc/adeck/nhc/2009/aal062009.dat']
 
A=AD.Adeck(adps,dtgopt=dtg)
 
 
 
##adps=['/w21/dat/tc/adeck/local/2009/wxmap.ecm2*']
##A=AD.Adeck(adps)
##stmopt='12e'
##tstms=AD.MakeStmList(stmopt)
 
##for tstm in tstms:
##    at=A.GetAidTrk('nec2',tstm)
 
##adps=['/w21/dat/tc/adeck/esrl/2009/rtfim/*F8C','/w21/dat/tc/adeck/esrl/2009/tacc/*F?EM','/w21/dat/tc/adeck/esrl/2009/tacc/*FIM9']
##A=AD.Adeck(adps)
 
##dopickle=1
##stmopt='12e'
##stmopt='08-11w'
##stmopt='e'
##tstmids=AD.MakeStmList(stmopt)
 
##taids=['f8em','f8c','f9em','f0em','fim9']
##taids=['f8em','f8c']
 
##bt={}
##at={}
 
##stmids=[]
##aids=[]
 
 
##for tstmid in tstmids:
##    bt[tstmid]=A.GetBestTrk(tstmid)
##    if(bt[tstmid] == None): continue
##    stmids.append(tstmid)
 
##    for taid in taids:
##        for aid in A.aids:
##            if(aid == taid):
##                aids.append(aid)
##                at[aid,tstmid]=A.GetAidTrk(aid,tstmid)
 
 
 
##        continue
 
 
##if(dopickle):
##    pyppath='/tmp/vd.pyp'
##    PW=open(pyppath,'wb')
 
##aids=mf.uniq(aids)
##stmids=mf.uniq(stmids)
 
##def loaddic(dic,key,list):
##    try:
##        dic[key]=dic[key]+list
##    except:
##        dic[key]=list
 
 
##vds={}
##fes={}
##for aid in aids:
##    for stmid in stmids:
##        atkey=(aid,stmid)
##        print 'vvvvvvvvvv ',atkey
##        vd=AD.MakeVdeck(bt[stmid],at[atkey])
##        #vd.printVar(vd.veriYN)
##        #vd.printVar(vd.pod)
##        #vd.printVar(vd.blf,format='%6.2f')
##        #vd.printVar(vd.flf,format='%6.2f')
##        #vd.printVar(vd.cpacYN)
##        #vd.printVar(vd.fe,format='%7.0f')
##        #vd.printVar(vd.vme,format='%6.1f')
##        #vd.printVar(vd.btvmax,format='%4.0f')
##        #vd.printVar(vd.fcvmax,format='%4.0f')
##        #vd.printVar(vd.vme,format='%4.0f')
##        #vd.printVar(vd.ate,format='%4.0f')
##        if(dopickle): pickle.dump(vd,PW)
##        list=vd.getVDVarlist('vme',0)
##        loaddic(fes,aid,list)
 
##for aid in aids:
##    list=fes[aid]
##    (mean,amean,sigma,n)=AD.SimpleListStats(list,verb=0,undef=vd.undef)
##    print 'SSSSSSSSSS ',aid,n,mean,amean,sigma
 
##fes=AD.HomoVDdics(fes,aids)
##for aid in aids:
##    list=fes[aid]
##    (mean,amean,sigma,n)=AD.SimpleListStats(list,verb=0,undef=vd.undef)
##    print 'SSSSSHHHHH ',aid,n,mean,amean,sigma
 
##if(dopickle):
##    PW.close()
##    PW=open(pyppath)
##    vdp=pickle.load(PW)
 
 
 
 
## ropt=''
## a=AD.AtcfAdeckPaths(dtgopt='ops6',verb=0)
 
## dtgs=a.dtgs
## stmids=a.stm1ids
 
## dtgs.sort()
## stmids.sort()
 
## adpaths=[]
## for dtg in dtgs:
 
##     for stmid in stmids:
##         adeck=a.adeckpaths[dtg,stmid]
##         adpaths.append(adeck)
##         (dir,file)=os.path.split(adeck)
##         print 'dtg: ',dtg,' stmid: ',stmid,' adeck: ',adeck
 
## A1=AD.Adeck(adpaths)
 
## for aid in A1.aids:
##     print aid