Add surface wetting module

[chrisjonesBSU]

This is still a work in progress, but begins a surface wetting module that we've discussed before. We can use this as a place to discuss the design choices we'll have to make.

[codecov]

Codecov Report

Merging #96 (0c8d139) into main (531cd82) will increase coverage by 0.15%.
The diff coverage is 95.23%.

❗ Current head 0c8d139 differs from pull request most recent head c8105bc. Consider uploading reports for the commit c8105bc to get more accurate results

Additional details and impacted files

@@            Coverage Diff             @@
##             main      #96      +/-   ##
==========================================
+ Coverage   93.87%   94.02%   +0.15%     
==========================================
  Files          20       23       +3     
  Lines        1550     1724     +174     
==========================================
+ Hits         1455     1621     +166     
- Misses         95      103       +8     

Files	Coverage Δ
flowermd/base/molecule.py	96.58% <ø> (ø)
flowermd/base/simulation.py	91.44% <ø> (ø)
flowermd/base/system.py	91.90% <100.00%> (+0.38%)	⬆️
flowermd/library/surfaces.py	100.00% <100.00%> (ø)
flowermd/modules/surface_wetting/__init__.py	100.00% <100.00%> (ø)
flowermd/modules/welding/__init__.py	100.00% <ø> (ø)
flowermd/modules/welding/utils.py	100.00% <ø> (ø)
flowermd/modules/welding/welding.py	100.00% <ø> (ø)
flowermd/modules/surface_wetting/utils.py	95.31% <95.31%> (ø)
...lowermd/modules/surface_wetting/surface_wetting.py	94.11% <94.11%> (ø)

[review-notebook-app]

Check out this pull request on 

See visual diffs & provide feedback on Jupyter Notebooks.

---

Powered by ReviewNB

[chrisjonesBSU]

The tutorial looks good! Just one suggestion, I think we should try to use the same formatting we have in pre-commit in the tutorial code snippets.

For example:

graphene = Graphene(x_repeat=20, y_repeat=20, n_layers=1,
                    base_units=droplet_system.reference_values)
graphene.apply_forcefield(force_field=OPLS_AA(), remove_charges=True,
                          auto_scale=False, r_cut=2.5)

Should become

graphene = Graphene(
    x_repeat=20,
    y_repeat=20,
    n_layers=1,
    base_units=droplet_system.reference_values
)

graphene.apply_forcefield(
    force_field=OPLS_AA(),
    remove_charges=True,
    auto_scale=False,
    r_cut=2.5
)

I think it is easier to read and following along this way.