Index: /issm/trunk-jpl/test/NightlyRun/test2084.m =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test2084.m (revision 25297) +++ /issm/trunk-jpl/test/NightlyRun/test2084.m (revision 25298) @@ -38,5 +38,5 @@ %loading love numbers field_names ={'LoveH_loading_elastic','LoveK_loading_elastic','LoveL_loading_elastic'}; -field_tolerances={4.3e-9,4.3e-9,4.3e-9}; +field_tolerances={2.0e-8,2.0e-8,2.0e-8}; field_values={... (md.results.LoveSolution.LoveHr(:,1)),... Index: /issm/trunk-jpl/test/NightlyRun/test2085.m =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test2085.m (revision 25297) +++ /issm/trunk-jpl/test/NightlyRun/test2085.m (revision 25298) @@ -1,45 +1,177 @@ +%Test Name: FourierLoveKernels +% Homogeneous Earth, for which analytic solutions exist. +% Love kernels for degree 2 (tested against analytic solns). -%Test Name: LovenumberstAtDepth. -%Same as test #1 of test2084.m +% for volumetric potential + md=model(); + + md.materials=materials('litho'); -md=model(); -md.cluster=generic('name',oshostname(),'np',1); + md.materials.numlayers = 10; + md.love.forcing_type = 9; -md.materials=materials('litho'); -md.miscellaneous.name='FourierLoveTest'; -md.groundingline.migration='None'; + md.materials.density= (1:md.materials.numlayers)'*0+5511; + md.materials.lame_mu= (1:md.materials.numlayers)'*0+0.75e11; + md.materials.viscosity=(1:md.materials.numlayers)'*0+1e21; + md.materials.lame_lambda=md.materials.lame_mu*0+5e17; + md.materials.issolid=ones(md.materials.numlayers,1); + md.materials.isburgers=zeros(md.materials.numlayers,1); + md.materials.burgers_mu=md.materials.lame_mu/3; + md.materials.burgers_viscosity=md.materials.viscosity/10; -md.verbose=verbose('111111101'); -cst=365.25*24*3600*1000; + md.materials.radius = linspace(10e3,6371e3,md.materials.numlayers+1)'; + md.love.g0 = 9.8134357285509388; % directly grabbed from fourierlovesolver for this particular case. - md.materials.numlayers=6; - md.materials.radius = [10 1222.5 3.4800e+03 5.7010e+03 5.9510e+03 6.3010e+03 6.3710e+03]'*1e3; - md.materials.density= [1.0750e4 1.0750e+04 4.9780e+03 3.8710e+03 3.4380e+03 3.0370e+03]'; - md.materials.lame_mu= [1e-5 0 2.2834e+00 1.0549e+00 7.0363e-01 5.0605e-01]'*1e11; - md.materials.viscosity=[0 0 2.0000e+00 1.0000e+00 1.0000e+00 1.0000e+25]'*1e21; - md.materials.lame_lambda=md.materials.lame_mu*0+5e14; - md.materials.issolid=[1 0 1 1 1 1]'; - md.materials.isburgers=zeros(md.materials.numlayers,1); + md.love.allow_layer_deletion=1; + md.love.frequencies=0; + md.love.nfreq=length(md.love.frequencies); -md.love.love_kernels=1; -md.love.allow_layer_deletion=1; -md.love.frequencies=([0]*2*pi)'/cst; -md.love.nfreq=length(md.love.frequencies); -md.love.sh_nmax=2; + md.love.sh_nmin = 2; + md.love.sh_nmax = 2; + md.love.love_kernels=1; -md.materials.burgers_mu=md.materials.lame_mu; -md.materials.burgers_viscosity=md.materials.viscosity; + md.miscellaneous.name='kernels'; + md.cluster=generic('name',oshostname(),'np',1); + md.verbose=verbose('111111101'); + + md=solve(md,'lv'); -md=solve(md,'lv'); + % extract love kernels {{{ + y1=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,1))); + y2=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,2))); + y3=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,3))); + y4=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,4))); + y5=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,5))); + y6=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,6))); -%Fields and tolerances to track changes -%loading love numbers -field_names = {'LoveH_loading_elastic','LoveK_loading_elastic','LoveL_loading_elastic','LoveKernels_degree1','LoveKernels_degree2'}; -field_tolerances={1e-10,1e-10,1e-10,4e-10,1e-10}; + y1_tidal_degree2_interface1 = y1(3,1); + y1_tidal_degree2_interface2 = y1(3,2); + y1_tidal_degree2_interface7 = y1(3,7); + y1_tidal_degree2_interface10 = y1(3,10); + y1_tidal_degree2_interface11 = y1(3,11); + + y2_tidal_degree2_interface1 = y2(3,1); + y2_tidal_degree2_interface2 = y2(3,2); + y2_tidal_degree2_interface7 = y2(3,7); + y2_tidal_degree2_interface10 = y2(3,10); + y2_tidal_degree2_interface11 = y2(3,11); + + y3_tidal_degree2_interface1 = y3(3,1); + y3_tidal_degree2_interface2 = y3(3,2); + y3_tidal_degree2_interface7 = y3(3,7); + y3_tidal_degree2_interface10 = y3(3,10); + y3_tidal_degree2_interface11 = y3(3,11); + + y4_tidal_degree2_interface1 = y4(3,1); + y4_tidal_degree2_interface2 = y4(3,2); + y4_tidal_degree2_interface7 = y4(3,7); + y4_tidal_degree2_interface10 = y4(3,10); + y4_tidal_degree2_interface11 = y4(3,11); + + y5_tidal_degree2_interface1 = y5(3,1); + y5_tidal_degree2_interface2 = y5(3,2); + y5_tidal_degree2_interface7 = y5(3,7); + y5_tidal_degree2_interface10 = y5(3,10); + y5_tidal_degree2_interface11 = y5(3,11); + + y6_tidal_degree2_interface1 = y6(3,1); + y6_tidal_degree2_interface2 = y6(3,2); + y6_tidal_degree2_interface7 = y6(3,7); + y6_tidal_degree2_interface10 = y6(3,10); + y6_tidal_degree2_interface11 = y6(3,11); + % }}} + +% for surface load. + + md.love.forcing_type = 11; + + md=solve(md,'lv'); + + % extract love kernels {{{ + y1=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,1))); + y2=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,2))); + y3=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,3))); + y4=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,4))); + y5=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,5))); + y6=squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(:,1,:,6))); + + y1_load_degree2_interface1 = y1(3,1); + y1_load_degree2_interface2 = y1(3,2); + y1_load_degree2_interface7 = y1(3,7); + y1_load_degree2_interface10 = y1(3,10); + y1_load_degree2_interface11 = y1(3,11); + + y2_load_degree2_interface1 = y2(3,1); + y2_load_degree2_interface2 = y2(3,2); + y2_load_degree2_interface7 = y2(3,7); + y2_load_degree2_interface10 = y2(3,10); + y2_load_degree2_interface11 = y2(3,11); + + y3_load_degree2_interface1 = y3(3,1); + y3_load_degree2_interface2 = y3(3,2); + y3_load_degree2_interface7 = y3(3,7); + y3_load_degree2_interface10 = y3(3,10); + y3_load_degree2_interface11 = y3(3,11); + + y4_load_degree2_interface1 = y4(3,1); + y4_load_degree2_interface2 = y4(3,2); + y4_load_degree2_interface7 = y4(3,7); + y4_load_degree2_interface10 = y4(3,10); + y4_load_degree2_interface11 = y4(3,11); + + y5_load_degree2_interface1 = y5(3,1); + y5_load_degree2_interface2 = y5(3,2); + y5_load_degree2_interface7 = y5(3,7); + y5_load_degree2_interface10 = y5(3,10); + y5_load_degree2_interface11 = y5(3,11); + + y6_load_degree2_interface1 = y6(3,1); + y6_load_degree2_interface2 = y6(3,2); + y6_load_degree2_interface7 = y6(3,7); + y6_load_degree2_interface10 = y6(3,10); + y6_load_degree2_interface11 = y6(3,11); + % }}} + +field_names = {... + 'y1_tidal_degree2_interface1','y1_tidal_degree2_interface2','y1_tidal_degree2_interface7','y1_tidal_degree2_interface10','y1_tidal_degree2_interface11',... + 'y2_tidal_degree2_interface1','y2_tidal_degree2_interface2','y2_tidal_degree2_interface7','y2_tidal_degree2_interface10','y2_tidal_degree2_interface11',... + 'y3_tidal_degree2_interface1','y3_tidal_degree2_interface2','y3_tidal_degree2_interface7','y3_tidal_degree2_interface10','y3_tidal_degree2_interface11',... + 'y4_tidal_degree2_interface1','y4_tidal_degree2_interface2','y4_tidal_degree2_interface7','y4_tidal_degree2_interface10','y4_tidal_degree2_interface11',... + 'y5_tidal_degree2_interface1','y5_tidal_degree2_interface2','y5_tidal_degree2_interface7','y5_tidal_degree2_interface10','y5_tidal_degree2_interface11',... + 'y6_tidal_degree2_interface1','y6_tidal_degree2_interface2','y6_tidal_degree2_interface7','y6_tidal_degree2_interface10','y6_tidal_degree2_interface11',... + 'y1_load_degree2_interface1','y1_load_degree2_interface2','y1_load_degree2_interface7','y1_load_degree2_interface10','y1_load_degree2_interface11',... + 'y2_load_degree2_interface1','y2_load_degree2_interface2','y2_load_degree2_interface7','y2_load_degree2_interface10','y2_load_degree2_interface11',... + 'y3_load_degree2_interface1','y3_load_degree2_interface2','y3_load_degree2_interface7','y3_load_degree2_interface10','y3_load_degree2_interface11',... + 'y4_load_degree2_interface1','y4_load_degree2_interface2','y4_load_degree2_interface7','y4_load_degree2_interface10','y4_load_degree2_interface11',... + 'y5_load_degree2_interface1','y5_load_degree2_interface2','y5_load_degree2_interface7','y5_load_degree2_interface10','y5_load_degree2_interface11',... + 'y6_load_degree2_interface1','y6_load_degree2_interface2','y6_load_degree2_interface7','y6_load_degree2_interface10','y6_load_degree2_interface11',... + }; +field_tolerances={... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-10,1e-10,1e-10,1e-10,1e-10,... + }; field_values={... - (md.results.LoveSolution.LoveHr(:,1)),... - (md.results.LoveSolution.LoveKr(:,1)),... - (md.results.LoveSolution.LoveLr(:,1)),... - squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(2,1,:,:))),... - squeeze(cell2mat(md.results.LoveSolution.LoveKernelsReal(3,1,:,:))),... + y1_tidal_degree2_interface1, y1_tidal_degree2_interface2, y1_tidal_degree2_interface7, y1_tidal_degree2_interface10,y1_tidal_degree2_interface11,... + y2_tidal_degree2_interface1, y2_tidal_degree2_interface2, y2_tidal_degree2_interface7, y2_tidal_degree2_interface10,y2_tidal_degree2_interface11,... + y3_tidal_degree2_interface1, y3_tidal_degree2_interface2, y3_tidal_degree2_interface7, y3_tidal_degree2_interface10,y3_tidal_degree2_interface11,... + y4_tidal_degree2_interface1, y4_tidal_degree2_interface2, y4_tidal_degree2_interface7, y4_tidal_degree2_interface10,y4_tidal_degree2_interface11,... + y5_tidal_degree2_interface1, y5_tidal_degree2_interface2, y5_tidal_degree2_interface7, y5_tidal_degree2_interface10,y5_tidal_degree2_interface11,... + y6_tidal_degree2_interface1, y6_tidal_degree2_interface2, y6_tidal_degree2_interface7, y6_tidal_degree2_interface10,y6_tidal_degree2_interface11,... + y1_load_degree2_interface1, y1_load_degree2_interface2, y1_load_degree2_interface7, y1_load_degree2_interface10,y1_load_degree2_interface11,... + y2_load_degree2_interface1, y2_load_degree2_interface2, y2_load_degree2_interface7, y2_load_degree2_interface10,y2_load_degree2_interface11,... + y3_load_degree2_interface1, y3_load_degree2_interface2, y3_load_degree2_interface7, y3_load_degree2_interface10,y3_load_degree2_interface11,... + y4_load_degree2_interface1, y4_load_degree2_interface2, y4_load_degree2_interface7, y4_load_degree2_interface10,y4_load_degree2_interface11,... + y5_load_degree2_interface1, y5_load_degree2_interface2, y5_load_degree2_interface7, y5_load_degree2_interface10,y5_load_degree2_interface11,... + y6_load_degree2_interface1, y6_load_degree2_interface2, y6_load_degree2_interface7, y6_load_degree2_interface10,y6_load_degree2_interface11,... }; +