for vtype in lvtype: #*** init **** lplainv = [] dfrontv = {} for ftype in lftype: dfrontv[ftype] = [] #************* for year in lyear: for mon in lmon: print year,mon #*************************************** # mask #------------------------ region = "ASAS" nx_fig,ny_fig = chart_para.ret_nxnyfig(region, year, mon) paradate = datetime.date(year,mon,1) xydatadir = "/media/disk2/out/chart/ASAS/const" if (paradate < datetime.date(2006,1,1)): name_x_corres = xydatadir + "/stereo.xfort.fig2saone.ASAS.2000.01.bn" name_y_corres = xydatadir + "/stereo.yfort.fig2saone.ASAS.2000.01.bn" name_domain_mask = xydatadir + "/domainmask_saone.%s.2000.01.bn"%(region) if ( datetime.date(2006,1,1)<=paradate<datetime.date(2006,3,1)): name_x_corres = xydatadir + "/stereo.xfort.fig2saone.ASAS.2006.01.bn" name_y_corres = xydatadir + "/stereo.yfort.fig2saone.ASAS.2006.01.bn" name_domain_mask = xydatadir + "/domainmask_saone.%s.2006.01.bn"%(region) if ( datetime.date(2006,3,1)<=paradate): name_x_corres = xydatadir + "/stereo.xfort.fig2saone.ASAS.2006.03.bn" name_y_corres = xydatadir + "/stereo.yfort.fig2saone.ASAS.2006.03.bn"
#ydom_last = 834 miss = -9999.0 #lperiod = [[2000,1],[2006,1],[2006,3]] lperiod = [[2000,1]] ny_dec, nx_dec = (700,1600) #------------------- #odir = "/home/utsumi/bin/dtanl/ctrack/chart" odir = "/media/disk2/out/chart/ASAS/const" for period in lperiod: #------------------ year, mon = period #------------------ nx,ny = chart_para.ret_nxnyfig(region, year, mon) lon0,lat0 = chart_para.ret_lonlat_center(region) latts = chart_para.ret_latts(region, year, mon) lllon, lllat, urlon, urlat \ = chart_para.ret_domain_corner(region, year, mon) xdom_first, xdom_last, ydom_first, ydom_last \ = chart_para.ret_xydom_first_last(region, year, mon) nxdom = xdom_last - xdom_first +1 nydom = ydom_last - ydom_first +1 #------------------ oname_lon = odir + "/stereo.lon.ASAS.%04d.%02d.bn"%(year,mon) oname_lat = odir + "/stereo.lat.ASAS.%04d.%02d.bn"%(year,mon) figname_lon = odir + "/stereo.lon.ASAS.%04d.%02d.png"%(year,mon) #**********************