################################################################################ #### UNIX Script Documentation Block # . . # Script name: cfs_chgres.sh # Script description: Changes resolution of cfs restart files. # # Author: Mark Iredell Org: NP23 Date: 1999-03-01 # # Abstract: This script changes the resolution of the cfs restart files, # namely the sigma file or the surface file or both. The resolution # of the output files is given in the argument list or as imported # environment variables. The resolution of the input files are taken # from the header records in the respective files. Resolution is # given as spectral truncation, number of levels, number of longitudes # and number of latitudes. The names of the input and output restart files # are also given in the argument list or as imported environment variables. # Other control variables may be specified in the imported environment, too. # For example, an alternate orography and sigma level structure may be # specified for the output sigma file, while an alternate sea-land mask # may be specified for the surface file. # # Script history log: # 1999-03-01 Mark Iredell # 2005-08-01 Hua-Lu Pan, Julia Zhu # # Usage: cfs_chgres.sh JCAP LEVS LONB LATB SIGINP SFCINP SIGOUT SFCOUT # # Input script positional parameters: # 1 Input sigma file # defaults to $SIGINP; skip sigma conversion if missing # 2 Input surface file # defaults to $SFCINP; skip surface conversion if missing # 3 Output sigma file # defaults to $SIGOUT, then to sigout # 4 Output sigma file # defaults to $SFCOUT, then to sfcout # 5 New spectral truncation # defaults to $JCAP; one or the other is required # 6 New number of levels # defaults to $LEVS; one or the other is required # 7 New number of longitudes # defaults to $LONB; one or the other is required # 8 New number of latitudes # defaults to $LATB; one or the other is required # # Imported Shell Variables: # SIGINP Input sigma file # overridden by $1; skip sigma conversion if missing # SFCINP Input surface file # overridden by $2; skip surface conversion if missing # SIGOUT Output sigma file # overridden by $3; defaults to sigout # SFCOUT Output surface file # overridden by $4; defaults to sfcout # JCAP New spectral truncation # overridden by $5; one or the other is required # LEVS New number of levels # overridden by $6; one or the other is required # LONB New number of longitudes # overridden by $7; one or the other is required # LATB New number of latitudes # overridden by $8; one or the other is required # NTRAC New number of tracers # defaults to input sigma file value # IDVC New vertical coordinate id (1 for sigma, 2 for hybrid) # defaults to input sigma file value # IDSL New midlayer pressure id (1 for phillips, 2 for mean) # defaults to input sigma file value # FIXGLOBAL Directory for global fixed files # defaults to /nwprod/fix # EXECGLOBAL Directory for global executables # defaults to /nwprod/exec # DATA working directory # (if nonexistent will be made, used and deleted) # defaults to current working directory # XC Suffix to add to executables # defaults to none # OROGRAPHY New orography ("NULL" to use from input sigma file) # defaults to ${FIXGLOBAL}/global_orography.t${JCAP}.grb # SIGLEVEL New sigma levels ("NULL" to use from input sigma file) # defaults to ${FIXGLOBAL}/global_siglevel.l${LEVS}.txt # O3CLIM New ozone climatology (required if adding ozone) # defaults to ${FIXGLOBAL}/global_o3clim.txt # SLMASK New sea-land mask ("NULL" to use from input surface file) # defaults to ${FIXGLOBAL}/global_slmask.t${JCAP}.grb # LONSPERLAT New lonsperlat ("NULL" to use from input surface file) # defaults to ${FIXGLOBAL}/global_lonsperlat.t${JCAP}.txt # CHGRESEXEC Change resolution executable # defaults to ${EXECGLOBAL}/cfs_chgres # INISCRIPT Preprocessing script # defaults to none # LOGSCRIPT Log posting script # defaults to none # ERRSCRIPT Error processing script # defaults to 'eval [[ $err = 0 ]]' # ENDSCRIPT Postprocessing script # defaults to none # CHGRESVARS Other namelist inputs to the change resolution executable # such as LONBI,LATBI,IGEN,MGG,JC,MQUICK # defaults to none set # NTHREADS Number of threads # defaults to 1 # NTHSTACK Size of stack per thread # defaults to 64000000 # PGMOUT Executable standard output # defaults to $pgmout, then to '&1' # PGMERR Executable standard error # defaults to $pgmerr, then to '&1' # pgmout Executable standard output default # pgmerr Executable standard error default # REDOUT standard output redirect ('1>' or '1>>') # defaults to '1>', or to '1>>' to append if $PGMOUT is a file # REDERR standard error redirect ('2>' or '2>>') # defaults to '2>', or to '2>>' to append if $PGMERR is a file # VERBOSE Verbose flag (YES or NO) # defaults to NO # # Exported Shell Variables: # PGM Current program name # pgm # ERR Last return code # err # # Modules and files referenced: # scripts : $INISCRIPT # $LOGSCRIPT # $ERRSCRIPT # $ENDSCRIPT # # programs : $CHGRESEXEC # # input data : $1 or $SIGINP # $2 or $SFCINP # $OROGRAPHY # $SIGLEVEL # $O3CLIM # $SLMASK # $LONSPERLAT # # output data: $3 or $SIGOUT # $4 or $SFCOUT # $PGMOUT # $PGMERR # # scratch : ${DATA}/chgres.inp.sig # ${DATA}/chgres.inp.orogb # ${DATA}/chgres.inp.siglevel # ${DATA}/chgres.inp.sfc # ${DATA}/chgres.inp.slmgb # ${DATA}/chgres.out.sig # ${DATA}/chgres.out.sfc # ${DATA}/NULL # # Remarks: # # Condition codes # 0 - no problem encountered # >0 - some problem encountered # # Control variable resolution priority # 1 Command line argument. # 2 Environment variable. # 3 Inline default. # # Attributes: # Language: POSIX shell # Machine: IBM SP # #### ################################################################################ # Set environment. export VERBOSE=${VERBOSE:-"NO"} if [[ "$VERBOSE" = "YES" ]] then echo $(date) EXECUTING $0 $* >&2 set -x fi # Command line arguments. export SIGINP=${1:-${SIGINP:-NULL}} export SFCINP=${2:-${SFCINP:-NULL}} export SIGOUT=${3:-${SIGOUT:-sigout}} export SFCOUT=${4:-${SFCOUT:-sfcout}} export JCAP=${5:-${JCAP:?}} export LEVS=${6:-${LEVS:?}} export LONB=${7:-${LONB:?}} export LATB=${8:-${LATB:?}} # Directories. export FIXGLOBAL=${FIXGLOBAL:-/nwprod/fix} export EXECGLOBAL=${EXECGLOBAL:-/nwprod/exec} export DATA=${DATA:-$(pwd)} # Filenames. export XC=${XC} export CHGRESEXEC=${CHGRESEXEC:-${EXECGLOBAL}/cfs_chgres$XC} export OROGRAPHY=${OROGRAPHY:-${FIXGLOBAL}/cfs_orography.t${JCAP}.grb} export SIGLEVEL=${SIGLEVEL:-${FIXGLOBAL}/cfs_siglevel.l${LEVS}.txt} export O3CLIM=${O3CLIM:-${FIXGLOBAL}/cfs_o3clim.txt} export SLMASK=${SLMASK:-${FIXGLOBAL}/cfs_slmask.t${JCAP}.grb} export LONSPERLAT=${LONSPERLAT:-${FIXGLOBAL}/global_lonsperlat.t${JCAP}.txt} export INISCRIPT=${INISCRIPT} export ERRSCRIPT=${ERRSCRIPT:-'eval [[ $err = 0 ]]'} export LOGSCRIPT=${LOGSCRIPT} export ENDSCRIPT=${ENDSCRIPT} # Other variables. export NTRAC=${NTRAC:-0} export IDVC=${IDVC:-0} export IDSL=${IDSL:-0} export CHGRESVARS=${CHGRESVARS} export NTHREADS=${NTHREADS:-1} export NTHSTACK=${NTHSTACK:-64000000} export XLSMPOPTS=${XLSMPOPTS:-"parthds=$NTHREADS:stack=NTHSTACK"} export PGMOUT=${PGMOUT:-${pgmout:-'&1'}} export PGMERR=${PGMERR:-${pgmerr:-'&2'}} typeset -L1 l=$PGMOUT [[ $l = '&' ]]&&a=''||a='>' export REDOUT=${REDOUT:-'1>'$a} typeset -L1 l=$PGMERR [[ $l = '&' ]]&&a=''||a='>' export REDERR=${REDERR:-'2>'$a} ################################################################################ # Preprocessing $INISCRIPT pwd=$(pwd) if [[ -d $DATA ]] then mkdata=NO else mkdir -p $DATA mkdata=YES fi cd $DATA||exit 99 ################################################################################ # Change resolution export XLSMPOPTS="parthds=$NTHREADS:stack=$NTHSTACK" export PGM=$CHGRESEXEC export pgm=$PGM $LOGSCRIPT rm -f NULL ln -sf $SIGINP chgres.inp.sig ln -sf $OROGRAPHY chgres.inp.orogb ln -sf $SIGLEVEL chgres.inp.siglevel ln -sf $O3CLIM chgres.inp.o3clim ln -sf $SFCINP chgres.inp.sfc ln -sf $SLMASK chgres.inp.slmgb ln -sf $LONSPERLAT chgres.inp.lonsperlat ln -sf $SIGOUT chgres.out.sig ln -sf $SFCOUT chgres.out.sfc eval $CHGRESEXEC <&2 fi exit $err