Python ret_thrvort Examples

Example #1

File: find.tc.py Project: nbykutsumi/dtanl

lmon   = range(1,12+1)
iday   = 1
lhour  = [0,6,12,18]
miss   = -9999.0

ny     = 180
nx     = 360
model  = "org"
dpgradrange   = ctrack_para.ret_dpgradrange()
thpgrad        = dpgradrange[1][0] 
thdura   = 36
#thwind   = 15 #m/s 
#thrvort  = 3.5e-5
#thwcore  = 1.5  # (K)
thwind   = tc_para.ret_thwind()
thrvort  = tc_para.ret_thrvort()
thwcore  = tc_para.ret_thwcore()
thsst    = tc_para.ret_thsst()
#
#thdura   = -9999.0
#thwind   = -9999.0 #m/s 
#thrvort  = -9999.0
#thwcore  = -9999.0  # (K)


plev_low = 850*100.0 # (Pa)
plev_mid = 500*100.0 # (Pa)
plev_up  = 250*100.0 # (Pa)

psldir_root     = "/media/disk2/data/JRA25/sa.one.%s/6hr/PRMSL"%(model)
pgraddir_root   = "/media/disk2/out/JRA25/sa.one.%s/6hr/pgrad"%(model)

Example #2

File: mk.tclines.cmip.py Project: nbykutsumi/dtanl

urlat = 80.0
lllon = 0.0
urlon = 360.0
# ----------------------------
for expr, model in [[expr, model] for expr in lexpr for model in lmodel]:
    # ----
    ens = cmip_para.ret_ens(model, expr, "psl")
    sunit, scalendar = cmip_para.ret_unit_calendar(model, expr)
    # ----
    # thsst    = 273.15 + 25.0
    # thwind   = 0.0 #m/s
    # thrvort  = 7.0e-5
    # thwcore  = 0.5  # (K)
    thsst = tc_para.ret_thsst()
    thwind = tc_para.ret_thwind()
    thrvort = tc_para.ret_thrvort(model)
    thwcore = tc_para.ret_thwcore(model)

    plev_low = 850 * 100.0  # (Pa)
    plev_mid = 500 * 100.0  # (Pa)
    plev_up = 250 * 100.0  # (Pa)

    tplev_low = 850 * 100.0

    psldir_root = "/media/disk2/data/CMIP5/sa.one.%s.%s/psl" % (model, expr)
    pgraddir_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/pgrad" % (model, expr)
    lifedir_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/life" % (model, expr)
    lastposdir_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/lastpos" % (model, expr)
    iposdir_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/ipos" % (model, expr)
    nextposdir_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/nextpos" % (model, expr)
    #

Example #3

File: freq.exc.jra.rad.py Project: nbykutsumi/dtanl

  pgradmin = ctrack_para.ret_dpgradrange()[2][0]  # Pa/1000km
  #-- orog ------------------------
  orogname = "/media/disk2/data/JRA25/sa.one.125/const/topo/topo.sa.one"
  a2orog   = fromfile(orogname, float32).reshape(ny,nx)

  a2shade  = ma.masked_where(a2orog >thorog, a2orog).filled(miss)


  #----------------------------
  for sresol in lsresol:
    #*******************
    # tc params
    #-------------------
    thsst    = tc_para.ret_thsst()
    thwind   = tc_para.ret_thwind()
    thrvort  = tc_para.ret_thrvort(sresol)
    thwcore  = tc_para.ret_thwcore(sresol)
    #-------------------
    psldir_root     = "/media/disk2/data/JRA25/sa.one.%s/6hr/PRMSL"%(sresol)
    pgraddir_root   = "/media/disk2/out/JRA25/sa.one.%s/6hr/pgrad"%(sresol)
    lifedir_root    = "/media/disk2/out/JRA25/sa.one.%s/6hr/life"%(sresol)
    #************************************
    # init
    #--------------------
    a2count = zeros([ny,nx],float32)

    #------------------------------------
    if calcflag ==True: 
      for year in range(iyear, eyear+1):
        # TC ----
        dbstxy   = tc_func.ret_ibtracs_dpyxy_saone(year)