Exemplo n.º 1
0
  def mkMaskRegrid(self, ltag, DTime, miss=0.0):
    """
    ltag = ["tag1", "tag2", ...], without "ot"
    """
    dictMask = {}
    for tag in ltag:
      #aMaskFrac     = biIntp(\
      aMaskFrac     = Regrid.biIntp(\
                      self.LatIn, self.LonIn   \
                    , self.mkMask(tag, DTime, miss=0.0) \
                    , self.LatOut, self.LonOut \
                    )[0]

      dictMask[tag] = ma.masked_greater_equal(\
                        ma.masked_less(aMaskFrac, 0.5).filled(0.0) \
                      , 0.5 \
                     ).filled(1.0)

    return dictMask
Exemplo n.º 2
0
Year = 2010
Mon = 3
Day = 1
Hour = 0
DTime = datetime(Year, Mon, Day, Hour)
model = "JRA55"
res = "bn"

T = Tag.Tag(model, res)
instMon = T.init_cyclone(Year, Mon)

a2t = T.mkMask("front.t", DTime, miss=0.0)
a2q = T.mkMask("front.q", DTime, miss=0.0)
a2c = T.mkMask("c", DTime)
a2cf = T.mkMask("cf", DTime)

LatIn = T.Front.Lat
LonIn = T.Front.Lon
LatOut = arange(-89.75, 89.75 + 0.01, 0.5)
LonOut = arange(0.25, 359.75 + 0.01, 0.5)

a2T = Regrid.biIntp(LatIn, LonIn, a2t, LatOut, LonOut)
# a2T = cf.regrids(LatIn, LonIn, a2t, LatOut, LonOut)
print shape(a2t)
# print shape(LatIn)
# print shape(LonIn)
print shape(a2T)
# print a2t
# print a2T