ACME surface dataset including land cover information for year 1850-2015 at ne120np4 resolution
This page borrows liberally from Gautam's excellent CLM4.5 dataset for year 1850 at ne120np4 resolution. It adds some more information on how to incorporate
land cover information while creating a surface dataset.
Date: 15-Nov-2016
Machine: Edison
1. Log on edison and clone ACME repository
ssh user-name@edison.nersc.gov
cd $SCRATCH
git clone git@github.com:ACME-Climate/ACME.git
cd ACME
export ACME_DIR=$PWD
We assume that mapping files have already been generated. If not, follow instructions from here: CLM4.5 dataset for year 1850 at ne120np4 resolution
Mapping files are available here: /global/project/projectdirs/acme/inputdata/lnd/clm2/mappingdata/maps/
2. Compile surface dataset source code
module load cray-netcdf
module load cray-hdf5
export LIB_NETCDF=$NETCDF_DIR/lib
export INC_NETCDF=$NETCDF_DIR/include
export USER_FC=ifort
export USER_CC=icc
export USER_LDFLAGS="-L$NETCDF_DIR/lib -lnetcdf -lnetcdff -lnetcdf_intel"
export USER_LDFLAGS=$USER_LDFLAGS" -L$HDF5_DIR/lib -lhdf5 -lhdf5_fortran -lhdf5_cpp -lhdf5_fortran_intel -lhdf5_hl_intel -lhdf5hl_fortran_intel"
cd $SCRATCH/ACME/components/clm/tools/clm4_5/mksurfdata_map/src/
gmake
3.Run the mksurfdata.pl script in "debug" mode to generate the namelist
cd $SCRATCH/ACME/components/clm/tools/clm4_5/mksurfdata_map/
For 1950 high res:
./mksurfdata.pl -res ne120np4 -d -dinlc /global/project/projectdirs/acme/inputdata -usr_mapdir /global/project/projectdirs/acme/inputdata/lnd/clm2/mappingdata/maps/ne120np4
For 1850 - 2015 low-res:
./mksurfdata.pl -res ne30np4 -years 1850,2015 -d -dinlc /global/project/projectdirs/acme/inputdata -usr_mapdir /global/project/projectdirs/acme/inputdata/lnd/clm2/mappingdata/maps/ne30np4
4. Modify namelist file if land cover information should be used
mksrf_fdynuse should contain name of file that contains all the land cover file paths
mksrf_fdynuse = 'landuse_1850-1855.txt'
mksrf_fvegtyp contains paths of all land cover files indicated in mksrf_fdynuse
mksrf_fvegtyp = '/scratch2/scratchdirs/ritvik/ACME/components/clm/tools/clm4_5/mksurfdata_map/luh2/mksrf_landuse_rc_1850_10202016.nc'
5. Make
cd $SCRATCH/ACME/components/clm/tools/clm4_5/mksurfdata_map/src
>> gmake clean
>> gmake
6. Create surface dataset
./mksurfdata_map < namelist