diff options
Diffstat (limited to 'src/gauss2d.x')
| -rw-r--r-- | src/gauss2d.x | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/src/gauss2d.x b/src/gauss2d.x new file mode 100644 index 0000000..19af15f --- /dev/null +++ b/src/gauss2d.x @@ -0,0 +1,5 @@ + gtil = amp*(exp(-((etagrd(i) - eta0)**2/sigma**2)) + exp(-((etagrd(i) + eta0 + & )**2/sigma**2))) + dngtil = amp*((-2.00000000000000d0*(etagrd(i) - eta0)*exp(-((etagrd(i) - eta + & 0)**2/sigma**2)))/sigma**2 - (2.00000000000000d0*(etagrd(i) + eta0)*ex + & p(-((etagrd(i) + eta0)**2/sigma**2)))/sigma**2) |