import os, sys
import ncClasses.ncField as ncField
from ncClasses.subdomains import setSSI
import matplotlib.pyplot as plt
from datetime import datetime, timedelta
import numpy as np
import multiprocessing as mp
os.chdir('00_newScripts/')
models = ['RAW1', 'SM1']
i_subdomain = 1
#var = 'FQV'
#var = 'FQVx'
var = 'FQVy'
ssI, domainName = setSSI(i_subdomain, {'4.4': {}, '2.2': {}, '1.1': {}})
dt0 = datetime(2006, 7, 11, 0)
#dt1 = datetime(2006,7,20,1)
dt1 = datetime(2006, 7, 13, 0)
dts = np.arange(dt0, dt1, timedelta(hours=1))
inpPath = '../01_rawData/topocut/'
if len(sys.argv) > 1:
njobs = int(sys.argv[1])
print('number of jobs is ' + str(njobs))
else:
print('Number of Jobs not given. Assume 1')
njobs = 1
import ncClasses.analysis as analysis
from datetime import datetime
from functions import *
from ncClasses.subdomains import setSSI
####################### NAMELIST INPUTS FILES #######################
# directory of input model folders
#inpPath = '../02_fields/subDomDiur'
inpPath = '../02_fields/diurnal'
#inpPath = '../02_fields/topocut'
fieldNames = ['nTOT_PREC', 'cHSURF']
#####################################################################
####################### NAMELIST DIMENSIONS #######################
i_subDomain = 1 # 0: full domain, 1: alpine region
ssI, domainName = setSSI(i_subDomain, {'4': {}, '2': {}, '1': {}})
ssI['diurnal'] = [20, 21, 22, 23, 0, 1, 2, 3, 4, 5, 6, 7] # list values
ssI['diurnal'] = [8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19] # list values
#ssI['diurnal'] = [12,13,14,15,16,17] # list values
#ssI['diurnal'] = [18,19,20,21,22,23] # list values
#ssI['diurnal'] = [0 ,1 ,2 ,3 ,4 ,5 ] # list values
#ssI['diurnal'] = [6 ,7 ,8 ,9 ,10,11] # list values
####################### NAMELIST AGGREGATE #######################
# Options: MEAN, SUM
ag_commnds = {}
#ag_commnds['rlat'] = 'MEAN'
#ag_commnds['rlon'] = 'MEAN'
#ag_commnds['time'] = 'SUM'
import matplotlib
if i_plot == 2:
matplotlib.use('Agg')
import matplotlib.pyplot as plt
import ncClasses.analysis as analysis
from ncClasses.subdomains import setSSI
from datetime import datetime
from functions import *
####################### NAMELIST INPUTS FILES #######################
# directory of input model folders
inpPath = '../02_fields/diurnal'
fieldNames = ['cHSURF']
i_subdomain = 1
#####################################################################
ssI = setSSI(i_subdomain, {'4.4': {}, '2.2': {}, '1.1': {}})
# DOMAIN 4 Greater Northern Italy Plains (use subDomDiur)
#lon0 = [45]
#lon1 = [118]
#lat0 = [12]
#lat1 = [62]
#plotName = 'locations_domain4_Greater_Northern_Italy_Plains.png'
## DOMAIN 3 Northern Italy Plains (use subDomDiur)
#lon0 = [48]
#lon1 = [100]
#lat0 = [25]
#lat1 = [56]
#plotName = 'locations_domain3_Northern_Italy_Plains.png'
###################### NAMELIST DIMENSIONS #######################
# SUBSPACE
subSpaceIndsIN = {}
import ncClasses.analysis as analysis
from datetime import datetime
from functions import *
from ncClasses.subdomains import setSSI
####################### NAMELIST INPUTS FILES #######################
# directory of input model folders
#inpPath = '../02_fields/topocut'
inpPath = '../02_fields/diurnal'
#fieldNames = ['zQC', 'nHPBL', 'cHSURF']
fieldNames = ['cHSURF', 'nTOT_PREC']
#####################################################################
####################### NAMELIST DIMENSIONS #######################
i_subDomain = 0 # 0: full domain, 1: alpine region
ssI, domainName = setSSI(i_subDomain, {'4.4': {}, '2.2': {}, '1.1': {}})
ssI_AR, domainName = setSSI(1, {'4.4': {}, '2.2': {}, '1.1': {}})
#####################################################################
####################### NAMELIST AGGREGATE #######################
# Options: MEAN, SUM, DIURNAL
ag_commnds = {}
#ag_commnds['rlat'] = 'MEAN'
#ag_commnds['rlon'] = 'MEAN'
ag_commnds['time'] = 'MEAN'
#ag_commnds['altitude'] = 'MEAN'
#####################################################################
####################### NAMELIST PLOT #######################
nDPlot = 2 # How many dimensions should plot have (1 or 2)
plotOutDir = '../00_plots'
import os
import ncClasses.ncField as ncField
from ncClasses.subdomains import setSSI
import matplotlib.pyplot as plt
import numpy as np
os.chdir('00_newScripts/')
i_resolutions = 3
i_subdomain = 0
ssI, domainName = setSSI(i_subdomain, {})
inpPath = '../01_rawData/topocut/'
#inpPath = '../01_rawData/'
res = '4.4'
mode = ''
ncFileName = 'lffd2006071515z.nc'
#ncFileName = 'constantParams.nc'
#inpPath = '../02_fields/topocut/'
#res = '4.4'
#mode = ''
#ncFileName = 'zQV.nc'
srcNCPath = inpPath + res + mode + '/' + ncFileName
print(srcNCPath)
qc = ncField.ncField('QC', srcNCPath, ssI)
#qc.saveToNewNC('../uTest.nc')
qc.addVarToExistingNC('../uTest.nc')
#plt.contourf(qc.vals[0,10,:,:].squeeze())
#plt.colorbar()
i_info = 2 # output some information [from 0 (off) to 5 (all you can read)]
import matplotlib
if i_plot == 2:
matplotlib.use('Agg')
import matplotlib.pyplot as plt
import ncClasses.analysis as analysis
from ncClasses.subdomains import setSSI
from datetime import datetime
from functions import *
####################### NAMELIST INPUTS FILES #######################
# directory of input model folders
inpPath = '../02_fields/topocut'
fieldNames = ['cHSURF']
#####################################################################
ssI_an, domainName = setSSI(0, {'4.4': {}, '2.2': {}, '1.1': {}})
i_subDomain = 5 # 0: full domain, 1: alpine region
ssI, domainName = setSSI(i_subDomain, {'4.4': {}, '2.2': {}, '1.1': {}})
#startHght = 10
#endHght = 10
#altInds = list(range(startHght,endHght+1))
#ssI['altitude'] = altInds
#startTime = datetime(2006,7,11,00)
#endTime = datetime(2006,7,11,23)
#ssI['time'] = [startTime,endTime] # border values (one value if only one time step desired)
#ssI['diurnal'] = [10] # list values
#####################################################################
####################### NAMELIST AGGREGATE #######################
# Options: MEAN, SUM, DIURNAL
