ACME surface dataset including land cover information for year 1850-2015 at ne120np4 resolution

Owned by Ritvik Sahajpal (Unlicensed)
Last updated: May 26, 2020 by Tanya Reshel

  NERSC Directory Change Notice  

Due to project's name change at NERSC from 'ACME' to 'E3SM' and NERSC’s file system update, the directory  '/project/projectdirs/acme/' is now '/cfs/cdirs/e3sm'

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