Community GSIv3.5 Online Tutorial

Prepare run script and anavinfo file for basic cases

[Back to the Index]

Setting up the Run Script

The GSI run script run/run_gsi_regional.ksh, which is provided with the source code, must be customized to the local environment. These changes include things such as the analysis time and date, the path to files, and MPI specific information, which are necessary for running all of the basic tutorial cases.

For this tutorial, start by copying the original run script to a working copy
cd ${PATH}/comGSIv3.5_EnKFv1.1/run/
cp run_gsi_regional.ksh run_gsi_regional.ksh_basic
This will be used as the template for all the later scripts used in the regional practice.

Make the following modifications to the script run_gsi_regional.ksh_basic

  • Set the variable GSIPROC to 4
  • Set the variable ARCH to LINUX_LSF, as an example value appropriate for Linux with LSF.
  • In the case set up section of the script, set the environment variables to values appropriate for your environment on your machine. For example:
    • ANAL_TIME=2014061700
    • WORK_ROOT=${PATH}/comGSIv3.5_EnKFv1.1/run/arw_${ANAL_TIME}
    • OBS_ROOT=${DATA_ROOT}/data/20140617/2014061700/obs
    • PREPBUFR=${OBS_ROOT}/nam.t00z.prepbufr.tm00.nr
    • BK_ROOT and BK_FILE
      • For ARW:
        • BK_ROOT=${DATA_ROOT}/data/20140617/2014061700/arw
        • BK_FILE=/${BK_ROOT}/wrfinput_d01.2014061700
      • For NMM:
        • BK_ROOT=${DATA_ROOT}/data/20140617/2014061700/nmm
        • BK_FILE=/${BK_ROOT}/wrfinput_d01.2014061700
    • CRTM_ROOT=${PATH_CRTM}/CRTM_2.2.3
    • GSI_ROOT=${PATH}/comGSIv3.5_EnKFv1.1
    • FIX_ROOT=${GSI_ROOT}/fix
    • GSI_EXE=${GSI_ROOT}/run/gsi.exe
    • GSI_NAMELIST=${GSI_ROOT}/run/comgsi_namelist.sh
    The difference between the ARW core and NMM core is the directory name (BK_ROOT) that holds the background files.

  • On top of the scripts, add an appropriate batch script. For example on Yellowstone:
  • 
    #BSUB -W 00:15                 # wall-clock time (hrs:mins)
    #BSUB -n 4                   # number of tasks in job         
    #BSUB -R "span[ptile=16]"      # run 16 MPI tasks per node
    #BSUB -J gsi                   # job name
    #BSUB -o gsi.%J.out          # output file name 
    #BSUB -e gsi.%J.err          # error file name 
    #BSUB -q regular             # queue
    

  • Still in case set up setion, we didn't change the setups in the last three variables:
    • bk_core=ARW
    • bkcv_option=NAM
    • if_clean=clean
  • Search Section # Link to the radiance data, comment out the links to the BUFR file under this section:
  • # ln -s ${OBS_ROOT}/gdas1.t${HH}z.1bamua.tm00.bufr_d amsuabufr
    # ln -s ${OBS_ROOT}/gdas1.t${HH}z.1bamub.tm00.bufr_d amsubbufr
    # ln -s ${OBS_ROOT}/gdas1.t${HH}z.1bhrs4.tm00.bufr_d hirs4bufr
    # ln -s ${OBS_ROOT}/gdas1.t${HH}z.1bmhs.tm00.bufr_d mhsbufr
    # ln -s ${OBS_ROOT}/gdas1.t${HH}z.gpsro.tm00.bufr_d gpsrobufr