From 2d02e155c465636fc30e7fb986654641991070ca Mon Sep 17 00:00:00 2001
From: Jean Braun <Jean.Braun@ujf-grenoble.fr>
Date: Sun, 12 Jun 2011 12:10:56 +0000
Subject: [PATCH] Added phase to sine/cosine wave surface geometries
---
src/create_surfaces.f90 | 10 ++++++----
1 file changed, 6 insertions(+), 4 deletions(-)
diff --git a/src/create_surfaces.f90 b/src/create_surfaces.f90
index 3ffd5e63..1f4833b5 100644
--- a/src/create_surfaces.f90
+++ b/src/create_surfaces.f90
@@ -557,9 +557,10 @@ select case(surface_type)
! a sinus,
! sp01 is the z level
! sp02 is the wavelength
- ! sp03 is the amplitude
+ ! sp03 is the amplitude
+ ! sp04 is the phase
do i=1,ns
- z(i)=sp01+sp03*sin(x(i)*2.d0*pi/sp02)
+ z(i)=sp01+sp03*sin(x(i)*2.d0*pi/sp02+sp04)
end do
case (7)
! a noisy surface,
@@ -583,9 +584,10 @@ select case(surface_type)
! a cosinus,
! sp01 is the z level
! sp02 is the wavelength
- ! sp03 is the amplitude
+ ! sp03 is the amplitude
+ ! sp04 is the phase
do i=1,ns
- z(i)=sp01+sp03*cos(x(i)*2.d0*pi/sp02)
+ z(i)=sp01+sp03*cos(x(i)*2.d0*pi/sp02+sp04)
end do
case (10)
! a 2D embankment
--
GitLab