Added in test code for new Perlin noise routine (commented out for now)

MintonGroup · Mar 3, 2020 · 02ee261 · 02ee261
1 parent 3d0cf1e
commit 02ee261
Showing 1 changed file with 3 additions and 0 deletions.
diff --git 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