Skip navigation
NASA Logo, National Aeronautics and Space Administration
Currently Being Moderated

Gocart2wrf Case Study

VERSION 6  Click to view document history
Created on: Jul 9, 2010 9:46 AM by Philip Hayes - Last Modified:  Sep 1, 2010 11:12 AM by Rahman Syed

INTRODUCTION

 

The purpose of this document is to describe the process of building and running the Gocart2wrf case study (#13) under the NU-WRF project.


This project is designed to extract chemical species data produced by the stand-alone Gocart model and append it to the WRF input and boundary condition files.  A case study was chosen that had high concentrations of aerosols throughout the domain in order to determine that the utility is properly working and WRF is reading in the data correctly.  The simulation contains one domain, located over the Eastern US on January 01, 2006. 


For more information on the Gocart2wrf utility, please check out the README file included with the code or visit the following link:

https://modelingguru.nasa.gov/clearspace/docs/DOC-1830

 

 

Running Case Study

 

The data that is needed to run Case 13 is located in:

 

                           /discover/nobackup/projects/nu-wrf/cases/13/usecase

 

First set up the environment variables for the workdir, WPS directory, WRFV3 directory, and the directory where your Emissions data set will be located:

For csh:
setenv workDir path_to_myWorkDir
setenv wrfDir $workDir/WRFV3
setenv wpsDir $workDir/WPS
setenv emissionDir path_to_emissions_data
setenv gocart2wrfDir path_to_gocart2wrf

For bash:
export workDir=path_to_myWorkDir
export wrfDir=$workDir/WRFV3
export wpsDir=$workDir/WPS
export emissionDir=path_to_emissions_data
export gocart2wrfDir=path_to_gocart2wrf

 

Make sure that when WRF is compiled that WRF-Chem is turned on.

 

First load the modules:

module purge
module load intel/11.1 mpi/openmpi-1.4

 

Step 1: Generate the input fields

 

* Using the GFS data, namelist.wps, and namelist.input in the case13 usecase directory, generate the wrfinput_d01 and wrfbdy_d01 for the case study. (Make sure that you change the path of the geography files in namelist.wps!)

 

Step 2: Set up Emissions Data set

 

* Change to the directory where your Emissions Data will be located:

cd emissionDir

 

* Two tar files are needed to generate the necessary emissions files used by WRF-Chem.  They are located at:

ftp://aftp.fsl.noaa.gov/divisions/taq/global_emissions/

 

* Download the files global_emissions_v3_04aug09.tar.gz and prep_sources_chem_13apr10.tar.gz and place them into $emissionDir.  *Note* These files may be updated and in that case grab the most recent versions.

 

* Untar/Unzip the two packages

tar xvfz global_emissions_v3_04aug09.tar.gz
tar xvfz prep_sources_chem_13apr10.tar.gz

 

* Go into the utils/bin directory

cd Prep_sources_chem_cptec_wrf/utils/bin

 

* Build the utility

make -f Make_utils

 

* Change directories to:

cd ../../bin

 

* Build the prep_chem_sources utility

make

 

*** You will get error message once the make is complete! ***  You must type make again to link all the libraries.

make

 

* Now edit the prep_chem_sources.inp.  The following variables must be changed:

 

            -     grid_type= 'lambert'

            -     ihour=00

            -     iday=01

            -     imon=01

            -     iyear=2006

            -     use_retro=0

            -     use_edgar=0

            -     gocart_data_dir='place full path to emissions data here'

                                  i.e., $emissionDir/Global_emissions_v3/Emission_data/GOCART/emissions

            -     use_bbem=0

            -     use_bbem_plumerise=0

            -     use_gocart_bg=1

            -     gocart_bg_data_dir='place full path to emissions data here'

                                  i.e., $emissionDir/Global_emissions_v3/Emission_data/GOCART

            -     NGRIDS = 1,

            -     NNXP     = 149,

            -     NNYP     = 119,

            -     NXTNEST = 0,

            -     DELTAAX = 18000.,

            -     DELTAAY = 18000.,

            -     NSTRATX = 1,

            -     NSTRATY = 1,

            -     POLELAT = 90.,

            -     POLELON = 0.,

            -     STDLAT1 = 37.,

            -     STDLAT2 = 37.,

            -     CENTLAT = 37.,

            -     CENTLON = -86.0,

 

 

* Change to main directory and copy prep_chem_sources.inp to that directory

cd ..
cp bin/prep_chem_sources.inp .

 

* Execute the prep_chem_sources in the bin directory

cd bin
./prep_chem_sources

 

* You should end up with a set of Emissions files:

FIRE-GOCART-T-2006-01-01-000000-g1-ab.bin
FIRE-GOCART-T-2006-01-01-000000-g1-bb.bin
FIRE-GOCART-T-2006-01-01-000000-g1-gocartBG.bin

 

Step 3: Implement Emission Data Set

 

* Move into the following directory:

cd $workDir/test/em_real

 

* Link the emissions files

ln -s $emissionDir/Prep_sources_chem_cptec_wrf/bin/FIRE-GOCART-T-2006-01-01-000000-g1-ab.bin emissopt3_d01
ln -s
$emissionDir/Prep_sources_chem_cptec_wrf/bin/FIRE-GOCART-T-2006-01-01-000000-g1-bb.bin emissfire_d01
ln -s $emissionDir/Prep_sources_chem_cptec_wrf/bin/FIRE-GOCART-T-2006-01-01-000000-g1-gocartBG.bin wrf_gocart_backg

 

* Move into WRF-Chem directory

cd $wrfDir/chem

 

* Edit convert_emiss.F so that each "ide" is "ide-1" and each "jde" is "jde-1".  There are many instances, so keeping a copy once it has been accomplished once is recommended.

 

* Move to WRF directory

cd $wrfDir/

 

* Compile convert_emiss program:

./compile emi_conv

 

* Edit namelist.input for the following variables:

 interval_seconds                    = 3600
frames_per_outfile                  = 1, 1000, 1000,
frames_per_emissfile                = 24,
io_style_emissions                  = 2,

 

* Link the convert_emiss executable:

ln -s chem/convert_emiss.exe .

 

* Execute convert_emiss:

./convert_emiss.exe

 

* Should produce two files:

              -  wrfchemi_d01

              -  wrfchemi_gocart_bg_d01

 

* Rename wrfchemi_d01 to wrfchemi_d01_2006-01-01_00:00:00

mv wrfchemi_d01 wrfchemi_d01_2006-01-01_00:00:00

 

* Change the namelist.input back to it's original form

 interval_seconds                    = 10800
frames_per_outfile                  = 1, 1000, 1000,
frames_per_emissfile                = 24,
io_style_emissions                  = 2,

 

Step 4: Run Gocart2wrf utility

 

* Move to Gocart2wrf directory

cd $gocart2wrfDir

 

* Copy over the Gocart data files to the Gocart2wrf directory

cp /discover/nobackup/projects/nu-wrf/cases/13/usecase/g4p0e530q0t_* . 

 

* Edit the gocart2wrf_namelist

 

* Build the Gocart2wrf utility

make

 

* Execute Gocart2wrf

./gocart2wrf

 

Step 5: Run the WRF model

 

* Submit wrf.exe

 

 

Example Output

 

The animation of hydrophyllic black carbon through the lifecycle of the simulation is attatched to the document.

Attachments:
Comments (0)

Bookmarked By (0)

More Like This

  • Retrieving data ...
USAGov logo NASA Logo - nasa.gov