#!/usr/bin/env python3

import os,sys,gc
import glob, errno
import numpy as np
import time
import calendar

# ========== SETTINGS ==========

#dirstat = 'B.Rosatti'
#dirstat = 'China_Cui'
#dirstat = 'CANADA'
#dirstat = 'Siberian_snow_2018'
#stat = 'FPOUT'
dirstat = 'NO0002R'
#stat = 'Kaifeng_2023_0'
#stat = 'ALT_2013'
#stat = 'DRY'
stat = 'OUT_202306_0'
precip = False # If precipitation file from eso is available, then calculate concentration of snow in ng/g instead of deposition in ng/m2
bb = 'CAMS' # GFED or CAMS or FINN

def imp():
    #fdir = '/xnilu_wrk/users/ne/BACKWARD_RUNS/'+dirstat+'/'+stat+ '/OUT'
    #fdir = '/xnilu_wrk/users/ne/BACKWARD_RUNS/'+dirstat+'/'+stat
    #fdir = '/xnilu_wrk/projects/ATMO-ACCESS/OUTDIR/RANDOM/'+dirstat+'/'+stat
    fdir = '/xnilu_wrk/projects/ATMO-ACCESS/OUTDIR/'+dirstat+'/'+stat[4:8]+'/'+stat
    if precip:
        precfile = "/xnilu_wrk/users/ne/BACKWARD_RUNS/Siberian_snow_2018/run2022.log"
        print("We have added a precipitation calculation from ECMWF and concentrations in ng/g will be calculated")
        print("Precipitation was taken from "+precfile)
        return fdir, dirstat, stat, precip, precfile, bb
    else:
        return fdir, dirstat, stat, precip, bb