From 3cfcf343434e285c28a171039b727a85228f4286 Mon Sep 17 00:00:00 2001 From: Do Ngoc Khanh <14163153+ngockhanhdo1997@users.noreply.github.com> Date: Wed, 24 Jan 2024 14:48:34 +0900 Subject: [PATCH 1/5] Add variables needed for slucm with 2d urban parameters --- geogrid/GEOGRID.TBL.ARW | 78 +++++++++++++++++++++++++++++++++++++++++ 1 file changed, 78 insertions(+) diff --git a/geogrid/GEOGRID.TBL.ARW b/geogrid/GEOGRID.TBL.ARW index a8a2681a..dd0dea61 100644 --- a/geogrid/GEOGRID.TBL.ARW +++ b/geogrid/GEOGRID.TBL.ARW @@ -39,6 +39,8 @@ name=LANDUSEF landmask_water = usgs_30s:16 # Calculate a landmask from this field landmask_water = usgs_lakes:16,28 # Calculate a landmask from this field landmask_water = modis_lakes:17,21 # Calculate a landmask from this field + landmask_water = modis_2010:17 # Calculate a landmask from this field + landmask_water = modis_2050:17 # Calculate a landmask from this field landmask_water = usgs_2m:16 # Calculate a landmask from this field landmask_water = usgs_5m:16 # Calculate a landmask from this field landmask_water = usgs_10m:16 # Calculate a landmask from this field @@ -57,6 +59,8 @@ name=LANDUSEF interp_option = usgs_30s:nearest_neighbor interp_option = usgs_lakes:nearest_neighbor interp_option = modis_lakes:nearest_neighbor + interp_option = modis_2010:nearest_neighbor + interp_option = modis_2050:nearest_neighbor interp_option = usgs_2m:four_pt interp_option = usgs_5m:four_pt interp_option = usgs_10m:four_pt @@ -75,6 +79,8 @@ name=LANDUSEF rel_path = usgs_30s:landuse_30s/ rel_path = usgs_lakes:landuse_30s_with_lakes/ rel_path = modis_lakes:modis_landuse_21class_30s/ + rel_path = modis_2010:modis_landuse_20class_15s_2010s/ + rel_path = modis_2050:modis_landuse_20class_15s_2050s/ rel_path = usgs_2m:landuse_2m/ rel_path = usgs_5m:landuse_5m/ rel_path = usgs_10m:landuse_10m/ @@ -966,3 +972,75 @@ name = BATHYMETRY rel_path = default:topobath_30s/ flag_in_output=FLAG_BATHYMETRY =============================== +name = MH_URB2D + priority=1 + optional=yes + dest_type=continuous + fill_missing = 0. + interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt + interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt + rel_path = y2010_urb:urban_params/2010/H_avg/ + rel_path = y2050_urb:urban_params/2050/H_avg/ +=============================== +name = ZD_URB2D + priority=1 + optional=yes + dest_type=continuous + fill_missing = 0. + interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt + interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt + rel_path = y2010_urb:urban_params/2010/d/ + rel_path = y2050_urb:urban_params/2050/d/ +=============================== +name = Z0_URB2D + priority=1 + optional=yes + dest_type=continuous + fill_missing = 0. + interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt + interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt + rel_path = y2010_urb:urban_params/2010/z_0/ + rel_path = y2050_urb:urban_params/2050/z_0/ +=============================== +name = BUILD_AREA_FRACTION + priority=1 + optional=yes + dest_type=continuous + fill_missing = 0. + interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt + interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt + rel_path = y2010_urb:urban_params/2010/lambda_p/ + rel_path = y2050_urb:urban_params/2050/lambda_p/ +=============================== +name = LF_URB2D_S + priority=1 + optional=yes + dest_type=continuous + fill_missing = 0. + interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt + interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt + rel_path = y2010_urb:urban_params/2010/lambda_f/ + rel_path = y2050_urb:urban_params/2050/lambda_f/ +=============================== +name = AHE + priority=1 + optional=yes + dest_type=continuous + fill_missing = 0. + interp_option = y2010_ahe:average_gcell(4.0)+four_pt+average_4pt + interp_option = y2050_ahe:average_gcell(4.0)+four_pt+average_4pt + rel_path = y2010_ahe:AHE_2010s/2_bytes/ + rel_path = y2050_ahe:AHE_2050s/2_bytes/ + z_dim_name = month_hour +=============================== +name = AHE + priority=2 + optional=yes + dest_type=continuous + fill_missing = 0. + interp_option = y2010_ahe:average_gcell(4.0)+four_pt+average_4pt + interp_option = y2050_ahe:average_gcell(4.0)+four_pt+average_4pt + rel_path = y2010_ahe:AHE_2010s/3_bytes/ + rel_path = y2050_ahe:AHE_2050s/3_bytes/ + z_dim_name = month_hour +=============================== From ac0b1886641680d0e263ab6b841cdc6c1b88b344 Mon Sep 17 00:00:00 2001 From: Do Ngoc Khanh <14163153+ngockhanhdo1997@users.noreply.github.com> Date: Wed, 14 Feb 2024 11:17:34 +0900 Subject: [PATCH 2/5] Add lakes --- geogrid/GEOGRID.TBL.ARW | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/geogrid/GEOGRID.TBL.ARW b/geogrid/GEOGRID.TBL.ARW index dd0dea61..ad52b3a2 100644 --- a/geogrid/GEOGRID.TBL.ARW +++ b/geogrid/GEOGRID.TBL.ARW @@ -40,7 +40,9 @@ name=LANDUSEF landmask_water = usgs_lakes:16,28 # Calculate a landmask from this field landmask_water = modis_lakes:17,21 # Calculate a landmask from this field landmask_water = modis_2010:17 # Calculate a landmask from this field + landmask_water = modis_2010_lake:17,21 # Calculate a landmask from this field landmask_water = modis_2050:17 # Calculate a landmask from this field + landmask_water = modis_2050_lake:17,21 # Calculate a landmask from this field landmask_water = usgs_2m:16 # Calculate a landmask from this field landmask_water = usgs_5m:16 # Calculate a landmask from this field landmask_water = usgs_10m:16 # Calculate a landmask from this field @@ -60,7 +62,9 @@ name=LANDUSEF interp_option = usgs_lakes:nearest_neighbor interp_option = modis_lakes:nearest_neighbor interp_option = modis_2010:nearest_neighbor + interp_option = modis_2010_lake:nearest_neighbor interp_option = modis_2050:nearest_neighbor + interp_option = modis_2050_lake:nearest_neighbor interp_option = usgs_2m:four_pt interp_option = usgs_5m:four_pt interp_option = usgs_10m:four_pt @@ -80,7 +84,9 @@ name=LANDUSEF rel_path = usgs_lakes:landuse_30s_with_lakes/ rel_path = modis_lakes:modis_landuse_21class_30s/ rel_path = modis_2010:modis_landuse_20class_15s_2010s/ + rel_path = modis_2010_lake:modis_landuse_20class_15s_2010s_with_lakes/ rel_path = modis_2050:modis_landuse_20class_15s_2050s/ + rel_path = modis_2050_lake:modis_landuse_20class_15s_2050s_with_lakes/ rel_path = usgs_2m:landuse_2m/ rel_path = usgs_5m:landuse_5m/ rel_path = usgs_10m:landuse_10m/ From 429d7da8ee226b39bbc4fb1aeea9b3f3c35e370d Mon Sep 17 00:00:00 2001 From: Do Ngoc Khanh <14163153+ngockhanhdo1997@users.noreply.github.com> Date: Fri, 16 Feb 2024 13:46:45 +0900 Subject: [PATCH 3/5] Add CGLS vegetation --- geogrid/GEOGRID.TBL.ARW | 8 ++++++++ 1 file changed, 8 insertions(+) diff --git a/geogrid/GEOGRID.TBL.ARW b/geogrid/GEOGRID.TBL.ARW index ad52b3a2..900c2016 100644 --- a/geogrid/GEOGRID.TBL.ARW +++ b/geogrid/GEOGRID.TBL.ARW @@ -161,6 +161,8 @@ name=ALBEDO12M name=GREENFRAC priority=1 dest_type=continuous + interp_option=cgls_veg:average_gcell(4.0)+four_pt+average_4pt+average_16pt + interp_option=modis_veg_no_search:average_gcell(4.0)+four_pt+average_4pt+average_16pt interp_option=modis_fpar:average_gcell(4.0)+four_pt+average_4pt+average_16pt+search interp_option=nesdis_greenfrac:four_pt+average_4pt+average_16pt+search interp_option= lowres:average_gcell(4.0)+four_pt+wt_average_4pt+wt_average_16pt+search @@ -168,6 +170,8 @@ name=GREENFRAC z_dim_name=month masked = water fill_missing = 0. + rel_path=cgls_veg:greenfrac_fcover_cgls/ + rel_path=modis_veg_no_search:greenfrac_fpar_modis/ rel_path=modis_fpar:greenfrac_fpar_modis/ rel_path=nesdis_greenfrac:greenfrac/ rel_path= lowres:greenfrac_fpar_modis_5m/ @@ -176,12 +180,16 @@ name=GREENFRAC name=LAI12M priority=1 dest_type=continuous + interp_option=cgls_veg:average_gcell(4.0)+four_pt+average_4pt+average_16pt + interp_option=modis_veg_no_search:average_gcell(4.0)+four_pt+average_4pt+average_16pt interp_option=modis_lai:average_gcell(4.0)+four_pt+average_4pt+average_16pt+search interp_option= lowres:average_gcell(4.0)+four_pt+wt_average_4pt+wt_average_16pt+search interp_option=default:average_gcell(4.0)+four_pt+average_4pt+average_16pt+search z_dim_name=month masked = water fill_missing = 0. + rel_path=cgls_veg:lai_cgls/ + rel_path=modis_veg_no_search:lai_modis_30s/ rel_path=modis_lai:lai_modis_30s/ rel_path= lowres:lai_modis_10m/ rel_path=default:lai_modis_10m/ From f26ae0cf6ef21ee90c1e08a4c59f106377e60538 Mon Sep 17 00:00:00 2001 From: Do Ngoc Khanh <14163153+epn09@users.noreply.github.com> Date: Tue, 7 May 2024 09:58:04 +0900 Subject: [PATCH 4/5] Move to slucm_kvnk subdirectory --- geogrid/GEOGRID.TBL.ARW | 40 ++++++++++++++++++++-------------------- 1 file changed, 20 insertions(+), 20 deletions(-) diff --git a/geogrid/GEOGRID.TBL.ARW b/geogrid/GEOGRID.TBL.ARW index 900c2016..0503365d 100644 --- a/geogrid/GEOGRID.TBL.ARW +++ b/geogrid/GEOGRID.TBL.ARW @@ -83,10 +83,10 @@ name=LANDUSEF rel_path = usgs_30s:landuse_30s/ rel_path = usgs_lakes:landuse_30s_with_lakes/ rel_path = modis_lakes:modis_landuse_21class_30s/ - rel_path = modis_2010:modis_landuse_20class_15s_2010s/ - rel_path = modis_2010_lake:modis_landuse_20class_15s_2010s_with_lakes/ - rel_path = modis_2050:modis_landuse_20class_15s_2050s/ - rel_path = modis_2050_lake:modis_landuse_20class_15s_2050s_with_lakes/ + rel_path = modis_2010:slucm_kvnk/modis_landuse_20class_15s_2010s/ + rel_path = modis_2010_lake:slucm_kvnk/modis_landuse_20class_15s_2010s_with_lakes/ + rel_path = modis_2050:slucm_kvnk/modis_landuse_20class_15s_2050s/ + rel_path = modis_2050_lake:slucm_kvnk/modis_landuse_20class_15s_2050s_with_lakes/ rel_path = usgs_2m:landuse_2m/ rel_path = usgs_5m:landuse_5m/ rel_path = usgs_10m:landuse_10m/ @@ -170,7 +170,7 @@ name=GREENFRAC z_dim_name=month masked = water fill_missing = 0. - rel_path=cgls_veg:greenfrac_fcover_cgls/ + rel_path=cgls_veg:slucm_kvnk/greenfrac_fcover_cgls/ rel_path=modis_veg_no_search:greenfrac_fpar_modis/ rel_path=modis_fpar:greenfrac_fpar_modis/ rel_path=nesdis_greenfrac:greenfrac/ @@ -188,7 +188,7 @@ name=LAI12M z_dim_name=month masked = water fill_missing = 0. - rel_path=cgls_veg:lai_cgls/ + rel_path=cgls_veg:slucm_kvnk/lai_cgls/ rel_path=modis_veg_no_search:lai_modis_30s/ rel_path=modis_lai:lai_modis_30s/ rel_path= lowres:lai_modis_10m/ @@ -993,8 +993,8 @@ name = MH_URB2D fill_missing = 0. interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_urb:urban_params/2010/H_avg/ - rel_path = y2050_urb:urban_params/2050/H_avg/ + rel_path = y2010_urb:slucm_kvnk/urban_params/2010/H_avg/ + rel_path = y2050_urb:slucm_kvnk/urban_params/2050/H_avg/ =============================== name = ZD_URB2D priority=1 @@ -1003,8 +1003,8 @@ name = ZD_URB2D fill_missing = 0. interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_urb:urban_params/2010/d/ - rel_path = y2050_urb:urban_params/2050/d/ + rel_path = y2010_urb:slucm_kvnk/urban_params/2010/d/ + rel_path = y2050_urb:slucm_kvnk/urban_params/2050/d/ =============================== name = Z0_URB2D priority=1 @@ -1013,8 +1013,8 @@ name = Z0_URB2D fill_missing = 0. interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_urb:urban_params/2010/z_0/ - rel_path = y2050_urb:urban_params/2050/z_0/ + rel_path = y2010_urb:slucm_kvnk/urban_params/2010/z_0/ + rel_path = y2050_urb:slucm_kvnk/urban_params/2050/z_0/ =============================== name = BUILD_AREA_FRACTION priority=1 @@ -1023,8 +1023,8 @@ name = BUILD_AREA_FRACTION fill_missing = 0. interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_urb:urban_params/2010/lambda_p/ - rel_path = y2050_urb:urban_params/2050/lambda_p/ + rel_path = y2010_urb:slucm_kvnk/urban_params/2010/lambda_p/ + rel_path = y2050_urb:slucm_kvnk/urban_params/2050/lambda_p/ =============================== name = LF_URB2D_S priority=1 @@ -1033,8 +1033,8 @@ name = LF_URB2D_S fill_missing = 0. interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_urb:urban_params/2010/lambda_f/ - rel_path = y2050_urb:urban_params/2050/lambda_f/ + rel_path = y2010_urb:slucm_kvnk/urban_params/2010/lambda_f/ + rel_path = y2050_urb:slucm_kvnk/urban_params/2050/lambda_f/ =============================== name = AHE priority=1 @@ -1043,8 +1043,8 @@ name = AHE fill_missing = 0. interp_option = y2010_ahe:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_ahe:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_ahe:AHE_2010s/2_bytes/ - rel_path = y2050_ahe:AHE_2050s/2_bytes/ + rel_path = y2010_ahe:slucm_kvnk/AHE_2010s/2_bytes/ + rel_path = y2050_ahe:slucm_kvnk/AHE_2050s/2_bytes/ z_dim_name = month_hour =============================== name = AHE @@ -1054,7 +1054,7 @@ name = AHE fill_missing = 0. interp_option = y2010_ahe:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_ahe:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_ahe:AHE_2010s/3_bytes/ - rel_path = y2050_ahe:AHE_2050s/3_bytes/ + rel_path = y2010_ahe:slucm_kvnk/AHE_2010s/3_bytes/ + rel_path = y2050_ahe:slucm_kvnk/AHE_2050s/3_bytes/ z_dim_name = month_hour =============================== From 72dafba605b6839ce137aba200705638f78049b2 Mon Sep 17 00:00:00 2001 From: Do Ngoc Khanh <14163153+epn09@users.noreply.github.com> Date: Thu, 9 May 2024 09:36:54 +0900 Subject: [PATCH 5/5] Rename directory --- geogrid/GEOGRID.TBL.ARW | 40 ++++++++++++++++++++-------------------- 1 file changed, 20 insertions(+), 20 deletions(-) diff --git a/geogrid/GEOGRID.TBL.ARW b/geogrid/GEOGRID.TBL.ARW index 0503365d..b3f4b83d 100644 --- a/geogrid/GEOGRID.TBL.ARW +++ b/geogrid/GEOGRID.TBL.ARW @@ -83,10 +83,10 @@ name=LANDUSEF rel_path = usgs_30s:landuse_30s/ rel_path = usgs_lakes:landuse_30s_with_lakes/ rel_path = modis_lakes:modis_landuse_21class_30s/ - rel_path = modis_2010:slucm_kvnk/modis_landuse_20class_15s_2010s/ - rel_path = modis_2010_lake:slucm_kvnk/modis_landuse_20class_15s_2010s_with_lakes/ - rel_path = modis_2050:slucm_kvnk/modis_landuse_20class_15s_2050s/ - rel_path = modis_2050_lake:slucm_kvnk/modis_landuse_20class_15s_2050s_with_lakes/ + rel_path = modis_2010:slucm/modis_landuse_20class_15s_2010s/ + rel_path = modis_2010_lake:slucm/modis_landuse_20class_15s_2010s_with_lakes/ + rel_path = modis_2050:slucm/modis_landuse_20class_15s_2050s/ + rel_path = modis_2050_lake:slucm/modis_landuse_20class_15s_2050s_with_lakes/ rel_path = usgs_2m:landuse_2m/ rel_path = usgs_5m:landuse_5m/ rel_path = usgs_10m:landuse_10m/ @@ -170,7 +170,7 @@ name=GREENFRAC z_dim_name=month masked = water fill_missing = 0. - rel_path=cgls_veg:slucm_kvnk/greenfrac_fcover_cgls/ + rel_path=cgls_veg:slucm/greenfrac_fcover_cgls/ rel_path=modis_veg_no_search:greenfrac_fpar_modis/ rel_path=modis_fpar:greenfrac_fpar_modis/ rel_path=nesdis_greenfrac:greenfrac/ @@ -188,7 +188,7 @@ name=LAI12M z_dim_name=month masked = water fill_missing = 0. - rel_path=cgls_veg:slucm_kvnk/lai_cgls/ + rel_path=cgls_veg:slucm/lai_cgls/ rel_path=modis_veg_no_search:lai_modis_30s/ rel_path=modis_lai:lai_modis_30s/ rel_path= lowres:lai_modis_10m/ @@ -993,8 +993,8 @@ name = MH_URB2D fill_missing = 0. interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_urb:slucm_kvnk/urban_params/2010/H_avg/ - rel_path = y2050_urb:slucm_kvnk/urban_params/2050/H_avg/ + rel_path = y2010_urb:slucm/urban_params/2010/H_avg/ + rel_path = y2050_urb:slucm/urban_params/2050/H_avg/ =============================== name = ZD_URB2D priority=1 @@ -1003,8 +1003,8 @@ name = ZD_URB2D fill_missing = 0. interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_urb:slucm_kvnk/urban_params/2010/d/ - rel_path = y2050_urb:slucm_kvnk/urban_params/2050/d/ + rel_path = y2010_urb:slucm/urban_params/2010/d/ + rel_path = y2050_urb:slucm/urban_params/2050/d/ =============================== name = Z0_URB2D priority=1 @@ -1013,8 +1013,8 @@ name = Z0_URB2D fill_missing = 0. interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_urb:slucm_kvnk/urban_params/2010/z_0/ - rel_path = y2050_urb:slucm_kvnk/urban_params/2050/z_0/ + rel_path = y2010_urb:slucm/urban_params/2010/z_0/ + rel_path = y2050_urb:slucm/urban_params/2050/z_0/ =============================== name = BUILD_AREA_FRACTION priority=1 @@ -1023,8 +1023,8 @@ name = BUILD_AREA_FRACTION fill_missing = 0. interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_urb:slucm_kvnk/urban_params/2010/lambda_p/ - rel_path = y2050_urb:slucm_kvnk/urban_params/2050/lambda_p/ + rel_path = y2010_urb:slucm/urban_params/2010/lambda_p/ + rel_path = y2050_urb:slucm/urban_params/2050/lambda_p/ =============================== name = LF_URB2D_S priority=1 @@ -1033,8 +1033,8 @@ name = LF_URB2D_S fill_missing = 0. interp_option = y2010_urb:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_urb:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_urb:slucm_kvnk/urban_params/2010/lambda_f/ - rel_path = y2050_urb:slucm_kvnk/urban_params/2050/lambda_f/ + rel_path = y2010_urb:slucm/urban_params/2010/lambda_f/ + rel_path = y2050_urb:slucm/urban_params/2050/lambda_f/ =============================== name = AHE priority=1 @@ -1043,8 +1043,8 @@ name = AHE fill_missing = 0. interp_option = y2010_ahe:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_ahe:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_ahe:slucm_kvnk/AHE_2010s/2_bytes/ - rel_path = y2050_ahe:slucm_kvnk/AHE_2050s/2_bytes/ + rel_path = y2010_ahe:slucm/AHE_2010s/2_bytes/ + rel_path = y2050_ahe:slucm/AHE_2050s/2_bytes/ z_dim_name = month_hour =============================== name = AHE @@ -1054,7 +1054,7 @@ name = AHE fill_missing = 0. interp_option = y2010_ahe:average_gcell(4.0)+four_pt+average_4pt interp_option = y2050_ahe:average_gcell(4.0)+four_pt+average_4pt - rel_path = y2010_ahe:slucm_kvnk/AHE_2010s/3_bytes/ - rel_path = y2050_ahe:slucm_kvnk/AHE_2050s/3_bytes/ + rel_path = y2010_ahe:slucm/AHE_2010s/3_bytes/ + rel_path = y2050_ahe:slucm/AHE_2050s/3_bytes/ z_dim_name = month_hour ===============================