def __init__(self, col='SW', subcol='FTOA'): ''' col is 'LW' or 'SW' subcol for SW: 'CWV', 'FABS', 'FALL', 'FCLR', 'FTOA', 'SALB', 'TCLD' for LW: 'sfc_down', 'sfc_up', 'toa_up', 'clr_sfc_down', 'clr_sfc_up', 'clr_toa_up' ''' name = 'SRB_' + col undef = -999 dir = '/discover/nobackup/projects/gmao/share/dao_ops/verification/SRB' rr = rrule.rrule(rrule.MONTHLY, dtstart=datetime.date(1983, 7, 1), until=datetime.date(2000, 6, 1)) flist = [ dir + '/' + col + '/' + subcol + '_monthly_' + str(date.year) + str(date.month).zfill(2) + '.binary' for date in rr ] lon = sp.arange(0.5, 360.1, 1) lat = sp.arange(-89.5, 90, 1.) lev = sp.zeros(1) grid = dset.Grid(lon, lat, lev) vlist = [(subcol, '>f4', grid.dims)] time = sp.array(rr[:], dtype='|O') super(Ctl, self).__init__(flist, vlist, grid, time, undef, name)
def __init__(self): name = 'GPCP' undef = -9.99990000e+04 syr = 1982 eyr = 2006 yrs = [str(i) for i in range(syr, eyr)] mon = [str(i).zfill(2) for i in range(1, 13)] dates = [i + j for i in yrs for j in mon] dir = '/discover/nobackup/projects/gmao/share/dao_ops/verification/GPCP_Version_2' flist = sp.array([dir + '/gpcp_v2_psg.' + i + '.data' for i in dates]) lon = sp.arange(1.25, 361, 2.5) lat = sp.arange(-88.75, 91, 2.5) lev = sp.zeros(1) grid = dset.Grid(lon, lat, lev) vlist = [('precip', '>f4', grid.dims)] r = rrule.rrule(rrule.MONTHLY, dtstart=datetime.date(syr, 1, 1), count=len(flist)) time = sp.array(r[:], dtype='|O') super(Ctl, self).__init__(flist, vlist, grid, time, undef, name)
def __init__(self): name = 'QSCAT' undef = -9999. oceanval = os.environ.get( 'OCEANVAL', '/discover/nobackup/projects/gmao/oceanval/verification') flist = [oceanval + '/QSCAT/qscat_clim.dat'] lon = sp.arange(-179.5, 180, 1) lat = sp.arange(-89.5, 90, 1) lev = sp.zeros(1) grid = dset.Grid(lon, lat, lev) vlist = [] vlist.append(('taux', 'f4', grid.dims)) vlist.append(('tauy', 'f4', grid.dims)) t = rrule.rrule(rrule.MONTHLY, dtstart=datetime.date(0001, 1, 1), count=12) time = sp.array(t[:], dtype='|O') super(Ctl, self).__init__(flist, vlist, grid, time, undef, name)
def __init__(self): name='Reynolds MERRA' undef=-999.0 flist=['/discover/nobackup/projects/gmao/share/dao_ops/fvInput/g5gcm/bcs/realtime/SST/360x180/dataoceanfile_MERRA_sst_1971-current.360x180.LE'] lon=sp.arange(-179.5,180,1); lat=sp.arange(-89.5,90,1); lev=sp.zeros(1) grid=dset.Grid(lon,lat,lev) vlist=[('','i4',1), ('head','f4',14), ('','i4',1), ('','i4',1), ('sst','f4',grid.dims), ('','i4',1)] time=sp.repeat(datetime.date(1,1,1),1732) a=sp.memmap(flist[0],dtype=vlist,mode='r') for i,hh in enumerate(a['head']): time[i]=datetime.date(hh[0],hh[1],hh[2]) a.close() super(Ctl,self).__init__(flist,vlist,grid,time,undef,name)
def __init__(self): name = 'ECMWF_Interim' undef = -1e-15 r = rrule.rrule(rrule.YEARLY, dtstart=datetime.date(1980, 1, 1), until=datetime.date(2011, 1, 1)) dir = '/discover/nobackup/projects/gmao/share/dao_ops/verification/ECMWF_Interim_FullRes/surface' flist = sp.array([ dir + '/Interim_monthly_FLX_' + str(date.year) + '.bin' for date in r[:] ]) lon = sp.arange(0, 360.5, 0.75) lat = sp.arange(-90, 90.5, 0.75) lev = sp.zeros(1) grid = dset.Grid(lon, lat, lev) vlist = [ ('cp', 'f4', grid.dims), ('e', 'f4', grid.dims), ('slhf', 'f4', grid.dims), ('sshf', 'f4', grid.dims), ('ssr', 'f4', grid.dims), ('ssrc', 'f4', grid.dims), ('ssrd', 'f4', grid.dims), ('str', 'f4', grid.dims), ('strc', 'f4', grid.dims), ('strd', 'f4', grid.dims), ('tisr', 'f4', grid.dims), ('tp', 'f4', grid.dims), ('tsr', 'f4', grid.dims), ('tsrc', 'f4', grid.dims), ('ttr', 'f4', grid.dims), ('ttrc', 'f4', grid.dims), ('csf', 'f4', grid.dims), ('es', 'f4', grid.dims), ('par', 'f4', grid.dims), ('ro', 'f4', grid.dims), ('sf', 'f4', grid.dims), ('smlt', 'f4', grid.dims), ('lsf', 'f4', grid.dims), ('lsp', 'f4', grid.dims), ('parcs', 'f4', grid.dims), ('uvb', 'f4', grid.dims), ] r = rrule.rrule(rrule.MONTHLY, dtstart=datetime.date(1980, 1, 1), until=datetime.date(2011, 12, 1)) time = sp.array(r[:], dtype='|O') super(Ctl, self).__init__(flist, vlist, grid, time, undef, name)
def __init__(self): name = 'lwmask' undef = 1e15 dir = '/discover/nobackup/projects/gmao/share/dao_ops/verification/lwmask' flist = [dir + '/lwmask1440721.data'] lon = sp.arange(-180, 180, 0.25) lat = sp.arange(-90, 90.1, 0.25) lev = sp.zeros(1) grid = dset.Grid(lon, lat, lev) vlist = [('mask', '>f4', grid.dims)] time = sp.array([date(2000, 1, 1)], dtype='|O') super(Ctl, self).__init__(flist, vlist, grid, time, undef, name)
def __init__(self): name = 'COADS' undef = -1e10 flist = [ os.environ['NOBACKUP'] + '/verification/coads/coads2x25_climate.dat' ] lon = sp.arange(0.0, 360., 2.5) lat = sp.arange(-90., 91, 2.) lev = sp.zeros(1) grid = dset.Grid(lon, lat, lev) rr = rrule.rrule(rrule.MONTHLY, dtstart=datetime.date(1980, 1, 15), count=12) time = sp.array(rr[:], dtype='|O') super(Ctl, self).__init__(flist, vlist, grid, time, undef, name)
def __init__(self): name = 'QSCAT' undef = -9999. path = os.environ['HOME'] + '/verification/stress_mon_clim' flist = [path + '/qscat_clim.dat'] lon = sp.arange(-179.5, 180, 1) lat = sp.arange(-89.5, 90, 1) lev = sp.zeros(1) grid = dset.Grid(lon, lat, lev) vlist = [] vlist.append(('taux', 'f4', grid.dims)) vlist.append(('tauy', 'f4', grid.dims)) t = rrule.rrule(rrule.MONTHLY, dtstart=datetime.date(0001, 1, 1), count=12) time = sp.array(t[:], dtype='|O') super(Ctl, self).__init__(flist, vlist, grid, time, undef, name)
def __init__(self): name='ERBE' undef=1e15 dir='/discover/nobackup/projects/gmao/share/dao_ops/verification/Clouds_radiation/erbe' flist=sp.array([dir+'/erbe2x25.data']) lon=sp.arange(0.,360,2.5); lat=sp.arange(-90.,91,2.); lev=sp.zeros(1) grid=dset.Grid(lon,lat,lev) print grid.dims vlist=[('OLR','>f4',grid.dims)] vlist.append(('OLRCLR','>f4',grid.dims)) vlist.append(('NSR','>f4',grid.dims)) vlist.append(('NSRCLR','>f4',grid.dims)) vlist.append(('OSR','>f4',grid.dims)) vlist.append(('OSRCLR','>f4',grid.dims)) r=rrule.rrule(rrule.MONTHLY,dtstart=datetime.date(1985,1,1),count=60) time=sp.array(r[:],dtype='|O') super(Ctl,self).__init__(flist,vlist,grid,time,undef,name)
def __init__(self, collection='lh'): expid = 'OAFlux' undef = 10e15 expdir = oceanval + '/OAfluxes/turbulence_1958-2006/' + collection + '/grads' r = rrule.rrule(rrule.MONTHLY, dtstart=datetime.date(1958, 1, 1), count=588) flist = sp.array([ expdir + '/a.' + collection + '_oaflux.' + str(year) + '.grads' for year in xrange(1958, 2007) ]) lon = sp.arange(0.5, 360, 1) lat = sp.arange(-89.5, 90, 1) lev = sp.zeros(1) grid = dset.Grid(lon, lat, lev) vlist = [(collection, '<f4', grid.dims), ('err', '<f4', grid.dims)] time = sp.array(r[:], dtype='|O') super(CtlTurb, self).__init__(flist, vlist, grid, time, undef, name=expid)
def __init__(self): expid = 'OAFlux' undef = 10e15 expdir = oceanval + '/OAfluxes/netheat_1983-2004/grads' r = rrule.rrule(rrule.MONTHLY, dtstart=datetime.date(1983, 7, 1), count=258) flist = sp.array([ expdir + '/a.qnet.' + str(year) + '.grads' for year in xrange(1983, 2005) ]) lon = sp.arange(0.5, 360, 1) lat = sp.arange(-89.5, 90, 1) lev = sp.zeros(1) grid = dset.Grid(lon, lat, lev) vlist = [('qnet', '<f4', grid.dims)] time = sp.array(r[:], dtype='|O') super(CtlQnet, self).__init__(flist, vlist, grid, time, undef, name=expid)