Working version of code for NERSC