From Ning Wang: this commit expands lake fraction and depth calculati…

…on to both global and regional domain (gtype=regional_gfdl). The resulting lake fraction and depth are added to oro data.
ufs-community · Oct 26, 2020 · 67be759 · 67be759
1 parent 5186965
commit 67be759
Show file tree

Hide file tree

Showing 5 changed files with 76 additions and 78 deletions.
diff --git a/driver_scripts/driver_grid.cray.sh b/driver_scripts/driver_grid.cray.sh
@@ -65,32 +65,32 @@ module list
 # Set grid specs here.
 #-----------------------------------------------------------------------
 
-export gtype=uniform       # 'uniform', 'stretch', 'nest'
-                           # 'regional_gfdl', 'regional_esg'
+export gtype=uniform    # 'uniform', 'stretch', 'nest', 
+                        # 'regional_gfdl', 'regional_esg'
+export lake_cutoff=0.20 # lake frac < lake_cutoff ignored when add_lake=T
 
 if [ $gtype = uniform ]; then
   export res=96
-  export add_lake=false   # Add lake frac and depth to orography data.
-                          # Uniform grids only.
-  export lake_cutoff=0.20 # lake frac less than lake_cutoff is ignored
+  export add_lake=false        # Add lake frac and depth to orography data.
 elif [ $gtype = stretch ]; then
   export res=96
   export stretch_fac=1.5       # Stretching factor for the grid
   export target_lon=-97.5      # Center longitude of the highest resolution tile
   export target_lat=35.5       # Center latitude of the highest resolution tile
 elif [ $gtype = nest ] || [ $gtype = regional_gfdl ]; then
-  export res=96
+  export add_lake=false        # Add lake frac and depth to orography data.
+  export res=768
   export stretch_fac=1.5       # Stretching factor for the grid
   export target_lon=-97.5      # Center longitude of the highest resolution tile
-  export target_lat=35.5       # Center latitude of the highest resolution tile
+  export target_lat=38.5       # Center latitude of the highest resolution tile
   export refine_ratio=3        # The refinement ratio
-  export istart_nest=27        # Starting i-direction index of nest grid in parent tile supergrid
-  export jstart_nest=37        # Starting j-direction index of nest grid in parent tile supergrid
-  export iend_nest=166         # Ending i-direction index of nest grid in parent tile supergrid
-  export jend_nest=164         # Ending j-direction index of nest grid in parent tile supergrid
-  export halo=3
+  export istart_nest=123       # Starting i-direction index of nest grid in parent tile supergrid
+  export jstart_nest=331       # Starting j-direction index of nest grid in parent tile supergrid
+  export iend_nest=1402        # Ending i-direction index of nest grid in parent tile supergrid
+  export jend_nest=1194        # Ending j-direction index of nest grid in parent tile supergrid
+  export halo=3                # Lateral boundary halo
 elif [ $gtype = regional_esg ] ; then
-  export res=-999              # equivalent res is computed.
+  export res=-999              # equivalent resolution is computed
   export target_lon=-97.5      # Center longitude of grid
   export target_lat=35.5       # Center latitude of grid
   export idim=301              # Dimension of grid in 'i' direction

diff --git a/driver_scripts/driver_grid.dell.sh b/driver_scripts/driver_grid.dell.sh
@@ -49,7 +49,7 @@
 #   6) For "regional_gfdl" grids, set the "halo".  Default is three
 #      rows/columns.
 #   7) For "regional_esg" grids, set center lat/lon of grid,
-#      - "target_lat/lon" - the i/j dimensions - "i/jdim", the 
+#      - "target_lat/lon" - the i/j dimensions - "i/jdim", the
 #      x/y grid spacing - "delx/y", and halo.
 #   8) Set working directory - TEMP_DIR - and path to the repository
 #      clone - home_dir.
@@ -67,32 +67,32 @@ module list
 # Set grid specs here.
 #-----------------------------------------------------------------------
 
-export gtype=uniform  # 'uniform', 'stretch', 'nest', 
-                      # 'regional_gfdl', 'regional_esg'
+export gtype=uniform    # 'uniform', 'stretch', 'nest', 
+                        # 'regional_gfdl', 'regional_esg'
+export lake_cutoff=0.20 # lake frac < lake_cutoff ignored when add_lake=T
 
 if [ $gtype = uniform ]; then
   export res=96
-  export add_lake=false   # Add lake frac and depth to orography data.
-                          # Uniform grids only.
-  export lake_cutoff=0.20 # lake frac less than lake_cutoff is ignored
+  export add_lake=false        # Add lake frac and depth to orography data.
 elif [ $gtype = stretch ]; then
   export res=96
   export stretch_fac=1.5       # Stretching factor for the grid
   export target_lon=-97.5      # Center longitude of the highest resolution tile
   export target_lat=35.5       # Center latitude of the highest resolution tile
 elif [ $gtype = nest ] || [ $gtype = regional_gfdl ]; then
-  export res=96
+  export add_lake=false        # Add lake frac and depth to orography data.
+  export res=768
   export stretch_fac=1.5       # Stretching factor for the grid
   export target_lon=-97.5      # Center longitude of the highest resolution tile
-  export target_lat=35.5       # Center latitude of the highest resolution tile
+  export target_lat=38.5       # Center latitude of the highest resolution tile
   export refine_ratio=3        # The refinement ratio
-  export istart_nest=27        # Starting i-direction index of nest grid in parent tile supergrid
-  export jstart_nest=37        # Starting j-direction index of nest grid in parent tile supergrid
-  export iend_nest=166         # Ending i-direction index of nest grid in parent tile supergrid
-  export jend_nest=164         # Ending j-direction index of nest grid in parent tile supergrid
-  export halo=3
+  export istart_nest=123       # Starting i-direction index of nest grid in parent tile supergrid
+  export jstart_nest=331       # Starting j-direction index of nest grid in parent tile supergrid
+  export iend_nest=1402        # Ending i-direction index of nest grid in parent tile supergrid
+  export jend_nest=1194        # Ending j-direction index of nest grid in parent tile supergrid
+  export halo=3                # Lateral boundary halo
 elif [ $gtype = regional_esg ] ; then
-  export res=-999              # equivalent res is computed.
+  export res=-999              # equivalent resolution is computed
   export target_lon=-97.5      # Center longitude of grid
   export target_lat=35.5       # Center latitude of grid
   export idim=301              # Dimension of grid in 'i' direction
@@ -109,7 +109,7 @@ fi
 #-----------------------------------------------------------------------
 # Check paths.
 #   home_dir - location of repository.
-#   TEMP_DIR   - working directory.
+#   TEMP_DIR - working directory.
 #   out_dir  - where files will be placed upon completion.
 #-----------------------------------------------------------------------
 

diff --git a/driver_scripts/driver_grid.hera.sh b/driver_scripts/driver_grid.hera.sh
@@ -23,7 +23,7 @@
 # Note: The sfc_climo_gen program only runs with an
 #       mpi task count that is a multiple of six.  This is
 #       an ESMF library requirement.  Large grids may require
-#       tasks spread across multiple nodes.  The orography code
+#       tasks spread across multiple nodes. The orography code
 #       benefits from threads.
 #
 # To run, do the following:
@@ -68,34 +68,32 @@ module list
 # Set grid specs here.
 #-----------------------------------------------------------------------
 
-export gtype=uniform       # 'uniform', 'stretch', 'nest'
-                           # 'regional_gfdl', 'regional_esg'
+export gtype=uniform    # 'uniform', 'stretch', 'nest', 
+                        # 'regional_gfdl', 'regional_esg'
+export lake_cutoff=0.20 # lake frac < lake_cutoff ignored when add_lake=T
 
 if [ $gtype = uniform ]; then
   export res=96
-  export add_lake=false   # Add lake frac and depth to orography data.
-                          # Uniform grids only.
-  export lake_cutoff=0.20 # lake frac less than lake_cutoff is ignored
+  export add_lake=false        # Add lake frac and depth to orography data.
 elif [ $gtype = stretch ]; then
   export res=96
   export stretch_fac=1.5       # Stretching factor for the grid
   export target_lon=-97.5      # Center longitude of the highest resolution tile
   export target_lat=35.5       # Center latitude of the highest resolution tile
 elif [ $gtype = nest ] || [ $gtype = regional_gfdl ]; then
-  export res=96
-  export add_lake=true   # Add lake frac and depth to orography data.
+  export add_lake=false        # Add lake frac and depth to orography data.
   export res=768
   export stretch_fac=1.5       # Stretching factor for the grid
   export target_lon=-97.5      # Center longitude of the highest resolution tile
   export target_lat=38.5       # Center latitude of the highest resolution tile
   export refine_ratio=3        # The refinement ratio
-  export istart_nest=123        # Starting i-direction index of nest grid in parent tile supergrid
-  export jstart_nest=331        # Starting j-direction index of nest grid in parent tile supergrid
-  export iend_nest=1402         # Ending i-direction index of nest grid in parent tile supergrid
-  export jend_nest=1194         # Ending j-direction index of nest grid in parent tile supergrid
+  export istart_nest=123       # Starting i-direction index of nest grid in parent tile supergrid
+  export jstart_nest=331       # Starting j-direction index of nest grid in parent tile supergrid
+  export iend_nest=1402        # Ending i-direction index of nest grid in parent tile supergrid
+  export jend_nest=1194        # Ending j-direction index of nest grid in parent tile supergrid
   export halo=3                # Lateral boundary halo
 elif [ $gtype = regional_esg ] ; then
-  export res=-999              # equivalent res is computed.
+  export res=-999              # equivalent resolution is computed
   export target_lon=-97.5      # Center longitude of grid
   export target_lat=35.5       # Center latitude of grid
   export idim=301              # Dimension of grid in 'i' direction
@@ -110,10 +108,10 @@ elif [ $gtype = regional_esg ] ; then
 fi
 
 #-----------------------------------------------------------------------
-# Check paths.  
+# Check paths.
 #   home_dir - location of repository.
-#   TEMP_DIR   - working directory.
-#   out_dir  - where files will be placed upon completion. 
+#   TEMP_DIR - working directory.
+#   out_dir  - where files will be placed upon completion.
 #-----------------------------------------------------------------------
 
 export home_dir=$SLURM_SUBMIT_DIR/..

diff --git a/driver_scripts/driver_grid.jet.sh b/driver_scripts/driver_grid.jet.sh
@@ -24,8 +24,8 @@
 # Note: The sfc_climo_gen program only runs with an
 #       mpi task count that is a multiple of six.  This is
 #       an ESMF library requirement.  Large grids may require
-#       tasks spread across multiple nodes.  The orography
-#       code benefits from threads.
+#       tasks spread across multiple nodes. The orography code
+#       benefits from threads.
 #
 # To run, do the following:
 #
@@ -68,29 +68,29 @@ module list
 # Set grid specs here.
 #-----------------------------------------------------------------------
 
-export gtype=uniform       # 'uniform', 'stretch', 'nest'
-                           # 'regional_gfdl', 'regional_esg'
+export gtype=uniform    # 'uniform', 'stretch', 'nest', 
+                        # 'regional_gfdl', 'regional_esg'
+export lake_cutoff=0.20 # lake frac < lake_cutoff ignored when add_lake=T
 
 if [ $gtype = uniform ]; then
   export res=96
-  export add_lake=false   # Add lake frac and depth to orography data.
-                          # Uniform grids only.
-  export lake_cutoff=0.20 # lake frac less than lake_cutoff is ignored
+  export add_lake=false        # Add lake frac and depth to orography data.
 elif [ $gtype = stretch ]; then
   export res=96
   export stretch_fac=1.5       # Stretching factor for the grid
   export target_lon=-97.5      # Center longitude of the highest resolution tile
   export target_lat=35.5       # Center latitude of the highest resolution tile
 elif [ $gtype = nest ] || [ $gtype = regional_gfdl ]; then
-  export res=96
+  export add_lake=false        # Add lake frac and depth to orography data.
+  export res=768
   export stretch_fac=1.5       # Stretching factor for the grid
   export target_lon=-97.5      # Center longitude of the highest resolution tile
-  export target_lat=35.5       # Center latitude of the highest resolution tile
+  export target_lat=38.5       # Center latitude of the highest resolution tile
   export refine_ratio=3        # The refinement ratio
-  export istart_nest=27        # Starting i-direction index of nest grid in parent tile supergrid
-  export jstart_nest=37        # Starting j-direction index of nest grid in parent tile supergrid
-  export iend_nest=166         # Ending i-direction index of nest grid in parent tile supergrid
-  export jend_nest=164         # Ending j-direction index of nest grid in parent tile supergrid
+  export istart_nest=123       # Starting i-direction index of nest grid in parent tile supergrid
+  export jstart_nest=331       # Starting j-direction index of nest grid in parent tile supergrid
+  export iend_nest=1402        # Ending i-direction index of nest grid in parent tile supergrid
+  export jend_nest=1194        # Ending j-direction index of nest grid in parent tile supergrid
   export halo=3                # Lateral boundary halo
 elif [ $gtype = regional_esg ] ; then
   export res=-999              # equivalent resolution is computed
@@ -108,10 +108,10 @@ elif [ $gtype = regional_esg ] ; then
 fi
 
 #-----------------------------------------------------------------------
-# Check paths.  
+# Check paths.
 #   home_dir - location of repository.
-#   TEMP_DIR   - working directory.
-#   out_dir  - where files will be placed upon completion. 
+#   TEMP_DIR - working directory.
+#   out_dir  - where files will be placed upon completion.
 #-----------------------------------------------------------------------
 
 export home_dir=$SLURM_SUBMIT_DIR/..

diff --git a/driver_scripts/driver_grid.orion.sh b/driver_scripts/driver_grid.orion.sh
@@ -23,7 +23,7 @@
 # Note: The sfc_climo_gen program only runs with an
 #       mpi task count that is a multiple of six.  This is
 #       an ESMF library requirement.  Large grids may require
-#       tasks spread across multiple nodes.  The orography code
+#       tasks spread across multiple nodes. The orography code
 #       benefits from threads.
 #
 # To run, do the following:
@@ -34,8 +34,8 @@
 #         "stretch"  - global stretched grid
 #         "nest"     - global stretched grid with nest
 #         "regional_gfdl" - stand-alone gfdl regional grid
-#         "regional_esg"  - stand-alone extended Schmidt
-#                           gnomonic (esg) regional grid
+#         "regional_esg"  - stand-alone extended Schmidt gnomonic
+#                           (esg) regional grid
 #   3) For "uniform" grids - to include lake fraction and
 #      depth, set "add_lake" to true, and the "lake_cutoff" value.
 #   4) For "stretch" and "nest" grids, set the stretching factor -
@@ -67,32 +67,32 @@ module list
 # Set grid specs here.
 #-----------------------------------------------------------------------
 
-export gtype=uniform   # 'uniform', 'stretch', 'nest', 
-                       # 'regional_gfdl', 'regional_esg'
+export gtype=uniform    # 'uniform', 'stretch', 'nest', 
+                        # 'regional_gfdl', 'regional_esg'
+export lake_cutoff=0.20 # lake frac < lake_cutoff ignored when add_lake=T
 
 if [ $gtype = uniform ]; then
   export res=96
-  export add_lake=false   # Add lake frac and depth to orography data.
-                          # Uniform grids only.
-  export lake_cutoff=0.20 # lake frac less than lake_cutoff is ignored
+  export add_lake=false        # Add lake frac and depth to orography data.
 elif [ $gtype = stretch ]; then
   export res=96
   export stretch_fac=1.5       # Stretching factor for the grid
   export target_lon=-97.5      # Center longitude of the highest resolution tile
   export target_lat=35.5       # Center latitude of the highest resolution tile
 elif [ $gtype = nest ] || [ $gtype = regional_gfdl ]; then
-  export res=96
+  export add_lake=false        # Add lake frac and depth to orography data.
+  export res=768
   export stretch_fac=1.5       # Stretching factor for the grid
   export target_lon=-97.5      # Center longitude of the highest resolution tile
-  export target_lat=35.5       # Center latitude of the highest resolution tile
+  export target_lat=38.5       # Center latitude of the highest resolution tile
   export refine_ratio=3        # The refinement ratio
-  export istart_nest=27        # Starting i-direction index of nest grid in parent tile supergrid
-  export jstart_nest=37        # Starting j-direction index of nest grid in parent tile supergrid
-  export iend_nest=166         # Ending i-direction index of nest grid in parent tile supergrid
-  export jend_nest=164         # Ending j-direction index of nest grid in parent tile supergrid
+  export istart_nest=123       # Starting i-direction index of nest grid in parent tile supergrid
+  export jstart_nest=331       # Starting j-direction index of nest grid in parent tile supergrid
+  export iend_nest=1402        # Ending i-direction index of nest grid in parent tile supergrid
+  export jend_nest=1194        # Ending j-direction index of nest grid in parent tile supergrid
   export halo=3                # Lateral boundary halo
 elif [ $gtype = regional_esg ] ; then
-  export res=-999              # equivalent res is computed.
+  export res=-999              # equivalent resolution is computed
   export target_lon=-97.5      # Center longitude of grid
   export target_lat=35.5       # Center latitude of grid
   export idim=301              # Dimension of grid in 'i' direction
@@ -107,10 +107,10 @@ elif [ $gtype = regional_esg ] ; then
 fi
 
 #-----------------------------------------------------------------------
-# Check paths.  
+# Check paths.
 #   home_dir - location of repository.
-#   TEMP_DIR   - working directory.
-#   out_dir  - where files will be placed upon completion. 
+#   TEMP_DIR - working directory.
+#   out_dir  - where files will be placed upon completion.
 #-----------------------------------------------------------------------
 
 export home_dir=$SLURM_SUBMIT_DIR/..