Commit c858279f authored by Christoph Heim's avatar Christoph Heim
Browse files

During extraction of DYAMOND_2

parent 539f3f3f
......@@ -7,7 +7,7 @@ date created: 05.07.2019
date changed: 05.09.2019
usage: arguments:
1st: n jobs for multiprocessing pool
MPAS_3.75 : 3D var: 2 jobs
MPAS_3.75 : 3D var: 2 jobs (sometimes 3 possible)
MPAS_7.50 : 3D var: 5 jobs
"""
###############################################################################
......@@ -131,29 +131,41 @@ if __name__ == '__main__':
# variables to extract
var_namess = {
'3D':['T', 'QV', 'QC', 'W'],
'2D':['LWUTOA', 'T2M', 'U10M', 'V10M', 'SWDSFC', 'TQC', 'PPCONV', 'PPGRID'],
'3D':['QV', 'QC'],
#'2D':['LWUTOA', 'T2M', 'U10M', 'V10M', 'SWDSFC',
# 'TQC', 'PPCONV', 'PPGRID'],
'2D':['U10M', 'V10M', 'SWDSFC',
'TQC', 'PPCONV', 'PPGRID'],
}
run_var_type = '3D'
run_var_type = '2D'
#run_var_type = '2D'
var_names = var_namess[run_var_type]
#var_names = ['T']
# working on
# 7.5 3D and 2D
# 3.75 2D and 3D
#ress = [7.5, 3.75]
ress = [3.75]
#ress = [7.5]
i_bash_output = 0
i_bash_output = 1
# date range
first_date = datetime(2016,8,11)
last_date = datetime(2016,9,9)
# position 1 mistpp1
# position 2 mistpp2
## position 3 mistpp3
## position 4 mistpp4
# position 5 mistpp5 3D
first_date = datetime(2016,8,11)
last_date = datetime(2016,9,9)
# options for computation
options = {}
options['recompute'] = 0
......
......@@ -93,12 +93,11 @@ if __name__ == '__main__':
var_names = ['QV', 'QC', 'T', 'W',
'U10M', 'V10M', 'T2M', 'LWUTOA', 'SWDSFC',
'SLHFLX', 'SSHFLX', 'TQC', 'PPCONV', 'PPGRID']
var_names = ['V10M', 'T2M', 'LWUTOA', 'SWDSFC',
'SLHFLX', 'SSHFLX', 'TQC', 'PPCONV', 'PPGRID']
#var_names = ['SSHFLX', 'TQC', 'PPCONV', 'PPGRID']
# model resolutions [km] of simulations
ress = [9,4]
ress = [4]
ress = [9]
# date range
first_date = datetime(2016,8,11)
......
......@@ -4,7 +4,7 @@
description: Extract lat-lon box of data from model ARPEGE-NH.
author: Christoph Heim
date created: 09.07.2019
date changed: 09.09.2019
date changed: 10.09.2019
usage: arguments:
1st: n jobs for multiprocessing pool
python: 3.5.2
......@@ -163,10 +163,12 @@ def sellatlon_ARPEGE(inp_file, out_file, dt, box, options, var_name, var_dicts,
search = filter_levels(search, box, var_dicts[var_name])
keep_files.extend(search)
search = glob.glob('{}/*.gp'.format(tmp_dir))
for file in search:
if file not in keep_files:
os.remove(file)
## For the moment do not remove files from other variables
## since they may be used later..
#search = glob.glob('{}/*.gp'.format(tmp_dir))
#for file in search:
# if file not in keep_files:
# os.remove(file)
search = glob.glob('{}/*.spectral'.format(tmp_dir))
for file in search:
os.remove(file)
......@@ -252,8 +254,9 @@ if __name__ == '__main__':
# variables to extract
var_namess = {
'3D':['T', 'QV', 'QC', 'W', 'H'],
'2D':['LWUTOA', 'T2M', 'U10M', 'V10M', 'SWDSFC',
'SLHFLX', 'SSHFLX', 'TQC' 'PP'],
#'2D':['LWUTOA', 'T2M', 'U10M', 'V10M', 'SWDSFC',
# 'SLHFLX', 'SSHFLX', 'TQC', 'PP'],
'2D':['SSHFLX', 'TQC', 'PP'],
}
# first variable that is read from grib file
main_vars = {'3D':'T','2D':'LWUTOA'}
......@@ -266,14 +269,22 @@ if __name__ == '__main__':
# model resolutions [km] of simulations
ress = [2.5]
# date range
# complete 3D -
# complete 2D -
# 3D 10 - 19
# 2D -
# 2D: done
first_date = datetime(2016,8,11)
last_date = datetime(2016,9,9)
## 3D: done
#first_date = datetime(2016,8,11)
#last_date = datetime(2016,9,12)
## 2D: laptop postition 3 mistpp4
#first_date = datetime(2016,8,13)
#last_date = datetime(2016,8,25)
## 3D: laptopt postition 4 mistpp5
#first_date = datetime(2016,8,26)
#last_date = datetime(2016,9,9)
# options for computation
options = {}
......@@ -300,7 +311,7 @@ if __name__ == '__main__':
'SLHFLX':{'file':'0.0.10', 'grb_srf':'t1', 'vdim':'2D',},
'SSHFLX':{'file':'0.0.11', 'grb_srf':'t1', 'vdim':'2D',},
'TQC' :{'file':'192.128.78', 'grb_srf':'t1', 'vdim':'2D',},
'PP' :{'file':'8.1.0', 'grb_srf':'t1', 'vdim':'2D',},
'PP' :{'file':'0.1.8', 'grb_srf':'t1', 'vdim':'2D',},
}
inc_min = {'3D':180, '2D':15}
remap_res = 2.5
......
2.5km_W_201608101500
2.5km_QV_201608110000
2.5km_T_201608110000
2.5km_QV_201608111200
2.5km_QV_201608111800
2.5km_QV_201608110600
2.5km_QV_201608110300
2.5km_QV_201608120000
2.5km_QV_201608112100
2.5km_QV_201608110000
2.5km_QV_201608111500
2.5km_QV_201608110600
2.5km_QV_201608110900
2.5km_QV_201608111800
2.5km_QV_201608112100
2.5km_QV_201608110300
2.5km_QV_201608111200
2.5km_QC_201608110000
2.5km_W_201608110000
2.5km_H_201608110000
2.5km_T2M_201608200000
2.5km_T2M_201608201200
2.5km_U10M_201608200000
2.5km_U10M_201608201200
2.5km_V10M_201608200000
2.5km_V10M_201608201200
2.5km_SWDSFC_201608200000
2.5km_SWDSFC_201608201200
2.5km_SWDSFC_201608200300
2.5km_SWDSFC_201608201500
2.5km_SLHFLX_201608200000
2.5km_SLHFLX_201608201200
2.5km_SLHFLX_201608200300
2.5km_SLHFLX_201608201500
2.5km_SSHFLX_201608200000
2.5km_SSHFLX_201608201200
2.5km_SSHFLX_201608201500
2.5km_SSHFLX_201608200300
2.5km_TQC_201608200000
2.5km_TQC_201608201200
2.5km_TQC_201608201500
2.5km_TQC_201608200300
2.5km_QV_201608120000
......@@ -19,8 +19,8 @@ models=(GEOS_3 ICON_10 ICON_2.5 IFS_4 IFS_9 MPAS_3.75 MPAS_7.5 \
NICAM_3.5 NICAM_7 SAM_4 UM_5 ARPEGE-NH_2.5 FV3_3.25)
models=(NICAM_7 NICAM_3.5)
#models=(SAM_4)
#models=(ICON_10 ICON_2.5)
models=(SAM_4)
models=(ICON_10 ICON_2.5)
#models=(UM_5)
#models=(MPAS_7.5 MPAS_3.75)
#models=(IFS_4 IFS_9)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment