Add extended Schmidt gnomonic grid capability

driver_scripts/driver_grid.cray.sh

@@ -33,19 +33,24 @@
 #         "uniform"  - global uniform grid
 #         "stretch"  - global stretched grid
 #         "nest"     - global stretched grid with nest
-#         "regional" - stand-alone regional grid
+#         "regional_gfdl" - stand-alone gfdl regional grid
+#         "regional_esg"  - stand-alone extended Schmidt gnomonic
+#                           (esg) regional grid
 #   3) For "stretch" and "nest" grids, set the stretching factor -
 #       "stretch_fac", and center lat/lon of highest resolution
 #      tile - "target_lat" and "target_lon".
 #   4) For "nest" grids, set the refinement ratio - "refine_ratio",
 #      the starting/ending i/j index location within the parent
 #      tile - "istart_nest", "jstart_nest", "iend_nest", "jend_nest"
-#   5) For "regional" grids, set the "halo".  Default is three
+#   5) For "regional_gfdl" grids, set the "halo".  Default is three
 #      rows/columns.
-#   6) Set working directory - TMPDIR - and path to the repository
+#   6) For "regional_esg" grids, set center lat/lon of grid,
+#      - "target_lat/lon" - the i/j dimensions - "i/jdim", the
+#      x/y grid spacing - "delx/y", and halo.
+#   7) Set working directory - TMPDIR - and path to the repository
 #      clone - home_dir.
-#   7) Submit script: "cat $script | bsub".
-#   8) All files will be placed in "out_dir".
+#   8) Submit script: "cat $script | bsub".
+#   9) All files will be placed in "out_dir".
 #
 #-----------------------------------------------------------------------
 
@@ -57,14 +62,18 @@ module list
 # Set grid specs here.
 #-----------------------------------------------------------------------
 
-export res=96
-export gtype=regional   # 'uniform', 'stretch', 'nest', or 'regional'
+export gtype=uniform        # 'uniform', 'stretch', 'nest'
+                           # 'regional_gfdl', 'regional_esg'
 
-if [ $gtype = stretch ]; then
+if [ $gtype = uniform ]; then
+  export res=96
+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 ]; then
+elif [ $gtype = nest ] || [ $gtype = regional_gfdl ]; 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
@@ -74,6 +83,19 @@ elif [ $gtype = nest ] || [ $gtype = regional ]; then
   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
+elif [ $gtype = regional_esg ] ; then
+  export res=-999              # equivalent res 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
+  export jdim=200              # Dimension of grid in 'j' direction
+  export delx=0.0585           # Grid spacing (in degrees) in the 'i' direction
+                               # on the SUPERGRID (which has twice the resolution of
+                               # the model grid).  The physical grid spacing in the 'i'
+                               # direction is related to delx as follows:
+                               #    distance = 2*delx*(circumf_Earth/360 deg)
+  export dely=0.0585           # Grid spacing (in degrees) in the 'j' direction.
+  export halo=3                # number of row/cols for halo
 fi
 
 #-----------------------------------------------------------------------
@@ -85,7 +107,7 @@ fi
 
 export home_dir=$LS_SUBCWD/..
 export TMPDIR=/gpfs/hps3/stmp/$LOGNAME/fv3_grid.$gtype
-export out_dir=/gpfs/hps3/stmp/$LOGNAME/C${res}
+export out_dir=/gpfs/hps3/stmp/$LOGNAME/my_grids
 
 export NODES=1
 export APRUN="aprun -n 1 -N 1 -j 1 -d 1 -cc depth"

driver_scripts/driver_grid.dell.sh

@@ -35,19 +35,24 @@
 #         "uniform"  - global uniform grid
 #         "stretch"  - global stretched grid
 #         "nest"     - global stretched grid with nest
-#         "regional" - stand-alone regional grid
+#         "regional_gfdl" - stand-alone gfdl regional grid
+#         "regional_esg"  - stand-alone extended Schmidt gnomonic
+#                           (esg) regional grid
 #   3) For "stretch" and "nest" grids, set the stretching factor -
 #       "stretch_fac", and center lat/lon of highest resolution
 #      tile - "target_lat" and "target_lon".
 #   4) For "nest" grids, set the refinement ratio - "refine_ratio",
 #      the starting/ending i/j index location within the parent
 #      tile - "istart_nest", "jstart_nest", "iend_nest", "jend_nest"
-#   5) For "regional" grids, set the "halo".  Default is three
+#   5) For "regional_gfdl" grids, set the "halo".  Default is three
 #      rows/columns.
-#   6) Set working directory - TMPDIR - and path to the repository
+#   6) For "regional_esg" grids, set center lat/lon of grid,
+#      - "target_lat/lon" - the i/j dimensions - "i/jdim", the 
+#      x/y grid spacing - "delx/y", and halo.
+#   7) Set working directory - TMPDIR - and path to the repository
 #      clone - home_dir.
-#   7) Submit script: "cat $script | bsub".
-#   8) All files will be placed in "out_dir".
+#   8) Submit script: "cat $script | bsub".
+#   9) All files will be placed in "out_dir".
 #
 #-----------------------------------------------------------------------
 
@@ -59,14 +64,18 @@ module list
 # Set grid specs here.
 #-----------------------------------------------------------------------
 
-export res=96
-export gtype=regional  # 'uniform', 'stretch', 'nest', or 'regional'
+export gtype=uniform  # 'uniform', 'stretch', 'nest', 
+                           # 'regional_gfdl', 'regional_esg'
 
-if [ $gtype = stretch ]; then
+if [ $gtype = uniform ]; then
+  export res=96
+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 ]; then
+elif [ $gtype = nest ] || [ $gtype = regional_gfdl ]; 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
@@ -76,6 +85,19 @@ elif [ $gtype = nest ] || [ $gtype = regional ]; then
   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
+elif [ $gtype = regional_esg ] ; then
+  export res=-999              # equivalent res 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
+  export jdim=200              # Dimension of grid in 'j' direction
+  export delx=0.0585           # Grid spacing (in degrees) in the 'i' direction
+                               # on the SUPERGRID (which has twice the resolution of
+                               # the model grid).  The physical grid spacing in the 'i'
+                               # direction is related to delx as follows:
+                               #    distance = 2*delx*(circumf_Earth/360 deg)
+  export dely=0.0585           # Grid spacing (in degrees) in the 'j' direction.
+  export halo=3                # number of row/cols for halo
 fi
 
 #-----------------------------------------------------------------------
@@ -87,7 +109,7 @@ fi
 
 export home_dir=$LS_SUBCWD/..
 export TMPDIR=/gpfs/dell1/stmp/$LOGNAME/fv3_grid.$gtype
-export out_dir=/gpfs/dell1/stmp/$LOGNAME/C${res}
+export out_dir=/gpfs/dell1/stmp/$LOGNAME/my_grids
 
 #-----------------------------------------------------------------------
 # Should not need to change anything below here.

driver_scripts/driver_grid.hera.sh

@@ -33,19 +33,24 @@
 #         "uniform"  - global uniform grid
 #         "stretch"  - global stretched grid
 #         "nest"     - global stretched grid with nest
-#         "regional" - stand-alone regional grid
+#         "regional_gfdl" - stand-alone gfdl regional grid
+#         "regional_esg"  - stand-alone extended Schmidt gnomonic
+#                           (esg) regional grid
 #   3) For "stretch" and "nest" grids, set the stretching factor -
 #       "stretch_fac", and center lat/lon of highest resolution
 #      tile - "target_lat" and "target_lon".
 #   4) For "nest" grids, set the refinement ratio - "refine_ratio",
 #      the starting/ending i/j index location within the parent
 #      tile - "istart_nest", "jstart_nest", "iend_nest", "jend_nest"
-#   5) For "regional" grids, set the "halo".  Default is three
+#   5) For "regional_gfdl" grids, set the "halo".  Default is three
 #      rows/columns.
-#   6) Set working directory - TMPDIR - and path to the repository
+#   6) For "regional_esg" grids, set center lat/lon of grid,
+#      - "target_lat/lon" - the i/j dimensions - "i/jdim", the
+#      x/y grid spacing - "delx/y", and halo.
+#   7) Set working directory - TMPDIR - and path to the repository
 #      clone - home_dir.
-#   7) Submit script: "sbatch $script".
-#   8) All files will be placed in "out_dir".
+#   8) Submit script: "sbatch $script".
+#   9) All files will be placed in "out_dir".
 #
 #-----------------------------------------------------------------------
 
@@ -59,14 +64,18 @@ module list
 # Set grid specs here.
 #-----------------------------------------------------------------------
 
-export res=96
-export gtype=uniform  # 'uniform', 'stretch', 'nest', or 'regional'
+export gtype=uniform       # 'uniform', 'stretch', 'nest'
+                           # 'regional_gfdl', 'regional_esg'
 
-if [ $gtype = stretch ]; then
+if [ $gtype = uniform ]; then
+  export res=96
+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 ]; then
+elif [ $gtype = nest ] || [ $gtype = regional_gfdl ]; 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
@@ -76,6 +85,19 @@ elif [ $gtype = nest ] || [ $gtype = regional ]; then
   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                # Lateral boundary halo
+elif [ $gtype = regional_esg ] ; then
+  export res=-999              # equivalent res 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
+  export jdim=200              # Dimension of grid in 'j' direction
+  export delx=0.0585           # Grid spacing (in degrees) in the 'i' direction
+                               # on the SUPERGRID (which has twice the resolution of
+                               # the model grid).  The physical grid spacing in the 'i'
+                               # direction is related to delx as follows:
+                               #    distance = 2*delx*(circumf_Earth/360 deg)
+  export dely=0.0585           # Grid spacing (in degrees) in the 'j' direction.
+  export halo=3                # number of row/cols for halo
 fi
 
 #-----------------------------------------------------------------------
@@ -87,7 +109,7 @@ fi
 
 export home_dir=$SLURM_SUBMIT_DIR/..
 export TMPDIR=/scratch2/NCEPDEV/stmp1/$LOGNAME/fv3_grid.$gtype
-export out_dir=/scratch2/NCEPDEV/stmp1/$LOGNAME/C${res}
+export out_dir=/scratch2/NCEPDEV/stmp1/$LOGNAME/my_grids
 
 #-----------------------------------------------------------------------
 # Should not need to change anything below here.

driver_scripts/driver_grid.jet.sh

@@ -1,13 +1,13 @@
 #!/bin/bash
 
 #SBATCH -J fv3_grid_driver
-#SBATCH -A emcda
+#SBATCH -A hfv3gfs
 #SBATCH --open-mode=truncate
 #SBATCH -o log.fv3_grid_driver
 #SBATCH -e log.fv3_grid_driver
 #SBATCH --nodes=1 --ntasks-per-node=24
 #SBATCH --partition=xjet
-#SBATCH -q windfall
+#SBATCH -q batch
 #SBATCH -t 00:10:00
 
 #-----------------------------------------------------------------------
@@ -34,19 +34,24 @@
 #         "uniform"  - global uniform grid
 #         "stretch"  - global stretched grid
 #         "nest"     - global stretched grid with nest
-#         "regional" - stand-alone regional grid
+#         "regional_gfdl" - stand-alone gfdl regional grid
+#         "regional_esg"  - stand-alone extended Schmidt gnomonic
+#                           (esg) regional grid
 #   3) For "stretch" and "nest" grids, set the stretching factor -
 #       "stretch_fac", and center lat/lon of highest resolution
 #      tile - "target_lat" and "target_lon".
 #   4) For "nest" grids, set the refinement ratio - "refine_ratio",
 #      the starting/ending i/j index location within the parent
 #      tile - "istart_nest", "jstart_nest", "iend_nest", "jend_nest"
-#   5) For "regional" grids, set the "halo".  Default is three
+#   5) For "regional_gfdl" grids, set the "halo".  Default is three
 #      rows/columns.
-#   6) Set working directory - TMPDIR - and path to the repository
+#   6) For "regional_esg" grids, set center lat/lon of grid,
+#      - "target_lat/lon" - the i/j dimensions - "i/jdim", the
+#      x/y grid spacing - "delx/y", and halo.
+#   7) Set working directory - TMPDIR - and path to the repository
 #      clone - home_dir.
-#   7) Submit script: "sbatch $script".
-#   8) All files will be placed in "out_dir".
+#   8) Submit script: "sbatch $script".
+#   9) All files will be placed in "out_dir".
 #
 #-----------------------------------------------------------------------
 
@@ -60,14 +65,18 @@ module list
 # Set grid specs here.
 #-----------------------------------------------------------------------
 
-export res=96
-export gtype=regional  # 'uniform', 'stretch', 'nest', or 'regional'
+export gtype=uniform       # 'uniform', 'stretch', 'nest'
+                           # 'regional_gfdl', 'regional_esg'
 
-if [ $gtype = stretch ]; then
+if [ $gtype = uniform ]; then
+  export res=96
+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 ]; then
+elif [ $gtype = nest ] || [ $gtype = regional_gfdl ]; 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
@@ -77,6 +86,19 @@ elif [ $gtype = nest ] || [ $gtype = regional ]; then
   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                # Lateral boundary halo
+elif [ $gtype = regional_esg ] ; then
+  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
+  export jdim=200              # Dimension of grid in 'j' direction
+  export delx=0.0585           # Grid spacing (in degrees) in the 'i' direction
+                               # on the SUPERGRID (which has twice the resolution of
+                               # the model grid).  The physical grid spacing in the 'i'
+                               # direction is related to delx as follows:
+                               #    distance = 2*delx*(circumf_Earth/360 deg)
+  export dely=0.0585           # Grid spacing (in degrees) in the 'j' direction.
+  export halo=3                # number of row/cols for halo
 fi
 
 #-----------------------------------------------------------------------
@@ -88,7 +110,7 @@ fi
 
 export home_dir=$SLURM_SUBMIT_DIR/..
 export TMPDIR=/lfs4/HFIP/emcda/$LOGNAME/stmp/fv3_grid.$gtype
-export out_dir=/lfs4/HFIP/emcda/$LOGNAME/stmp/C${res}
+export out_dir=/lfs4/HFIP/emcda/$LOGNAME/stmp/my_grids
 
 #-----------------------------------------------------------------------
 # Should not need to change anything below here.

reg_tests/grid_gen/c96.uniform.sh

@@ -35,7 +35,7 @@ echo "Ending at: " `date`
 # Compare output to baseline set of data.
 #-----------------------------------------------------------------------------
 
-cd $out_dir
+cd $out_dir/C96
 
 test_failed=0
 for files in *tile*.nc ./fix_sfc/*tile*.nc

reg_tests/grid_gen/driver.cray.sh

@@ -44,6 +44,7 @@ export machine=WCOSS_C
 export KMP_AFFINITY=disabled
 export NCCMP=/gpfs/hps3/emc/global/noscrub/George.Gayno/util/netcdf/nccmp
 export HOMEreg=/gpfs/hps3/emc/global/noscrub/George.Gayno/ufs_utils.git/reg_tests/grid_gen/baseline_data
+export NCDUMP=/gpfs/hps/usrx/local/prod/NetCDF/4.2/intel/sandybridge/bin/ncdump
 
 rm -fr $WORK_DIR
 
@@ -58,16 +59,16 @@ bsub -e $LOG_FILE -o $LOG_FILE -q $QUEUE -P $PROJECT_CODE -J c96.uniform -W 0:15
         -extsched 'CRAYLINUX[]' "export NODES=1; $PWD/c96.uniform.sh"
 
 #-----------------------------------------------------------------------------
-# C96 regional grid
+# gfdl regional grid
 #-----------------------------------------------------------------------------
 
-bsub -e $LOG_FILE -o $LOG_FILE -q $QUEUE -P $PROJECT_CODE -J c96.regional -W 0:10 -M 2400 \
-        -w 'ended(c96.uniform)' -extsched 'CRAYLINUX[]' "export NODES=1; $PWD/c96.regional.sh"
+bsub -e $LOG_FILE -o $LOG_FILE -q $QUEUE -P $PROJECT_CODE -J gfdl.regional -W 0:10 -M 2400 \
+        -w 'ended(c96.uniform)' -extsched 'CRAYLINUX[]' "export NODES=1; $PWD/gfdl.regional.sh"
 
 #-----------------------------------------------------------------------------
 # Create summary log.
 #-----------------------------------------------------------------------------
 
-bsub -o $LOG_FILE -q $QUEUE -P $PROJECT_CODE -J summary -R "rusage[mem=100]" -W 0:01 -w 'ended(c96.regional)' "grep -a '<<<' $LOG_FILE >> $SUM_FILE"
+bsub -o $LOG_FILE -q $QUEUE -P $PROJECT_CODE -J summary -R "rusage[mem=100]" -W 0:01 -w 'ended(gfdl.regional)' "grep -a '<<<' $LOG_FILE >> $SUM_FILE"
 
 exit