Changeset 699
- Timestamp:
- 04/11/16 18:36:51 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/file_Tbcorr_daily.py
r698 r699 6 6 import argparse 7 7 8 from pylab import * 8 import inspect 9 from datetime import datetime 10 import logging 11 12 #from pylab import * 9 13 import glob 10 from datetime import date14 #from datetime import date 11 15 from useCorrTb_daily import * 12 16 import itertools as it … … 28 32 parser = argparse.ArgumentParser(description=__doc__, 29 33 formatter_class=argparse.RawDescriptionHelpFormatter) 34 35 parser.add_argument("--verbose", 36 action='store_true', 37 default=False, 38 help=("run in verbose mode.")) 30 39 31 40 ### Satellite, Instrument, Produit Level, Channel … … 97 106 98 107 opts = vars(parser.parse_args()) 108 is_verbose = opts['verbose'] 99 109 lvl = opts['lvl'] 100 110 sat = opts['sat'] … … 106 116 latvec = arange(-34.5,34.6,resol) 107 117 lonvec = arange(-180.,179.76,resol) 108 118 119 # define log file 120 command = '{0}'.format(inspect.stack()[0][3]) 121 log_date = datetime.datetime.utcnow().strftime('%Y%m%dT%H%M%SZ') 122 log = '{0}/{1}.log.{2}'.format(os.environ['PROJECT_LOG'], command, log_date) 123 124 logger = logging.getLogger() 125 handler = logging.FileHandler(log) 126 logger.addHandler(handler) 127 if is_verbose: 128 logger.setLevel(logging.DEBUG) 129 else: 130 logger.setLevel(logging.INFO) 131 132 logger = logging.getLogger() 133 logger.info('[Context]') 134 logger.info('PID={0}'.format(os.getpid())) 135 logger.info('command={0}'.format(command)) 136 import socket 137 logger.info('hostname={0}'.format(socket.getfqdn())) 138 logger.info('runtime={0}'.format(log_date)) 139 logger.info('log={0}'.format(log)) 140 141 logger.info('[Parameters]') 142 logger.info('lvl={0}'.format(lvl)) 143 logger.info('sat={0}'.format(sat)) 144 logger.info('inst={0}'.format(inst)) 145 logger.info('ch={0}'.format(ch)) 146 logger.info('dmin={0}'.format(dmin)) 147 logger.info('dmax={0}'.format(dmax)) 148 109 149 # number of processes 110 150 npr = 10 111 151 112 152 for day in rrule(DAILY,dtstart=dmin,until=dmax): 113 print day153 logger.info('{0}'.format(day)) 114 154 ### File list 115 155 repday = [day.strftime('%Y_%m_%d')]#[ff.strftime('%Y_%m_%d') for ff in rrule(DAILY,dtstart=day+relativedelta(days=-1),until=day+relativedelta(days=+1))] … … 123 163 pool.join() 124 164 125 print 'Join' 165 logger.info('Join') 166 126 167 tbcorr, tobs, satvec = [dstack([proc[var] for proc in res if size(proc[var])]) for var in xrange(3)] 127 168 del res … … 141 182 142 183 ### Save in netcdf file 143 print 'Save in nc'184 logger.info('Save in nc') 144 185 nsat = satvec.size 145 186 nlon = lonvec.size … … 191 232 # Fermeture du fichier 192 233 ncfile.close() 193 print output+' created'234 logger.info('{0} created'.format(output)) 194 235 195 236
Note: See TracChangeset
for help on using the changeset viewer.