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)
#**********************