MCIP_TO_HYSPLIT - 30 September 2006


 This new software tool has been developed to link the Eulerian-based 
 Community Multiscale Air Quality (CMAQ) modeling system with the 
 Lagrangian-based HYSPLIT (HYbrid Single-Particle Lagrangian Integrated 
 Trajectory) model.  Both models require many of the same hourly 
 meteorological parameter fields in their applications, however, their 
 meteorological input data files exhibit different formats and data 
 structures.  MCIP2ARL reads four different MCIP data files to retrieve the 
 2-D and 3-D meteorological parameter fields needed by HYSPLIT and it 
 generates a packed-binary (in ARL format) meteorological data file for 
 direct input into the HYSPLIT modeling system. The MCIP2ARL interface 
 program source code is in Fortran free form and executables have been 
 successfully created with the Fortran 90 compiler on various computer 
 platforms.  The Input/Output Applications Programming Interface (I/O-API) 
 and netCDF libraries are also required to successfully compile this 
 program.

 Thus, the HYSPLIT trajectory model or other modeling components, can now 
 be driven with the same wind fields applied in CMAQ simulations. The 
 linkage of CMAQ and HYSPLIT through the meteorology permits innovative 
 analyses of various CMAQ output files to be performed within the HYSPLIT 
 Lagrangian trajectory framework.

 Further information about the MCIP2ARL program and an example application 
 are provided in an Extended Abstract by J. Godowitch and R. Draxler for 
 the 5th Annual CMAS conference on the web site: www.cmascenter.org. 
 Another utility program,which is under development, converts trajectory 
 coordinates into grid cell indices (row, column, layer) and retrieves user 
 selected species from a CMAQ 3-D concentration file.  The example 
 application in the paper shows CMAQ species concentration differences 
 along HYSPLIT forward and backward trajectory paths in downwind and upwind 
 directions, respectively, from selected locations based on model 
 simulation results with various emission scenarios.

 A sample Unix run script to apply MCIP2ARL is provided below. Although 
 MCIP files typically contain the hourly fields for a single day, by 
 specifying the particular "days" to be processed, multiple days could be 
 processed during a run and each output file ready for HYSPLIT also 
 contains the hourly meteorological files for a single day. The MCIP files 
 required to be read by MCIP2ARL are; GRIDDOT2D, METCRO2D, METCRO3D, and 
 METDOT3D.

      #!
      # C-shell script to run MCIP2ARL
      date
      set timestamp
      set echo
       setenv IOAPI_LOG_WRITE F
      # These are the parameters for this run
      set BASE     = /home/...
      set XBASE    = /home/...
      set EXEC     = MCIP2ARL

      # Define various directory paths.
       set METDIR   = /home/.../mcipfiles
       set OUTDIR   = /home/.../output

      #  Specify the year of the MCIP data files.
      set YR = 95

      ##### for a MONTH to be processed
      foreach MONTH (07)
      # time step parameters through the days
      if ( $MONTH == 07 ) then
       set days  = "7"
        else
       echo " Month is not correct "
       exit (1)
      endif

       foreach day ($days)
       echo $day
       set STDAT  = `echo $day | awk '{printf("%2.2d",$1)}'`

      cd $BASE; set timestamp; echo " "; set echo
      #
      ##############################################################
      #                  Set up for the MCIP2ARL test run          #
      ##############################################################
      #  Define the date for the data files.
       set extn = ${YR}${MONTH}${STDAT}

      # Define the output file.
       set OUT = MCIP.ARL${extn}
       setenv OUTFILE $OUTDIR/$OUT
      if (-e $OUTDIR/$OUT) rm $OUTDIR/$OUT
      #
      ln -s  $METDIR/GRIDDOT2D_${extn} GRIDDOT2D_${extn}
      ln -s  $METDIR/METCRO2D_${extn}  METCRO2D_${extn}
      ln -s  $METDIR/METCRO3D_${extn}  METCRO3D_${extn}
      ln -s  $METDIR/METDOT3D_${extn}  METDOT3D_${extn}
      #
       if ($status != 0) exit (2)

       ls -l $XBASE/$EXEC; size $XBASE/$EXEC

       $XBASE/$EXEC ${extn} ${OUTFILE}

       end
       end