source: CMIP6dreqbuild/trunk/src/framework/audit.py @ 967

Subversion URL: http://proj.badc.rl.ac.uk/svn/exarch/CMIP6dreqbuild/trunk/src/framework/audit.py@967
Revision 967, 21.5 KB checked in by mjuckes, 3 years ago (diff)

misc

Line 
1
2from dreqPy import dreq
3from dreqPy import misc_utils
4from ingest import mipTab
5import cmchecl
6import collections, sys, os
7
8python2 = True
9if sys.version_info[0] == 3:
10  python2 = False
11
12if len(sys.argv) > 1:
13  lname = sys.argv[1]
14else:
15  lname = 'audit'
16logFarm = misc_utils.dreqLog(dir='logs')
17log = logFarm.getLog(lname)
18
19class c1(object):
20  def __init__(self):
21    self.a = collections.defaultdict( list )
22
23class checkUnits(object):
24
25  def __init__(self, cc, extra=None):
26    assert python2, 'The checkUnits module requires cf, which is not available in python 3'
27    import cf
28    if extra != None:
29      for k in sorted( extra.keys() ):
30        try:
31          a = cf.units.Units( k )
32        except:
33          log.warn( 'WARN.001.0004: bad units %s (%s)' % (k,str(extra[k])) )
34
35    for k in sorted( cc.keys() ):
36      try:
37        a = cf.units.Units( k )
38      except:
39        log.warn( 'Bad CF units: %s' % k )
40        a = None
41
42      if a != None:
43        for k2 in cc[k].a:
44          try:
45            b = cf.units.Units( k2 )
46          except:
47            b = None
48            log.warn( 'WARN.001.0003: unit invalid : %s .... %s: %s' % (k,k2,str( cc[k].a[k2]) ) )
49          if b != None and not a.equivalent(b):
50            log.warn( 'WARN.001.0002: unit mismatch: %s .... %s: %s' % (k,k2,str( cc[k].a[k2]) ) )
51
52def splitStr( s ):
53  if s == '':
54    return ['']
55  if s.find( ' ') == -1:
56    return [s,]
57  return s.split()
58
59class checkDims(object):
60
61  def __init__(self, dq):
62    """Generate list of used dimensions in ss;
63       Check against defined dimenions [todo]"""
64    ss = set()
65    self.cellm = set()
66    cc = collections.defaultdict(list)
67    duds = ('', '****', '?',)
68    for i in dq.coll['spatialShape'].items:
69      if i.dimensions != '':
70        for x in i.dimensions.split( '|' ):
71          if x in duds:
72            cc[x].append( i.uid )
73          else:
74            ss.add(x)
75    for i in dq.coll['temporalShape'].items:
76      if i.dimensions != '':
77        for x in i.dimensions.split( '|' ):
78          if x in duds:
79            cc[x].append( i.uid )
80          else:
81            ss.add(x)
82    for i in dq.coll['structure'].items:
83      self.cellm.add( i.cell_methods )
84      if i.odims != '':
85        for x in i.odims.split( '|' ):
86          if x in duds:
87            cc[x].append( i.uid )
88          else:
89            ss.add(x)
90      if i.coords != '':
91        for x in i.coords.split( '|' ):
92          if x in duds:
93            cc[x].append( i.uid )
94          else:
95            ss.add(x)
96
97    cellm = sorted( list( self.cellm ) )
98    for c in cellm:
99      log.info( 'cell_methods\t %s' % c )
100    log.info( str( sorted( list(ss) ) ) )
101    ll = set( [i.label for i in dq.coll['grids'].items] )
102    for x in ss:
103      if x not in ll:
104        log.error( 'ERROR: dimension not found: %s' % x )
105
106class auditor(object):
107  def __init__(self,dq):
108    self.dq = dq
109    self.run()
110
111  def run(self):
112    self.loadDraftNewVar()
113    self.audit_var()
114    self.audit_sn()
115    self.audit_cmv()
116    self.audit_rqv()
117    self.audit_rql()
118    self.audit_rqg()
119    self.audit_str()
120    self.audit_cm()
121
122  def audit_cm(self):
123    log.info( 'Running Cell Methods Audit' )
124    cc = cmchecl.check_cm()
125    nf = 0
126    for i in self.dq.coll['cellMethods'].items:
127      res = cc.test( i.cell_methods )
128      if not res:
129         log.error( 'ERROR.cm.00001: cannot parse cell methods: "%s" [%s]' % (i.cell_methods,i.label) )
130         nf += 1
131    if nf == 0:
132       log.info( 'INFO.cm.00001: all cell methods strings parsed by cmchecl' )
133
134    l1 = [i for i in self.dq.coll['structure'].items if i.cell_methods == ""]
135    ss = set()
136    if len( l1 ) > 0:
137      log.error( 'ERROR.cm.00010: %s blank cell methods in structure records: ' % len(l1) )
138      for i in l1:
139        log.error( 'ERROR.cm.00011: %s: %s' % ( i.label, i.title )  )
140    for i in l1:
141      if 'CMORvar' in self.dq.inx.iref_by_sect[i.uid].a:
142        for u in dq.inx.iref_by_sect[i.uid].a['CMORvar']:
143           ss.add(u)
144    if len(ss) > 0:
145      log.error( 'ERROR.cm.00020: %s CMOR vars affected ' % len(ss) )
146
147  def audit_var(self):
148    cc = collections.defaultdict( list )
149    ss = collections.defaultdict( c1 )
150    ssx = collections.defaultdict( set )
151    umap = { 'Wm-2':'W m-2', 'string':''}
152    for i in self.dq.coll['var'].items:
153      cc[i.label].append(i)
154      isn = self.dq.inx.uid[i.sn]
155      if isn._h.label != 'remarks':
156        u = umap.get( isn.units, isn.units )
157        ss[u].a[i.units].append( i.label )
158      else:
159        ssx[i.units].add( i.label )
160
161    ii = [k for k in cc if len(cc[k]) > 1]
162    log.warn( 'var.0001: %48s [%s]: %s' % ('Duplicate variable names',len( self.dq.coll['var'].items ),len(ii)) )
163    showDupVar=False
164    showDupVar=True
165    showDupVarBrief=False
166    if showDupVar:
167      for i in ii:
168        log.info(  '----- %s -----' % i )
169        for x in cc[i]:
170          log.info( '%s, %s, %s, %s' % ( x.label,x.title,x.sn,x.prov ))
171    elif showDupVarBrief:
172      lg.info( str( ii ))
173    self.cc = cc
174
175    if not python2:
176      log.warn( 'Skipping the units check .. not available in python3' )
177    else:
178      checkUnits( ss, extra=ssx )
179   
180  def audit_str(self):
181    nms = 0
182    ss = set()
183    cc = collections.defaultdict( list )
184    for i in self.dq.coll['structure'].items:
185      cc[ (i.spid, i.tmid,i.odims,i.coords,i.cmid) ].append( i )
186      if i.cmid != '__unset__' and self.dq.inx.uid[i.cmid]._h.label == 'remarks' and i.cell_methods != '':
187        ss.add(i.cmid)
188        nms += 1
189    if nms > 0:
190      log.error( 'str.00010: missing cell methods records for %s records (%s)' % (nms, len(ss) ) )
191      log.error( 'str.00011: %s' % str( sorted(list(ss)) ) )
192    else:
193      log.info( 'str.00010: all structure methods have valid cmid' )
194
195    ks = [k for k in cc if len(cc[k]) > 1]
196    log.info( 'str.00020: count of repeated space/time structure: %s' % len(ks) )
197    if len(ks) > 0:
198      oo = open( 'structureRepeats.csv', 'w' )
199      for k in ks:
200        oo.write( '\t'.join( ['NEXT', k[0], k[1], k[2], k[3], self.dq.inx.uid[k[0]].label, self.dq.inx.uid[k[1]].label] ) + '\n' )
201        for i in cc[k]:
202          oo.write( '\t'.join( ['+',i.uid, i.label, i.title, i.odims, i.coords, i.cell_methods, i.cell_measures, i.flag_meanings, i.spid ] ) + '\n' )
203      oo.close()
204
205  def saveRevisedNewVar(self,ll):
206    oo = open( 'revised_newvars.csv', 'w' )
207    for v in sorted( list( ll ) ):
208      oo.write( self.draftNewVar[v] )
209    oo.close()
210
211  def loadDraftNewVar(self):
212    self.draftNewVar = {}
213    for f in ['draft_newvars.csv','draft_newvars02.csv']:
214      ii = open( f, 'r' ).readlines()
215      ss = set()
216      for l in ii:
217        r = l.split( '\t' )
218        vn = r[0].strip()
219        assert vn not in ss, 'Duplicate variable name in %s: %s' % (f,vn)
220        ss.add(vn)
221        if vn in self.draftNewVar:
222          log.info( 'INFO.newvar.0001:  overwriting new var: %s' % vn )
223        self.draftNewVar[vn] = l
224   
225  def audit_cmv(self):
226    zmtabs = ['AERmonZ', 'E6hrZ', 'EdayZ', 'EmonZ']
227    mipt = mipTab.mipt()
228    zmtabsi = [mipt.invmap[x] for x in zmtabs]
229    ccmv = collections.defaultdict( list )
230    ccmv2 = collections.defaultdict( list )
231    ccmv4 = collections.defaultdict( list )
232    ccmv6 = collections.defaultdict( list )
233    ccp = collections.defaultdict( set )
234
235    cmlk = {}
236    zmstr = set()
237    for i in self.dq.coll['structure'].items:
238      cmlk[i.uid] = i.cell_methods
239      if i.cell_methods.find( 'longitude: mean' ) != -1 or i.cell_methods.find( 'global: mean' ):
240        zmstr.add( i.uid )
241    print 'ZMSTR', zmstr
242     
243##
244## check that vars in AERmonZ, E6hrZ, EdayZ, EmonZ have cell_methods with longitude mean.
245##
246    nms = 0
247    ss = set()
248    badRealm = set()
249    badzm = set()
250    badzmt = set()
251    for i in self.dq.coll['CMORvar'].items:
252      if i.mipTable in zmtabsi and (i.stid not in zmstr):
253        badzm.add( i.uid )
254        badzmt.add( i.mipTable )
255        log.error( 'cmv.00121: bad zonal mean variable: %s.%s, %s, %s, %s, %s' % (i.mipTable, i.label,i.uid, i.prov, i.stid, cmlk[i.stid]) )
256
257      ccp[ i.label ].add( i.positive )
258      ccmv[i.label].append(i)
259      ccmv4[ (i.frequency,i.label,i.stid)].append(i)
260      ccmv2[ (i.mipTable,i.label)].append(i)
261      if self.dq.inx.uid[i.stid]._h.label == 'remarks':
262        ss.add( i.stid )
263        nms += 1
264      else:
265        st = self.dq.inx.uid[i.stid]
266        ccmv6[ (i.frequency,i.label,st.spid,st.tmid)].append((i,st))
267      if i.modeling_realm in ['','__unset__']:
268        badRealm.add( (i.mipTable, i.label ) )
269    kkp = [k for k in ccp if len(ccp[k]) > 1]
270    if len(kkp) > 0:
271       log.error( 'cmv.00170: bad positive attributes: %s, %s' % (len( kkp ),str(kkp)) )
272    else:
273       log.info( 'cmv.00170: positive attributes OK' )
274
275    if len( badzm ) > 0:
276       log.error( 'cmv.00120: bad zonal mean variables: %s, %s' % (len( badzm ),str(badzmt)) )
277    if len( badRealm ) > 0:
278      log.error( 'cmv.00030: bad modeling realms set in %s cmv records' % len( badRealm ) )
279      cc1 = collections.defaultdict( list )
280      for tab, lab in badRealm:
281        cc1[tab].append(lab)
282      for tab in sorted( cc1.keys() ):
283        log.error( 'cmv.00031: table %s: %s' % (tab,str(sorted(cc1[tab]))) )
284    else:
285      log.info( 'cmv.00030: modeling realms set ok cmv records' )
286
287    if nms > 0:
288      log.error( 'cmv.00010: missing variable structures for %s CMORvar records, bad stids: %s' % (nms,len(ss)) )
289      log.error( 'cmv.00011: bad stids: %s' % str( sorted( list( ss ) ) ) )
290    else:
291      log.info( 'cmv.00010: variable structures for all CMORvar records present' )
292
293    ii0 = [i for i in self.dq.coll['CMORvar'].items if i.mipTable == '']
294    if len(ii0) > 0:
295      log.error( 'cmv.00020: %48s [%s]: %s' % ('Empty MIP table in CMORvar',len(self.dq.coll['CMORvar'].items),len(ii0)) )
296    else:
297      log.info( 'cmv.00020: all CMORvar records have mip table defined' )
298
299    ii = [i for i in self.dq.coll['CMORvar'].items if self.dq.inx.uid[i.vid]._h.label == 'remarks']
300    if len(ii) == 0:
301      log.info( 'cmv.00200: all CMORvar records have valid vid' )
302    else:
303      ii1 = [i for i in ii if 'requestVar' in self.dq.inx.iref_by_sect[i.uid].a]
304      log.error( 'cmv.00200: %48s [%s]: %s (requested: %s)' % ('CMORvar records with bad vid',len(self.dq.coll['CMORvar'].items),len(ii),len(ii1)) )
305
306      iix2 = []
307      iix3 = []
308      for i in ii1:
309        ok = False
310        ok3 = False
311        for u in self.dq.inx.iref_by_sect[i.uid].a['requestVar']:
312          vg = self.dq.inx.uid[ self.dq.inx.uid[u].vgid ]
313          if vg._h.label != 'remarks':
314           ok = True
315           iix2.append(i)
316           if 'requestLink' in self.dq.inx.iref_by_sect[vg.uid].a:
317             iix3.append(i)
318      log.error( 'cmv.0021: of these: valid requestVarGroup: %s; valid requestLink(s): %s' %  (len(iix2),len(iix3)) )
319      rnv = set()
320      for i in iix3:
321        if i.label in self.draftNewVar:
322          log.error( 'cmv.0125: %s, %s, %s' % (i.uid,i.label,i.title) )
323          rnv.add( i.label )
324        else:
325          log.error( 'cmv.0122: %s, %s, %s' % (i.uid,i.label,i.title) )
326
327      if len( rnv ) > 0:
328        self.saveRevisedNewVar( rnv )
329       
330         
331    ii2 = [k for k in ccmv2 if len(ccmv2[k]) > 1]
332    if len(ii2) == 0:
333      log.info( 'cmv.0030: unique CMOR variable names OK' )
334    else:
335      log.error( 'cmv.0030: non-unique CMOR variable naming for %s records' % len(ii2) )
336      log.error( 'cmv.0031: %s' % str(ii2) )
337
338    ii4 = [k for k in ccmv4 if len(ccmv4[k]) > 1]
339    ii6 = [k for k in ccmv6 if len(ccmv6[k]) > 1]
340    ii4b = list()
341    ii6b = list()
342    for k in ii4:
343      if len(ccmv4[k]) == 2 and set( [i.mipTable[-3:] for i in ccmv4[k]] ) == set(['Ant','Gre']):
344        pass
345      else:
346        ii4b.append(k)
347
348    for k in ii6:
349      if len(ccmv6[k]) == 2 and set( [i.mipTable[-3:] for i,s in ccmv6[k]] ) == set(['Ant','Gre']):
350        pass
351      else:
352        ii6b.append(k)
353
354    if len(ii4b) == 0:
355      log.info( 'cmv.0040: unique CMOR variable usage [4] OK' )
356    else:
357      log.error( 'cmv.0040: non-unique CMOR variable usage possible for %s records' % len(ii4b) )
358      log.error( 'cmv.0041: %s' % str(ii4b) )
359      kk=0
360      oo = open( 'cmv_probDuplication.csv', 'w' )
361      for k in ii4b:
362        kk += 1
363        rec = ('NEXT\t%s\t' % kk ) + str( '\t'.join( k ) )
364        print 'INFO.cmvdup.00001: ',rec
365        oo.write( rec + '\n' )
366        vid = None
367        for i in ccmv4[k]:
368          assert vid == None or i.vid == vid, 'Unexpected vid variation ...: %s%s %s' % (i.label, i.vid, vid)
369          vid = i.vid
370          rec = '+\t%s\t%s\t%s\t%s\t%s' % (i.uid,i.mipTable,i.title,i.description,i.prov)
371          oo.write( rec + '\n' )
372          print 'INFO.cmvdup.00001: ',rec
373      oo.close()
374
375    if len(ii6b) == 0:
376      log.info( 'cmv.0060: unique CMOR variable usage [6] OK' )
377    else:
378      log.info( 'cmv.0060: unique CMOR variable issues [6] %s ' % len(ii6b) )
379      kk=0
380      ccx = collections.defaultdict( set )
381      for k in ii6b:
382        ss = set()
383        freq,label,spid,tmid = k
384        for i,st in ccmv6[k]:
385          ss.add( (st.uid,st.label,st.title) )
386        tt = tuple( sorted( list( ss ) ) )
387        ccx[tt].add( (freq,label) )
388      log.info( 'cmv.0061: unique CMOR variable issues [6.1] %s ' % len(ccx.keys()) )
389      for tt in ccx:
390        ll = [t[2] for t in tt]
391        log.info( 'cmv.0062: unique CMOR variable issues [6.2] %s: %s ' % ('; '.join(ll),str(ccx[tt]) ) )
392       
393      oo = open( 'cmv_probDuplication6.csv', 'w' )
394      for k in ii6b:
395        kk += 1
396        rec = ('NEXT\t%s\t' % kk ) + str( '\t'.join( k ) )
397        oo.write( rec + '\n' )
398        vid = None
399        for i,st in ccmv6[k]:
400          assert vid == None or i.vid == vid, 'Unexpected vid variation ...: %s%s %s' % (i.label, i.vid, vid)
401          vid = i.vid
402          rec = '+\t%s\t%s\t%s\t%s\t%s\t%s' % (i.uid,i.mipTable,i.title,st.title,i.description,i.prov)
403          oo.write( rec + '\n' )
404      oo.close()
405
406    ccmv3 = collections.defaultdict( list )
407    for t in ccmv2:
408      s = set()
409      for i in ccmv2[t]:
410        id = i.uid
411        if 'varRelLnk' in self.dq.inx.iref_by_sect[id].a:
412          for x in self.dq.inx.iref_by_sect[id].a['varRelLnk']:
413            s.add( x )
414      s = list(s)
415      if len(s) > 0:
416        ccmv3[t] = s
417
418    ii5 = [k for k in ii2 if k not in ccmv3 or len(ccmv2[k]) > len(ccmv3[k])]
419    cc4 = collections.defaultdict( list )
420    for t,v in ii2:
421      cc4[t].append( v )
422
423    log.warn( '%48s [%s]: %s [%s tables]' % ('Duplicate variables in CMORvar',len(self.dq.coll['CMORvar'].items),len(ii2), len( cc4.keys() ) ) )
424    showAllCmvErrors=False
425    showAllCmvErrors=True
426    if showAllCmvErrors:
427       rats = ['uid','label','title','description','frequency','mipTable','stid','prov','provNote']
428       oocm = open( 'cmvDup.csv', 'w' )
429       c1 = collections.defaultdict( list )
430       ##for i in ii2:
431         ##c1[i.prov].append( i.label )
432       ##log.warn( 'No MIP var: %s ... %s' % (str( sorted( [i.label for i in ii if i.label in self.cc] ) ), str(c1)))
433       for k in sorted( cc4.keys() ):
434         log.warn(  '%16s::  %s' % (k,str(sorted(cc4[k]))) )
435         for v in sorted(cc4[k]):
436             sv = set()
437             for i in ccmv2[ (k,v) ]:
438                sv.add(i.vid)
439
440             if len(list(sv)) == 1:
441               vun = True
442               vv = self.dq.inx.uid[ list(sv)[0] ]
443               oocm.write( '@VAR:\t\t\t%s\t%s\t%s\n' % (vv.label,vv.title,vv.description) )
444             else:
445               vun = False
446
447             for i in ccmv2[ (k,v) ]:
448                stt = self.dq.inx.uid[i.stid].title
449                if 'Temporal mean' in stt:
450                   if i.frequency == '3hr':
451                     i.title += ' [3 hour mean]'
452                xtra = ''
453                oocm.write( '\t\t' + '\t'.join( [str(i.__dict__[x]) for x in rats] ) + '\t%s%s\n' % (stt,xtra) )
454                if not vun:
455                  vv = self.dq.inx.uid[i.vid]
456                  xtra = '\t\t\t\t%s\t%s\t%s\n' % (vv.label,vv.title,vv.description)
457                  oocm.write(xtra)
458             oocm.write( '####\n' )
459       oocm.close()
460    t2_ii = ii[:]
461    self.ccmv = ccmv
462
463  def audit_sn(self):
464    cc = collections.defaultdict( list )
465    cc2 = collections.defaultdict( int )
466    nrem = 0
467    snm = set()
468    snmm = collections.defaultdict( int )
469    for i in self.dq.coll['var'].items:
470      if 'CMORvar' in self.dq.inx.iref_by_sect[i.uid].a:
471        srq = set()
472        for u in self.dq.inx.iref_by_sect[i.uid].a['CMORvar']:
473          if 'requestVar' in self.dq.inx.iref_by_sect[u].a:
474            for uu in self.dq.inx.iref_by_sect[u].a['requestVar']:
475              rqv = self.dq.inx.uid[uu]
476              srq.add( rqv.priority )
477        if len(srq) > 0:
478          pm = min( srq )
479           
480          kk = i.sn
481          if self.dq.inx.uid[kk]._h.label == 'remarks':
482            nrem += 1
483            snm.add(kk)
484           ## if pm == 1:
485            print 'MISSING STANDARD NAME: ',i.label, i.title, i.sn, i.provmip
486            snmm[pm] += 1
487            cc2[i.provmip] += 1
488          if i.procnote != []:
489            ll = sorted( i.procnote )
490            kk += ':' + '|'.join(ll )
491          cc[kk].append( i.uid )
492
493    log.warn( 'sn.001: Number of missing standard names: %s [for %s variables] {%s}' % (len(snm),nrem,str(snmm)) )
494    log.warn( 'sn.002: Number of missing standard names by MIP: %s ' % (str(cc2)) )
495    ll = []
496    for k in cc:
497      if len( cc[k] ) > 1:
498        ll.append( k )
499
500    log.warn( 'sn.003: Duplicate use of standard names: %s' % len(ll) )
501    oo = open( 'audit_sn_repeats.csv', 'w' )
502    for k in sorted( ll ):
503      orc = [k,]
504      l1 = cc[k][:]
505      l2 = []
506      l1b = []
507      for u in cc[k]:
508        i = self.dq.inx.uid[u]
509        if i.prov == "CMIP6 endorsement [SIMIP]":
510          l2.append(u)
511        else:
512          l1b.append(u)
513      if len( l2) > 1:
514        log.error( 'WARN.simip.00001: repeat use of standard name in SIMIP: %s:: %s' % (k,str(l2)) )
515      if len(l2) == 1:
516        orc.append( '*' )
517      else:
518        orc.append( '' )
519
520      nnn = 0
521      for u in l2 +l1b:
522        i = self.dq.inx.uid[u]
523        if "CMORvar" not in self.dq.inx.iref_by_sect[u].a:
524           orc += [i.label +'**', i.title]
525        else:
526          nn = 0
527          for uu in self.dq.inx.iref_by_sect[u].a['CMORvar']:
528            cmv = self.dq.inx.uid[uu]
529            if 'requestVar' in self.dq.inx.iref_by_sect[uu].a:
530              nn += 1
531          if nn == 0:
532             orc += [i.label + '*', i.title]
533          else:
534             nnn += 1
535             orc += [i.label, i.title]
536         
537      if nnn > 1:
538        oo.write( '\t'.join( orc ) + '\n' )
539    oo.close()
540     
541  def audit_rqv(self):
542    ixx = [i for i in self.dq.coll['requestVar'].items if i.vid in self.dq.inx.uid]
543    if len(ixx) < len( self.dq.inx.uid ):
544       log.error( 'SEVERE.var.00301: request vars with bad links ....' )
545   
546    ii = [i for i in ixx if self.dq.inx.uid[i.vid]._h.label == 'remarks']
547    ii2 = [i for i in ii if i.label not in self.cc]
548    ii3 = [i.vgid for i in ii if self.dq.inx.uid[i.vgid]._h.label != 'remarks']
549##ii4: set of invalid requestVarGroup records
550    ii4 = [i for i in ii3 if 'requestLink' not in self.dq.inx.iref_by_sect[i].]
551######
552    ii5Info = 'Set of bad requestVar records linking to a valid requestVarGroup records.'
553    ii5 = [i for i in ii if i.vgid not in ii4]
554##ii6: set of bad requestVar records linking to a valid requestVarGroup records and no variable name match.
555    ii6 = [i for i in ii5 if i.label not in self.cc]
556    log.error( '%48s [%s]: %s [%s, %s, %s; %s]' % ('Bad variable links in requestVar',len(self.dq.coll['requestVar'].items),len(ii),len(ii2), len(ii5), len(ii6), len(ii4)) )
557    log.info ( str( ii5Info ))
558    svg = set()
559    for i in ii5:
560      log.info( 'ii5: %s' % str([i.label, i.title, i.uid, i.mip, self.dq.inx.uid[i.vgid].title] ) )
561      svg.add( i.vgid )
562    for i in ii2:
563      if i.vgid in self.dq.inx.uid and self.dq.inx.uid[i.vgid]._h.label != 'remarks':
564        xx = '%s [%s]' % (self.dq.inx.uid[i.vgid].title, i.vgid)
565      else:
566        xx = i.vgid
567      log.info( 'ii2: %s' % str([i.label, i.title, i.uid, i.mip, xx] ) )
568
569    log.warn( 'rvg.0010: requestVarGroups with bad request vars:' )
570    for u in sorted( list(svg) ):
571      log.warn( 'rvg.0011: %s: %s' % (u, self.dq.inx.uid[u].title ))
572
573    log.info( '    [# records]: # broken [var name not known, link to valid group, valid group and no var name; valid request groups' )
574    showAllRqvErrors=False
575    if showAllRqvErrors:
576      log.info( '%s\n%s' % (str( [i.label for i in ii if i.label in self.cc]), str( [i.label for i in ii2] ) ) )
577
578    ii1 = [i for i in ii if i.label in self.ccmv]
579    listBadRequestVar=False
580    if listBadRequestVar:
581      for i in ii1:
582        log.info( str( i.label,i.mip,i.table,[x.frequency for x in self.ccmv[i.label]] ) )
583
584  def audit_rql(self):
585    ii = [i for i in self.dq.coll['requestLink'].items if self.dq.inx.uid[i.refid]._h.label == 'remarks']
586    log.info( '%48s: %s (from %s)' % ('ERROR.rql.001: Bad request group links in requestLink',len(ii), len(self.dq.coll['requestLink'].items)) )
587    rql_detail=True
588    if rql_detail:
589      for i in ii:
590         log.info( '%s: %s, %s [%s]' % ('ERROR.rql.002: ',i.mip,i.label,i.uid) )
591
592  def audit_rqg(self):
593    ii = [i for i in self.dq.coll['requestVarGroup'].items if len( dq.inx.iref_by_sect[i.uid].a['requestVar'])  == 0]
594    ii1 = [i for i in ii if len( dq.inx.iref_by_sect[i.uid].a['requestLink'])  != 0]
595    log.info( '%48s: %s/%s (from %s)' % ('Request groups with no request links:',len(ii) - len(ii1),len(ii1), len(self.dq.coll['requestVarGroup'].items)) )
596    listEmptyRequestVarGroups=True
597    if listEmptyRequestVarGroups:
598      for i in ii1:
599        log.info ( 'INFO.rqgempty.0001: %s: %s, %s [%s]' % (i.label, i.title, i.mip, i.uid) )
600
601dq = dreq.loadDreq(manifest=None)
602a = auditor( dq )
603cd = checkDims( dq )
Note: See TracBrowser for help on using the repository browser.