Difference between revisions of "FlexAOD"

From Geos-chem
Jump to: navigation, search
Line 1: Line 1:
[http://pumpkin.aquila.infn.it/flexaod/ FlexAOD] (Flexible Aerosol Optical Depth) is a post-processing tool of GEOS-Chem, aimed at extending the model on-line capabilities concerning calculations of [[Aerosol optical properties]].
+
[http://pumpkin.aquila.infn.it/flexaod/ FlexAOD] (Flexible Aerosol Optical Depth) is a post-processing tool of GEOS-Chem, aimed at extending the model on-line capabilities concerning calculations of [[Aerosol optical properties]]. For a brief overview, check out the [http://acmg.seas.harvard.edu/geos/meetings/2011/pdf/Day1/MonD_AerClim_curci_gabriele_1_pc.pdf presentation] at the [http://acmg.seas.harvard.edu/geos/geos_meeting_2011.html 5th International GEOS-Chem Meeting]!
  
 
'''THIS PAGE IS STILL UNDER CONTRUCTION! PLEASE COME BACK LATER!'''
 
'''THIS PAGE IS STILL UNDER CONTRUCTION! PLEASE COME BACK LATER!'''
Line 5: Line 5:
 
== User's Guide ==
 
== User's Guide ==
  
=== Download code & data ===
+
 
Source code, input & test data are at this link:
+
=== Download code ===
http://pumpkin.aquila.infn.it/flexaod/download
+
 
 +
Source code & input data are available at this link:
 +
http://pumpkin.aquila.infn.it/flexaod
  
 
From a Linux/Unix terminal, you may follow this procedure:
 
From a Linux/Unix terminal, you may follow this procedure:
  
  mkdir -p geos-chem/flexaod
+
  mkdir -p ~/geos-chem/flexaod
  cd geos-chem/flexaod
+
  cd ~/geos-chem/flexaod
 
  wget http://pumpkin.aquila.infn.it/flexaod/download/flexaod_v1.0.tgz
 
  wget http://pumpkin.aquila.infn.it/flexaod/download/flexaod_v1.0.tgz
 
  tar zxvf flexaod_v1.0.tgz
 
  tar zxvf flexaod_v1.0.tgz
Line 66: Line 68:
  
 
=== Set up GEOS-Chem timeseries & data dir ===
 
=== Set up GEOS-Chem timeseries & data dir ===
 +
 +
==== GEOS-Chem tracers ====
  
 
FlexAOD needs GEOS-Chem timeseries in input (instantaneous (ND49) or satellite (ND51)). There is a minimum set of tracers to be saved in output timeseries:
 
FlexAOD needs GEOS-Chem timeseries in input (instantaneous (ND49) or satellite (ND51)). There is a minimum set of tracers to be saved in output timeseries:
 +
  
 
{| border=1 cellspacing=0 cellpadding=5
 
{| border=1 cellspacing=0 cellpadding=5
 
|- bgcolor="#CCCCCC"
 
|- bgcolor="#CCCCCC"
 
!Aerosol concentrations
 
!Aerosol concentrations
!Tracer #              
+
!Tracer #<br>(v8-02-02*)       
 
!Description
 
!Description
 
|-  
 
|-  
Line 165: Line 170:
 
|-
 
|-
 
|}
 
|}
 +
<nowiki>*</nowiki>NOTE: tracer # valid for v8-02-02, please check for next versions!
  
 +
==== GEOS-Chem rundir ====
  
 
By default, FlexAOD will look for GEOS-Chem output into a subdirectory defined at line 168 of <tt>flexaod.sh</tt>:
 
By default, FlexAOD will look for GEOS-Chem output into a subdirectory defined at line 168 of <tt>flexaod.sh</tt>:
Line 171: Line 178:
 
  rundir=${HOME}/geos-chem/GEOS-Chem.${ver}.stdrun/runs/${met}.${res}/output/${lab}
 
  rundir=${HOME}/geos-chem/GEOS-Chem.${ver}.stdrun/runs/${met}.${res}/output/${lab}
  
where $ver, $met, $res, and $lab are all defined in [[FlexAOD#input.dat|input.dat]]. A GEOS-Chem user is familiar with the first three variables, the latter is introduced for further organization of the run directory, but its use is optional. Tailor the naming conventions of the run directory to your needs.
+
where <tt>$ver</tt>, <tt>$met</tt>, <tt>$res</tt>, and <tt>$lab</tt> are all defined in <tt>[[FlexAOD#input.dat|input.dat]]</tt>. A GEOS-Chem user is familiar with the first three variables, the latter is introduced for further organization of the run directory, but its use is optional. Tailor the naming conventions of the run directory to your needs.
  
In addition to time series, FlexAOD expects to find under rundir:
+
In addition to time series, FlexAOD expects to find under <tt>rundir</tt>:
  
 
  tracerinfo.dat.${tlab}
 
  tracerinfo.dat.${tlab}
 
  diaginfo.dat.${tlab}
 
  diaginfo.dat.${tlab}
  
where $tlab is an additional "time label" also defined in [[FlexAOD#input.dat|input.dat]].
+
where <tt>$tlab</tt> is an additional "time label" also defined in <tt>[[FlexAOD#input.dat|input.dat]]</tt>. Please rename your <tt>*info.dat</tt> accordingly.
  
Sample GEOS-Chem timeseries output and info files are provided for testing. Go under your "rundir" and download data:
+
Sample GEOS-Chem timeseries output and info files are provided for testing under: http://pumpkin.aquila.infn.it/flexaod/download/testdata
 +
 
 +
From shell, go under your "rundir" and you may download all data with:
 
  wget -r -nd -np -R*html* http://pumpkin.aquila.infn.it/flexaod/download/testdata
 
  wget -r -nd -np -R*html* http://pumpkin.aquila.infn.it/flexaod/download/testdata
 +
 +
 +
==== Optional GEOS-Chem launch script ====
 +
 +
A launch script for GEOS-Chem model is provided with FlexAOD in order to:
 +
*launch the simulation within your <tt>rundir</tt>
 +
*move/rename the output files after run in an output directory suitable for FlexAOD (see [[FlexAOD#GEOS-Chem rundir|GEOS-Chem rundir]])
 +
You may download and launch the script with:
 +
wget http://pumpkin.aquila.infn.it/flexaod/download/run_geos.sh
 +
./run_geos.sh <lab> <tlab>
 +
where <tt><lab></tt> and <tt><tlab></tt> are the simulation and time labels to subsequently set into <tt>[[FlexAOD#input.dat|input.dat]]</tt> for use with FlexAOD.
  
  
Line 314: Line 334:
 
  ./flexaod.sh &> log &
 
  ./flexaod.sh &> log &
 
  tail -f log
 
  tail -f log
 +
 +
The utility <tt>loop.sh</tt> allows automatic FlexAOD launch for several user selected wavelengths and time periods. The <tt>loop_wave_rtm.sh</tt> is for specific use with <tt>point_rtm=T</tt> option: it will paste in a unique file the output of several wavelengths.
 +
 +
 +
=== FlexAOD output ===
 +
 +
Similarly to [[Aerosol optical properties|GEOS-Chem on-line AOD calculations]], FlexAOD performs calculations for 6 "optical" species, assumed in external mixing state:
 +
 +
{| border=1 cellspacing=0 cellpadding=5
 +
|- bgcolor="#CCCCCC"
 +
!"Optical" species
 +
![[FlexAOD#GEOS-Chem tracers|GEOS-Chem tracers]] included
 +
|-
 +
|SULF
 +
|SO4, NH4, NIT
 +
|-
 +
|OC
 +
|OCPI, OCPO, SOA1-4
 +
|-
 +
|BC
 +
|BCPI, BCPO
 +
|-
 +
|SSa
 +
|SSa
 +
|-
 +
|SSc
 +
|SSc
 +
|-
 +
|DUST
 +
|DST1-4
 +
|-
 +
|}
 +
 +
 +
FlexAOD has the following types of user selected output:
 +
 +
{| border=1 cellspacing=0 cellpadding=5
 +
|- bgcolor="#CCCCCC"
 +
!Description
 +
!<tt>[[FlexAOD#input.dat|input.dat]]</tt> option
 +
!default output dir
 +
!sample file name
 +
!output variables
 +
|-
 +
|Full 3-D gridded output
 +
|<tt>gridout=T</tt>
 +
|<tt>[[FlexAOD#GEOS-Chem rundir|<rundir>]]/../GrADS</tt>
 +
|<tt>ts_satellite_<tlab>_<lab>_<dini>-<dend>_<wl>nm_<aodlab>.ctl/dat</tt>
 +
|AOD of "optical" species (same as GEOS-Chem on-line output)
 +
[[FlexAOD#Gridded (GrADS) output|[details]]]
 +
|-
 +
|Profile at selected points
 +
|<tt>point_prof=T</tt><br><tt>point_file=<filename></tt>
 +
|<tt>[[FlexAOD#GEOS-Chem rundir|<rundir>]]/../PROFILES</tt>
 +
|<tt>ts_satellite_<tlab>_<lab>_<dini>-<dend>_<wl>nm_<aodlab>_000012</tt>
 +
|Profile of mass concentrations, number concentration; extinction, scattering and backscattering coefficients; single-scattering albedo; asymmetry factor; phase function coefficients
 +
[[FlexAOD#Profile output|[details]]]
 +
|-
 +
|Profile at selected points for use with RTM (libRadtran)
 +
|<tt>point_rtm=T</tt><br><tt>point_file=<filename></tt>
 +
|<tt>[[FlexAOD#GEOS-Chem rundir|<rundir>]]/../libRadtran</tt>
 +
|<tt>ts_satellite_<tlab>_<lab>_<dini>-<dend>_<wl>nm_<aodlab>_000012_FILES</tt> (list of files with model level info)<br>
 +
<tt>ts_satellite_<tlab>_<lab>_<dini>-<dend>_<wl>nm_<aodlab>_CONC</tt> ("optical" species mass column, g/m^2)<br>
 +
|extinction coefficient, single-scattering albedo, phase function coefficients
 +
|-
 +
|}
 +
 +
Any number of output types may be requested (at least one!).
 +
 +
 +
==== Gridded (GrADS) output ====
 +
 +
The <tt>.ctl/.dat</tt> pairs generated for gridded GrADS binary output may be used for direct visualization and analysis in [http://www.iges.org/grads/ GrADS].
 +
 +
If you don't like GrADS, you may read data with other computer languages, such as Fortran:
 +
 +
program read_grads
 +
 +
  ! Parameters
 +
  integer,paramter :: iou = 10  ! I/O unit
 +
  integer,paramter :: ni = 144  ! n. lons
 +
  integer,paramter :: nj = 91    ! n. lats
 +
  integer,paramter :: nl = 38    ! n. levs
 +
  integer,paramter :: nspecs = 6 ! n. specs
 +
  integer,paramter :: ntimes = 5 ! n. times
 +
 +
  ! Arrays
 +
  real, dimension(ni,nj,nl) :: aottot, boxh
 +
  real, dimension(ni,nj,nl,nspecs) :: aot
 +
 +
 +
  ! Open file in "binary" mode
 +
  open(iou, file="<filename>", status="old", form="binary")
 +
 +
  ! Loop over times
 +
  do it = 1, ntimes
 +
 +
    ! AOT
 +
    read(iou) aottot
 +
 +
    ! Optical Depths of single species
 +
    do k=1,nspecs+1
 +
      do l=1,nl
 +
        read(iou) aot(:,:,l,k)
 +
      enddo
 +
    enddo
 +
 +
    ! Box Height (m)
 +
    do l=1,nl
 +
      read(iou) boxh(:,:,l)
 +
    enddo
 +
 +
  enddo  ! times
 +
 +
  ! Close file
 +
  close(iou)
 +
 +
end program read_grads
 +
 +
If your Fortran compiler does not understand the "binary" form, then just change it to "unformatted" at line 547 of <tt>src/flexaod.f90</tt>. In that case, however, you will NOT be able to use the output with GrADS.
 +
 +
 +
==== Profile output ====
 +
 +
Profile output (<tt>point_prof=T</tt>) is much easier to handle, because it is in ASCII format. One file per point specified in <tt>[[FlexAOD#pointout.dat|pointout.dat]]</tt> list is generated. Each file thus contains the calculated profile for a given time/place.
 +
 +
Each output file has <tt>nl+1</tt> rows, where <tt>nl</tt> is the number of levels in GEOS-Chem timeseries. Here is the content of rows:
 +
* '''1st row:''' column average/integrated values
 +
* '''2nd to (<tt>nl+1</tt>)th row:''' values at model levels
 +
 +
Let <tt>nspecs = 6</tt> the number of "optical" species and <tt>noutpc</tt> is the number of output coefficients (default is 16, set at line <tt>1431</tt> of <tt>src/flexaod.f90</tt>).
 +
 +
Output variables for '''1st row''' are:
 +
* '''column 1:''' 0.000
 +
* '''column 2 to <tt>nspecs+1</tt>:''' column concentrations of "optical" species (g/m^2)
 +
* '''column <tt>nspecs+2</tt>:''' column number concentration (#/cm^2)
 +
* '''column <tt>nspecs+3</tt>:''' aerosol optical depth (unitless)
 +
* '''column <tt>nspecs+4</tt>:''' scattering aerosol optical depth (unitless)
 +
* '''column <tt>nspecs+5</tt>:''' back-scattering aerosol optical depth (unitless)
 +
* '''column <tt>nspecs+6</tt>:''' column-weighted single-scattering albedo (unitless)
 +
* '''column <tt>nspecs+7</tt>:''' column-weighted asymmetry factor (unitless)
 +
* '''column <tt>nspecs+8</tt> to <tt>nspecs+8+noutpc-1</tt>:''' column-weighted phase function coefficients (unitless)
 +
 +
Output variables for '''other rows''' are:
 +
* '''column 1:''' GEOS-Chem level height (km)
 +
* '''column 2 to <tt>nspecs+1</tt>:''' mass concentrations of "optical" species (ug/m^3)
 +
* '''column <tt>nspecs+2</tt>:''' number concentration (#/cm^3)
 +
* '''column <tt>nspecs+3</tt>:''' extinction coefficient (1/km)
 +
* '''column <tt>nspecs+4</tt>:''' scattering coefficient (1/km)
 +
* '''column <tt>nspecs+5</tt>:''' back-scattering coefficient (1/km)
 +
* '''column <tt>nspecs+6</tt>:''' single-scattering albedo (unitless)
 +
* '''column <tt>nspecs+7</tt>:''' asymmetry factor (unitless)
 +
* '''column <tt>nspecs+8</tt> to <tt>nspecs+8+noutpc-1</tt>:''' phase function coefficients (unitless)
  
  

Revision as of 14:57, 25 May 2011

FlexAOD (Flexible Aerosol Optical Depth) is a post-processing tool of GEOS-Chem, aimed at extending the model on-line capabilities concerning calculations of Aerosol optical properties. For a brief overview, check out the presentation at the 5th International GEOS-Chem Meeting!

THIS PAGE IS STILL UNDER CONTRUCTION! PLEASE COME BACK LATER!

User's Guide

Download code

Source code & input data are available at this link: http://pumpkin.aquila.infn.it/flexaod

From a Linux/Unix terminal, you may follow this procedure:

mkdir -p ~/geos-chem/flexaod
cd ~/geos-chem/flexaod
wget http://pumpkin.aquila.infn.it/flexaod/download/flexaod_v1.0.tgz
tar zxvf flexaod_v1.0.tgz

This will create the subdirectory "flexaod_v1.0", here is a quick look at the content:

Filename/Dirname Description
flexaod.sh Launch script
input.dat Configuration file
pointout.dat Point output configuration file (the actual name is set with point_file variable)
random_pointout.sh Random pointout.dat generator
loop.sh Loop script for automatic launch at several wavelengths and/or date intervals
loop_wave_rtm.sh Loop script for automatic launch and post-processing of libRadtran output at several wavelengths
src/ Fortran source code
src/Makefile.hdr.* Makefile headers for Fortran compiler configuration
data/ Optical input data (i.e. all but GEOS-Chem timeseries)
data/gads/optdat GADS/OPAC input database (Hess et al., 1998)
data/lsmask Land/sea mask calculated from gridded national IDs (CIESIN and CIAT, 2005), only needed if point_watonly=T
data/size_dist Parameters of size distributions and hygroscopic growth factors
data/water Water spectral refractive index (Segelstein, 1981)


Set up GEOS-Chem timeseries & data dir

GEOS-Chem tracers

FlexAOD needs GEOS-Chem timeseries in input (instantaneous (ND49) or satellite (ND51)). There is a minimum set of tracers to be saved in output timeseries:


Aerosol concentrations Tracer #
(v8-02-02*)
Description
SO4 27 Sulfate
NH4 31 Ammonium
NIT 32 Nitrate
BCPI 34 Black carbon hydrophillic
OCPI 35 Primary organic carbon hydrophillic
BCPO 36 Black carbon hydrophobic
OCPO 37 Primary organic carbon hydrophobic
SOA1-4 45-48 Secondary organic aerosols classes 1-4
DST1-4 49-52 Mineral dust aerosols (size bins 1-4)
SSa 53 Sea salt aerosols accumulation mode
SSc 54 Sea salt aerosols coarse mode
Meteo fields Tracer # Description
AIRD 78 Air density (molecules/cm^3)
BOXH 93 Box Height (m)
RELH 94 Relative Humidity (%)
Optional (for comparison with on-line AOD) Tracer # Description
OPSO4 82 Sulfate aerosols AOD
OPBC 83 Black carbon AOD
OPOC 84 Organic carbon AOD
OPSSa 85 SSa AOD
OPSSc 86 SSc AOD
OPD 87 Total dust AOD

*NOTE: tracer # valid for v8-02-02, please check for next versions!

GEOS-Chem rundir

By default, FlexAOD will look for GEOS-Chem output into a subdirectory defined at line 168 of flexaod.sh:

rundir=${HOME}/geos-chem/GEOS-Chem.${ver}.stdrun/runs/${met}.${res}/output/${lab}

where $ver, $met, $res, and $lab are all defined in input.dat. A GEOS-Chem user is familiar with the first three variables, the latter is introduced for further organization of the run directory, but its use is optional. Tailor the naming conventions of the run directory to your needs.

In addition to time series, FlexAOD expects to find under rundir:

tracerinfo.dat.${tlab}
diaginfo.dat.${tlab}

where $tlab is an additional "time label" also defined in input.dat. Please rename your *info.dat accordingly.

Sample GEOS-Chem timeseries output and info files are provided for testing under: http://pumpkin.aquila.infn.it/flexaod/download/testdata

From shell, go under your "rundir" and you may download all data with:

wget -r -nd -np -R*html* http://pumpkin.aquila.infn.it/flexaod/download/testdata


Optional GEOS-Chem launch script

A launch script for GEOS-Chem model is provided with FlexAOD in order to:

  • launch the simulation within your rundir
  • move/rename the output files after run in an output directory suitable for FlexAOD (see GEOS-Chem rundir)

You may download and launch the script with:

wget http://pumpkin.aquila.infn.it/flexaod/download/run_geos.sh
./run_geos.sh <lab> <tlab>

where <lab> and <tlab> are the simulation and time labels to subsequently set into input.dat for use with FlexAOD.


input.dat

Here is a description of the FlexAOD configuration file input.dat:

Parameter Description Examples/Possible values
ver GEOS-Chem version v8-02-02
met GEOS-Chem met-fields version geos5
res GEOS-Chem horizontal resolution 2x25
compiler Fortran compiler option IFORT (used to link to corresponding src/Makefile.hdr.<compiler>)
lab GEOS-Chem simulation label CTRL (optional, see Set up GEOS-Chem timeseries & data_dir)
tlab GEOS-Chem time label 200701 (mandatory, see Set up GEOS-Chem timeseries & data_dir)
typ GEOS-Chem timeseries type inst (ND49) or
sat (ND51)
dini Initial date of FlexAOD run 20070101 (YYYYMMDD)
dend Final date of FlexAOD run 20070131 (YYYYMMDD)
wl Wavelegth (nm) 443 (range limited by optical input)
optsrc Optical properties option OPAC (OPAC/GADS, only option currently)
distsrc Size distribution option STD (associates to data/size_dist/dist_<distsrc>.dat)
hgfsrc Hygroscopic growth option OPAC (from OPAC/GADS) or
CHIN (from Chin et al, 2002)
watsrc Water refractive index option SEGEL (Segelstein, 1981, only option currently)
fullmie Mie calculations option F (create Mie look-up table) or
T (Mie calculation for every grid-cell, NOT YET IMPLEMENTED, VERY CPU INTENSIVE!)
aodlab FlexAOD output label CTRL
gridout Switch for full 3-D output T or F
point_prof Switch for point profile output T or F
point_rtm Switch for point RTM output (libRadtran) T or F
point_file Point output list filename `pwd`/pointout.dat
point_watonly Point output only over water pixels T or F


pointout.dat

The pointout.dat contains the list of points (time and place) for FlexAOD point output. The actual name of the file may vary and is set with variable point_file in input.dat.

The format (free ASCII format) of the file is the following:

column 1: point number
column 2: time (YYYYMMDDHH)
column 3: latitude (degrees north)
column 4: longitude (degrees east)

Here is an example:

#   N     Date/hour     Lat      Lon
#       (YYYYMMDDHH)  (deg N)  (deg E)
1        2006010109    42.0      13.0
2        2006010110    42.0     -33.0

All the lines beginning with "#" will be removed by FlexAOD.

The utility random_pointout.sh generates a random list of points, between user defined dates and coordinates, that may be used for testing or random data model extraction.


Program Launch

Once you set up input.dat and, optionally, pointout.dat, launch FlexAOD from a shell:

./flexaod.sh

or, if you want a background job and log file (bash syntax):

./flexaod.sh &> log &
tail -f log

The utility loop.sh allows automatic FlexAOD launch for several user selected wavelengths and time periods. The loop_wave_rtm.sh is for specific use with point_rtm=T option: it will paste in a unique file the output of several wavelengths.


FlexAOD output

Similarly to GEOS-Chem on-line AOD calculations, FlexAOD performs calculations for 6 "optical" species, assumed in external mixing state:

"Optical" species GEOS-Chem tracers included
SULF SO4, NH4, NIT
OC OCPI, OCPO, SOA1-4
BC BCPI, BCPO
SSa SSa
SSc SSc
DUST DST1-4


FlexAOD has the following types of user selected output:

Description input.dat option default output dir sample file name output variables
Full 3-D gridded output gridout=T <rundir>/../GrADS ts_satellite_<tlab>_<lab>_<dini>-<dend>_<wl>nm_<aodlab>.ctl/dat AOD of "optical" species (same as GEOS-Chem on-line output)

[details]

Profile at selected points point_prof=T
point_file=<filename>
<rundir>/../PROFILES ts_satellite_<tlab>_<lab>_<dini>-<dend>_<wl>nm_<aodlab>_000012 Profile of mass concentrations, number concentration; extinction, scattering and backscattering coefficients; single-scattering albedo; asymmetry factor; phase function coefficients

[details]

Profile at selected points for use with RTM (libRadtran) point_rtm=T
point_file=<filename>
<rundir>/../libRadtran ts_satellite_<tlab>_<lab>_<dini>-<dend>_<wl>nm_<aodlab>_000012_FILES (list of files with model level info)

ts_satellite_<tlab>_<lab>_<dini>-<dend>_<wl>nm_<aodlab>_CONC ("optical" species mass column, g/m^2)

extinction coefficient, single-scattering albedo, phase function coefficients

Any number of output types may be requested (at least one!).


Gridded (GrADS) output

The .ctl/.dat pairs generated for gridded GrADS binary output may be used for direct visualization and analysis in GrADS.

If you don't like GrADS, you may read data with other computer languages, such as Fortran:

program read_grads

 ! Parameters
 integer,paramter :: iou = 10   ! I/O unit
 integer,paramter :: ni = 144   ! n. lons
 integer,paramter :: nj = 91    ! n. lats
 integer,paramter :: nl = 38    ! n. levs
 integer,paramter :: nspecs = 6 ! n. specs
 integer,paramter :: ntimes = 5 ! n. times

 ! Arrays
 real, dimension(ni,nj,nl) :: aottot, boxh
 real, dimension(ni,nj,nl,nspecs) :: aot


 ! Open file in "binary" mode
 open(iou, file="<filename>", status="old", form="binary")

 ! Loop over times
 do it = 1, ntimes

   ! AOT
   read(iou) aottot

   ! Optical Depths of single species
   do k=1,nspecs+1
     do l=1,nl
       read(iou) aot(:,:,l,k)
     enddo
   enddo

   ! Box Height (m)
   do l=1,nl
     read(iou) boxh(:,:,l)
   enddo

 enddo   ! times

 ! Close file
 close(iou)

end program read_grads

If your Fortran compiler does not understand the "binary" form, then just change it to "unformatted" at line 547 of src/flexaod.f90. In that case, however, you will NOT be able to use the output with GrADS.


Profile output

Profile output (point_prof=T) is much easier to handle, because it is in ASCII format. One file per point specified in pointout.dat list is generated. Each file thus contains the calculated profile for a given time/place.

Each output file has nl+1 rows, where nl is the number of levels in GEOS-Chem timeseries. Here is the content of rows:

  • 1st row: column average/integrated values
  • 2nd to (nl+1)th row: values at model levels

Let nspecs = 6 the number of "optical" species and noutpc is the number of output coefficients (default is 16, set at line 1431 of src/flexaod.f90).

Output variables for 1st row are:

  • column 1: 0.000
  • column 2 to nspecs+1: column concentrations of "optical" species (g/m^2)
  • column nspecs+2: column number concentration (#/cm^2)
  • column nspecs+3: aerosol optical depth (unitless)
  • column nspecs+4: scattering aerosol optical depth (unitless)
  • column nspecs+5: back-scattering aerosol optical depth (unitless)
  • column nspecs+6: column-weighted single-scattering albedo (unitless)
  • column nspecs+7: column-weighted asymmetry factor (unitless)
  • column nspecs+8 to nspecs+8+noutpc-1: column-weighted phase function coefficients (unitless)

Output variables for other rows are:

  • column 1: GEOS-Chem level height (km)
  • column 2 to nspecs+1: mass concentrations of "optical" species (ug/m^3)
  • column nspecs+2: number concentration (#/cm^3)
  • column nspecs+3: extinction coefficient (1/km)
  • column nspecs+4: scattering coefficient (1/km)
  • column nspecs+5: back-scattering coefficient (1/km)
  • column nspecs+6: single-scattering albedo (unitless)
  • column nspecs+7: asymmetry factor (unitless)
  • column nspecs+8 to nspecs+8+noutpc-1: phase function coefficients (unitless)


References

  1. Center for International Earth Science Information Network (CIESIN), Columbia University; and Centro Internacional de Agricultura Tropical (CIAT) (2005), Gridded Population of the World Version 3 (GPWv3): National Identifier Grid. Palisades, NY: Socioeconomic Data and Applications Center (SEDAC), Columbia University. http://sedac.ciesin.columbia.edu/gpw/
  2. Chin M, Ginoux P, Kinne S, Torres O, Holben BN, Duncan BN, Martin RV, Logan JA, Higurashi A, Nakajima T (2002), Tropospheric aerosol optical thickness from the GOCART model and comparisons with satellite and Sun photometer measurements, J. Atmos. Sci., 59(3), 461-483. (Table 3 of the paper)
  3. Hess, M., P. Koepke, and I. Schult (1998), Optical Properties of Aerosols and clouds: The software package OPAC, Bull. Am. Met. Soc., 79, 831-844. http://www.lrz.de/~uh234an/www/radaer/opac.html
  4. Mishchenko, M. I., J. M. Dlugach, E. G. Yanovitskij, and N. T. Zakharova, (1999), Bidirectional reflectance of flat, optically thick particulate laters: an efficient radiative transfer solution and applications to snow and soil surfaces, J. Quant. Spectrosc. Radiat. Transfer, Vol. 63, 409-432.
  5. Mishchenko, M. I., L. D. Travis, and A. A. Lacis (2002), Scattering, Absorption, and Emission of Light by Small Particles, Cambridge University Press, Cambridge. http://www.giss.nasa.gov/staff/mmishchenko/brf/
  6. Segelstein, D. (1981), The Complex Refractive Index of Water, M.S. Thesis, University of Missouri, Kansas City. http://www.philiplaven.com/p20.html