output.py 文件源码

def create_diag(self,diag):
         nc = Dataset(self.diagfile,'w',format='NETCDF4')
         nc.createDimension('t',None)   # define dimensions ('None' is record dim)

         d = nc.createVariable('t','f',('t',))
         d.long_name = 'model time'

         d = nc.createVariable('kt','i',('t',))
         d.long_name = 'model iteration'

         for v in self.list_diag:
             d = nc.createVariable(v,'f',('t',))
             d.long_name = v
         nc.close()
         self.kdiag = 0

         # set up internal buffer to avoid too frequent disk access
         self.ndiags = len(self.list_diag)+2
         self.buffersize = 10
         self.buffer=zeros((self.buffersize,self.ndiags))