#! /bin/csh -f

# RCS file, release, date & time of last delta, author, state, [and locker]
# $Header$ 

# what(1) key, module and SID; SCCS file; date and time of last delta:
# %W% %P% %G% %U%

#> method: run.icon >&! icon_d1b.log &

### setenv M3DATA /project/air4/yoj/rel/data

 if ( ! -e $M3DATA ) then
    echo "   $M3DATA path does not exist"
    exit 1
    endif
 echo " "; echo " Input data path, M3DATA set to $M3DATA"; echo " "

 set APPL     = us36
#set CFG      = d1b
 set CFG      = $APPL
 set EXEC     = ICON_${CFG}
 setenv NPCOL_NPROW "1 1"

#set BASE     = /tmp/you/your_directory
 set BASE     = $cwd

 cd $BASE; date; set timestamp; cat $BASE/cfg.${CFG}; echo " "; set echo

#> horizontal grid defn
#setenv GRIDDESC /project/cmaq/yoj/GRIDDESC1
 setenv GRIDDESC ../GRIDDESC1
 setenv GRID_NAME US36_132X90
#setenv GRID_NAME M_08_99TUT02

#> set the path for the output file containing the IC's
 set OUTDIR = $M3DATA/icon
 if ( ! -d "$OUTDIR" ) mkdir -p $OUTDIR

##############################################################################
# This script to run the ICON preprocessor has two major sections;           #
# 1) use default profile inputs, or 2) use CMAQ CTM concentration files.     #
# Follow these steps:                                                        #
#   1) Set IC equal to "profile" or "m3conc"                                 #
#   2) set the remainder of the environment variables for the section being  #
#      used (see below)                                                      #
##############################################################################

 set IC = profile      # use default profile file
#set IC = m3conc       # use CMAQ CTM concentration files (nested runs)

#> date or profile
 if ( $IC == profile ) then
    setenv INIT_CONC_1    "$OUTDIR/ICON_${APPL}_${GRID_NAME}_profile -v"     
    endif
 if ( $IC == m3conc ) then 
    set DATE = 1999183  # July 02, 1999                                       
    setenv INIT_CONC_1    "$OUTDIR/ICON_${APPL}_${GRID_NAME}_${DATE} -v"     
    endif  

##############################################################################
# Default profile environment variables                                      #
#----------------------------------------------------------------------------#
# IC_PROFILE     = the file containing the default IC profiles               #
#                                                                            #
# MECH_CONV_FL   = file containing mechanism conversion rules (optional) -   #
#                  must be set if a user-defined mechanism conversion file   #
#                  is being used                                             #
#                                                                            #
# NOTE: SDATE and STIME do not need to be set                                #
##############################################################################

 if ( $IC == profile ) then

    setenv IC_PROFILE      $M3DATA/raw/icon/ic_profile_v6.dat

#   setenv MECH_CONV_FILE  /home/ggb/models3/icbc/mech_conv_fl.v1

    endif
 
##############################################################################
# CMAQ CTM concentration file environment variables (e.g., for nesting )     #
#----------------------------------------------------------------------------#
# CTM_CONC_1     = the CTM concentration file for the coarse domain          #
#                                                                            #
# ICON_PING      = flag controlling PinG IC's creation - must be set to YES  #
#                  if PinG is included, otherwise set ICON_PING to NO        #
#                                                                            #
# CTM_PING_1     = the PinG concentration file for the coarse domain         #
#                                                                            #
# PING_PDM_1     = the PinG plume dynamics file for the coarse domain        #
#                                                                            #
# MET_CRO_3D_CRS = the MET_CRO_3D met file for the coarse domain - must be   #
#                  set if PinG is included or if the vertical grid type is   #
#                  changed between nests                                     #
#                                                                            #
# MET_CRO_3D_FIN = the MET_CRO_3D met file for the inner, nested, domain -   #
#                  must be set if the vertical grid type is changed between  #
#                  nests                                                     #
#                                                                            #
# MECH_CONV_FL   = file containing mechanism conversion rules (optional) -   #
#                  must be set if a user-defined mechanism conversion file   #
#                  is being used                                             #
#                                                                            #
# NOTE: SDATE (yyyyddd) and STIME (hhmmss) must always be set                #
##############################################################################
 
 if ( $IC == m3conc ) then 

    setenv CTM_CONC_1 $M3DATA/cctm/CCTM_d1bCONC.d1b

    setenv ICON_PING

    setenv CTM_PING_1
    setenv PING_PDM_1
    setenv MET_CRO_3D_CRS

    setenv MET_CRO_3D_FIN

    setenv MECH_CONV_FILE

    setenv SDATE           ${DATE}
    setenv STIME           000000

    endif

#/opt/SUNWspro/bin/workshop -D "$BASE/$EXEC"; exit
#/usr/pgi/linux86/bin/Xpgdbg "$BASE/$EXEC"; exit

 /usr/bin/time $BASE/$EXEC

 exit()