#!/usr/bin/env python

##########################################################################
#var = ["151","146","167"]
#var = ["167"]
##var = ["78"]
##var = ["137"]
##var = ["174"] albedo marche pas.
#var = ["134","151"] 
#fieldtype = "2d"
##var = ["130"]
#lev = ["700.","850."]
#tim = ["00"]
#date = ['10','08','2010','10','08','2010']
#date = ['01','09','2009','01','09','2009']
##date = ['01','09','2010','01','09','2010']
#area = "Whole"
###########################################################################


##########################################################################
######### FIGURES DE LA PRESSION
var = ["134","151"]
fieldtype = "2d"
lev = ["9999."]
tim = ["00"]
date = ['01','09','2009','01','09','2009']
area = "Whole"
##########################################################################
var = ["130"]
fieldtype = "3d"
lev = ["1000.","900.","800.","700.","600."]
tim = ["00"]
date = ['01','09','2009','01','09','2009']
area = "Whole"

var=["189","228"]
fieldtype = "2d"
lev = ["9999."]
tim = ["00","12"]
date = ['01','09','2009','10','09','2009']
area = "Whole"



##########################################################################
from myplot import latinterv
from myecmwf import get_ecmwf
from os import system
##########################################################################
[wlon,wlat] = latinterv(area)
nc = get_ecmwf (var, fieldtype, wlat, wlon, lev, date, tim) 
system("cp -f create_ecmwf.py output.py")