Index: ../trunk-jpl/test/NightlyRun/test2085.m =================================================================== --- ../trunk-jpl/test/NightlyRun/test2085.m (nonexistent) +++ ../trunk-jpl/test/NightlyRun/test2085.m (revision 22384) @@ -0,0 +1,44 @@ + +%Test Name: LovenumberstAtDepth. +%Same as test #1 of test2084.m + +md=model(); +md.cluster=generic('name',oshostname(),'np',1); + +md.materials=materials('litho'); +md.miscellaneous.name='FourierLoveTest'; + +md.verbose=verbose('111111101'); +cst=365.25*24*3600*1000; + + 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.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.materials.burgers_mu=md.materials.lame_mu; +md.materials.burgers_viscosity=md.materials.viscosity; + +md=solve(md,'lv'); + +%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,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,:,:))),... + }; Index: ../trunk-jpl/test/Archives/Archive2085.arch =================================================================== Cannot display: file marked as a binary type. svn:mime-type = application/octet-stream Index: ../trunk-jpl/test/Archives/Archive2085.arch =================================================================== --- ../trunk-jpl/test/Archives/Archive2085.arch (nonexistent) +++ ../trunk-jpl/test/Archives/Archive2085.arch (revision 22384) Property changes on: ../trunk-jpl/test/Archives/Archive2085.arch ___________________________________________________________________ Added: svn:mime-type ## -0,0 +1 ## +application/octet-stream \ No newline at end of property