Esempio n. 1
0
    datdir     = odir
    datname    = oname

    if region == "GLOB":
      bnd        = [-40,-20,-10,-5,-2,2,5,10,20,40]
      figdir     = datdir
      dotstep    = 5
    if region == "JPN":
      bnd        = [-40,-20,-10,-5,-2,2,5,10,20,40]
      figdir     = datdir + ".JPN"
      dotstep    = 1
      ctrack_func.mk_dir(figdir)
    #----

    figname    = figdir + "/" + datname.split("/")[-1][:-7] + ".png"
    cbarname   = figdir + "/dpr.dec.cbar.png"

    mons     = len(ctrack_para.ret_lmon(season))
    a2figdat = fromfile(datname,float32).reshape(ny,nx) *60*60*24.*totaltimes /mons  # (mm/month)

    ##-------------------------------
    ##--- test ---
    #a2shade = ones([ny,nx],float32)
    ##------------
    ctrack_fig.mk_pict_saone_reg_symm_dotshade(a2figdat, lllat=lllat, lllon=lllon, urlat=urlat, urlon=urlon, bnd=bnd, mycm=mycm, soname=figname, stitle=stitle, miss=miss, a2shade=a2shade, cbarname=cbarname, dotcolor="k",dotstep=dotstep )
    #ctrack_fig.mk_pict_saone_reg_symm(a2figdat, lllat=lllat, lllon=lllon, urlat=urlat, urlon=urlon, bnd=bnd, mycm=mycm, soname=figname, stitle=stitle, miss=miss, a2shade=a2shade, cbarname=cbarname)

    print figname


Esempio n. 2
0
# at the center of original grid box

lllat   = 25.
urlat   = 60.
lllon   = 110.
urlon   = 180.

#lllat   = -89.5
#urlat   = 89.5
#lllon   = 0.5
#urlon   = 359.5

thdura  = 72

#----------------------------
lmon = ctrack_para.ret_lmon(season)
#lmon = [1]
#----------------------------
dlat    = 1.0
dlon    = 1.0
a1lat   = arange(-89.5, 89.5 +dlat*0.5,  dlat)
a1lon   = arange(0.5,   359.5 +dlon*0.5, dlon)
#----------------------------
#----------------------------

psldir_root     = "/media/disk2/out/chart/ASAS/exc"
lifedir_root    = "/media/disk2/out/chart/ASAS/exc/life"
nextposdir_root = "/media/disk2/out/chart/ASAS/exc/nextpos"

#************************************
ltrack     = []
Esempio n. 3
0
  if region == "GLOB":
    bnd        = [-40,-20,-10,-5,-2,2,5,10,20,40]
    bnd        = [-40,-20,-10,-5,-2,2,5,10,20,40]
    figdir     = datdir
    dotstep    = 5
  if region == "JPN":
    bnd        = [-40,-20,-10,-5,-2,2,5,10,20,40]
    figdir     = datdir + ".JPN"
    dotstep    = 1
    ctrack_func.mk_dir(figdir)
  #----

  figname    = figdir + "/" + datname.split("/")[-1][:-7] + ".png"
  cbarname   = figdir + "/dpr.dec.cbar.png"

  mons     = len(lyear) * len(ctrack_para.ret_lmon(season))

  a2figdat = fromfile(datname,float32).reshape(ny,nx) *60*60*24.*totaldays /mons  # (mm/month)

  ##-------------------------------
  ##--- test ---
  #a2shade = ones([ny,nx],float32)
  ##------------
  ctrack_fig.mk_pict_saone_reg_symm_dotshade(a2figdat, lllat=lllat, lllon=lllon, urlat=urlat, urlon=urlon, bnd=bnd, mycm=mycm, soname=figname, stitle=stitle, miss=miss, a2shade=a2shade, cbarname=cbarname, dotcolor="k",dotstep=dotstep )
  #ctrack_fig.mk_pict_saone_reg_symm(a2figdat, lllat=lllat, lllon=lllon, urlat=urlat, urlon=urlon, bnd=bnd, mycm=mycm, soname=figname, stitle=stitle, miss=miss, a2shade=a2shade, cbarname=cbarname)

  print figname