for i in range(len(penna_wetz_list)):
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,
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)),
ambres = False
cache_path = '/dev/shm'#'/scratch/bogdanm/cache'
tropNom_input = 'trop'
IONEX_products = '/scratch/bogdanm/Products/IONEX_Products'
rate = 300
eterna_path='/scratch/bogdanm/Products/otl/eterna'
hardisp_path = '/scratch/bogdanm/Products/otl/hardisp/hardisp'
tree_options_code = 'trees_options.rw_otl'
tqdm=False
ElDepWeight = 'SqrtSin'
pbs_base = _os.path.join('/scratch/bogdanm/pbs',project_name) #break down by project folders as gets slow on hpc with multiple files
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)
#generating tree files that won't be overwritten as crc32 will be the same
gen_trees( ionex_type=ionex_type,tmp_dir=tmp_dir,tree_options=tree_options,blq_file=blq_file,mode = 'GPS+GLONASS',ElDepWeight=ElDepWeight,
ElMin = ElMin,pos_s = pos_s,wetz_s = wetz_s,PPPtype = PPPtype,years_list=years_list,cache_path = cache_path,
VMF1_dir = VMF1_dir,project_name = project_name_construct,static_clk = static_clk,ambres = ambres)#the GNSS_class single project name
staDb_path = gen_staDb(tmp_dir = tmp_dir, project_name = project_name, stations_list = stations_list, IGS_logs_dir = IGS_logs_dir)
stations_list_arrays = _np.array_split(stations_list,num_nodes)
for i in range(len(stations_list_arrays)):
code = gen_code(stations_list = list(stations_list_arrays[i]), 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,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='gps.gd2e();kinematic_project.gps.envs(dump=True)')
qsub_python_code(code,name='{}{}{}{}'.format(project_name,str(i),'AR' if ambres else '',ElMin if ElMin!=7 else ''),email='*****@*****.**',cleanup=False,pbs_base = pbs_base,walltime='06:00:00')
