README.aerosols.public


Contact Nan Rosenbloom  or Christine Shields to download a tar file of the scripts described on this page.


current aerosol method
(developed by n.mahowald, and c. shields)

overview
--------
1. run CAM/CCSM4 with bulk aerosols,("bam", not default prescribed), i.e. trop_bam
2. this gives us dust and sea salt (primary aerosols)
3. for bam: create zonal distributions of oxid, dms, bc, and oc and soil erodiablity based on land
   i.e., avg modern grid points over ocean (where appropriate), avg modern grid pts over land, give
   avg values to paleo 2d data forcing files over respective lat and land
4. create prescribed aerosol (radiative) set based on BAM run
5. create prescribed aerosol (deposition) based on run BAM run
6. run CAM with prescribed aerosols


directory structure
-------------------
bulk_aero  -> scripts, programs, and readmes necessary for bulk_aerosol run
docs       -> paleo ccsm3 tech note + ccsm4 updates
presc_aero -> scripts, programs, and readmes necessary for prescribed aerosol run
README.aerosols.public -> this readme file with generic paleo instructions


details
-------

===== BAM run ============

--(1) create forcing data w/ncl scripts in bulk_aero subdir

notes:
see cam.input_data_list for bulk aerosol forcing files:
don't need deposition velocity stuff unless using chemistry
1) ignore depvel_monthly.nc, clim_soilw.nc, and season_wes.nc;
   legacy files, not needed
2) ignore regrid_vegetation.nc -- not needed

need new datasets for:
1) oxid: oxid_1.9x2.5_L26_1850clim_c091123.nc
2) soil erod:  dst_48x96_c090203.nc
3) so2:  aerocom_SO2_vertical_1750.c080807.nc
         aerocom_SO2_surface_1750.c080807.nc  
4) so4:  aerocom_SO4_vertical_1750.c080807.nc
         aerocom_SO4_surface_1750.c080807.nc
5) CB1:  aerocom_CB1_1750.c080807.nc 
6) DMS:  aerocom_DMS_2000.c080417.nc
7) OC1:  aerocom_OC1_1750.nosoa.c080807.nc

modify ncl scripts to read in template (default) file, and modify accordingly
details: see README.whatididforPT

--(2) run CAM w/bulk aerosols, i.e. in env_conf

      ***see step (4) if you want to include extra output variables for deposition
      in this run. RECOMMENDED.(space/time resource dependent).

<!--"CAM configure options, see CAM configure utility for details (char) " -->
<entry id="CAM_CONFIG_OPTS"   value="-phys cam4 -chem trop_bam"  />

<!--"CAM namelist use_case (char) " -->
<entry id="CAM_NML_USE_CASE"   value="1850_cam4_trop_bam"  />

**tracer_cnst date spec for oxide file changes from
tracer_cnst_ymd =  = 18500101  (cesm_1_0_beta13)
 to
tracer_cnst_cycle_yr = 1850   (cesm_1_0_beta23)
 


--(3) take default output from BAM and create prescribed forcing aerosol (raditative)
        climatology (from jf lamarque)


        (a) for each month:
ncra -h -v
lat,lon,lev,time,date,datesec,hyam,hybm,hyai,hybi,PS,P0,T,CB1,CB2,OC1,OC2,SO4,
DST01,DST02,DST03,DST04,SSLT01,SSLT02,SSLT03,SSLT04

       (b) concatenate all into a single file (beware that sometime p0 becomes 0)

       (c) you probably need to overwrite date and time as those also get messed
          up with ncra.

        SCRIPT written: presc_aero/aero_climo.csh

        (d) NOTE: does NOT need to be interpolated to fv2deg res (as default). model will accept
            any resolution, (eaton), and interp, on the fly, if necessary. normally all model
            configs use the fv2deg, but you can supply a file at the correct config resolution to
            begin with...

                ** although (d) is true, older released code is buggy and the aero file may not
                   work,  use aero_climo_interp.ncl to put onto fv2deg grid. (march 4 2011). bug
                   reported to b.eaton and f.vitt.


--(4) branch basic default (BAM) run (or change in original) to include vars necessary for deposition
       to create aero dep file for prescribed aero run.
        (from jf lamarque)

        (a) monthly avg hist needed:
        DF_CB1   (need to include in finc1)
        DF_CB2   (need to include in finc1)
        CB2SFWET  (may need to include depending on release: bam)
        DF_OC1    (need to include in finc1)
        DF_OC2    (need to include in finc1)
        OC2SFWET  (may need to include depending on release:  bam)
        DST01DD   (default)
        DST01PP   ("")
        DST02DD  ("")
        DST02PP  ("")
        DST03DD  ("")
        DST03PP  ("")
        DST04DD  ("")
        DST04PP  ("")

        i.e. in cam.buildnml.csh, 
        fincl1                               = 'DF_CB1:A','DF_CB2:A','DF_OC1:A','DF_OC2:A'

        (b)  run dep_snow.f90 (from jf), in presc_aero/dep  subdir:

        (c)  overwrite time and date with 1850 default values, and add source run to global history
              script to do this:                aerodep.nco.csh

ncks --append -v date,time date_time_1850clim.nc aerosoldep_monthly_12-21_mean_b40.t31x3.new.nc
ncatted -a history,global,a,c,"b40.t31x3.new.cam2.h0* files years 12-21 cshields yymmdd" aerosoldep_monthly_12-21_mean_b40.t31x3.pt.new.nc


        (d) this seems to work fine w/T31.


--- notes on oc1/2 and bc1/2:
from nick heavens (cornell):
OC1,2 and BC1,2 are internal mixtures of hydrophilic and
hydrophobic carbon species. 1 is 80:20 hydrophobic:hydrophilic, 2 is 50:50.