plotFields.py

Go to the documentation of this file.

#!/usr/local/bin/python
# (C) Copyright 2009-2016 ECMWF.
# 
# This software is licensed under the terms of the Apache Licence Version 2.0
# which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 
# In applying this licence, ECMWF does not waive the privileges and immunities 
# granted to it by virtue of its status as an intergovernmental organisation nor
# does it submit to any jurisdiction.


import QgPlot
import QgRead

#fdbdir="/home/rd/dai/Data/FDB"
fdbdir="/home/da/dasg/Data/FDB"
attributes = {}
attributes["expver"]        = "example"
attributes["validitydate"]  = "2010-01-01T00:00:00Z"

#truth
#attributes["referencedate"] = "2009-12-15T00:00:00Z"
#attributes["type"]          = "truth"

#bg
#attributes["type"]          = "fc"
#attributes["referencedate"] = "2009-12-31T00:00:00Z"

#an
attributes["type"]          = "an"
attributes["referencedate"] = "2010-01-02T00:00:00Z"

# Read the fields

qgReader = QgRead.QgRead()

pv, psi = qgReader.readFields(fdbdir,attributes)

# Plot the fields

qgplt = QgPlot.QgPlot()

qgplt.plot(pv,psi)