p2api 1.0.6__py2.py3-none-any.whl

p2api/utils/__init__.py

@@ -0,0 +1 @@
+# $Id: __init__.py 244843 2019-09-12 07:31:09Z jpritcha $

p2api/utils/add_absolute_time_intervals.py

@@ -0,0 +1,399 @@
+#! /usr/bin/env python
+
+# $Id: generate_finding_charts.py 267033 2021-09-16 11:08:57Z jpritcha $
+from __future__ import print_function
+
+_VERSION_ = "$Revision: 267033 $"
+
+import sys                            # interaction with Python interpreter
+from optparse import OptionParser     # for parsing command line arguments
+from optparse import OptionGroup
+
+import logging
+
+import keyring
+import datetime
+import re
+
+import p2api
+import p2api.utils.p2api_utils as p2api_utils
+
+from astropy.time import Time
+import numpy as np
+from astroplan import Observer
+import astropy.units as u
+
+################################################################################
+default_eso_period=p2api_utils.current_ESO_period()+1
+min_ATI_length=1.*u.hr
+################################################################################
+class get_options:
+  "get command line options"
+
+  def __init__(self, vers="0.0"):
+    # define options and arguments
+    parser = OptionParser(version="%prog "+vers, formatter=p2api_utils.IndentedHelpFormatterWithNL())
+
+    group = OptionGroup(parser, "Input OB(s)")
+    group.add_option(
+        "-i", "--ids",
+        dest="ids",
+        help="""\
+List of run, container, folder, and/or OB IDs. \
+Note that the ATCs will be applied to all OBs recursively when a \
+run, container or folder is specified.
+[required]\
+""",
+        default=None
+    )
+    parser.add_option_group(group)
+
+    group = OptionGroup(parser, "Specify the Absolute Time Intervals explicity.")
+    group.add_option(
+        "--ATCs", dest="ATCs",
+        help="""\
+Underscore separated list of comma separated Absolute Time Constraints \
+e.g. yyyy-mm-ddThh:mm,yyyy-mm-ddThh:mm_yyyy-mm-ddThh:mm,yyyy-mm-ddThh:mm{_...}
+[optional]\
+""",
+        default=None
+    )
+    group.add_option("-f", "--file", dest="file", help="A file with one ATC per line [optional]", default=None)
+    parser.add_option_group(group)
+
+    group = OptionGroup(parser, "Specify the Absolute Time Intervals using an orbital ephemeris. This functionality is EXPERIMENTAL, please check carefully!!!")
+    group.add_option("-E", "--ephem", dest="ephem", metavar="T0,PERIOD", help="Comma separated list of orbital T0[JD] and period[days] [optional]", default=None)
+    group.add_option("--phase_intervals", dest="phase_intervals",  help="Comma separated list of phase intervals e.g. 0.1-0.2,0.305-0.375 [optional]", default=None)
+    group.add_option("--phase_target", dest="phase_target",  help="Comma separated list of target phase, start_before [mins], min_start_before [mins] e.g. 0.01,120,30 [optional]", default=None)
+    group.add_option("--periods", dest="periods", help=f"Comma separated list of periods, default={default_eso_period} [optional]", default=str(default_eso_period))
+    parser.add_option_group(group)
+
+    group = OptionGroup(parser, "p2 server selection, options and password management")
+    group.add_option("--env", dest="env", help="Use the %s environment [optional], default = production" %("|".join(p2api.API_URL.keys())), default='production')
+    group.add_option("-u", "--user", dest="user", metavar="USER", help="the User Portal account of the user to use, e.g. ASMITHSONIAN [required for production (not required for demo)]", default=None)
+    group.add_option("-S", "--PWstore", dest="pwstore", metavar="PWSTORE", help="Store the password in the keyring, default=False", action="store_true", default=False)
+    group.add_option("-R", "--PWreenter", dest="pwreenter", metavar="PWREENTER", help="Reenter the password, thus overriding any PW in the keyring, default=False", action="store_true", default=False)
+    group.add_option("-Z", "--PWremove", dest="pwremove", metavar="PWREMOVE", help="Remove any password from the keyring, default=False", action="store_true", default=False)
+    parser.add_option_group(group)
+
+    group = OptionGroup(parser, "Verbosity...",
+                    "All optional")
+    group.add_option("-n", "--dryRun", dest="dryrun", help="Do everything but make any changes to tickets [optional]", action="store_true", default=False)
+    group.add_option("-D", "--debug", dest="debug", help="Write out debugging info [optional]", action="store_true", default=False)
+    #group.add_option("-v", "--verbose", dest="vbose", help="Write out verbose info [optional]", action="store_true", default=False)
+    group.add_option("-q", "--quiet", dest="quiet", help="Write out minimal info [optional]", action="store_true", default=False)
+    parser.add_option_group(group)
+
+    # parse arguments
+    (options, args) = parser.parse_args()
+
+    if options.ids is None or (
+        (
+            options.ATCs is None
+            and
+            options.file is None
+            and
+            (
+                options.ephem is None
+                or
+                options.periods is None
+                or
+                options.phase_intervals is None
+            )
+            and
+            (
+                options.ephem is None
+                or
+                options.periods is None
+                or
+                options.phase_target is None
+            )
+        )
+    ) :
+        parser.print_help()
+        sys.exit(2)
+
+    self.ids = options.ids
+    self.periods = options.periods
+    self.ATCs = options.ATCs
+    self.file = options.file
+    self.ephem = options.ephem
+    self.phase_intervals = options.phase_intervals
+    self.phase_target = options.phase_target
+
+    self.user = options.user
+    self.env = options.env
+    self.pwstore = options.pwstore
+    self.pwreenter = options.pwreenter
+    self.pwremove = options.pwremove
+
+    self.debug = options.debug
+    self.dryrun = options.dryrun
+    self.quiet = options.quiet
+
+    logFMT = "%(asctime)s %(module)15s[%(process)5d] [%(levelname)s] %(message)s"
+    logLVL=logging.INFO
+    if options.debug:
+      logLVL=logging.DEBUG
+    if options.quiet:
+      logLVL=logging.WARN
+    logging.basicConfig(level=logLVL, format=logFMT)
+
+options = get_options(_VERSION_)
+
+# ------------------------------------------------------------------------------
+def get_ob_ATCs(telescope=None, ob_info=None, run_info=None) :
+    ob_ATCs=[]
+    ATCs=None
+    check_ATCs=False
+    if options.ATCs is not None :
+        ATCs=[
+            [p2api_utils.decode_datetime(x.split(',')[0]),p2api_utils.decode_datetime(x.split(',')[1])]
+            for x in options.ATCs.splt('_')
+        ]
+
+    elif options.file is not None :
+        with open(options.file) as file:
+            lines = file.readlines()
+            lines = [line.rstrip().lstrip().split(',') for line in lines]
+        ATCs=[
+            [
+                p2api_utils.decode_datetime(x[0])+datetime.timedelta(seconds=30), # rounding to 1min
+                p2api_utils.decode_datetime(x[1])+datetime.timedelta(seconds=30), # rounding to 1min
+            ] for x in lines
+        ]
+
+    elif options.ephem is not None :
+        ob_exec_time=ob_info['executionTime']
+        T0,P=(float(x) for x in options.ephem.split(','))
+        periods=[
+            item for sublist in 
+            [[y for y in range(int(x.split('-')[0]),int(x.split('-')[-1])+1)] for x in options.periods.split(',')]
+            for item in sublist
+        ]
+        periods.sort()
+        period_range=[periods[0],]
+        period_ranges=[]
+        _pp=periods[0]
+        sf_jds=[]
+        for _p in periods[1:] :
+            if _p > _pp+1 :
+                period_range+=[_pp,]
+                sf_jd=Time([p2api_utils.ESO_period_start_end_datetimes(ESO_P) for ESO_P in period_range]).jd
+                sf_jds+=[[sf_jd[0][0],sf_jd[1][1]],]
+                period_ranges+=[period_range,]
+                period_range=[_p,]
+            _pp=_p
+        period_range+=[_pp,]
+        sf_jd=Time([p2api_utils.ESO_period_start_end_datetimes(ESO_P) for ESO_P in period_range]).jd
+        sf_jds+=[[sf_jd[0][0],sf_jd[1][1]],]
+        period_ranges+=[period_range,]
+
+        post_ob_exec_time=0.
+        if options.phase_intervals is not None :
+            phase_intervals=[[float(y) for y in x.split('-')] for x in options.phase_intervals.split(',')]
+            min_start_before=0.
+        elif options.phase_target is not None :
+            phase_target,start_before,min_start_before=(float(x) for x in options.phase_target.split(','))
+            phase_intervals=[[phase_target-start_before/P/24./60.,phase_target],]
+            post_ob_exec_time=ob_exec_time
+            ob_exec_time=0.
+        else :
+            logging.error('No means to calculate phase intervals specified.')
+
+        checked_phase_intervals=[]
+        for p_int in phase_intervals :
+            if p_int[1] < p_int[0] :
+                p_int[0]=p_int[0]-1.
+            p_int_in_sec=(p_int[1]-p_int[0])*P*24.*60.*60.
+            if (
+                ob_exec_time < p_int_in_sec
+                and
+                p_int_in_sec > min_ATI_length.to_value(u.second)
+            ) :
+                checked_phase_intervals+=[p_int,]
+            else :
+                if ob_exec_time < p_int_in_sec :
+                    logging.warning(f'Phase interval {p_int[0]}--{p_int[1]} == {p_int_in_sec:3.1f}sec is too short to start and finish this OB with execution time of {ob_exec_time}secs.')
+                if p_int_in_sec > min_ATI_length.to_value(u.second) :
+                    logging.warning(f'Phase interval {p_int[0]}--{p_int[1]} == {p_int_in_sec:3.1f}sec is shorter than the minimum ATI length time of {min_ATI_length.to_value(u.second)}secs.')
+
+        sf_ints=[]
+        for sf_jd in sf_jds :
+            orbits=(np.array(sf_jd)-T0)//P
+            for o in range(int(orbits[0])-1,int(orbits[1])+1) :
+                for p_int in checked_phase_intervals :
+                    sf_int=[
+                        T0+(o+p_int[0])*P,
+                        T0+(o+p_int[1])*P-ob_exec_time/24./60./60.,
+                    ]
+                    if (
+                        ( sf_int[0] > sf_jd[0] and sf_int[0] < sf_jd[1] )
+                        or
+                        ( sf_int[1] > sf_jd[0] and sf_int[1] < sf_jd[1] )
+                    ) :
+                        # Interval is within the time range of the specified period(s)...
+                        original_sf_int=sf_int
+                        int_spans_nighttime=True
+                        if sf_int[1] - sf_int[0] < 1. :
+                            sssr=[]
+                            # Check that some part of the interval occurs during nighttime...
+                            int_spans_nighttime=False
+                            for which in ['previous','next',] :
+                                sun_set_time=telescope.sun_set_time(Time(sf_int[0], format='jd'), which=which)
+                                sun_rise_time=telescope.sun_rise_time(sun_set_time, which='next')
+                                sssr+=[[sun_set_time,sun_rise_time],]
+                                if (
+                                    (sf_int[0] > sun_set_time.jd and sf_int[0] < sun_rise_time.jd-ob_exec_time/24./60./60. )
+                                    or
+                                    (sf_int[1] > sun_set_time.jd and sf_int[1] < sun_rise_time.jd-ob_exec_time/24./60./60. )
+                                    or
+                                    (sf_int[0] < sun_set_time.jd and sf_int[1] > sun_rise_time.jd-ob_exec_time/24./60./60. )
+                                ) :
+                                    sf_int=[
+                                        np.max([sf_int[0],sun_set_time.jd]),
+                                        np.min([sf_int[1],sun_rise_time.jd-ob_exec_time/24./60./60.,sun_rise_time.jd-post_ob_exec_time/24./60./60.]),
+                                    ]
+                                    if sf_int[1] > sf_int[0]+min_start_before/24./60. :
+                                        int_spans_nighttime=True                                    
+                                        break
+                        if (
+                            int_spans_nighttime
+                            and
+                            ((sf_int[1]-sf_int[0])*24*60.*60. > ob_exec_time)
+                            and
+                            (sf_int[1]-sf_int[0])*24*60.*60. > min_ATI_length.to_value(u.second)
+                        ) :
+                            sf_ints+=[sf_int,]
+                        else :
+                            if not int_spans_nighttime :
+                                logging.debug(f'Interval does not span any part of the night, original interval = {original_sf_int}, SSSR of previous and next nights {sssr[0][0].jd}--{sssr[0][1].jd} and {sssr[1][0].jd}--{sssr[1][1].jd}')
+                            if ((sf_int[1]-sf_int[0])*24*60.*60. < ob_exec_time) :
+                                logging.debug(f'Checked Interval shorter [{(sf_int[1]-sf_int[0])*24*60.*60.:3.1f}sec] than OB execution time [{ob_exec_time:3.1f}sec], original interval = {original_sf_int}, checked interval = {sf_int}')
+                            if ((sf_int[1]-sf_int[0])*24*60.*60. < min_ATI_length.to_value(u.second)) :
+                                logging.debug(f'Checked Interval shorter [{(sf_int[1]-sf_int[0])*24*60.*60.:3.1f}sec] than minimum ATI length time [{min_ATI_length.to_value(u.second):3.1f}sec], original interval = {original_sf_int}, checked interval = {sf_int}')
+                pass
+        ATCs=[Time(sf_int, format='jd') for sf_int in sf_ints]
+        check_ATCs=True
+    else :
+        raise "Don't know how to compute Absolute Time Constraints."
+ 
+    ob_ATCs=[]
+    for i,ATC in enumerate(ATCs) :
+        t=Time(ATC)
+        t.format = 'isot'
+        t.out_subfmt = 'date_hm'
+        ob_ATCs+=[{
+            'from': t[0].isot,
+            'to':   t[1].isot,
+        },]
+        if check_ATCs :
+            phase_range=''
+            if options.ephem is not None :
+                phase_range=f"  [Phases: {(((t[0]-Time(T0,format='jd'))/P).to_value(u.day)%1):5.3f} -- {(((t[1]-Time(T0,format='jd'))/P).to_value(u.day)%1):5.3f}]"
+            print(f"{i+1:3d} {ob_ATCs[-1]['from']} -- {ob_ATCs[-1]['to']}  [duration={((t[1]-t[0]).to_datetime())}]{phase_range}")
+    if check_ATCs :
+        '''
+        logging.warning("****************************************************************************************")
+        logging.warning("****************************************************************************************")
+        logging.warning("****************************************************************************************")
+        logging.warning("Intervals via Orbital Ephemeris functionality is EXPERIMENTAL")
+        logging.warning("Please check the above intervals carefully and only continue if you are sure they are what YOU require.")
+        logging.warning("YOU, and YOU ALONE are responsible for the accuracy of these inetervals.")
+        '''
+        print("****************************************************************************************")
+        print("****************************************************************************************")
+        print("****************************************************************************************")
+        print("Intervals via Orbital Ephemeris functionality is EXPERIMENTAL")
+        print("Please check the above intervals carefully and only continue if you are sure they are what YOU require.")
+        print("YOU, and YOU ALONE are responsible for the accuracy of these intervals.")
+        yorn = input("Continue [y/N] ")
+        if yorn.lower() not in ['y',] :
+            print("Please contact ESO User Support via a ticket at https://support.eso.org/tickets")
+            print("to report/discuss the errors/problems in the above calculations.")
+            print("Please include in the description the EXACT command line used to run the script.")
+            sys.exit(0)
+            #raise Exception(f"Aborted.")
+    return ob_ATCs
+# ------------------------------------------------------------------------------
+def oneOB( api, obId ) :
+    try :
+        OB, _  = api.getOB(obId)
+    except p2api.P2Error as e :
+        ## Looks like this ID is not an OB...
+        logging.error('OB %d could not be found...' %(obId))
+        return
+    except:
+        raise
+    ob_info = api.executionTime(obId)
+    run_info,_=api.getRun(OB['runId'])
+
+    telescope = Observer(location=p2api_utils.telescope_locations[run_info['telescope']])
+
+    if OB['itemType'] not in ['OB',] :
+        logging.warning('OB %d is not a science OB' %(obId))
+        return
+
+    if OB['obStatus'] not in ['P','-'] :
+        logging.warning('OB status is %s, OB not modifiable' %(OB['obstatus']))
+        return
+
+    runId, _ = api.getRun(OB['runId'])
+    ob_url=api.apiUrl+'/obsBlocks/'+str(obId)
+    oneOB_sdt=datetime.datetime.now()
+    ## For the moment only setting is supported,
+    ## ToDo: implement support for adding to...
+    logging.info("Setting Absolute Time Intervals for OBid = %d, OBname = '%s' @ %s..." % (OB['obId'],OB['name'],OB['instrument']))
+    absTCs, atcVersion = api.getAbsoluteTimeConstraints(obId)
+
+    ob_ATCs=get_ob_ATCs(
+        telescope=telescope,
+        ob_info=ob_info,
+        run_info=run_info,
+    )
+
+    if not options.dryrun :
+        api.saveAbsoluteTimeConstraints(
+            obId,
+            ob_ATCs,
+            atcVersion,
+        )
+    oneOB_edt=datetime.datetime.now()
+    logging.info(
+        "Added %d Absolute Time Intervals for OBid = %d, OBname = '%s' @ %s done in %s."
+        %(
+            len(ob_ATCs), obId, OB['name'], OB['instrument'], str(oneOB_edt-oneOB_sdt)
+        )
+    )
+# ------------------------------------------------------------------------------
+
+# ------------------------------------------------------------------------------
+# ------------------------------------------------------------------------------
+# ------------------------------------------------------------------------------
+
+def main():
+
+    ################################################################################
+    # remove password from keyring...
+    if options.pwremove and options.user is not None :
+        keyring.delete_password('usd_p2api_interactive:www.eso.org', options.user)
+        sys.exit(0)
+
+    api=p2api_utils.get_api(options)
+    obIds=p2api_utils.r_get_all_OBIds( api, options.ids.split(',') )
+
+    if len(obIds) > 0 :
+        if len(obIds) > 1 :
+            sdt=datetime.datetime.now()
+            logging.warning('Adding Absolute Time Constraints to %d OBs...' %( len(obIds) ))
+        for obId in obIds :
+            oneOB(api, obId, )
+        if len(obIds) > 1 :
+            edt=datetime.datetime.now()
+            logging.info('Total elapsed time to process %d finding charts %s...' %(len(obIds),str(edt-sdt)))
+    else :
+        logging.error('No OBs found.')
+
+
+if __name__ == '__main__' :
+    main()

p2api/utils/generate_finding_charts.py

@@ -0,0 +1,260 @@
+#! /usr/bin/env python
+
+# $Id: generate_finding_charts.py 267033 2021-09-16 11:08:57Z jpritcha $
+from __future__ import print_function
+
+_VERSION_ = "$Revision: 267033 $"
+
+import sys                            # interaction with Python interpreter
+from optparse import OptionParser     # for parsing command line arguments
+from optparse import OptionGroup
+
+import logging
+
+import keyring
+import datetime
+import json
+import requests
+
+import p2api
+import p2api.utils.p2api_utils as p2api_utils
+################################################################################
+
+################################################################################
+class get_options:
+  "get command line options"
+
+  def __init__(self, vers="0.0"):
+    # define options and arguments
+    parser = OptionParser(version="%prog "+vers, formatter=p2api_utils.IndentedHelpFormatterWithNL())
+
+    group = OptionGroup(parser, "Input OB(s)")
+    group.add_option("-i", "--ids", dest="ids", metavar="Ids", help="List of OB IDs [required]", default=None)
+    group.add_option("--obsdate", dest="obsdate", metavar="OBS_DATE", help="Set the obsdate [optional]", default=None)
+    parser.add_option_group(group)
+
+    group = OptionGroup(parser, "Image display parameters",
+"""\
+The following parameters affect the way the 'default' images are displayed in the FCs.
+They correspond directly to the corresponding parameters of the APLPY FITSFigure.show_grayscale method, see:
+  https://aplpy.readthedocs.io/en/stable/api/aplpy.FITSFigure.html#aplpy.FITSFigure.show_grayscale
+for details.
+"""
+    )
+    group.add_option("--pmax", dest="pmax", metavar="PMAX", help="pmax to apply to the display of first FC [optional]", type='float', default=None)
+    group.add_option("--pmin", dest="pmin", metavar="PMIN", help="pmin to apply to the display of first FC [optional]", type='float', default=None)
+    group.add_option("--stretch", dest="stretch", metavar="STRETCH", help="Stretch to apply to the display of first FC [optional]\n", default=None)
+    group.add_option("--exponent", dest="exponent", metavar="EXPONENT", help="Exponent if stretch='power' [optional]\n", default=None)
+    parser.add_option_group(group)
+
+    group = OptionGroup(parser, "Custom Image/Survey",
+"""\
+The following parameters affect the way the 'default' images are displayed in the FCs.
+They correspond directly to the corresponding parameters without 'c_' prefixes of the APLPY FITSFigure.show_grayscale method, see:
+  https://aplpy.readthedocs.io/en/stable/api/aplpy.FITSFigure.html#aplpy.FITSFigure.show_grayscale
+for details.
+"""
+    )
+    group.add_option("--survey", dest="survey", metavar="SURVEY", help="An alternative SkyView survey to use. Use --supported_surveys to see which surveys are supported. [optional]", default=None)
+    group.add_option("--bkg_image", dest="bkg_image", metavar="BK_IMAGE", help="A custom background image to use [optional]", default=None)
+    group.add_option("--bkg_lam", dest="bkg_lam", metavar="BK_LAM", help="Description of the WL coverage of the custom background image/survey [optional]", default=None)
+    group.add_option("--c_pmax", dest="c_pmax", metavar="PMAX", help="pmax to apply to the display of first FC [optional]", type='float', default=None)
+    group.add_option("--c_pmin", dest="c_pmin", metavar="PMIN", help="pmin to apply to the display of first FC [optional]", type='float', default=None)
+    group.add_option("--c_stretch", dest="c_stretch", metavar="STRETCH", help="Stretch to apply to the display of first FC [optional]\n", default=None)
+    group.add_option("--c_exponent", dest="c_exponent", metavar="EXPONENT", help="Exponent if stretch='power' [optional]\n", default=None)
+    parser.add_option_group(group)
+
+    group = OptionGroup(parser, "Gaia pseudo images parameters",
+        "The following parameters affect the way the Gaia pseudo images (if any) are constructed."
+    )
+    group.add_option("--GAIA_im_IQ", dest="GAIA_im_IQ", metavar="FWHM",
+        help="FWHM to use for any Gaia images that will be generated [optional]",
+        type='float',
+        default=None
+    )
+    group.add_option("--GAIA_im_noise", dest="GAIA_im_noise", metavar="NOISE",
+        help=
+"""\
+GAIA_im_noise is a dimensionless factor, with the noise computed via the numpy method random.normal(loc,scale) with loc=0.0 and:
+    scale=GAIA_im_noise*<brightest_flux_value>
+with <brightest_flux_value> being the maximum pixel value in the noiseless pseudo image.
+The default value is (as at 2020-01): 1e-6
+Together with pmin and pmax this parameter gives (some) control over the relative appearance of stars of different brightness, and can thus be helpful to (for example) highlight faint, but not insignificant, field star(s).
+[optional] default=1e-6
+""",
+        type='float', default=None
+    )
+    group.add_option("--c_GAIA_im_IQ", dest="c_GAIA_im_IQ", metavar="FWHM",
+        help="Same as --GAIA_im_IQ but for any Gaia pseudo images associated with the custom FC.",
+        type='float',
+        default=None
+    )
+    group.add_option("--c_GAIA_im_noise", dest="c_GAIA_im_noise", metavar="NOISE",
+        help="Same as --GAIA_im_noise but for any Gaia pseudo images associated with the custom FC.",
+        type='float',
+        default=None
+    )
+    parser.add_option_group(group)
+
+    group = OptionGroup(parser, "p2 server selection, options and password management")
+    group.add_option("--env", dest="env", help="Use the %s environment [optional], default = production" %("|".join(p2api.API_URL.keys())), default='production')
+    group.add_option("-u", "--user", dest="user", metavar="USER", help="the User Portal account of the user to use, e.g. ASMITHSONIAN [required for production (not required for demo)]", default=None)
+    group.add_option("-S", "--PWstore", dest="pwstore", metavar="PWSTORE", help="Store the password in the keyring, default=False", action="store_true", default=False)
+    group.add_option("-R", "--PWreenter", dest="pwreenter", metavar="PWREENTER", help="Reenter the password, thus overriding any PW in the keyring, default=False", action="store_true", default=False)
+    group.add_option("-Z", "--PWremove", dest="pwremove", metavar="PWREMOVE", help="Remove any password from the keyring, default=False", action="store_true", default=False)
+    parser.add_option_group(group)
+
+    group = OptionGroup(parser, "help...")
+    group.add_option("--supported_surveys", dest="supSurs", metavar="Ids", help="Report the list of supported surveys",  action="store_true", default=False)
+    parser.add_option_group(group)
+
+    group = OptionGroup(parser, "Verbosity...",
+                    "All optional")
+    group.add_option("-D", "--debug", dest="debug", help="Write out debugging info [optional]", action="store_true", default=False)
+    #group.add_option("-v", "--verbose", dest="vbose", help="Write out verbose info [optional]", action="store_true", default=False)
+    group.add_option("-q", "--quiet", dest="quiet", help="Write out minimal info [optional]", action="store_true", default=False)
+    parser.add_option_group(group)
+
+    # parse arguments
+    (options, args) = parser.parse_args()
+
+    if not options.ids and not options.supSurs :
+        parser.print_help()
+        sys.exit(2)
+
+    self.ids = options.ids
+    self.obsdate = options.obsdate
+
+    self.pmin = options.pmin
+    self.pmax = options.pmax
+    self.stretch = options.stretch
+    self.exponent = options.exponent
+    self.GAIA_im_IQ = options.GAIA_im_IQ
+    self.GAIA_im_noise = options.GAIA_im_noise
+
+    self.survey = options.survey
+    self.bkg_image = options.bkg_image
+    self.bkg_lam = options.bkg_lam
+    self.c_pmin = options.c_pmin
+    self.c_pmax = options.c_pmax
+    self.c_stretch = options.c_stretch
+    self.c_exponent = options.c_exponent
+    self.c_GAIA_im_IQ = options.c_GAIA_im_IQ
+    self.c_GAIA_im_noise = options.c_GAIA_im_noise
+
+    self.user = options.user
+    self.env = options.env
+    self.pwstore = options.pwstore
+    self.pwreenter = options.pwreenter
+    self.pwremove = options.pwremove
+
+    self.supSurs = options.supSurs
+
+    self.debug = options.debug
+    #self.vbose = options.vbose
+    self.quiet = options.quiet
+
+    logFMT = "%(asctime)s %(module)15s[%(process)5d] [%(levelname)s] %(message)s"
+    logLVL=logging.INFO
+    #if options.vbose:
+    #  logLVL=logging.INFO
+    if options.debug:
+      logLVL=logging.DEBUG
+    if options.quiet:
+      logLVL=logging.WARN
+    logging.basicConfig(level=logLVL, format=logFMT)
+
+options = get_options(_VERSION_)
+# ------------------------------------------------------------------------------
+def oneOB( api, obId ) :
+    try :
+        OB, _  = api.getOB(obId)
+    except p2api.P2Error as e :
+        ## Looks like this ID is not an OB...
+        logging.error('OB %d could not be found...' %(obId))
+        return
+    except:
+        raise
+    
+    if OB['itemType'] not in ['OB',] :
+        logging.warning('OB %d is not a science OB' %(obId))
+        return
+
+    if OB['obStatus'] not in ['P','-'] :
+        logging.warning('OB status is %s, OB not modifiable' %(OB['obstatus']))
+        return
+
+    runId, _ = api.getRun(OB['runId'])
+    ob_url=api.apiUrl+'/obsBlocks/'+str(obId)
+    oneOB_sdt=datetime.datetime.now()
+    logging.info("Creating finding charts for OBid = %d, OBname = '%s' @ %s..." % (OB['obId'],OB['name'],OB['instrument']))
+    data = api.generateFindingChart(
+        obId,
+        obs_date=options.obsdate,
+        stretch=options.stretch,
+        exponent=options.exponent,
+        pmin=options.pmin,
+        pmax=options.pmax,
+        GAIA_im_IQ=options.GAIA_im_IQ,
+        GAIA_im_noise=options.GAIA_im_noise,
+        survey=options.survey,
+        bkg_image=options.bkg_image,
+        bkg_lam=options.bkg_lam,
+        c_stretch=options.c_stretch,
+        c_exponent=options.c_exponent,
+        c_pmin=options.c_pmin,
+        c_pmax=options.c_pmax,
+        c_GAIA_im_IQ=options.c_GAIA_im_IQ,
+        c_GAIA_im_noise=options.c_GAIA_im_noise,
+    )
+    oneOB_edt=datetime.datetime.now()
+    logging.info("Finding chart for OBid = %d, OBname = '%s' @ %s done in %s." % (obId, OB['name'], OB['instrument'], str(oneOB_edt-oneOB_sdt)))
+    if data is not None :
+        try :
+            jwarn=json.loads(r.content.decode('utf-8'))
+            if 'warnings' in jwarn.keys() :
+                logging.warning(jwarn['warnings'])
+        except :
+            pass
+# ------------------------------------------------------------------------------
+
+
+# ------------------------------------------------------------------------------
+# ------------------------------------------------------------------------------
+# ------------------------------------------------------------------------------
+
+def main():
+
+    ################################################################################
+    # remove passowrd from keyring...
+    if options.pwremove and options.user is not None :
+        keyring.delete_password('usd_p2api_interactive:www.eso.org', options.user)
+        sys.exit(0)
+    
+    ################################################################################
+    if options.supSurs :
+        r = requests.get('%s/supported_surveys' %(p2api.P2FC_URL[options.env]))
+        if r.status_code in [200,201,204] :
+            print(r.text)
+        else:
+            raise(r.text)
+        sys.exit(0)
+
+    #login to the API
+    api=p2api_utils.get_api(options)
+    obIds=p2api_utils.r_get_all_OBIds( api, options.ids.split(',') )
+
+    if len(obIds) > 0 :
+        if len(obIds) > 1 :
+            sdt=datetime.datetime.now()
+            logging.warning('Submitting %d OBs for FC generation...' %( len(obIds) ))
+        for obId in obIds :
+            oneOB(api, obId)
+        if len(obIds) > 1 :
+            edt=datetime.datetime.now()
+            logging.info('Total elapsed time to process %d finding charts %s...' %(len(obIds),str(edt-sdt)))
+    else :
+        logging.error('No OBs found.')
+
+if __name__ == '__main__' :
+    main()