FetchSyn

#!/usr/bin/env python
#
# FetchSyn
#
# For help type:  FetchSyn --u
#
# FetchSyn is a python command line client to download synthetic
# seismograms from the IRIS Synthetics Engine (Syngine).
# Find the most current version at http://service.iris.edu/clients/
#
# SYNGINE OVERVIEW:           http://ds.iris.edu/ds/products/syngine
# SYNGINE MODELS:             http://ds.iris.edu/ds/products/syngine/#models
# SYNGINE PARAMETERS:         http://service.iris.edu/irisws/syngine/1
#
# Python Dependencies:
#    *must have the requests package
#    *currently written for python2 and python3
#
# FetchSyn usage: A source and receiver is required, everything else has a default.
#
# RECEIVER: one of the following is required
#  A) rlat     Receiver longitude
#     rlon     Receiver latitude 
#  B) sta,S    Station name.  use * to wildcard.
#     net,N    Network name
#                  Notes: *if only a net is given, then all sta will be included
#                         *virtual networks are valid networks, e.g. _GSN
#  C) recfile  Textfile with a list of receivers as (NET STA) or
#                  (RLAT RLON [netcode stacode loccode])
#     Valid examples in a recfile:
#     IU ANMO 
#     II *
#     -89 -179
#     -89 -178 S0002
#     -89 -177 N2 S0003 L2
#     -89 -176 netcode=N2 stacode=S0004 loccode=L3
#
# SOURCES:  one of the following is required
#  A) evid     Event id [catalog]:[eventid] 
#  B) slat     source latitude
#     slon     source longitude
#     sdepm    source depth in meters
#   b1) mt     source moment tensor in NM: mrr,mtt,mpp,mrt,mrp,mtp
#   b2) dc     source double couple:  strike,dip,rake [,scalar moment in NM]
#   b3) force  source force in NM: Fr,Ft,Fp
#
# ORIGIN/START/END TIME (optional)
#     o      Origin Time. Must be absolute (YYYY-MM-DD,HH:MM:SS) or (YYYY,MM,DD,HH,MM,SS)
#     s      Start Time. Either 1) absolute time (YYYY-MM-DD,HH:MM:SS)
#                               2) phase-relative offset e.g. P+10 or ScS-100
#                               3) an offset from origin in seconds
#     e      End Time. Either   1) absolute time (YYYY-MM-DD,HH:MM:SS)
#                               2) phase-relative offset e.g. P+10 or ScS-100
#                               3) an offset (duration) from start time in seconds
#
#
# OUTPUT:
#  Output files are either in miniSEED (mseed) or a zip archive of SAC
#      files (saczip).  SAC files will have the headers populated.
#
# QUICK EXAMPLES
#  FetchSyn -recfile rec.txt -slat 38.3 -slon 142.3 -sdepm 24400 -dc 203,10,88,5.3e22
#  FetchSyn -N IU -S ANMO,KIP -evid GCMT:M201103110546A
#  FetchSyn -N IU -S ANMO,KIP -evid GCMT:M201103110546A -sdepm 35000
#  FetchSyn -N _GSN -evid GCMT:M201103110546A -model iasp91_2s -dt 0.05 -label TOHOKU
#           -format mseed -C ZRT -units velocity -s P-10 -e 300
#
# ## Change history ##
#
#  2015.343:
#  -Initial release
#
#  2015:345:
#  -Removed modelname validation, since server already does that.
#  -added .strip() to slat & slon to remove whitespace when negative values added
#  -updated error messaging
#
#  2015.352:
#  -added support for the doublecouple, momenttensors & force sources.
#  -added .strip() to rlat,rlon,dc,mt, and force to remove whitespace when negative 
#   values added
#  -modified argparse flags to accept single dash arguments
#
#  2016.007
#  -added user-agent info sent with each request as well as user-added agent info.
#
#  Author: Alex Hutko, IRIS Data Management Center

from __future__ import print_function

FSversion = '2016.007'

import requests
import datetime
import re
import sys
import time
import argparse
import csv
try:
    import StringIO
    pythonversion = 2
except ImportError:
    from io import BytesIO
    pythonversion = 3

########################
# Used to check if an input number is an actual number and within bounds
def check_number(s,min,max):
    try:
        float(s)
        s=float(s)
        if ( s >= min and s <= max ):
            ierror = 0
        else:
            ierror = 1
    except:
        ierror = 2
    return ierror
########################

CurrentTime0 = time.time()

for i, arg in enumerate(sys.argv):
    if (arg[0] == '-') and arg[1].isdigit(): sys.argv[i] = ' ' + arg

parser = argparse.ArgumentParser()

parser.add_argument('-A','--A', action='store', dest='agent',help='UserAgent info')

parser.add_argument('-recfile','--recfile','--ReceiverFile', action='store', dest='ReceiverFile',help='receiver file name with NET STATION')
parser.add_argument('-u','--u', action='store_true', dest='helpmenu',help='extended HELP MENU')
parser.add_argument('-outfile','--outfile', action='store', dest='outfile',help='output file label' )

parser.add_argument('-model','--model', action='store', dest='model',help='model from list: http://ds.iris.edu/ds/products/syngine' )
parser.add_argument('-format','--format', action='store', dest='formatstring',help='miniseed or saczip')
parser.add_argument('-label','--label', action='store', dest='label',help='label')
parser.add_argument('-C','--C','--components', action='store', dest='components',help='components: any combination of ZNERT')
parser.add_argument('-units','--units', action='store', dest='units',help='displacement, velocity, or acceleration (or d,v,a)')
parser.add_argument('-scale','--scale', action='store', dest='scale',help='amplitude scaling factor')
parser.add_argument('-dt','--dt', action='store', dest='dt',help='sampling interval in sec')
parser.add_argument('-kernelwidth','--kernelwidth', action='store', dest='kernelwidth',help='width of sinc resampling kernel')

parser.add_argument('-origin','-origintime','--origin','--origintime', action='store', dest='originstring',help='end time YYYY-MM-DDTHH:MM:SS')
parser.add_argument('-s','-start','-starttime','--s','--start','--starttime', action='store', dest='startstring',help='source time YYYY-MM-DDTHH:MM:SS or PHASE+-offset')
parser.add_argument('-e','-end','-endtime','--e','--end','-endtime', action='store', dest='endstring',help='end time YYYY-MM-DDTHH:MM:SS or PHASE+-offset or offset_from_start')

parser.add_argument('-N','-network','--N','--network', action='store', dest='network',help='network')
parser.add_argument('-S','-station','--S','--station', action='store', dest='station',help='station')
parser.add_argument('-rlat','--rlat', action='store', dest='rlat',help='reciver lat')
parser.add_argument('-rlon','--rlon', action='store', dest='rlon',help='receiver lon')
parser.add_argument('-Ncode','-networkcode','--Ncode','--networkcode', action='store', dest='netcode',help='networkcode')
parser.add_argument('-Scode','-stationcode','--Scode','--stationcode', action='store', dest='stacode',help='stationcode')
parser.add_argument('-Lcode','-locationcode','--Lcode','--locationcode', action='store', dest='loccode',help='locationcode')

parser.add_argument('-evid','-eventid','--evid','--eventid', action='store', dest='evid',help='event id')
parser.add_argument('-slat','--slat', action='store', dest='slat',help='source lat')
parser.add_argument('-slon','--slon', action='store', dest='slon',help='source lon')
parser.add_argument('-sdepm','--sdepm', action='store', dest='sdepm',help='source depth in meters')

parser.add_argument('-mt','--mt','--sourcemomenttensor', action='store', dest='momentstring',help='moment: mrr,mtt,mpp,mrt,mrp,mtp in N-m (1 J = 1 N-m = 1e7 dyne-cm).  Comma delimited, no spaces or brackets.')
parser.add_argument('-dc','--dc','--sourcedoublecouple', action='store', dest='doublecouplestring',help='strike,dip,rake[,ScalarMoment_in_NM]')
parser.add_argument('-force','--force','-forces','--forces','--sourceforce', action='store', dest='forcestring',help='force source in NM: Fr,Ft,Fp')

results = parser.parse_args()
helpmenu = results.helpmenu
agent = results.agent
ReceiverFile = results.ReceiverFile
evid = results.evid
slat = results.slat
slon = results.slon
sdepm = results.sdepm
rlat = results.rlat
rlon = results.rlon
networkinput = results.network
stationinput = results.station
netcode = results.netcode
stacode = results.stacode
loccode = results.loccode
network = results.network
station = results.station
momentstring = results.momentstring
doublecouplestring = results.doublecouplestring
forcestring = results.forcestring
startstring = results.startstring
originstring = results.originstring
endstring = results.endstring
components = results.components
units = results.units
scale = results.scale
dt = results.dt
kernelwidth = results.kernelwidth
modelstring = results.model
formatstring = results.formatstring
outfile = results.outfile
label = results.label

#---------- help menu

if ( helpmenu is True  ):
    print ('')
    print ('FetchSyn: collect Syngine synthetic seismograms (version 2015.300)')
    print ('http://service.iris.edu/clients')
    print ('')
    print ('Usage: FetchSyn [options]')
    print ('')
    print (' EXAMPLES')
    print (' FetchSyn -recfile rec.txt -slat 38.3 -slon 142.3 -sdepm 24400 -dc 203,10,88,5.3e22')
    print (' FetchSyn -N IU -S ANMO,KIP -evid GCMT:M201103110546A')
    print (' FetchSyn -N _GSN -evid GCMT:M201103110546A -model iasp91_2s -dt 0.05 -label TOHOKU')
    print ('          -format mseed -C ZRT -units velocity -s P-10 -e 300')
    print (' ')
    print ('Options:         ')
    print ('------------- At least one receiver required -----')
    print ('-recfile            Text file with list of receivers (NET STA) or ')
    print ('                         (RLAT RLON [netcode stacode loccode])')
    print ('-N, -net            Network name.  Virtual nets valid e.g. : _GSN')
    print ('-S, -sta            Station name, or comma separated names')
    print ('-rlat               Receiver latitude')
    print ('-rlon               Receiver longitude')
    print ('-Ncode,-networkcode    Network name when rlat/rlon used')
    print ('-Scode,-stationcode    Station name when rlat/rlon used')
    print ('-Lcode,-locationcode   Location code for synthetic')
    print ('------------- A source is required ---------------')
    print ('-evid               Event id [catalog]:[eventid] e.g. GCMT:C201002270634A' )
    print ('-slat               Source latitude')
    print ('-slon               Source longitude ')
    print ('-sdepm              Source depth in meters ')
    print ('-mt                 Source moment tensor in NM (mrr,mtt,mpp,mrt,mrp,mtp)')
    print ('-dc                 Source double couple (strike,dip,rake [,scalar moment])')
    print ('-force              Source force in NM (Fr,Ft,Fp)')
    print ('------------- Optional --------------------------')
    print ('-model              Model name  http://ds.iris.edu/ds/products/syngine/#models')
    print ('-format             Either miniseed or saczip for output')
    print ('-label              Apply a label to be included in file names')
    print ('-C, -components     Orientation, any combination of ZNERT')
    print ('-units              Either (D)displacement, (V)velocity or (A)acceleration')
    print ('-dt                 Sampling interval in seconds')
    print ('-kernelwidth        Width of the sinc kernel used for resampling')
    print ('-o, -origintime     Source origin time (YYYY-MM-DD,HH:MM:SS)')
    print ('-s, -starttime      Trace start time (YYYY-MM-DD,HH:MM:SS) or phase-+offset e.g. P+10')
    print ('                         or offset from origin time in seconds')
    print ('-e, -endtime        Trace end time (YYYY-MM-DD,HH:MM:SS) or phase-+offset or duration_in_s')
    print ('-A                  Add a user-agent string for identification ')
    print (' ')

#---------- User-agent info

agentstatement = 'FetchSyn/' + FSversion + ' Python/' + sys.version.split(' ')[0] + ' Requests/' + requests.__version__
if ( agent is not None ):
    agentstatement = agentstatement + ' (' + str(agent) + ')'
user_agent = {'User-Agent': str(agentstatement) }

#---------- Make sure there is a source and a receiver

if ( rlat is None and rlon is None and network is None and station is None and ReceiverFile is None ):
    print ('Error: a receiver is required ')
    exit()
if ( slat is None and slon is None and evid is None ):
    print ('Error: a source is required ')
    exit()

#---------- validate format

if ( formatstring is not None ):
    format_list = ['mseed','miniseed','ms','sac','SAC','zip','saczip']
    if ( formatstring not in format_list) :
        print ('Error: invalid format type:',formatstring,' Must be miniseed or saczip ')
        exit()
    if ( formatstring[0] == 'm' ):
        formatstring = 'miniseed'
    else:
        formatstring = 'saczip'

#---------- validate components

if ( components is not None ):
    component_list = ['Z','N','E','R','T']
    while re.search(r'([A-Z])(.*)\1', components):
        components = re.sub(r'([A-Z])(.*)\1', r'\1\2', components)
    for i in range(len(components)):
        if ( components[i] not in component_list) :
            print ('Error: invalid components:',components,' Must be any combination of ZNERT ')
            exit()

#---------- validate units

if ( units is not None ):
    unit_list = ['displacement','velocity','acceleration','d','v','a','D','V','A']
    if ( units not in unit_list) :
        print ('Error: invalid unit type:',unit,' Must be displacement,velocity,acceleration,d,v,or a ')
        exit()
    if ( units == 'd' or units == 'D' ):
        units = 'displacement'
    if ( units == 'v' or units == 'V' ):
        units = 'velocity'
    if ( units == 'a' or units == 'A' ):
        units = 'acceleration'

#---------- validate scale

if ( scale is not None ):
    scalecheck = check_number(scale,-1e26,1e26)
    if ( scalecheck != 0 ):
        print ('Error: bad value for scale:',scale,' Must be between -1e26 and 1e26')
        exit()

#--------- validate dt

if ( dt is not None ):
    dtcheck = check_number(dt,0.001,5)
    if ( dtcheck != 0 ):
        print ('Error: bad value for dt:',dt,' Must be between 0.001 and 5')
        exit()

#---------- validate kernelwidth

if ( kernelwidth is not None ):
    kernelwidthcheck = check_number(kernelwidth,1,1000)
    if ( kernelwidthcheck != 0 ):
        if ( isinstance(kernelwidth,int) == False ):
            print ('Error: bad value for kernelwidth:',kernelwidth,' Must be positive integer')
            exit()
    kernelwidth = str(int(float(kernelwidth)))

#---------- validate source location: slat,slon,sdepm

if ( slat is not None ):
    slatcheck = check_number(slat,-90,90)
    if ( slatcheck != 0 ):
        print ('Error: bad value for slat',slat,' Must be between -90 and 90.')
        exit()

if ( slon is not None ):
    sloncheck = check_number(slon,-180,360)
    if ( sloncheck != 0 ):
        print ('Error: bad value for slon',slon,' Must be between -180 and 360.')
        exit()
    if ( float(slon) > 180 ):
       slonF = float(slon) - 360
       slon = str(slonF)

if ( sdepm is not None ):
    sdepmcheck = check_number(sdepm,0,700000)
    if ( sdepmcheck != 0 ):
        print ('Error: bad value for sdepm',sdepm,' Must be between 0 and 700000.')
        exit()

if ( slat is not None and slon is not None and sdepm is None ):
    print ('Error: sourcedepthinmeters (sdepm) is required when using slat & slon ')
    exit()

#---------- take care of the origntime if given

if ( originstring is not None ):
    origin = originstring.replace(':',',')
    origin = origin.replace('T',',')
    origin = origin.replace('-',',')
    origin = origin.replace('/',',')
    origin = origin.split(',')
    if ( len(origin) != 6):
        print ('Error: invalid origintime format. Use: 1900,1,31,23,59,59 or 1900-12-31T23:59:59')
        exit()
    else:
        origin = origin[0] + "-" + origin[1] + "-" + origin[2] + "T" + origin[3] + ":" + origin[4] + ":" + origin[5]

#---------- take care of the starttime if given

if ( startstring is not None ):
    start = startstring.replace(':',',')
    start = start.replace('T',',')
    start = start.replace('/',',')
    start = start.replace('-',',')
    start = start.replace('+',',')
    start = start.split(',')
    if ( len(start) == 6 ):
        start = start[0] + "-" + start[1] + "-" + start[2] + "T" + start[3] + ":" + start[4] + ":" + start[5]
    elif ( len(start) == 2 ):
        startcheck = check_number(start[1],0,18000)
        if ( startcheck !=0 ):
            print ('Error: invalid starttime:',startstring,' Valid: 1900,01,31,23,59,0 or 1900-01-31T23:59:0 or P-20 or ScS+200 or 100')
            exit()
        else:
            start=startstring
            start = start.replace('+','%2B')
    elif ( len(start) == 1):
        startcheck = check_number(start[0],0,18000)
        if ( startcheck !=0 ):
            print ('Error: invalid starttime:',startstring,' Valid: 1900,01,31,23,59,0 or 1900-01-31T23:59:0 or P-20 or ScS+200 or 100')
            exit()
        else:
            start=startstring
    else:
        print ('Error: invalid starttime:',startstring,' Valid: 1900,01,31,23,59,0 or 1900-01-31T23:59:0 or P-20 or ScS+200 or 100')
        exit()

#---------- take care of the endtime if given

if ( endstring is not None ):
    end = endstring.replace(':',',')
    end = end.replace('T',',')
    end = end.replace('/',',')
    end = end.replace('-',',')
    end = end.replace('+',',')
    end = end.split(',')
    if ( len(end) == 6 ):
        end = end[0] + "-" + end[1] + "-" + end[2] + "T" + end[3] + ":" + end[4] + ":" + end[5]
    elif ( len(end) == 2 ):
        endcheck = check_number(end[1],0,18000)
        if ( endcheck !=0 ):
            print ('Error: invalid endtime:',endstring,' Valid: 1900,01,31,23,59,0 or 1900-01-31T23:59:0 or P-20 or ScS+200 or 100')
            exit()
        else:
            end=endstring
            end = end.replace('+','%2B')
    elif ( len(end) == 1 ):
        endcheck = check_number(end[0],0,18000)
        if ( endcheck !=0 ):
            print ('Error: invalid endtime:',endstring,' Valid: 1900,01,31,23,59,0 or 1900-01-31T23:59:0 or P-20 or ScS+200 or 100')
            exit()
        else:
            end=endstring
    else:
        print ('Error: invalid endtime:',endstring,' Valid: 1900,01,31,23,59,0 or 1900-01-31T23:59:0 or P-20 or ScS+200 or 100')
        exit()

#---------- validate sourcemomenttensor

if ( momentstring is not None ):
    moment = momentstring.split(',')
    if ( len(moment) != 6 ):
        print ('Error: sourcemomenttensor must have 6 comma separated numbers')
        exit()
    else:
        for i in range(0,len(moment)):
            momentcheck = check_number(moment[i],-1e26,1e26)
            if ( momentcheck != 0 ):
                print ('Error: bad value for moment component ',i+1,':',moment[i],' Must be between -1e26 and 1e26 (Nm)')
                exit()

#---------- validate sourcedoublecouple

if ( doublecouplestring is not None ):
    doublecouple = doublecouplestring.split(',')
    if ( len(doublecouple) == 3 or len(doublecouple) == 4 ):
        strikecheck = check_number(doublecouple[0],-180,360)
        if ( strikecheck != 0 ):
            print ('Error: bad value for strike:',doublecouple[0],' Must be between -180 and 360')
            exit()
        dipcheck = check_number(doublecouple[1],0,90)
        if ( dipcheck != 0 ):
            print ('Error: bad value for dip:',doublecouple[1],' Must be between 0 and 90')
            exit()
        rakecheck = check_number(doublecouple[2],-180,360)
        if ( rakecheck != 0 ):
            print ('Error: bad value for rake:',doublecouple[2],' Must be between -180 and 180')
            exit()
        if ( len(doublecouple) == 4 ):
            Mocheck = check_number(doublecouple[3],0,1e26)
            if ( Mocheck != 0 ):
                print ('Error: bad value for scalar moment :',doublecouple[3],' Must be between 0 and 1e26 (Nm)')
                exit()
    else:
        print ('Error: sourcedoublecouple must have 3 or 4 comma separated numbers')
        exit()

#---------- validate sourceforce

if ( forcestring is not None ):
    force = forcestring.split(',')
    if ( len(force) != 3 ):
        print ('Error: sourceforce must have 3 comma separated numbers')
        exit()
    for i in range(0,len(force)):
        forcecheck = check_number(force[i],-1e26,1e26)
        if ( forcecheck != 0 ):
            print ('Error: bad value for force component ',i+1,':',force[i],' Must be between -1e26 and 1e26 (Nm)')
            exit()

#---------- validate rlat & rlon

if ( rlat is not None ):
    rlatcheck = check_number(rlat,-90,90)
    if ( rlatcheck != 0 ):
        print ('Error: bad value for rlat',rlat,' Must be between -90 and 90.')
        exit()

if ( rlon is not None ):
    rloncheck = check_number(rlon,-180,360)
    if ( rloncheck != 0 ):
        print ('Error: bad value for rlon',rlon,' Must be between -180 and 360.')
        exit()
    if ( float(rlon) > 180 ):
        rlonF = float(rlon) - 360
        rlon = str(rlonF)

if ( stacode is None ):
    stacode = 'S0001'

#---------- validate the network(s)

if ( network is not None ):
    network = network.split(',')
    Nnetworks = len(network)
    if ( Nnetworks == 1 ):
        netstring = 'http://service.iris.edu/fdsnws/station/1/query?level=network&format=text&network=' + network[0]
        r=requests.get(netstring,headers=user_agent)
        if ( r.status_code != 200 ):
            print ('Error: invalid network: ', network)
            exit()
    else:
        for i in range(0,Nnetworks):
            r=requests.get('http://service.iris.edu/fdsnws/station/1/query?level=network&format=text&network=' + network[i],headers=user_agent)
            if ( r.status_code != 200 ):
                print ('Error: invalid network: ',network[i])
                exit()

#---------- validate the station(s)

if ( station is not None ):
    station = station.split(',')
    Nstations = len(station)
    if ( Nstations == 1 ):
        stastring = 'http://service.iris.edu/fdsnws/station/1/query?level=station&format=text&network=' + network[0] + '&station=' + station[0]
        r=requests.get(stastring,headers=user_agent)
        if ( r.status_code != 200 ):
            print ('Error: invalid net station: ',network[0], station)
            exit()
    else:
        for i in range(0,Nstations):
            r=requests.get('http://service.iris.edu/fdsnws/station/1/query?level=station&format=text&network=' + network[0] + '&station=' + station[i],headers=user_agent)
            if ( r.status_code != 200 ):
                print ('Error: invalid net station: ',network[0], station[i])
                exit()
    if ( station[0] == "*" ):
        station = None;

#---------- read in list of stations

if ( ReceiverFile is not None ):
    f = open(ReceiverFile,'r')
    i = 0
    ListReceivers = []
    for line in f:
        i=i+1
        fields = line.split()
        if ( len(fields) < 2 ):
            print ('Error: invalid ReceiverFile line: ',line,' Must have at least (NET STAT) or (rlat rlon)')
            exit()
        else:
            LR1 = fields[0]
            LR2 = fields[1]
            if any( c in "ABCDEFGHIJKLMNOPQRSTUVWXYZ" for c in LR1 ):
                stastring = 'http://service.iris.edu/fdsnws/station/1/query?level=station&format=text&network=' + LR1 + '&station=' + LR2
                r=requests.get(stastring,headers=user_agent)
                if ( r.status_code != 200 ):
                    print ('Error: invalid net station: ',LR1, LR2,' on line ',i,' of ',ReceiverFile)
                    exit()
                ListReceivers.append(line)
            else:
                LR1check = check_number(LR1,-90,90)
                if ( LR1check != 0 ):
                    print ('Error: bad value for rlat:',LR1,' on line ', i,' of ',ReceiverFile, ' Must be between -90 and 90.')
                LR2check = check_number(LR2,-180,360)
                if ( LR2check != 0 ):
                    print ('Error: bad value for rlon:',LR2,' on line ', i,' of ',ReceiverFile, ' Must be between -180 and 360.')
                LR_Lon = float(LR2)
                if ( LR_Lon > 180 ):
                    LR_Lon = LR_Lon - 360
                Nfields = len(fields)
                NewLine = LR1 + " " + str(LR_Lon).strip()
                if any ( c in "code" for c in line ):
                    if ( Nfields == 3 ):
                        NewLine = NewLine + " " + fields[2]
                    if ( Nfields == 4 ):
                        NewLine = NewLine + " " + fields[2] + " " + fields[3]
                    if ( Nfields == 5 ):
                        NewLine = NewLine + " " + fields[2] + " " + fields[3] + " " + fields[4]
                else:
                    if ( Nfields == 3 ):
                         NewLine = NewLine + " stacode=" + fields[2]
                    if ( Nfields == 4 ):
                        NewLine = NewLine + " netcode=" + fields[2] + " stacode=" + fields[3]
                    if ( Nfields == 5 ):
                        NewLine = NewLine + " netcode=" + fields[2] + " stacode=" + fields[3] + " loccode=" + fields[4]
                if ( Nfields == 2 ):
                    NewLine = NewLine + " stacode=S" + str(i).zfill(4)
                ListReceivers.append(NewLine)
    NReceivers = i

#------------- Build the POST file for the webservice query -----------
#-------- Two versions here since python 2 vs 3 handle StrinIO differently -----

url = ""

if ( modelstring is not None):
    url = url + "\nmodel=" + modelstring

if ( formatstring is not None ):
    url = url + "\nformat=" + formatstring

if ( label is not None ):
    url = url + "\nlabel=" + label

if ( components is not None ):
    url = url + "\ncomponents=" + components

if ( units is not None ):
    url = url + "\nunits=" + units

if ( scale is not None ):
    url = url + "\nscale=" + scale

if ( dt is not None ):
    url = url + "\ndt=" + dt

if ( kernelwidth is not None ):
    url = url + "\nkernelwidth=" + kernelwidth

if ( originstring is not None ):
    url = url + "\norigintime=" + origin

if ( startstring is not None ):
    url = url + "\nstarttime=" + start

if ( endstring is not None ):
    url = url + "\nendtime=" + end

if ( evid is not None ):
    url = url + "\neventid=" + evid

if ( slat is not None ):
    url = url + "\nsourcelatitude=" + slat.strip()

if ( slon is not None ):
    url = url + "\nsourcelongitude=" + slon.strip()

if ( sdepm is not None ):
    url = url + "\nsourcedepthinmeters=" + sdepm

if ( momentstring is not None ):
   url = url + "\nsourcemomenttensor=" + momentstring.strip()

if ( doublecouplestring is not None ):
   url = url + "\nsourcedoublecouple=" + doublecouplestring.strip()

if ( forcestring is not None ):
   url = url + "\nsourceforce=" + forcestring.strip()

if ( network is not None ):
    if ( station is not None ):
        for i in range(0,Nstations):
            url = url + "\n" + network[0] + " " + station[i]
    else:
        for i in range(0,Nnetworks):
            url = url + "\n" + network[i] + " *" 

if ( rlat is not None and rlon is not None ):
    url = url + "\n" + rlat.strip() + " " + rlon.strip()
    if ( netcode is not None ):
        url = url + " netcode=" + netcode + " "
    if ( stacode is not None ):
        url = url + " stacode=" + stacode + " "
    if ( loccode is not None ):
        url = url + " loccode=" + loccode 

if ( ReceiverFile is not None ):
    for i in range(0,NReceivers):
        url = url + "\n" + ListReceivers[i]

print ('-------- your POST file ---------')
print (url)
print ('---------------------------------')

if ( pythonversion == 2 ):
    inputstring = StringIO.StringIO(url)
elif ( pythonversion == 3 ):
    url = bytes(url,encoding="ascii")
    inputstring = BytesIO(url)
else:
    print ('StringIO is incompatible with this version of python. Try v2.X or v3')
    exit()

#------------- Build the output file name ------------

if ( outfile is not None ):
    outputfile = outfile + ".zip"
    if ( formatstring == 'miniseed' ):
        outputfile = outfile + ".mseed"
elif ( label is not None ):
    outputfile = label + ".zip"
    if ( formatstring == 'miniseed' ):
        outputfile = label + ".mseed"
else:
    outputfile = 'Synthetics.zip'
    if ( formatstring == 'miniseed' ):
        outputfile = 'Synthetics.mseed'

#------------- Make the request --------------

with open(outputfile, 'wb') as output:
    r = requests.post(url="http://service.iris.edu/irisws/syngine/1/query", data=inputstring,headers=user_agent)
    if ( r.status_code == 200 ):
        print ('status = 200.  Successful request.  Writing output to: ',outputfile)
        output.write(r.content)
        if ( formatstring is not 'miniseed' ):
            print ('Logfile:  Synthetics.log')
    else:
        print ('Request was not served. ',r.status_code,r.content)

now = datetime.datetime.utcnow()
print (now.strftime("%Y-%m-%d %H:%M:%S"), 'UTC')

CurrentTime1 = time.time()
print ('Request took: ',CurrentTime1-CurrentTime0,'seconds')