[Ncep.list.emc.fv3gfs_tickets] FV3GFS Ticket #12: Changes for nests

[FV3GFS Trac Ticket]

#12: Changes for nests
--------------------------+----------------------------
  Reporter:  tom.black@…  |      Owner:  george.gayno@…
      Type:  enhancement  |     Status:  accepted
  Priority:  major        |  Milestone:
 Component:  component1   |    Version:
Resolution:               |   Keywords:  nest
--------------------------+----------------------------

Comment (by george.gayno@…):

 Next, the `global_chgres` program was tested on Luna for the 'uniform',
 'stretched' and 'nest' cases.  The grid and orography outputs from the C96
 tests created under this ticket were used to define the grid.  To do so, I
 had to modify the `global_chgres.sh` script to point to these files.  For
 the 'nest', I did:

 {{{
 /gpfs/hps/emc/global/noscrub/George.Gayno/fv3gfs/branches/fv3gfs.tom.black/global_shared.v15.0.0/ush\
 $ svn diff fv3gfs_chgres.sh
 Index: fv3gfs_chgres.sh
 ===================================================================
 --- fv3gfs_chgres.sh    (revision 90988)
 +++ fv3gfs_chgres.sh    (working copy)
 @@ -187,8 +187,8 @@
  $NLN NULL    chgres.inp.sfc

  for tile in `seq 1 $ntiles`; do
 -  $NLN $FIX_FV3/$CASE/${CASE}_grid.tile${tile}.nc
 chgres.fv3.grd.t$tile
 -  $NLN $FIX_FV3/$CASE/${CASE}_oro_data.tile${tile}.nc
 chgres.fv3.orog.t$tile
 +  $NLN
 /gpfs/hps/ptmp/George.Gayno/fv3_grid.nest/C96r15n3_c96s/grid/${CASE}_grid.tile${tile}.nc
 chgres.fv3.grd.t$tile
 +  $NLN
 /gpfs/hps/ptmp/George.Gayno/fv3_grid.nest/C96r15n3_c96s/orog/oro.${CASE}.tile${tile}.nc
 chgres.fv3.orog.t$tile
  done

  $APRUN $GLOBAL_CHGRES_EXEC << EOF '1>&1' '2>&2'
 @@ -215,8 +215,8 @@

  for tile in `seq 1 $ntiles`; do

 -  $NLN $FIX_FV3/$CASE/${CASE}_grid.tile$tile.nc     chgres$tile
 -  $NLN $FIX_FV3/$CASE/${CASE}_oro_data.tile$tile.nc
 chgres.fv3.orog.t$tile
 +  $NLN
 /gpfs/hps/ptmp/George.Gayno/fv3_grid.nest/C96r15n3_c96s/grid/${CASE}_grid.tile${tile}.nc
 chgres$tile
 +  $NLN
 /gpfs/hps/ptmp/George.Gayno/fv3_grid.nest/C96r15n3_c96s/orog/oro.${CASE}.tile${tile}.nc
 chgres.fv3.orog.t$tile

    $APRUN $GLOBAL_CHGRES_EXEC << EOF '1>&1' '2>&2'
 }}}

 And I made the following changes to the `fv3gfs_driver_chgres.sh` script
 to point to my branch and to define the stretch factor and refine ratio
 (which appear to only be used to name the output directory):
 {{{
 /gpfs/hps/emc/global/noscrub/George.Gayno/fv3gfs/branches/fv3gfs.tom.black/global_shared.v15.0.0/ush\
 $ svn diff fv3gfs_driver_chgres.sh
 Index: fv3gfs_driver_chgres.sh
 ===================================================================
 --- fv3gfs_driver_chgres.sh     (revision 90988)
 +++ fv3gfs_driver_chgres.sh     (working copy)
 @@ -5,8 +5,8 @@
  #BSUB -oo log.chgres
  #BSUB -eo log.chgres
  #BSUB -J chgres_fv3
 -#BSUB -q devonprod
 -#BSUB -W 06:00
 +#BSUB -q dev
 +#BSUB -W 01:00
  #BSUB -M 1024
  #BSUB -extsched 'CRAYLINUX[]'

 @@ -51,7 +51,7 @@
  ulimit -s unlimited

  export CASE=${CASE:-C96}                     # resolution of tile: 48,
 96, 192, 384, 768, 1152, 3072
 -export CDATE=${CDATE:-${cdate:-2017030800}}  # format yyyymmddhh
 yyyymmddhh ...
 +export CDATE=${CDATE:-${cdate:-2017040600}}  # format yyyymmddhh
 yyyymmddhh ...

  export machine=${machine:-WCOSS_C}
  export NODES=1
 @@ -82,10 +82,10 @@
   exit
  fi

 -export BASE_GSM=${BASE_GSM:-/gpfs/hps/ptmp/emc.glopara/EXP-
 cyc/global_shared.v15.0.0}
 +export
 BASE_GSM=${BASE_GSM:-/gpfs/hps/emc/global/noscrub/George.Gayno/fv3gfs/branches/fv3gfs.tom.black/global_shared.v15.0.0}
  export script_dir=$BASE_GSM/ush

 -export gtype=uniform             # grid type = uniform, stretch, or
 nested
 +export gtype=nest          # grid type = uniform, stretch, or nested

  if [ ! -s $out_dir ]; then mkdir -p $out_dir; fi
  if [ ! -s $tmp_dir ]; then mkdir -p $tmp_dir; fi
 @@ -98,15 +98,15 @@
    export name=${CASE}
    export ntiles=6
  elif [ $gtype = stretch ]; then
 -  export stetch_fac=                            # Stretching factor for
 the grid
 +  export stetch_fac=1.5                         # Stretching factor for
 the grid
    export rn=$( echo "$stetch_fac * 10" | bc | cut -c1-2 )
    export name=${CASE}r${rn}                     # identifier based on
 refined location (same as grid)
    export ntiles=6
    echo "creating stretched ICs"
  elif [ $gtype = nest ]; then
 -  export stetch_fac=                            # Stretching factor for
 the grid
 +  export stetch_fac=1.5                         # Stretching factor for
 the grid
    export rn=$( echo "$stetch_fac * 10" | bc | cut -c1-2 )
 -  export refine_ratio=                          # Specify the refinement
 ratio for nest grid
 +  export refine_ratio=3                         # Specify the refinement
 ratio for nest grid
    export name=${CASE}r${rn}n${refine_ratio}      # identifier based on
 nest location (same as grid)
    export ntiles=7
 }}}

 All three tests ran to completion.  The output was spot checked using
 Grads.  No problems were noted.

--
Ticket URL: <https://svnemc.ncep.noaa.gov/trac/fv3gfs/ticket/12#comment:16>
fv3gfs <https://svnemc.ncep.noaa.gov/trac/fv3gfs>
NGGPS FV3GFS Development