From fec24f551219a0808181ce042a4578138b45e64d Mon Sep 17 00:00:00 2001 From: chengzhuzhang Date: Sat, 12 Oct 2024 17:02:37 -0500 Subject: [PATCH 1/3] add H2OLNZ --- ...onal_mean_2d_stratosphere_model_vs_obs.cfg | 68 +++++++++++++++++++ 1 file changed, 68 insertions(+) diff --git a/e3sm_diags/driver/default_diags/zonal_mean_2d_stratosphere_model_vs_obs.cfg b/e3sm_diags/driver/default_diags/zonal_mean_2d_stratosphere_model_vs_obs.cfg index 674d80911..2e1113b3d 100644 --- a/e3sm_diags/driver/default_diags/zonal_mean_2d_stratosphere_model_vs_obs.cfg +++ b/e3sm_diags/driver/default_diags/zonal_mean_2d_stratosphere_model_vs_obs.cfg @@ -40,6 +40,40 @@ seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "1 contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] diff_levels = [-3, -2,-1.5,-1,-0.5,-0.1,0.1,0.5,1,1.5,2,3] +[#] +sets = ["zonal_mean_2d_stratosphere"] +case_id = "MERRA2" +variables = ["H2OLNZ"] +ref_name = "MERRA2" +reference_name = "MERRA2 Reanalysis" +seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"] +contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] +diff_levels = [-1.5, -1, -0.5, -0.2, -0.1,-0.05,-0.01, 0.01,0.05, 0.1, 0.2, 0.5, 1, 1.5] + +[#] +sets = ["zonal_mean_2d_stratosphere"] +case_id = "MERRA2 (relative difference)" +variables = ["Q"] +ref_name = "MERRA2" +reference_name = "MERRA2 Reanalysis" +seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"] +diff_type = "relative" +diff_title = "(Model - Obs.)/Obs. * 100" +contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] +diff_levels = [-100, -75, -50, -30, -15, -5, 5, 15, 30, 50, 75, 100] + +[#] +sets = ["zonal_mean_2d_stratosphere"] +case_id = "MERRA2 (relative difference)" +variables = ["H2OLNZ"] +ref_name = "MERRA2" +reference_name = "MERRA2 Reanalysis" +seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"] +diff_type = "relative" +diff_title = "(Model - Obs.)/Obs. * 100" +contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] +diff_levels = [-100, -75, -50, -30, -15, -5, 5, 15, 30, 50, 75, 100] + [#] sets = ["zonal_mean_2d_stratosphere"] case_id = "MERRA2" @@ -105,3 +139,37 @@ reference_name = "ERA5 Reanalysis" seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"] contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] diff_levels = [-3, -2,-1.5,-1,-0.5,-0.1,0.1,0.5,1,1.5,2,3] + +[#] +sets = ["zonal_mean_2d_stratosphere"] +case_id = "ERA5" +variables = ["H2OLNZ"] +ref_name = "ERA5" +reference_name = "ERA5 Reanalysis" +seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"] +contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] +diff_levels = [-1.5, -1, -0.5, -0.2, -0.1, -0.05,-0.01, 0.01,0.05, 0.1, 0.2, 0.5, 1, 1.5] + +[#] +sets = ["zonal_mean_2d_stratosphere"] +case_id = "ERA5 (relative difference)" +variables = ["Q"] +ref_name = "ERA5" +reference_name = "ERA5 Reanalysis" +seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"] +diff_type = "relative" +diff_title = "(Model - Obs.)/Obs. * 100" +contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] +diff_levels = [-100, -75, -50, -30, -15, -5, 5, 15, 30, 50, 75, 100] + +[#] +sets = ["zonal_mean_2d_stratosphere"] +case_id = "ERA5 (relative difference)" +variables = ["H2OLNZ"] +ref_name = "ERA5" +reference_name = "ERA5 Reanalysis" +seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"] +diff_type = "relative" +diff_title = "(Model - Obs.)/Obs. * 100" +contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] +diff_levels = [-100, -75, -50, -30, -15, -5, 5, 15, 30, 50, 75, 100] From ac352711ffa51ead6bb0da7f8fc5ed3ea131b135 Mon Sep 17 00:00:00 2001 From: chengzhuzhang Date: Mon, 14 Oct 2024 11:58:57 -0500 Subject: [PATCH 2/3] refine contour levels --- .../zonal_mean_2d_stratosphere_model_vs_obs.cfg | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/e3sm_diags/driver/default_diags/zonal_mean_2d_stratosphere_model_vs_obs.cfg b/e3sm_diags/driver/default_diags/zonal_mean_2d_stratosphere_model_vs_obs.cfg index 2e1113b3d..2346bec34 100644 --- a/e3sm_diags/driver/default_diags/zonal_mean_2d_stratosphere_model_vs_obs.cfg +++ b/e3sm_diags/driver/default_diags/zonal_mean_2d_stratosphere_model_vs_obs.cfg @@ -47,8 +47,8 @@ variables = ["H2OLNZ"] ref_name = "MERRA2" reference_name = "MERRA2 Reanalysis" seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"] -contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] -diff_levels = [-1.5, -1, -0.5, -0.2, -0.1,-0.05,-0.01, 0.01,0.05, 0.1, 0.2, 0.5, 1, 1.5] +contour_levels = [0.001, 0.0015, 0.002,0.0025,0.003,0.0035,0.004,0.0045,0.005] +diff_levels = [-0.001, -0.00075, -0.00050, -0.00025, -0.00005,0.00005,0.00025, 0.0005,0.00075,0.001] [#] sets = ["zonal_mean_2d_stratosphere"] @@ -71,7 +71,7 @@ reference_name = "MERRA2 Reanalysis" seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"] diff_type = "relative" diff_title = "(Model - Obs.)/Obs. * 100" -contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] +contour_levels = [0.001, 0.0015, 0.002,0.0025,0.003,0.0035,0.004,0.0045,0.005] diff_levels = [-100, -75, -50, -30, -15, -5, 5, 15, 30, 50, 75, 100] [#] @@ -147,8 +147,8 @@ variables = ["H2OLNZ"] ref_name = "ERA5" reference_name = "ERA5 Reanalysis" seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"] -contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] -diff_levels = [-1.5, -1, -0.5, -0.2, -0.1, -0.05,-0.01, 0.01,0.05, 0.1, 0.2, 0.5, 1, 1.5] +contour_levels = [0.001, 0.0015, 0.002,0.0025,0.003,0.0035,0.004,0.0045,0.005] +diff_levels = [-0.001, -0.00075, -0.00050, -0.00025, -0.00005,0.00005,0.00025, 0.0005,0.00075,0.001] [#] sets = ["zonal_mean_2d_stratosphere"] @@ -171,5 +171,5 @@ reference_name = "ERA5 Reanalysis" seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"] diff_type = "relative" diff_title = "(Model - Obs.)/Obs. * 100" -contour_levels = [0.6,0.8, 1, 1.2, 1.6, 2, 2.5, 3, 3.5, 4] +contour_levels = [0.001, 0.0015, 0.002,0.0025,0.003,0.0035,0.004,0.0045,0.005] diff_levels = [-100, -75, -50, -30, -15, -5, 5, 15, 30, 50, 75, 100] From df18f777edf640c7931fff4a920e1d164b85a77f Mon Sep 17 00:00:00 2001 From: chengzhuzhang Date: Mon, 14 Oct 2024 12:16:40 -0500 Subject: [PATCH 3/3] fix UnboundLocalError in testing --- e3sm_diags/plot/cartopy/arm_diags_plot.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/e3sm_diags/plot/cartopy/arm_diags_plot.py b/e3sm_diags/plot/cartopy/arm_diags_plot.py index cf485dd06..dc77b31b8 100644 --- a/e3sm_diags/plot/cartopy/arm_diags_plot.py +++ b/e3sm_diags/plot/cartopy/arm_diags_plot.py @@ -149,7 +149,7 @@ def plot_convection_onset_statistics( cwv_min = 28 bin_width = 2.0 sitename = "Darwin" - if region == "sgp": # sgp + if region == "sgpc1": # sgp cwv_max = 75 cwv_min = 20 bin_width = 2.0