diff --git a/src/realistic/realistic_crater_topography.f90 b/src/realistic/realistic_crater_topography.f90 index 81cb1dab..6768fa6b 100644 --- a/src/realistic/realistic_crater_topography.f90 +++ b/src/realistic/realistic_crater_topography.f90 @@ -164,11 +164,11 @@ subroutine make_realistic_crater(user,surf,crater,deltaMtot) psd_elevation%num_sine=psd_elevation%num_vertices/2 psd_elevation%input_x_max=PI*psd_elevation%diameter_in_km psd_elevation%diameter_in_km_trans=20.0_DP - psd_elevation%bp2_x_k_b_sigma=[0.0417732075218209_DP , 2.48787783435609e-18_DP , 0.496231041304548_DP , 0.00210183162260986_DP , 0.793427517984220_DP , 0.308862248697524_DP] - psd_elevation%bp2_y_k_b_sigma=[0.0593460236785858_DP , 0.408350328939867_DP , 1.53480069842064_DP , -0.00201958329120022_DP , 1.63566246833559_DP , 0.809757199553671_DP] - psd_elevation%bp3_y_k_b_sigma=[0.0561324763037810_DP , 2.46101909781936_DP , 0.537503435665947_DP , -0.00292059238533839_DP , 3.64208047160175_DP , 0.462627466875698_DP] - psd_elevation%bp4_y_k_b_sigma=[0.0159696581958859_DP , 3.25122353414224_DP , 0.660958440041299_DP , 0.00638368649794204_DP , 3.44294296810112_DP , 0.769783029462181_DP] - psd_elevation%slope12_k_b_sigma=[0.0279866462789263_DP , 2.53280320608848_DP , 0.492438940196540_DP , 0.000799959919780047_DP , 3.07653693327141_DP , 0.400693610967074_DP] + psd_elevation%bp2_x_k_b_sigma=[0.0417732075218209_DP , 2.48787783435609e-18_DP , 0.496231041304548_DP , 0.00210183162260986_DP , 0.793427517984220_DP , 0.308862248697524_DP] + psd_elevation%bp2_y_k_b_sigma=[0.0593460236785858_DP , 0.408350328939867_DP , 1.53480069842064_DP , -0.00201958329120022_DP , 1.63566246833559_DP , 0.809757199553671_DP] + psd_elevation%bp3_y_k_b_sigma=[0.0561324763037810_DP , 2.46101909781936_DP , 0.537503435665947_DP , -0.00292059238533839_DP , 3.64208047160175_DP , 0.462627466875698_DP] + psd_elevation%bp4_y_k_b_sigma=[0.0159696581958859_DP , 3.25122353414224_DP , 0.660958440041299_DP , 0.00638368649794204_DP , 3.44294296810112_DP , 0.769783029462181_DP] + psd_elevation%slope12_k_b_sigma=[0.0279866462789263_DP , 2.53280320608848_DP , 0.492438940196540_DP , 0.000799959919780047_DP , 3.07653693327141_DP , 0.400693610967074_DP] psd_elevation%misfit_k_b_sigma=0.000000_DP call Calculate_am_wl_phase_from_diameter(psd_elevation,amplitude_rim_elevation,wavelength_rim_elevation,phase_rim_elevation)