staDb_path=staDb_path,
years_list=years_list,
num_cores=num_cores,
tmp_dir=tmp_dir,
project_name=project_name,
IGS_logs_dir=IGS_logs_dir,
blq_file=blq_file,
VMF1_dir=VMF1_dir,
pos_s=pos_s,
wetz_s=penna_wetz_list[i],
PPPtype=PPPtype,
ionex_type=ionex_type,
IONEX_products=IONEX_products,
rate=rate,
gnss_products_dir=gnss_products_dir,
eterna_path=eterna_path,
hardisp_path=hardisp_path,
rnx_dir=rnx_dir,
tree_options=tree_options_code,
tqdm=False,
command='gd2e();kinematic_project.gather_mGNSS()')
qsub_python_code(code,
name='{}{}'.format(project_name, str(i)),
email='*****@*****.**',
cleanup=False,
pbs_base='/scratch/bogdanm/pbs',
walltime='02:00:00')
#gen_tropNom can not be run rhis way as we need all the stations to be present
# kinematic_project.get_drInfo() .gather_drInfo()
# 'gd2e()' was executed!!! ALL OK
num_cores=num_cores,
tmp_dir=tmp_dir,
project_name=project_name,
IGS_logs_dir=IGS_logs_dir,
blq_file=blq_file,
VMF1_dir=VMF1_dir,
pos_s=pos_s,
wetz_s=wetz_s,
PPPtype=PPPtype,
ionex_type=ionex_type,
IONEX_products=IONEX_products,
rate=rate,
gnss_products_dir=gnss_products_dir,
eterna_path=eterna_path,
hardisp_path=hardisp_path,
rnx_dir=rnx_dir,
hatanaka=hatanaka,
cddis=cddis,
tree_options=tree_options_code,
tqdm=False,
command=
'dr_merge();kinematic_project.gd2e();kinematic_project.gather_mGNSS()')
qsub_python_code(code,
name='{}{}{}'.format(project_name,
str(ElMin) if ElMin != 7 else '',
str(i)),
email='*****@*****.**',
cleanup=False,
pbs_base=pbs_base)
years_list=years_list,
num_cores=num_cores,
tmp_dir=tmp_dir,
project_name=project_name,
IGS_logs_dir=IGS_logs_dir,
blq_file=blq_file,
VMF1_dir=VMF1_dir,
pos_s=pos_s,
wetz_s=wetz_s,
PPPtype=PPPtype,
ionex_type=ionex_type,
IONEX_products=IONEX_products,
rate=rate,
cddis=cddis,
gnss_products_dir=gnss_products_dir,
eterna_path=eterna_path,
hardisp_path=hardisp_path,
rnx_dir=rnx_dir,
hatanaka=hatanaka,
tree_options=tree_options_code,
tqdm=False,
command=
'dr_merge();kinematic_project.gd2e();kinematic_project.gather_mGNSS()')
qsub_python_code(code,
name='{}{}{}'.format(project_name, str(i),
ElMin if ElMin != 7 else ''),
email='*****@*****.**',
cleanup=False,
pbs_base=pbs_base,
walltime='06:00:00')
staDb_path = gen_staDb(tmp_dir = tmp_dir, project_name = project_name, stations_list = stations_list, IGS_logs_dir = IGS_logs_dir)
pbs_base = _os.path.join('/scratch/bogdanm/pbs',project_name) #break down by project folders as gets slow on hpc with multiple files
for i in range(len(penna_pos_s_list)):
pos_s = penna_pos_s_list[i]
prepare_dir_struct_dr(begin_year=_np.min(years_list), end_year = _np.max(years_list),tmp_dir=tmp_dir) #prepare dir struct for dr files
project_name_construct = _project_name_construct(project_name=project_name,PPPtype=PPPtype,pos_s=pos_s,wetz_s=wetz_s,tropNom_input=tropNom_input,ElMin=ElMin,ambres=ambres)
prepare_dir_struct_gathers(tmp_dir=tmp_dir,project_name=project_name_construct)
code = gen_code(stations_list = stations_list, cache_path = cache_path,tropNom_input=tropNom_input, ambres = ambres,ElMin=ElMin, ElDepWeight=ElDepWeight,
staDb_path = staDb_path,years_list=years_list,num_cores=num_cores,tmp_dir=tmp_dir,project_name=project_name,IGS_logs_dir=IGS_logs_dir,blq_file=blq_file,VMF1_dir = VMF1_dir,
hatanaka=hatanaka,cddis=cddis,
pos_s = pos_s,
wetz_s = wetz_s,
PPPtype = PPPtype,ionex_type=ionex_type,IONEX_products = IONEX_products,rate = rate,
gnss_products_dir = gnss_products_dir,eterna_path=eterna_path,hardisp_path = hardisp_path,rnx_dir=rnx_dir,tree_options = tree_options_code,tqdm=False,
command='gd2e();kinematic_project.gather_mGNSS()')
qsub_python_code(code,name='{}pos_{}'.format(project_name,str(i)),email='*****@*****.**',cleanup=False,pbs_base = pbs_base, walltime = '02:00:00')
#single static run
code = gen_code(stations_list = stations_list, cache_path = cache_path,tropNom_input=tropNom_input, ambres = ambres,ElMin=ElMin, ElDepWeight=ElDepWeight,
staDb_path = staDb_path,years_list=years_list,num_cores=num_cores,tmp_dir=tmp_dir,project_name=project_name,IGS_logs_dir=IGS_logs_dir,blq_file=blq_file,VMF1_dir = VMF1_dir,
hatanaka=hatanaka,cddis=cddis,
pos_s = 0,
wetz_s = 0.05,
PPPtype = 'static',ionex_type=ionex_type,IONEX_products = IONEX_products,rate = rate,
gnss_products_dir = gnss_products_dir,eterna_path=eterna_path,hardisp_path = hardisp_path,rnx_dir=rnx_dir,tree_options = tree_options_code,tqdm=False,
command='gd2e();kinematic_project.gather_mGNSS()')
qsub_python_code(code,name='{}_static'.format(project_name),email='*****@*****.**',cleanup=False,pbs_base = pbs_base, walltime = '02:00:00')
