From 02ee261ab791c9f5422ac5abe7f28d5d39d8613a Mon Sep 17 00:00:00 2001 From: David Minton Date: Tue, 3 Mar 2020 12:55:10 -0500 Subject: [PATCH] Added in test code for new Perlin noise routine (commented out for now) --- src/crater/crater_form_interior.f90 | 3 +++ 1 file changed, 3 insertions(+) diff --git a/src/crater/crater_form_interior.f90 b/src/crater/crater_form_interior.f90 index c039c773..b9f6d4ae 100644 --- a/src/crater/crater_form_interior.f90 +++ b/src/crater/crater_form_interior.f90 @@ -56,6 +56,9 @@ subroutine crater_form_interior(user,surfi,crater,x_relative, y_relative ,newele pikeD = 1.044e3_DP * (crater%fcrat * 1e-3_DP)**(0.301_DP) ! Pike (1977) if ((crater%fcrat > crater%cxtran * 2) .and. newdem < (crater%melev - pikeD)) then newdem = crater%melev - pikeD ! Flatten out the bottom of the crater + !write(*,*) x_relative,y_relative,util_perlin_noise(x_relative,y_relative)*1e3_DP + !newdem = newdem + max(util_perlin_noise(x_relative/crater%fcrat,y_relative/crater%fcrat)*1e3_DP,0.0_DP) + !newdem = newdem + abs(util_perlin_noise(x_relative,y_relative)*1e3_DP) end if if (newdem < (crater%melev - user%deplimit)) then newdem = crater%melev - user%deplimit ! Flatten out the bottom of the crater