"""
This programme makes the planeflight*.dat files required to output for specific locations and times in the model.
NOTES:
- This programme can be used to produce files to output data for ship and aricraft campaigns
"""
# --- Packages
import numpy as np
from time import gmtime, strftime
import time
import glob
import AC_tools as AC
# --- Settings
try:
wd = AC.get_dir( 'dwd' )
except:
wd = './'
# the directory where files of required output locations are (one per UTC day)
camp = './' # Set location for output files here
# Set first and last year to look
start_year, end_year = 2013, 2015
# debug the oubtput
debug = True
# tag to (up to 4 characters)
tag= 'CON' #'TRB'#'MAL'#'ANT'
# do the altitude values need converting from input files ( must be in hPa)
convert_km_2_hPa, convert_m_2_hPa, convert_2_m = False, False, False
time_str ='%H:%M' # '%H%M' # '%H:%M:%S' # # '%h/%m/%s',
# Which (halogen) code version is being used?
#ver = '1.6' # Iodine simulation in v9-2
# --- Packages
import AC_tools as AC
import numpy as np
from time import gmtime, strftime
import time
import glob
# --- Settings
res='0.25x0.3125'
#res='0.5x0.666'
tpwd = AC.get_dir('tpwd')
#start_year, end_year = 2006,2007
#start_year, end_year = 2012,2013
#start_year, end_year = 2014,2017
start_year, end_year = 2015,2017
debug = False
#ver='1.7'
ver='3.0'
pf_locs_file = 'EU_GRID_{}.dat'.format( res )
#pf_locs_file ='ROW_GRID.dat'
# --- Read in site Detail
location = AC.readin_gaw_sites( pf_locs_file, all=True )
numbers, locs, lats, lons, pres = [ location[:,i] for i in range(5) ]
lats, lons, pres = [ np.float64(i) for i in lats, lons, pres ]
locs = np.array(locs)
print lats[0:4]
# --- Set Variables
#slist = pf_var( 'slist_v9_2_NREA_red_NOy', ver=ver )#'slist_v9_2_NREA_red' )
#slist = pf_var( 'slist_ClearFlo', ver=ver )
"""
This programme makes the planeflight*.dat files required to output for specific locations and times in the model.
NOTES:
- This programme can be used to produce files to output data for ship and aricraft campaigns
"""
# --- Packages
import numpy as np
from time import gmtime, strftime
import time
import glob
import AC_tools as AC
# --- Settings
try:
wd = AC.get_dir('dwd')
except:
wd = './'
# the directory where files of required output locations are (one per UTC day)
camp = './' # Set location for output files here
# Set first and last year to look
start_year, end_year = 2013, 2015
# debug the oubtput
debug = True
# tag to (up to 4 characters)
tag = 'CON' #'TRB'#'MAL'#'ANT'
# do the altitude values need converting from input files ( must be in hPa)
convert_km_2_hPa, convert_m_2_hPa, convert_2_m = False, False, False
time_str = '%H:%M' # '%H%M' # '%H:%M:%S' # # '%h/%m/%s',
# Which (halogen) code version is being used?
#ver = '1.6' # Iodine simulation in v9-2