starttime='20160301' end__time='20160308' INPUTDIR='/gpfs/work/IscrC_MYMEDBIO/COPERNICUS/online_validation_data/TMP/' #args.inputdir OUTDIR ='/gpfs/work/IscrC_MYMEDBIO/COPERNICUS/online_validation_data/TMP/' #args.outdir maskfile='/pico/home/usera07ogs/a07ogs00/OPA/V2C/etc/static-data/MED1672_cut/MASK/meshmask.nc' # args.maskfile TI=TimeInterval(starttime,end__time,'%Y%m%d') archive_dir='/pico/home/usera07ogs/a07ogs00/OPA/V4/archive' TheMask=Mask(maskfile) TS = TimeSeries(TI, archive_dir,postfix_dir='POSTPROC/AVE_FREQ_1/',glob_pattern="ave*gz") forecasts =TS.get_forecast_days(rundays=[2]) forecasts_sublist=TS.get_sublist(forecasts,[2,3,4]) #forecast tuesday and wed,thu sat_archive="/gss/gss_work/DRES_OGS_BiGe/Observations/TIME_RAW_DATA/ONLINE/SAT/MODIS/DAILY/CHECKED/" DAILY_SAT_LIST=TS.get_daily_sat(forecasts_sublist,sat_archive) # float aggregator already done by others day=0 surf_layer=Layer(0,10) for time,archived_file,satfile in DAILY_SAT_LIST: avefile=INPUTDIR + os.path.basename(archived_file)[:-3] day=day+1 outfile=OUTDIR + "misfit+%dh.nc" % (day*24) print avefile continue Sat16 = Sat.convertinV4format( Sat.readfromfile(satfile) ) De = DataExtractor(TheMask,filename=avefile, varname='P_i') Model = MapBuilder.get_layer_average(De, surf_layer) Misfit = Sat16-Model