else: parser.error("Invalid time parameter") tstep = int(tstep) coords = options.coords.split(",") if len(coords) != 4: parser.error("Invalid coords parameter: %s", options.coords) xmin, ymin, xmax, ymax = map(float, coords) dtfmt = "%Y%m%dT%H%M%S" mars = MARS3D() if options.cfgfile: mars.load_default_config(options.cfgfile, nested=True) mars.load_config(options.cfgfile, nested=True) for varname in variables: log.notice("Variable %s", varname) for itv in Intervals((tmin, tmax, tbb), (tstep, ustep)): log.notice("Interval %s", itv) select = dict(time=itv[:2]) if options.level: select["level"] = slice(*map(int, options.level.split(","))) if options.bbox: select.update(dict(latitude=(latmin, latmax), longitude=(lonmin, lonmax))) pylab.figure() try: sec = mars.plot_section( varname=varname, select=select, xmin=xmin, ymin=ymin,
dtfmt = '%Y%m%dT%H%M%S' mars = MARS3D() profiles = ProfilesDataset() if options.cfgfile: mars.load_default_config(options.cfgfile, nested=True) profiles.load_default_config(options.cfgfile, nested=True) mars.load_config(options.cfgfile, nested=True) profiles.load_config(options.cfgfile, nested=True) if not options.bbox: #mars_grid = mars.get_grid() lat, lon = mars.get_latitude(), mars.get_longitude() lonmin, latmin, lonmax, latmax = min(lon), min(lat), max(lon), max(lat) for varname in variables: log.notice('Variable %s', varname) for depth in depths: for itv in Intervals((tmin, tmax, tbb),(tstep, ustep)): log.notice('Interval %s', itv) plotkw = dict(plot_show=False) select = dict(time=itv[:2], latitude=(latmin,latmax,'ccb'), longitude=(lonmin,lonmax,'ccb')) if not options.coloc: pylab.figure() try: mp = mars.plot_layer(varname=varname, depth=depth, select=select, **plotkw) plotkw['map'] = mp
tstep = int(tstep) coords = options.coords.split(',') if len(coords) != 3: parser.error('Invalid coords parameter: %s', options.coords) xymin, xymax, xy = map(float, coords) dtfmt = '%Y%m%dT%H%M%S' mars = MARS3D() sat = Satellite() if options.cfgfile: mars.load_default_config(options.cfgfile, nested=True) sat.load_default_config(options.cfgfile, nested=True) mars.load_config(options.cfgfile, nested=True) sat.load_config(options.cfgfile, nested=True) for varname in variables: log.notice('Variable %s', varname) for itv in Intervals((tmin, tmax, tbb),(tstep, ustep)): log.notice('Interval %s', itv) select = dict(time=itv[:2]) if options.bbox: select.update(dict(latitude=(latmin,latmax), longitude=(lonmin,lonmax))) datakw = dict( varname=('temperature', 'TEMP'), select=select, xorymin=xymin, xorymax=xymax, xory=xy, meridonal=options.meridional, cur_axes_rect=[0.1, 0.1, 0.6, 0.7], map_axes_rect=[0.75, 0.1, 0.2, 0.3], plot_show=False) if options.meridional: lat = ([xymin], [xymax])
dtfmt = '%Y%m%dT%H%M%S' mars = MARS3D() profiles = ProfilesDataset() if options.cfgfile: mars.load_default_config(options.cfgfile, nested=True) profiles.load_default_config(options.cfgfile, nested=True) mars.load_config(options.cfgfile, nested=True) profiles.load_config(options.cfgfile, nested=True) if not options.bbox: #mars_grid = mars.get_grid() lat, lon = mars.get_latitude(), mars.get_longitude() lonmin, latmin, lonmax, latmax = min(lon), min(lat), max(lon), max(lat) for varname in variables: log.notice('Variable %s', varname) for depth in depths: for itv in Intervals((tmin, tmax, tbb), (tstep, ustep)): log.notice('Interval %s', itv) plotkw = dict(plot_show=False) select = dict(time=itv[:2], latitude=(latmin, latmax, 'ccb'), longitude=(lonmin, lonmax, 'ccb')) if not options.coloc: pylab.figure() try:
select = dict() mars = MARS3D() profiles = ProfilesDataset() if options.cfgfile: mars.load_default_config(options.cfgfile, nested=True) profiles.load_default_config(options.cfgfile, nested=True) mars.load_config(options.cfgfile, nested=True) profiles.load_config(options.cfgfile, nested=True) if not options.bbox: #mars_grid = mars.get_grid() lat, lon = mars.get_latitude(), mars.get_longitude() lonmin, latmin, lonmax, latmax = min(lon), min(lat), max(lon), max(lat) for itv in Intervals((tmin, tmax, tbb), (tstep, ustep)): log.notice('Interval %s', itv) select = dict(time=itv, latitude=(latmin, latmax, 'ccb'), longitude=(lonmin, lonmax, 'ccb')) plotkw = dict(plot_show=False) # ========== # Methode 1. Tracé des observations (points) sur fond 2D grillés de champs de modèles moyennés sur la période choisie. if not options.coloc: if 'mld' in options.plots: pylab.figure() kw = plotkw.copy() m = c = None try:
select=dict() mars = MARS3D() profiles = ProfilesDataset() if options.cfgfile: mars.load_default_config(options.cfgfile, nested=True) profiles.load_default_config(options.cfgfile, nested=True) mars.load_config(options.cfgfile, nested=True) profiles.load_config(options.cfgfile, nested=True) if not options.bbox: #mars_grid = mars.get_grid() lat, lon = mars.get_latitude(), mars.get_longitude() lonmin, latmin, lonmax, latmax = min(lon), min(lat), max(lon), max(lat) for itv in Intervals((tmin, tmax, tbb),(tstep, ustep)): log.notice('Interval %s', itv) select = dict(time=itv, latitude=(latmin,latmax,'ccb'), longitude=(lonmin,lonmax,'ccb')) plotkw = dict(plot_show=False) # ========== # Methode 1. Tracé des observations (points) sur fond 2D grillés de champs de modèles moyennés sur la période choisie. if not options.coloc: if 'mld' in options.plots: pylab.figure() kw = plotkw.copy() m = c = None try: m,c = mars.plot_mld(select, **kw) except: log.exception('Cannot plot mars mld') kw.update(map=m, cmap=c)