Ejemplo n.º 1
0
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"
Ejemplo n.º 2
0
#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)
  #**********************