skip to main content


Title: Data from: The role of Southeast Asian island topography on Indo-Pacific climate and silicate weathering
The modern configuration of the South East Asian Islands (SEAI) evolved over the last fifteen million years, as a result of subduction, arc magmatism, and arc-continent collisions, contributing to both increased land area and high topography.  The presence of the additional land area has been postulated to enhance convective rainfall, facilitating both increased silicate weathering and the development of the modern-day Walker circulation.  Using an Earth System Model in conjunction with a climate-silicate weathering model, we argue instead for a significant role of SEAI topography for both effects.  This dataset archives model output used in this investigation, including simulations using the Community Earth System Model version 1.2, and the climate-silicate weathering model GEOCLIM. All data are in Netcdf format, and were generated either by the Community Earth System Model 1.2 (Hurrell et al. 2013) or the climate-silicate weathering model GEOCLIM (Park et al. 2020).  Model output is organized into 4 tar files: 1) B1850C5.tar Contains model output for the fully coupled CESM1.2 runs, for 2D fields and for 3D pressure vertical velocity (W) between 10S-10N.  Monthly mean data for years 41-110 of the simulations.   Naming convention is No SEAI topography: B1850C5_noSEAItopo_y41-110.nc and B1850C5_noSEAItopo_W_y41-110.nc 50% SEAI topography: B1850C5_0.5SEAItopo_y41-110.nc and B1850C5_0.5SEAItopo_W_y41-110.nc 100% SEAI topography: B1850C5_y41-110.nc and B1850C5_W_y41-110.nc 150% SEAO topogaphy: B1850C5_1.5SEAItopo_y41-110.nc and B1850C5_1.5SEAItopo_W_y41-110.nc 2) E1850C5.tar Contains model output for the slab ocean CESM1.2 runs, for 2D fields and for 3D pressure vertical velocity (W) between 10S-10N.  Monthly mean data for years 21-50 of the simulations.  Naming convention is No SEAI topography: E1850C5_noSEAItopo_y21-50.nc and E1850C5_noSEAItopo_W_y21-50.nc 50% SEAI topography: E1850C5_0.5SEAItopo_y21-50.nc and E1850C5_0.5SEAItopo_W_y21-50.nc 100% SEAI topography: E1850C5_y21-50.nc and E1850C5_W_y21-50.nc 150% SEAO topogaphy:  E1850C5_1.5SEAItopo_y21-50.nc and E1850C5_1.5SEAItopo_W_y21-50.nc 3) GEOCLIM.tar Contains model output from the climate-silicate weathering model GEOCLIM.  Data is provided for all 573 parameter combinations.  All values are climatological annual means. All files contain these variables: GMST: global mean surface temperature (in K) atm_CO2_level: atmospheric pCO2 (in ppm) degassing: globally-integrated CO2 flux (in mol/yr) The files ending with 1xCO2.nc also contain these spatial fields: lithology fraction: fraction of land covered by a lithology class erosion: Regolith erosion rate (m/yr) weathering: Ca-Mg weathering rate (mol/m^2/yr) E1850C5_1xCO2.nc - GEOCLIM output using the Modern SEAI simulation as input, and for CO2 fixed to 286.7ppm.  E1850C5_noSEAI_1xCO2.nc - GEOCLIM output using the no SEAI simulation as input, and for CO2 fixed to 286.7ppm.  E1850C5_noSEAItopo_1xCO2.nc - GEOCLIM output using the flat SEAI simulation as input, and for CO2 fixed to 286.7ppm.  E1850C5_noSEAI_equil.nc - GEOCLIM output using the no SEAI simulation as input, and CO2 adjusted so that system is in carbon flux equilibrium.   E1850C5_noSEAItopo_flatSEAIslope_equil.nc - GEOCLIM output using the flat SEAI simulation as input, and CO2 adjusted so that system is in carbon flux equilibrium.   4) Surface.tar Contains land fraction and surface geopotential fields for the modern SEAI (Landfrac.nc) and no SEAI (Landfrac_noSEAI.nc) simulations References Hurrell, J.W., Holland, M.M., Gent, P.R., Ghan, S., Kay, J.E., Kushner, P.J., Lamarque, J.F., Large, W.G., Lawrence, D., Lindsay, K. and Lipscomb, W.H., 2013. The community earth system model: a framework for collaborative research. Bulletin of the American Meteorological Society, 94(9), pp.1339-1360. Park, Y., Maffre, P., Goddéris, Y., Macdonald, F.A., Anttila, E.S. and Swanson-Hysell, N.L., 2020. Emergence of the Southeast Asian islands as a driver for Neogene cooling. Proceedings of the National Academy of Sciences, 117(41), pp.25319-25326.  more » « less
Award ID(s):
1925990
NSF-PAR ID:
10431447
Author(s) / Creator(s):
;
Publisher / Repository:
Dryad
Date Published:
Edition / Version:
5
Subject(s) / Keyword(s):
["FOS: Earth and related environmental sciences"]
Format(s):
Medium: X Size: 22163909902 bytes
Size(s):
["22163909902 bytes"]
Sponsoring Org:
National Science Foundation
More Like this
  1. MCMC chains for the GWB analyses performed in the paper "The NANOGrav 15 yr Data Set: Search for Signals from New Physics". 

    The data is provided in pickle format. Each file contains a NumPy array with the MCMC chain (with burn-in already removed), and a dictionary with the model parameters' names as keys and their priors as values. You can load them as

    with open ('path/to/file.pkl', 'rb') as pick: temp = pickle.load(pick) params = temp[0] chain = temp[1]

    The naming convention for the files is the following:

    • igw: inflationary Gravitational Waves (GWs)
    • sigw: scalar-induced GWs
      • sigw_box: assumes a box-like feature in the primordial power spectrum.
      • sigw_delta: assumes a delta-like feature in the primordial power spectrum.
      • sigw_gauss: assumes a Gaussian peak feature in the primordial power spectrum.
    • pt: cosmological phase transitions
      • pt_bubble: assumes that the dominant contribution to the GW productions comes from bubble collisions.
      • pt_sound: assumes that the dominant contribution to the GW productions comes from sound waves.
    • stable: stable cosmic strings
      • stable-c: stable strings emitting GWs only in the form of GW bursts from cusps on closed loops.
      • stable-k: stable strings emitting GWs only in the form of GW bursts from kinks on closed loops.
      • stable-m: stable strings emitting monochromatic GW at the fundamental frequency.
      • stable-n: stable strings described by numerical simulations including GWs from cusps and kinks.
    • meta: metastable cosmic strings
      • meta-l: metastable strings with GW emission from loops only.
      • meta-ls metastable strings with GW emission from loops and segments.
    • super: cosmic superstrings.
    • dw: domain walls
      • dw-sm: domain walls decaying into Standard Model particles.
      • dw-dr: domain walls decaying into dark radiation.

    For each model, we provide four files. One for the run where the new-physics signal is assumed to be the only GWB source. One for the run where the new-physics signal is superimposed to the signal from Supermassive Black Hole Binaries (SMBHB), for these files "_bhb" will be appended to the model name. Then, for both these scenarios, in the "compare" folder we provide the files for the hypermodel runs that were used to derive the Bayes' factors.

    In addition to chains for the stochastic models, we also provide data for the two deterministic models considered in the paper (ULDM and DM substructures). For the ULDM model, the naming convention of the files is the following (all the ULDM signals are superimposed to the SMBHB signal, see the discussion in the paper for more details)

    • uldm_e: ULDM Earth signal.
    • uldm_p: ULDM pulsar signal
      • uldm_p_cor: correlated limit
      • uldm_p_unc: uncorrelated limit
    • uldm_c: ULDM combined Earth + pulsar signal direct coupling 
      • uldm_c_cor: correlated limit
      • uldm_c_unc: uncorrelated limit
    • uldm_vecB: vector ULDM coupled to the baryon number
      • uldm_vecB_cor: correlated limit
      • uldm_vecB_unc: uncorrelated limit 
    • uldm_vecBL: vector ULDM coupled to B-L
      • uldm_vecBL_cor: correlated limit
      • uldm_vecBL_unc: uncorrelated limit
    • uldm_c_grav: ULDM combined Earth + pulsar signal for gravitational-only coupling
      • uldm_c_grav_cor: correlated limit
        • uldm_c_cor_grav_low: low mass region  
        • uldm_c_cor_grav_mon: monopole region
        • uldm_c_cor_grav_low: high mass region
      • uldm_c_unc: uncorrelated limit
        • uldm_c_unc_grav_low: low mass region  
        • uldm_c_unc_grav_mon: monopole region
        • uldm_c_unc_grav_low: high mass region

    For the substructure (static) model, we provide the chain for the marginalized distribution (as for the ULDM signal, the substructure signal is always superimposed to the SMBHB signal)

     
    more » « less
  2. The intended use of this archive is to facilitate meta-analysis of the Data Observation Network for Earth (DataONE, [1]). 

    DataONE is a distributed infrastructure that provides information about earth observation data. This dataset was derived from the DataONE network using Preston [2] between 17 October 2018 and 6 November 2018, resolving 335,213 urls at an average retrieval rate of about 5 seconds per url, or 720 files per hour, resulting in a data gzip compressed tar archive of 837.3 MB .  

    The archive associates 325,757 unique metadata urls [3] to 202,063 unique ecological metadata files [4]. Also, the DataONE search index was captured to establish provenance of how the dataset descriptors were found and acquired. During the creation of the snapshot (or crawl), 15,389 urls [5], or 4.7% of urls, did not successfully resolve. 

    To facilitate discovery, the record of the Preston snapshot crawl is included in the preston-ls-* files . There files are derived from the rdf/nquad file with hash://sha256/8c67e0741d1c90db54740e08d2e39d91dfd73566ea69c1f2da0d9ab9780a9a9f . This file can also be found in the data.tar.gz at data/8c/67/e0/8c67e0741d1c90db54740e08d2e39d91dfd73566ea69c1f2da0d9ab9780a9a9f/data . For more information about concepts and format, please see [2]. 

    To extract all EML files from the included Preston archive, first extract the hashes assocated with EML files using:

    cat preston-ls.tsv.gz | gunzip | grep "Version" | grep -v "deeplinker" | grep -v "query/solr" | cut -f1,3 | tr '\t' '\n' | grep "hash://" | sort | uniq > eml-hashes.txt

    extract data.tar.gz using:

    ~/preston-archive$ tar xzf data.tar.gz 

    then use Preston to extract each hash using something like:

    ~/preston-archive$ preston get hash://sha256/00002d0fc9e35a9194da7dd3d8ce25eddee40740533f5af2397d6708542b9baa
    <eml:eml xmlns:eml="eml://ecoinformatics.org/eml-2.1.1" xmlns:xs="http://www.w3.org/2001/XMLSchema" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:stmml="http://www.xml-cml.org/schema/stmml_1.1" packageId="doi:10.18739/A24P9Q" system="https://arcticdata.io" scope="system" xsi:schemaLocation="eml://ecoinformatics.org/eml-2.1.1 ~/development/eml/eml.xsd">
      <dataset>
        <alternateIdentifier>urn:x-wmo:md:org.aoncadis.www::d76bc3b5-7b19-11e4-8526-00c0f03d5b7c</alternateIdentifier>
        <alternateIdentifier>d76bc3b5-7b19-11e4-8526-00c0f03d5b7c</alternateIdentifier>
        <title>Airglow Image Data 2011 4 of 5</title>
    ...

    Alternatively, without using Preston, you can extract the data using the naming convention:

    data/[x]/[y]/[z]/[hash]/data

    where x is the first 2 characters of the hash, y the second 2 characters, z the third 2 characters, and hash the full sha256 content hash of the EML file.

    For example, the hash hash://sha256/00002d0fc9e35a9194da7dd3d8ce25eddee40740533f5af2397d6708542b9baa can be found in the file: data/00/00/2d/00002d0fc9e35a9194da7dd3d8ce25eddee40740533f5af2397d6708542b9baa/data . For more information, see [2].

    The intended use of this archive is to facilitate meta-analysis of the DataONE dataset network. 

    [1] DataONE, https://www.dataone.org
    [2] https://preston.guoda.bio, https://doi.org/10.5281/zenodo.1410543 . DataONE was crawled via Preston with "preston update -u https://dataone.org".
    [3] cat preston-ls.tsv.gz | gunzip | grep "Version" | grep -v "deeplinker" | grep -v "query/solr" | cut -f1,3 | tr '\t' '\n' | grep -v "hash://" | sort | uniq | wc -l
    [4] cat preston-ls.tsv.gz | gunzip | grep "Version" | grep -v "deeplinker" | grep -v "query/solr" | cut -f1,3 | tr '\t' '\n' | grep "hash://" | sort | uniq | wc -l
    [5] cat preston-ls.tsv.gz | gunzip | grep "Version" | grep  "deeplinker" | grep -v "query/solr" | cut -f1,3 | tr '\t' '\n' | grep -v "hash://" | sort | uniq | wc -l

    This work is funded in part by grant NSF OAC 1839201 from the National Science Foundation.

     
    more » « less
  3. This dataset contains monthly average output files from the iCAM6 simulations used in the manuscript "Enhancing understanding of the hydrological cycle via pairing of process-oriented and isotope ratio tracers," in review at the Journal of Advances in Modeling Earth Systems. A file corresponding to each of the tagged and isotopic variables used in this manuscript is included. Files are at 0.9° latitude x 1.25° longitude, and are in NetCDF format. Data from two simulations are included: 1) a simulation where the atmospheric model was "nudged" to ERA5 wind and surface pressure fields, by adding an additional tendency (see section 3.1 of associated manuscript), and 2) a simulation where the atmospheric state was allowed to freely evolve, using only boundary conditions imposed at the surface and top of atmosphere. Specific information about each of the variables provided is located in the "usage notes" section below. Associated article abstract: The hydrologic cycle couples the Earth's energy and carbon budgets through evaporation, moisture transport, and precipitation. Despite a wealth of observations and models, fundamental limitations remain in our capacity to deduce even the most basic properties of the hydrological cycle, including the spatial pattern of the residence time (RT) of water in the atmosphere and the mean distance traveled from evaporation sources to precipitation sinks. Meanwhile, geochemical tracers such as stable water isotope ratios provide a tool to probe hydrological processes, yet their interpretation remains equivocal despite several decades of use. As a result, there is a need for new mechanistic tools that link variations in water isotope ratios to underlying hydrological processes. Here we present a new suite of “process-oriented tags,” which we use to explicitly trace hydrological processes within the isotopically enabled Community Atmosphere Model, version 6 (iCAM6). Using these tags, we test the hypotheses that precipitation isotope ratios respond to parcel rainout, variations in atmospheric RT, and preserve information regarding meteorological conditions during evaporation. We present results for a historical simulation from 1980 to 2004, forced with winds from the ERA5 reanalysis. We find strong evidence that precipitation isotope ratios record information about atmospheric rainout and meteorological conditions during evaporation, but little evidence that precipitation isotope ratios vary with water vapor RT. These new tracer methods will enable more robust linkages between observations of isotope ratios in the modern hydrologic cycle or proxies of past terrestrial environments and the environmental processes underlying these observations.   Details about the simulation setup can be found in section 3 of the associated open-source manuscript, "Enhancing understanding of the hydrological cycle via pairing of process‐oriented and isotope ratio tracers." In brief, we conducted two simulations of the atmosphere from 1980-2004 using the isotope-enabled version of the Community Atmosphere Model 6 (iCAM6) at 0.9x1.25° horizontal resolution, and with 30 vertical hybrid layers spanning from the surface to ~3 hPa. In the first simulation, wind and surface pressure fields were "nudged" toward the ERA5 reanalysis dataset by adding a nudging tendency, preventing the model from diverging from observed/reanalysis wind fields. In the second simulation, no additional nudging tendency was included, and the model was allowed to evolve 'freely' with only boundary conditions provided at the top (e.g., incoming solar radiation) and bottom (e.g., observed sea surface temperatures) of the model. In addition to the isotopic variables, our simulation included a suite of 'process-oriented tracers,' which we describe in section 2 of the manuscript. These variables are meant to track a property of water associated with evaporation, condensation, or atmospheric transport. Metadata are provided about each of the files below; moreover, since the attached files are NetCDF data - this information is also provided with the data files. NetCDF metadata can be accessed using standard tools (e.g., ncdump). Each file has 4 variables: the tagged quantity, and the associated coordinate variables (time, latitude, longitude). The latter three are identical across all files, only the tagged quantity changes. Twelve files are provided for the nudged simulation, and an additional three are provided for the free simulations: Nudged simulation files iCAM6_nudged_1980-2004_mon_RHevap: Mass-weighted mean evaporation source property: RH (%) with respect to surface temperature. iCAM6_nudged_1980-2004_mon_Tevap: Mass-weighted mean evaporation source property: surface temperature in Kelvin iCAM6_nudged_1980-2004_mon_Tcond: Mass-weighted mean condensation property: temperature (K) iCAM6_nudged_1980-2004_mon_columnQ: Total (vertically integrated) precipitable water (kg/m2).  Not a tagged quantity, but necessary to calculate depletion times in section 4.3 (e.g., Fig. 11 and 12). iCAM6_nudged_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_nudged_1980-2004_mon_d18Oevap_0: Mass-weighted mean evaporation source property - d18O of the evaporative flux (e.g., the 'initial' isotope ratio prior to condensation), (‰ VSMOW) iCAM6_nudged_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_nudged_1980-2004_mon_dexevap_0: Mass-weighted mean evaporation source property - deuterium excess of the evaporative flux iCAM6_nudged_1980-2004_mon_lnf: Integrated property - ln(f) calculated from the constant-fractionation d18O tracer (see section 3.2). iCAM6_nudged_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation. iCAM6_nudged_1980-2004_mon_residencetime: Mean atmospheric water residence time (in days). iCAM6_nudged_1980-2004_mon_transportdistance: Mean atmospheric water transport distance (in km). Free simulation files iCAM6_free_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_free_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_free_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation. 
    more » « less
  4. This data set for the manuscript entitled "Design of Peptides that Fold and Self-Assemble on Graphite" includes all files needed to run and analyze the simulations described in the this manuscript in the molecular dynamics software NAMD, as well as the output of the simulations. The files are organized into directories corresponding to the figures of the main text and supporting information. They include molecular model structure files (NAMD psf or Amber prmtop format), force field parameter files (in CHARMM format), initial atomic coordinates (pdb format), NAMD configuration files, Colvars configuration files, NAMD log files, and NAMD output including restart files (in binary NAMD format) and trajectories in dcd format (downsampled to 10 ns per frame). Analysis is controlled by shell scripts (Bash-compatible) that call VMD Tcl scripts or python scripts. These scripts and their output are also included.

    Version: 2.0

    Changes versus version 1.0 are the addition of the free energy of folding, adsorption, and pairing calculations (Sim_Figure-7) and shifting of the figure numbers to accommodate this addition.


    Conventions Used in These Files
    ===============================

    Structure Files
    ----------------
    - graph_*.psf or sol_*.psf (original NAMD (XPLOR?) format psf file including atom details (type, charge, mass), as well as definitions of bonds, angles, dihedrals, and impropers for each dipeptide.)

    - graph_*.pdb or sol_*.pdb (initial coordinates before equilibration)
    - repart_*.psf (same as the above psf files, but the masses of non-water hydrogen atoms have been repartitioned by VMD script repartitionMass.tcl)
    - freeTop_*.pdb (same as the above pdb files, but the carbons of the lower graphene layer have been placed at a single z value and marked for restraints in NAMD)
    - amber_*.prmtop (combined topology and parameter files for Amber force field simulations)
    - repart_amber_*.prmtop (same as the above prmtop files, but the masses of non-water hydrogen atoms have been repartitioned by ParmEd)

    Force Field Parameters
    ----------------------
    CHARMM format parameter files:
    - par_all36m_prot.prm (CHARMM36m FF for proteins)
    - par_all36_cgenff_no_nbfix.prm (CGenFF v4.4 for graphene) The NBFIX parameters are commented out since they are only needed for aromatic halogens and we use only the CG2R61 type for graphene.
    - toppar_water_ions_prot_cgenff.str (CHARMM water and ions with NBFIX parameters needed for protein and CGenFF included and others commented out)

    Template NAMD Configuration Files
    ---------------------------------
    These contain the most commonly used simulation parameters. They are called by the other NAMD configuration files (which are in the namd/ subdirectory):
    - template_min.namd (minimization)
    - template_eq.namd (NPT equilibration with lower graphene fixed)
    - template_abf.namd (for adaptive biasing force)

    Minimization
    -------------
    - namd/min_*.0.namd

    Equilibration
    -------------
    - namd/eq_*.0.namd

    Adaptive biasing force calculations
    -----------------------------------
    - namd/eabfZRest7_graph_chp1404.0.namd
    - namd/eabfZRest7_graph_chp1404.1.namd (continuation of eabfZRest7_graph_chp1404.0.namd)

    Log Files
    ---------
    For each NAMD configuration file given in the last two sections, there is a log file with the same prefix, which gives the text output of NAMD. For instance, the output of namd/eabfZRest7_graph_chp1404.0.namd is eabfZRest7_graph_chp1404.0.log.

    Simulation Output
    -----------------
    The simulation output files (which match the names of the NAMD configuration files) are in the output/ directory. Files with the extensions .coor, .vel, and .xsc are coordinates in NAMD binary format, velocities in NAMD binary format, and extended system information (including cell size) in text format. Files with the extension .dcd give the trajectory of the atomic coorinates over time (and also include system cell information). Due to storage limitations, large DCD files have been omitted or replaced with new DCD files having the prefix stride50_ including only every 50 frames. The time between frames in these files is 50 * 50000 steps/frame * 4 fs/step = 10 ns. The system cell trajectory is also included for the NPT runs are output/eq_*.xst.

    Scripts
    -------
    Files with the .sh extension can be found throughout. These usually provide the highest level control for submission of simulations and analysis. Look to these as a guide to what is happening. If there are scripts with step1_*.sh and step2_*.sh, they are intended to be run in order, with step1_*.sh first.


    CONTENTS
    ========

    The directory contents are as follows. The directories Sim_Figure-1 and Sim_Figure-8 include README.txt files that describe the files and naming conventions used throughout this data set.

    Sim_Figure-1: Simulations of N-acetylated C-amidated amino acids (Ac-X-NHMe) at the graphite–water interface.

    Sim_Figure-2: Simulations of different peptide designs (including acyclic, disulfide cyclized, and N-to-C cyclized) at the graphite–water interface.

    Sim_Figure-3: MM-GBSA calculations of different peptide sequences for a folded conformation and 5 misfolded/unfolded conformations.

    Sim_Figure-4: Simulation of four peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) at the graphite–water interface at 370 K.

    Sim_Figure-5: Simulation of four peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) at the graphite–water interface at 295 K.

    Sim_Figure-5_replica: Temperature replica exchange molecular dynamics simulations for the peptide cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) with 20 replicas for temperatures from 295 to 454 K.

    Sim_Figure-6: Simulation of the peptide molecule cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) in free solution (no graphite).

    Sim_Figure-7: Free energy calculations for folding, adsorption, and pairing for the peptide CHP1404 (sequence: cyc(GTGSGTG-GPGG-GCGTGTG-SGPG)). For folding, we calculate the PMF as function of RMSD by replica-exchange umbrella sampling (in the subdirectory Folding_CHP1404_Graphene/). We make the same calculation in solution, which required 3 seperate replica-exchange umbrella sampling calculations (in the subdirectory Folding_CHP1404_Solution/). Both PMF of RMSD calculations for the scrambled peptide are in Folding_scram1404/. For adsorption, calculation of the PMF for the orientational restraints and the calculation of the PMF along z (the distance between the graphene sheet and the center of mass of the peptide) are in Adsorption_CHP1404/ and Adsorption_scram1404/. The actual calculation of the free energy is done by a shell script ("doRestraintEnergyError.sh") in the 1_free_energy/ subsubdirectory. Processing of the PMFs must be done first in the 0_pmf/ subsubdirectory. Finally, files for free energy calculations of pair formation for CHP1404 are found in the Pair/ subdirectory.

    Sim_Figure-8: Simulation of four peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) where the peptides are far above the graphene–water interface in the initial configuration.

    Sim_Figure-9: Two replicates of a simulation of nine peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) at the graphite–water interface at 370 K.

    Sim_Figure-9_scrambled: Two replicates of a simulation of nine peptide molecules with the control sequence cyc(GGTPTTGGGGGGSGGPSGTGGC) at the graphite–water interface at 370 K.

    Sim_Figure-10: Adaptive biasing for calculation of the free energy of the folded peptide as a function of the angle between its long axis and the zigzag directions of the underlying graphene sheet.

     

    This material is based upon work supported by the US National Science Foundation under grant no. DMR-1945589. A majority of the computing for this project was performed on the Beocat Research Cluster at Kansas State University, which is funded in part by NSF grants CHE-1726332, CNS-1006860, EPS-1006860, and EPS-0919443. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1548562, through allocation BIO200030. 
    more » « less
  5. This archive contains COAWST model input, grids and initial conditions, and output used to produce the results in a submitted manuscript. The files are:

    model_input.zip: input files for simulations presented in this paper
      ocean_rip_current.in: ROMS ocean model input file
      swan_rip_current.in: SWAN wave model input file (example with Hs=1m)
      coupling_rip_current.in: model coupling file
      rip_current.h: model header file
      
    model_grids_forcing.zip: bathymetry and initial condition files
         hbeach_grid_isbathy_2m.nc: ROMS bathymetry input file
         hbeach_grid_isbathy_2m.bot: SWAN bathymetry input file
         hbeach_grid_isbathy_2m.grd: SWAN grid input file
         hbeach_init_isbathy_14_18_17.nc: Initial temperature, cool surf zone dT=-1C case
         hbeach_init_isbathy_14_18_19.nc: Initial temperature, warm surf zone dT=+1C case
         hbeach_init_isbathy_14_18_16.nc: Initial temperature, cool surf zone dT=-2C case
         hbeach_init_isbathy_14_18_20.nc: Initial temperature, warm surf zone dT=+2C case
         hbeach_init_isbathy_14_18_17p5.nc: Initial temperature, cool surf zone dT=-0.5C case
         hbeach_init_isbathy_14_18_18p5.nc: Initial temperature, warm surf zone dT=+0.5C case

    model_output files: model output used to produce the figures
         netcdf files, zipped
         variables included:
              x_rho (cross-shore coordinate, m)
              y_rho (alongshore coordinate, m)
              z_rho (vertical coordinate, m)
              ocean_time (time since initialization, s, output every 5 mins)
              h (bathymetry, m)
              temp (temperature, Celsius)
              dye_02 (surfzone-released dye)
              Hwave (wave height, m)
              Dissip_break (wave dissipation W/m2) 
              ubar (cross-shore depth-average velocity, m/s, interpolated to rho-points)
         Case_141817.nc: cool surf zone dT=-1C Hs=1m
         Case_141819.nc: warm surf zone dT=+1C Hs=1m
         Case_141816.nc: cool surf zone dT=-2C Hs=1m
         Case_141820.nc: warm surf zone dT=-2C Hs=1m
         Case_141817p5.nc: cool surf zone dT=-0.5C Hs=1m
         Case_141818p5.nc: warm surf zone dT=+0.5C Hs=1m
         Case_141817_Hp5.nc: cool surf zone dT=-1C Hs=0.5m
         Case_141819_Hp5.nc: warm surf zone dT=+1C Hs=0.5m
         Case_141817_Hp75.nc: cool surf zone dT=-1C Hs=0.75m
         Case_141819_Hp75.nc: warm surf zone dT=+1C Hs=0.75m

    COAWST is an open source code and can be download at https://coawstmodel-trac.sourcerepo.com/coawstmodel_COAWST/. Descriptions of the input and output files can be found in the manual distributed with the model code and in the glossary at the end of the ocean.in file.

    Corresponding author: Melissa Moulton, mmoulton@uw.edu

     
    more » « less