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