Commit 37392f0a authored by Kennedy, Joseph H's avatar Kennedy, Joseph H
Browse files

Provide data for dynamics extension MWE

The necessary data files have been added to the LEX data directory for the
dynamics extension to provide a full minimal working example on any machine.
Additionally, the config file and plotting modules have been updated to
reflect this change.

Fixes #18
parent 75f9cb28
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
...@@ -11,13 +11,12 @@ describe = """CISMA_js plot""" ...@@ -11,13 +11,12 @@ describe = """CISMA_js plot"""
# plots velnorm, zooms in on Jakobshavn 69.2N 51.1W # plots velnorm, zooms in on Jakobshavn 69.2N 51.1W
def make_plot(config=None, out_path='.', def make_plot(config, out_path='.'):
cism_path='/lustre/atlas1/cli115/world-shared/4ue/'):
# ---------------- Read Data ----------------------------------- # ---------------- Read Data -----------------------------------
input_file1 = os.path.join(cism_path, 'Gris_init/greenland_1km_2016_12_01.mcb_velnorm.nc') input_file1 = os.path.join(config['cism_data'], 'greenland_1km_2016_12_01.mcb_velnorm.nc')
input_file4 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_velnorm.nc') input_file4 = os.path.join(config['cism_data'], 'albany.out.tstep.ANN_velnorm.nc')
input_file5 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_thk.nc') input_file5 = os.path.join(config['cism_data'], 'albany.out.tstep.ANN_thk.nc')
input_file6 = os.path.join(cism_path, 'CISM-grids/climo_greenland_1km_2017_05_06_notime.nc') input_file6 = os.path.join(config['cism_data'], 'climo_greenland_1km_2017_05_06_notime.nc')
# -------------------------------------------------------------- # --------------------------------------------------------------
img_list = [] img_list = []
...@@ -165,12 +164,7 @@ def make_plot(config=None, out_path='.', ...@@ -165,12 +164,7 @@ def make_plot(config=None, out_path='.',
img_elem = el.image('CISMA_js', img_elem = el.image('CISMA_js',
' '.join(describe.split()), ' '.join(describe.split()),
img_link) img_link)
if config: img_elem['Height'] = config['image_height']
img_elem['Height'] = config['image_height']
img_list.append(img_elem) img_list.append(img_elem)
return img_list return img_list
if __name__ == '__main__':
make_plot()
...@@ -11,14 +11,12 @@ describe = """CISMA_pm plot""" ...@@ -11,14 +11,12 @@ describe = """CISMA_pm plot"""
# Plots velnorm, zooms in on Petermann glacier, 80.5N -59.5 # Plots velnorm, zooms in on Petermann glacier, 80.5N -59.5
def make_plot(config=None, out_path='.', def make_plot(config, out_path='.'):
cism_path='/lustre/atlas1/cli115/world-shared/4ue/'):
# ---------------- Read Data ----------------------------------- # ---------------- Read Data -----------------------------------
input_file1 = os.path.join(cism_path, 'Gris_init/greenland_1km_2016_12_01.mcb_velnorm.nc') input_file1 = os.path.join(config['cism_data'], 'greenland_1km_2016_12_01.mcb_velnorm.nc')
input_file4 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_velnorm.nc') input_file4 = os.path.join(config['cism_data'], 'albany.out.tstep.ANN_velnorm.nc')
input_file5 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_thk.nc') input_file5 = os.path.join(config['cism_data'], 'albany.out.tstep.ANN_thk.nc')
input_file6 = os.path.join(cism_path, 'CISM-grids/climo_greenland_1km_2017_05_06_notime.nc') input_file6 = os.path.join(config['cism_data'], 'climo_greenland_1km_2017_05_06_notime.nc')
# input_file7 = '/lustre/atlas1/cli115/world-shared/4ue/CISM-grids/zurich_mask.nc' # x1 y1 mask for simulated CISMA
# -------------------------------------------------------------- # --------------------------------------------------------------
img_list = [] img_list = []
...@@ -166,12 +164,7 @@ def make_plot(config=None, out_path='.', ...@@ -166,12 +164,7 @@ def make_plot(config=None, out_path='.',
img_elem = el.image('CISMA_pm', img_elem = el.image('CISMA_pm',
' '.join(describe.split()), ' '.join(describe.split()),
img_link) img_link)
if config: img_elem['Height'] = config['image_height']
img_elem['Height'] = config['image_height']
img_list.append(img_elem) img_list.append(img_elem)
return img_list return img_list
if __name__ == '__main__':
make_plot()
...@@ -11,21 +11,19 @@ from livvkit.util import elements as el ...@@ -11,21 +11,19 @@ from livvkit.util import elements as el
describe = """scatter_CISMA_RACMO23 plot.""" describe = """scatter_CISMA_RACMO23 plot."""
def make_plot(config=None, out_path='.', def make_plot(config, out_path='.'):
racmo_path='/lustre/atlas1/cli115/world-shared/4ue/racmo23_GRN_monthly/',
cism_path='/lustre/atlas1/cli115/world-shared/4ue/'):
img_list = []
# -------------------------------------------------------------- # --------------------------------------------------------------
# first use Remap_CISM2RACMO.sh to get remapped file and source data file # first use Remap_CISM2RACMO.sh to get remapped file and source data file
# which were saved in my local directory NCLplot # which were saved in my local directory NCLplot
# ---------------- Read Data ----------------------------------- # ---------------- Read Data -----------------------------------
input_file1 = os.path.join(racmo_path, 'racmo23.ANN_with_latlon_smb.nc') input_file2 = os.path.join(config['cism_data'], 'CISM-Albany.acab.remap2racmo.ANN.nc')
input_file2 = os.path.join(cism_path, 'CISM-Albany/postproc/remapped_cism/CISM-Albany.acab.remap2racmo.ANN.nc') input_file1 = os.path.join(config['racmo_data'], 'racmo23.ANN_with_latlon_smb.nc')
input_file3 = os.path.join(racmo_path, 'RACMO23_masks_ZGRN11.nc') input_file3 = os.path.join(config['racmo_data'], 'RACMO23_masks_ZGRN11.nc')
# -------------------------------------------------------------- # --------------------------------------------------------------
img_list = []
# input_file1 get the following, lat(rlat,rlon),lat=312,lon=306 # input_file1 get the following, lat(rlat,rlon),lat=312,lon=306
# smb(0,:,:), smb(time, lat, lon) with unit"mmWE", have some values= 0.005072316 # smb(0,:,:), smb(time, lat, lon) with unit"mmWE", have some values= 0.005072316
# input_file2 get y1(y1,x1), y1=312,x1=306, here lat_name is y1 # input_file2 get y1(y1,x1), y1=312,x1=306, here lat_name is y1
...@@ -88,12 +86,7 @@ def make_plot(config=None, out_path='.', ...@@ -88,12 +86,7 @@ def make_plot(config=None, out_path='.',
img_elem = el.image('scatter_CISMA_RACMO23', img_elem = el.image('scatter_CISMA_RACMO23',
' '.join(describe.split()), ' '.join(describe.split()),
img_link) img_link)
if config: img_elem['Height'] = config['image_height']
img_elem['Height'] = config['image_height']
img_list.append(img_elem) img_list.append(img_elem)
return img_list return img_list
if __name__ == '__main__':
make_plot()
...@@ -10,21 +10,19 @@ from livvkit.util import elements as el ...@@ -10,21 +10,19 @@ from livvkit.util import elements as el
describe = """CISMA_RACMO23 plot""" describe = """CISMA_RACMO23 plot"""
def make_plot(config=None, out_path='.', def make_plot(config, out_path='.'):
racmo_path='/lustre/atlas1/cli115/world-shared/4ue/racmo23_GRN_monthly/',
cism_path='/lustre/atlas1/cli115/world-shared/4ue/'):
# -------------------------------------------------------------- # --------------------------------------------------------------
# first use remap_CISMA_to_RACMO.sh to get remapped file and source data file # first use remap_CISMA_to_RACMO.sh to get remapped file and source data file
# which were saved in my local directory NCLplot # which were saved in my local directory NCLplot
# ---------------- Read Data ----------------------------------- # ---------------- Read Data -----------------------------------
input_file1 = os.path.join(racmo_path, 'racmo23.ANN_with_latlon_smb.nc') input_file2 = os.path.join(config['cism_data'], 'CISM-Albany.acab.remap2racmo.ANN.nc')
input_file2 = os.path.join(cism_path, 'CISM-Albany/postproc/remapped_cism/CISM-Albany.acab.remap2racmo.ANN.nc') input_file4 = os.path.join(config['cism_data'], 'albany.out.tstep.ANN_acab.nc')
input_file3 = os.path.join(racmo_path, 'RACMO23_masks_ZGRN11.nc') input_file5 = os.path.join(config['cism_data'], 'albany.out.tstep.ANN_thk.nc')
input_file4 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_acab.nc') input_file6 = os.path.join(config['cism_data'], 'climo_greenland_1km_2017_05_06_notime.nc')
input_file5 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_thk.nc') input_file1 = os.path.join(config['racmo_data'], 'racmo23.ANN_with_latlon_smb.nc')
input_file6 = os.path.join(cism_path, 'CISM-grids/climo_greenland_1km_2017_05_06_notime.nc') input_file3 = os.path.join(config['racmo_data'], 'RACMO23_masks_ZGRN11.nc')
# -------------------------------------------------------------- # --------------------------------------------------------------
img_list = [] img_list = []
...@@ -195,12 +193,7 @@ def make_plot(config=None, out_path='.', ...@@ -195,12 +193,7 @@ def make_plot(config=None, out_path='.',
img_elem = el.image('CISMA_RACMO23', img_elem = el.image('CISMA_RACMO23',
' '.join(describe.split()), ' '.join(describe.split()),
img_link) img_link)
if config: img_elem['Height'] = config['image_height']
img_elem['Height'] = config['image_height']
img_list.append(img_elem) img_list.append(img_elem)
return img_list return img_list
if __name__ == '__main__':
make_plot()
...@@ -10,13 +10,10 @@ from livvkit.util import elements as el ...@@ -10,13 +10,10 @@ from livvkit.util import elements as el
describe = """CISMA_temp plot""" describe = """CISMA_temp plot"""
def make_plot(config=None, out_path='.', def make_plot(config, out_path='.'):
cism_path='/lustre/atlas1/cli115/world-shared/4ue/'):
# ---------------- Read Data ----------------------------------- # ---------------- Read Data -----------------------------------
input_file4 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_tempstag.nc') input_file4 = os.path.join(config['cism_data'], 'albany.out.tstep.ANN_tempstag.nc')
# input_file5 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_thk.nc') input_file6 = os.path.join(config['cism_data'], 'climo_greenland_1km_2017_05_06_notime.nc')
input_file6 = os.path.join(cism_path, 'CISM-grids/climo_greenland_1km_2017_05_06_notime.nc')
# input_file7 = '/lustre/atlas1/cli115/world-shared/4ue/CISM-grids/zurich_mask.nc' # x1 y1 mask for simulated CISMA
# -------------------------------------------------------------- # --------------------------------------------------------------
img_list = [] img_list = []
...@@ -156,12 +153,7 @@ def make_plot(config=None, out_path='.', ...@@ -156,12 +153,7 @@ def make_plot(config=None, out_path='.',
img_elem = el.image('CISMA_temp', img_elem = el.image('CISMA_temp',
' '.join(describe.split()), ' '.join(describe.split()),
img_link) img_link)
if config: img_elem['Height'] = config['image_height']
img_elem['Height'] = config['image_height']
img_list.append(img_elem) img_list.append(img_elem)
return img_list return img_list
if __name__ == '__main__':
make_plot()
...@@ -11,14 +11,12 @@ from livvkit.util import elements as el ...@@ -11,14 +11,12 @@ from livvkit.util import elements as el
describe = """CISMA_velnorm plot""" describe = """CISMA_velnorm plot"""
def make_plot(config=None, out_path='.', def make_plot(config, out_path='.'):
cism_path='/lustre/atlas1/cli115/world-shared/4ue/'):
# ---------------- Read Data ----------------------------------- # ---------------- Read Data -----------------------------------
input_file1 = os.path.join(cism_path, 'Gris_init/greenland_1km_2016_12_01.mcb_velnorm.nc') input_file1 = os.path.join(config['cism_data'], 'greenland_1km_2016_12_01.mcb_velnorm.nc')
input_file4 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_velnorm.nc') input_file4 = os.path.join(config['cism_data'], 'albany.out.tstep.ANN_velnorm.nc')
input_file5 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_thk.nc') input_file5 = os.path.join(config['cism_data'], 'albany.out.tstep.ANN_thk.nc')
input_file6 = os.path.join(cism_path, 'CISM-grids/climo_greenland_1km_2017_05_06_notime.nc') input_file6 = os.path.join(config['cism_data'], 'climo_greenland_1km_2017_05_06_notime.nc')
# input_file7 = '/lustre/atlas1/cli115/world-shared/4ue/CISM-grids/zurich_mask.nc' # x1 y1 mask for simulated CISMA
# -------------------------------------------------------------- # --------------------------------------------------------------
img_list = [] img_list = []
...@@ -166,12 +164,7 @@ def make_plot(config=None, out_path='.', ...@@ -166,12 +164,7 @@ def make_plot(config=None, out_path='.',
img_elem = el.image('CISMA_velnorm', img_elem = el.image('CISMA_velnorm',
' '.join(describe.split()), ' '.join(describe.split()),
img_link) img_link)
if config: img_elem['Height'] = config['image_height']
img_elem['Height'] = config['image_height']
img_list.append(img_elem) img_list.append(img_elem)
return img_list return img_list
if __name__ == '__main__':
make_plot()
...@@ -12,14 +12,12 @@ describe = """CISMA_zach plot""" ...@@ -12,14 +12,12 @@ describe = """CISMA_zach plot"""
# zoomed in velnorm over Zachariae glacier at 78N, -30E # zoomed in velnorm over Zachariae glacier at 78N, -30E
def make_plot(config=None, out_path='.', def make_plot(config, out_path='.'):
cism_path='/lustre/atlas1/cli115/world-shared/4ue/'):
# ---------------- Read Data ----------------------------------- # ---------------- Read Data -----------------------------------
input_file1 = os.path.join(cism_path, 'Gris_init/greenland_1km_2016_12_01.mcb_velnorm.nc') input_file1 = os.path.join(config['cism_data'], 'greenland_1km_2016_12_01.mcb_velnorm.nc')
input_file4 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_velnorm.nc') input_file4 = os.path.join(config['cism_data'], 'albany.out.tstep.ANN_velnorm.nc')
input_file5 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_thk.nc') input_file5 = os.path.join(config['cism_data'], 'albany.out.tstep.ANN_thk.nc')
input_file6 = os.path.join(cism_path, 'CISM-grids/climo_greenland_1km_2017_05_06_notime.nc') input_file6 = os.path.join(config['cism_data'], 'climo_greenland_1km_2017_05_06_notime.nc')
# input_file7 = '/lustre/atlas1/cli115/world-shared/4ue/CISM-grids/zurich_mask.nc' # x1 y1 mask for simulated CISMA
# -------------------------------------------------------------- # --------------------------------------------------------------
img_list = [] img_list = []
...@@ -167,12 +165,7 @@ def make_plot(config=None, out_path='.', ...@@ -167,12 +165,7 @@ def make_plot(config=None, out_path='.',
img_elem = el.image('CISMA_zach', img_elem = el.image('CISMA_zach',
' '.join(describe.split()), ' '.join(describe.split()),
img_link) img_link)
if config: img_elem['Height'] = config['image_height']
img_elem['Height'] = config['image_height']
img_list.append(img_elem) img_list.append(img_elem)
return img_list return img_list
if __name__ == '__main__':
make_plot()
...@@ -2,6 +2,8 @@ ...@@ -2,6 +2,8 @@
"dynamics_cisma" : { "dynamics_cisma" : {
"module" : "dynamics/dynamics_cisma.py", "module" : "dynamics/dynamics_cisma.py",
"references" : "dynamics/dynamics_cisma.bib", "references" : "dynamics/dynamics_cisma.bib",
"cism_data": "data/cism/albany",
"racmo_data": "data/racmo",
"image_height" : 300 "image_height" : 300
} }
} }
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