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