source: issm/oecreview/Archive/21724-22754/ISSM-22059-22060.diff@ 22755

Last change on this file since 22755 was 22755, checked in by Mathieu Morlighem, 7 years ago

CHG: added 21724-22754
File size: 17.0 KB

  • ../trunk-jpl/test/NightlyRun/test2084.m

     
     1
     2%Test Name: GiaCaron.
     3%Forward Love number solution for a viscoelastic earth, model M3-L70-V01 from Spada, G., Barletta, V. R., Klemann, V., Riva, R. E. M., Martinec, Z., Gasperini, P., Lund, B., Wolf, D., Vermeersen, L. L. A. and King, M. A. (2011), A benchmark study for glacial isostatic adjustment codes. Geophysical Journal International, 185: 106–132. doi:10.1111/j.1365-246X.2011.04952.x
     4
     5md=model();
     6md.cluster=generic('name',oshostname(),'np',1);
     7
     8md.materials=materials('litho');
     9md.miscellaneous.name='FourierLoveTest';
     10
     11md.verbose=verbose('111111101');
     12cst=365.25*24*3600*1000;
     13
     14        md.materials.numlayers=6;
     15        md.materials.radius =  [10 1222.5 3.4800e+03   5.7010e+03   5.9510e+03   6.3010e+03   6.3710e+03]'*1e3;
     16        md.materials.density=  [1.0750e4 1.0750e+04   4.9780e+03   3.8710e+03   3.4380e+03   3.0370e+03]';
     17        md.materials.lame_mu=  [1e-5         0   2.2834e+00   1.0549e+00   7.0363e-01   5.0605e-01]'*1e11;
     18        md.materials.viscosity=[0             0   2.0000e+00   1.0000e+00   1.0000e+00   1.0000e+25]'*1e21;
     19        md.materials.lame_lambda=md.materials.lame_mu*0+5e14;
     20        md.materials.issolid=[1 0 1 1 1 1]';
     21        md.materials.isburgers=zeros(md.materials.numlayers,1);
     22
     23md.love.allow_layer_deletion=1;
     24md.love.frequencies=([0]*2*pi)'/cst;
     25md.love.nfreq=length(md.love.frequencies);
     26md.love.sh_nmax=256;
     27
     28md.materials.burgers_mu=md.materials.lame_mu;
     29md.materials.burgers_viscosity=md.materials.viscosity;
     30
     31md=solve(md,'lv');
     32
     33
     34
     35%Fields and tolerances to track changes
     36%loading love numbers
     37field_names     ={'LoveH_loading_elastic','LoveK_loading_elastic','LoveL_loading_elastic'};
     38field_tolerances={1e-13,1e-13,1e-13};
     39field_values={...
     40        (md.results.LoveSolution.LoveHr(:,1)),...
     41        (md.results.LoveSolution.LoveKr(:,1)),...
     42        (md.results.LoveSolution.LoveLr(:,1)),...
     43        };
     44
     45
     46md.love.frequencies=([1e-3 1e-2 1e-1 1 -1e-3 -1e-2 -1e-1 -1]*2*pi)'/cst;
     47md.love.nfreq=length(md.love.frequencies);
     48md.love.sh_nmax=256;
     49md.materials.burgers_mu=md.materials.lame_mu;
     50md.materials.burgers_viscosity=md.materials.viscosity;
     51
     52md=solve(md,'lv');
     53%Fields and tolerances to track changes
     54%loading love numbers
     55field_names     ={field_names{:},'LoveH_loading_realpart','LoveK_loading_realpart','LoveL_loading_realpart','LoveH_loading_imagpart','LoveK_loading_imagpart','LoveL_loading_imagpart'};
     56field_tolerances={field_tolerances{:},1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     57field_values={field_values{:},...
     58        (md.results.LoveSolution.LoveHr(:,:)),...
     59        (md.results.LoveSolution.LoveKr(:,:)),...
     60        (md.results.LoveSolution.LoveLr(:,:)),...
     61        (md.results.LoveSolution.LoveHi(:,:)),...
     62        (md.results.LoveSolution.LoveKi(:,:)),...
     63        (md.results.LoveSolution.LoveLi(:,:)),...
     64        };
     65
     66
     67return
     68md.love.forcing_type=9;
     69md.love.sh_nmin=2
     70md.love.frequencies=([0 1e-4 1e-3 1e-2 1e-1 1 -1e-4 -1e-3 -1e-2 -1e-1 -1]*2*pi)'/cst;
     71md.love.nfreq=length(md.love.frequencies);
     72md=solve(md,'lv');
     73%tidal love numbers
     74field_names     ={field_names{:},'LoveH_tidal_elastic','LoveK_tidal_elastic','LoveL_tidal_elastic','LoveH_tidal_realpart','LoveK_tidal_realpart','LoveL_tidal_realpart','LoveH_tidal_imagpart','LoveK_tidal_imagpart','LoveL_tidal_imagpart'};
     75field_tolerances={field_tolerances{:},1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     76field_values={field_values{:},...
     77        (md.results.LoveSolution.LoveHr(:,1)),...
     78        (md.results.LoveSolution.LoveKr(:,1)),...
     79        (md.results.LoveSolution.LoveLr(:,1)),...
     80        (md.results.LoveSolution.LoveHr(:,2:end)),...
     81        (md.results.LoveSolution.LoveKr(:,2:end)),...
     82        (md.results.LoveSolution.LoveLr(:,2:end)),...
     83        (md.results.LoveSolution.LoveHi(:,2:end)),...
     84        (md.results.LoveSolution.LoveKi(:,2:end)),...
     85        (md.results.LoveSolution.LoveLi(:,2:end)),...
     86        };
     87
     88
     89
     90 %Many layers PREM-based model
     91%data=load('../Data/PREM_500layers');
     92%md.love.sh_nmin=1;
     93%md.materials.radius=data(2:end-2,1);
     94%md.materials.density=data(3:end-2,2);
     95%md.materials.lame_lambda=data(3:end-2,3);
     96%md.materials.lame_mu=data(3:end-2,4);
     97%md.materials.issolid=data(3:end-2,4)>0;
     98%ind=find(md.materials.issolid==0);
     99%md.materials.density(ind(1))=sum((md.materials.radius(ind+1).^3-md.materials.radius(ind).^3).*md.materials.density(ind))/(md.materials.radius(ind(end)+1).^3-md.materials.radius(ind(1)+1).^3);
     100%md.materials.lame_lambda(ind(1))=sum((md.materials.radius(ind+1).^3-md.materials.radius(ind).^3).*md.materials.lame_lambda(ind))/(md.materials.radius(ind(end)+1).^3-md.materials.radius(ind(1)+1).^3);
     101%md.materials.lame_mu(ind(1))=sum((md.materials.radius(ind+1).^3-md.materials.radius(ind).^3).*md.materials.lame_mu(ind))/(md.materials.radius(ind(end)+1).^3-md.materials.radius(ind(1)).^3);
     102%md.materials.radius(ind(2:end)+1)=[];
     103%md.materials.density(ind(2:end))=[];
     104%md.materials.lame_lambda(ind(2:end))=[];
     105%md.materials.lame_mu(ind(2:end))=[];
     106%md.materials.issolid(ind(2:end))=[];
     107%md.materials.viscosity=10.^interp1([0 3479e3 3480e3 3680e3 5720e3 5800e3 6270e3 6371e3], log10([1e8 1e8 5e21 1e23 1e22 1e20 1e21 1e40]), md.materials.radius(2:end),'PCHIP');
     108%md.materials.viscosity=md.materials.viscosity.*md.materials.issolid;
     109%md.materials.burgers_mu=md.materials.lame_mu;
     110%md.materials.burgers_viscosity=md.materials.viscosity;
     111%md.materials.isburgers=md.materials.issolid*0;
     112%md.love.forcing_type=11;
     113%md.materials.numlayers=length(md.materials.viscosity);
     114%md=solve(md,'lv');
     115%
     116%field_names     ={field_names{:},'LoveH_loadingVSS96_elastic','LoveK_loadingVSS96_elastic','LoveL_loadingVSS96_elastic','LoveH_loadingVSS96_realpart','LoveK_loadingVSS96_realpart','LoveL_loadingVSS96_realpart','LoveH_loadingVSS96_imagpart','LoveK_loadingVSS96_imagpart','LoveL_loadingVSS96_imagpart'};
     117%field_tolerances={field_tolerances{:},1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     118%field_values={field_values{:},...
     119%       (md.results.LoveSolution.LoveHr(:,1)),...
     120%       (md.results.LoveSolution.LoveKr(:,1)),...
     121%       (md.results.LoveSolution.LoveLr(:,1)),...
     122%       (md.results.LoveSolution.LoveHr(:,2:end)),...
     123%       (md.results.LoveSolution.LoveKr(:,2:end)),...
     124%       (md.results.LoveSolution.LoveLr(:,2:end)),...
     125%       (md.results.LoveSolution.LoveHi(:,2:end)),...
     126%       (md.results.LoveSolution.LoveKi(:,2:end)),...
     127%       (md.results.LoveSolution.LoveLi(:,2:end)),...
     128%       };
     129
     130
     131%Model VSS96 from Vermeersen, L.L.A., Sabadini, R. & Spada, G., 1996a. Analytical visco-elastic relaxation models, Geophys. Res. Lett., 23, 697–700.
     132md.materials.radius=[10, 1222.5, 3480., 3600., 3630.5, 3700., 3900., 4000., 4200., 4300., 4500., 4600., 4800., 4900., 5100., 5200., 5400., 5500., 5600.5, 5650., 5701., 5736., 5771.5, 5821., 5951., 5970.5, 6016., 6061., 6150.5, 6151.5, 6251., 6371.]'*1e3;
     133md.materials.lame_mu=[1e-5, 0., 2.933, 2.8990002, 2.8550003, 2.7340002, 2.675, 2.559, 2.502, 2.388, 2.331, 2.215, 2.157, 2.039, 1.979, 1.8560001, 1.794, 1.73, 1.639, 1.2390001, 1.224, 1.21, 1.128, 0.97700006, 0.906, 0.79, 0.773, 0.741, 0.656, 0.665, 0.602]'*1e11;
     134md.materials.density=[10925., 10925., 5506.42, 5491.45, 5456.57, 5357.06, 5307.24, 5207.13, 5156.69, 5054.69, 5002.99, 4897.83, 4844.22, 4734.6, 4678.44, 4563.07, 4503.72, 4443.16, 4412.41, 3992.14, 3983.99, 3975.84, 3912.82, 3786.78, 3723.78, 3516.39, 3489.51, 3435.78, 3359.5, 3367.1, 3184.3]';
     135md.materials.viscosity=[0., 0., 8E+21, 8.5E+21, 9E+21, 3.E+22, 4.E+22, 5.E+22, 6.E+22, 5.E+22, 4.5E+22, 3.E+22, 2.5E+22, 1.8E+22, 1.3E+22, 8E+21, 7+21, 6.5E+21, 6.E+21, 5.5E+21, 5.E+21, 4.5E+21, 4E+21, 2.5E+21, 2.E+21, 1.5E+21, 10.E+20, 6.E+20, 5.5E+20, 2.E+20, 1.E40]';
     136md.materials.lame_lambda=md.materials.lame_mu*0+5e14;
     137md.materials.issolid=md.materials.lame_mu>0;
     138md.materials.numlayers=length(md.materials.lame_mu);
     139md.materials.burgers_mu=md.materials.lame_mu;
     140md.materials.burgers_viscosity=md.materials.viscosity;
     141md.materials.isburgers=md.materials.issolid*0;
     142md.love.forcing_type=11;
     143md.love.sh_nmin=1;
     144md.love.sh_nmax=100;
     145md=solve(md,'lv');
     146md.love.frequencies=([0 1e-4 1e-2 1 -1e-4 -1e-2 -1]*2*pi)'/cst;
     147md.love.nfreq=length(md.love.frequencies);
     148
     149
     150field_names     ={field_names{:},'LoveH_loadingVSS96_elastic','LoveK_loadingVSS96_elastic','LoveL_loadingVSS96_elastic','LoveH_loadingVSS96_realpart','LoveK_loadingVSS96_realpart','LoveL_loadingVSS96_realpart','LoveH_loadingVSS96_imagpart','LoveK_loadingVSS96_imagpart','LoveL_loadingVSS96_imagpart'};
     151field_tolerances={field_tolerances{:},1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     152field_values={field_values{:},...
     153        (md.results.LoveSolution.LoveHr(:,1)),...
     154        (md.results.LoveSolution.LoveKr(:,1)),...
     155        (md.results.LoveSolution.LoveLr(:,1)),...
     156        (md.results.LoveSolution.LoveHr(:,2:end)),...
     157        (md.results.LoveSolution.LoveKr(:,2:end)),...
     158        (md.results.LoveSolution.LoveLr(:,2:end)),...
     159        (md.results.LoveSolution.LoveHi(:,2:end)),...
     160        (md.results.LoveSolution.LoveKi(:,2:end)),...
     161        (md.results.LoveSolution.LoveLi(:,2:end)),...
     162        };

  • ../trunk-jpl/test/NightlyRun/test2084.py

     
     1
     2#Test Name: GiaCaron. Forward Love number solution for a viscoelastic earth, model M3-L70-V01 from Spada, G., Barletta, V. R., Klemann, V., Riva, R. E. M., Martinec, Z., Gasperini, P., Lund, B., Wolf, D., Vermeersen, L. L. A. and King, M. A. (2011), A benchmark study for glacial isostatic adjustment codes. Geophysical Journal International, 185: 106–132. doi:10.1111/j.1365-246X.2011.04952.x
     3from model import *
     4from socket import gethostname
     5from solve import *
     6from numpy import *
     7md=model()
     8md.cluster=generic('name',oshostname(),'np',1)
     9md.materials=materials('litho')
     10md.miscellaneous.name='FourierLoveTest'
     11md.verbose=verbose('111111101')
     12cst=365.25*24*3600*1000
     13md.materials.numlayers=6
     14md.materials.radius =[10, 1222.5, 3.4800e+03,5.7010e+03,5.9510e+03,6.3010e+03,6.3710e+03]*1e3
     15md.materials.density=[1.0750e4, 1.0750e+04,4.9780e+03,3.8710e+03,3.4380e+03,3.0370e+03]
     16md.materials.lame_mu=[1e-5,0,2.2834e+00,1.0549e+00,7.0363e-01,5.0605e-01]*1e11
     17md.materials.viscosity=[0,0,2.0000e+00,1.0000e+00,1.0000e+00,1.0000e+25]*1e21
     18md.materials.lame_lambda=md.materials.lame_mu*0+5e14
     19md.materials.issolid=[1, 0, 1, 1, 1, 1]
     20md.materials.isburgers=zeros(md.materials.numlayers,1)
     21md.love.allow_layer_deletion=1
     22
     23md.love.frequencies=([0]*2*pi)/cst
     24md.love.nfreq=len(md.love.frequencies)
     25md.love.sh_nmax=256
     26md.materials.burgers_mu=md.materials.lame_mu
     27md.materials.burgers_viscosity=md.materials.viscosity
     28md=solve(md,'lv')
     29#Fields and tolerances to track changes
     30#loading love numbers
     31field_names=['LoveHe','LoveKe','LoveLe']
     32field_tolerances=[1e-13,1e-13,1e-13]
     33field_values=[\
     34(md.results.LoveSolution.LoveHr[:][0]),\
     35(md.results.LoveSolution.LoveKr[:][0]),\
     36(md.results.LoveSolution.LoveLr[:][0]),\
     37]
     38
     39
     40md.love.frequencies=([1e-3, 1e-2, 1e-1, 1, -1e-3, -1e-2, -1e-1, -1]*2*pi)/cst
     41md.love.nfreq=len(md.love.frequencies)
     42md.love.sh_nmax=256
     43md.materials.burgers_mu=md.materials.lame_mu
     44md.materials.burgers_viscosity=md.materials.viscosity
     45md=solve(md,'lv')
     46#Fields and tolerances to track changes
     47#loading love numbers
     48field_names=concatenate([field_names,['LoveHlr','LoveKlr','LoveLlr','LoveHli','LoveKli','LoveLli']])
     49field_tolerances=concatenate([field_tolerances,[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]])
     50field_values=concatenate([field_values,[\
     51(md.results.LoveSolution.LoveHr[:][:]),\
     52(md.results.LoveSolution.LoveKr[:][:]),\
     53(md.results.LoveSolution.LoveLr[:][:])),\
     54(md.results.LoveSolution.LoveHi[:][:]),\
     55(md.results.LoveSolution.LoveKi[:][:]),\
     56(md.results.LoveSolution.LoveLi[:][:]),\
     57]])
     58
     59
     60
     61md.love.forcing_type=9
     62md.love.sh_nmin=2
     63md.love.frequencies=([0, 1e-4, 1e-3, 1e-2, 1e-1, 1, -1e-4, -1e-3, -1e-2, -1e-1, -1]*2*pi)/cst
     64md.love.nfreq=len(md.love.frequencies)
     65md=solve(md,'lv')
     66#tidal love numbers
     67field_names=concatenate([field_names,['LoveHtr','LoveKtr','LoveLtr','LoveHti','LoveKti','LoveLti']])
     68field_tolerances=concatenate([field_tolerances,[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]])
     69field_values=concatenate([field_values,[\
     70(md.results.LoveSolution.LoveHr[:][0]),\
     71(md.results.LoveSolution.LoveKr[:][0]),\
     72(md.results.LoveSolution.LoveLr[:][0]),\
     73(md.results.LoveSolution.LoveHr[:][1:]),\
     74(md.results.LoveSolution.LoveKr[:][1:]),\
     75(md.results.LoveSolution.LoveLr[:][1:]),\
     76(md.results.LoveSolution.LoveHi[:][1:]),\
     77(md.results.LoveSolution.LoveKi[:][1:]),\
     78(md.results.LoveSolution.LoveLi[:][1:]),\
     79]])
     80
     81
     82
     83#Many layers PREM-based model
     84#data=load('../Data/PREM_500layers')
     85#md.love.sh_nmin=1
     86#md.materials.radius=data(2:end-2,1)
     87#md.materials.density=data(3:end-2,2)
     88#md.materials.lame_lambda=data(3:end-2,3)
     89#md.materials.lame_mu=data(3:end-2,4)
     90#md.materials.issolid=data(3:end-2,4)>0
     91#ind=find(md.materials.issolid==0)
     92#md.materials.density(ind(1))=sum((md.materials.radius(ind+1).^3-md.materials.radius(ind).^3).*md.materials.density(ind))/(md.materials.radius(ind(end)+1).^3-md.materials.radius(ind(1)+1).^3)
     93#md.materials.lame_lambda(ind(1))=sum((md.materials.radius(ind+1).^3-md.materials.radius(ind).^3).*md.materials.lame_lambda(ind))/(md.materials.radius(ind(end)+1).^3-md.materials.radius(ind(1)+1).^3)
     94#md.materials.lame_mu(ind(1))=sum((md.materials.radius(ind+1).^3-md.materials.radius(ind).^3).*md.materials.lame_mu(ind))/(md.materials.radius(ind(end)+1).^3-md.materials.radius(ind(1)).^3)
     95#md.materials.radius(ind(2:end)+1)=[]
     96#md.materials.density(ind(2:end))=[]
     97#md.materials.lame_lambda(ind(2:end))=[]
     98#md.materials.lame_mu(ind(2:end))=[]
     99#md.materials.issolid(ind(2:end))=[]
     100#md.materials.viscosity=10.^interp1([0 3479e3 3480e3 3680e3 5720e3 5800e3 6270e3 6371e3], log10([1e8 1e8 5e21 1e23 1e22 1e20 1e21 1e40]), md.materials.radius(2:end),'PCHIP')
     101#md.materials.viscosity=md.materials.viscosity.*md.materials.issolid
     102#md.materials.burgers_mu=md.materials.lame_mu
     103#md.materials.burgers_viscosity=md.materials.viscosity
     104#md.materials.isburgers=md.materials.issolid*0
     105#md.love.forcing_type=11
     106#md.materials.numlayers=len(md.materials.viscosity)
     107#md=solve(md,'lv')
     108#
     109#field_names=[field_names,'LoveHmr','LoveKmr','LoveLmr','LoveHmi','LoveKmi','LoveLmi']
     110#field_tolerances=[field_tolerances,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
     111#field_values=[field_values,\
     112#       (md.results.LoveSolution.LoveHr[:][0]),\
     113#       (md.results.LoveSolution.LoveKr[:][0]),\
     114#       (md.results.LoveSolution.LoveLr[:][0]),\
     115#       (md.results.LoveSolution.LoveHr[:][1:]),\
     116#       (md.results.LoveSolution.LoveKr[:][1:]),\
     117#       (md.results.LoveSolution.LoveLr[:][1:]),\
     118#       (md.results.LoveSolution.LoveHi[:][1:]),\
     119#       (md.results.LoveSolution.LoveKi[:][1:]),\
     120#       (md.results.LoveSolution.LoveLi[:][1:]),\
     121#       ]
     122#Model VSS96 from Vermeersen, L.L.A., Sabadini, R. & Spada, G., 1996a. Analytical visco-elastic relaxation models, Geophys. Res. Lett., 23, 697–700.
     123md.materials.radius=[10, 1222.5, 3480., 3600., 3630.5, 3700., 3900., 4000., 4200., 4300., 4500., 4600., 4800., 4900., 5100., 5200., 5400., 5500., 5600.5, 5650., 5701., 5736., 5771.5, 5821., 5951., 5970.5, 6016., 6061., 6150.5, 6151.5, 6251., 6371.]*1e3
     124md.materials.lame_mu=[1e-5, 0., 2.933, 2.8990002, 2.8550003, 2.7340002, 2.675, 2.559, 2.502, 2.388, 2.331, 2.215, 2.157, 2.039, 1.979, 1.8560001, 1.794, 1.73, 1.639, 1.2390001, 1.224, 1.21, 1.128, 0.97700006, 0.906, 0.79, 0.773, 0.741, 0.656, 0.665, 0.602]*1e11
     125md.materials.density=[10925., 10925., 5506.42, 5491.45, 5456.57, 5357.06, 5307.24, 5207.13, 5156.69, 5054.69, 5002.99, 4897.83, 4844.22, 4734.6, 4678.44, 4563.07, 4503.72, 4443.16, 4412.41, 3992.14, 3983.99, 3975.84, 3912.82, 3786.78, 3723.78, 3516.39, 3489.51, 3435.78, 3359.5, 3367.1, 3184.3]
     126md.materials.viscosity=[0., 0., 8E+21, 8.5E+21, 9E+21, 3.E+22, 4.E+22, 5E+22, 6.E+22, 5E+22, 4.5E+22, 3.E+22, 2.5E+22, 1.8E+22, 1.3E+22, 8E+21, 7E+21, 6.5E+21, 6.E+21, 5.5E+21, 5.E+21, 55E+21, 45E+21, 2.5E+21, 27E+21, 1.5E+21, 10E+20, 6.E+20, 5.5E+20, 2.E+20, 1.E40]
     127md.materials.lame_lambda=md.materials.lame_mu*0+5e14
     128md.materials.issolid=ones(len(md.materials.mu))
     129md.materials.issolid[1]=0
     130md.materials.numlayers=len(md.materials.lame_mu)
     131md.materials.burgers_mu=md.materials.lame_mu
     132md.materials.burgers_viscosity=md.materials.viscosity
     133md.materials.isburgers=md.materials.issolid*0
     134md.love.forcing_type=11
     135md.love.sh_nmin=1
     136md.love.sh_nmax=100
     137md=solve(md,'lv')
     138md.love.frequencies=([0, 1e-4, 1e-2, 1, -1e-4, -1e-2, -1]*2*pi)/cst
     139md.love.nfreq=len(md.love.frequencies)
     140field_names=concatenate([field_names,['LoveHmr','LoveKmr','LoveLmr','LoveHmi','LoveKmi','LoveLmi']])
     141field_tolerances=concatenate([field_tolerances,[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]])
     142field_values=concatenate([field_values,[\
     143(md.results.LoveSolution.LoveHr[:][0]),\
     144(md.results.LoveSolution.LoveKr[:][0]),\
     145(md.results.LoveSolution.LoveLr[:][0]),\
     146(md.results.LoveSolution.LoveHr[:][1:]),\
     147(md.results.LoveSolution.LoveKr[:][1:]),\
     148(md.results.LoveSolution.LoveLr[:][1:]),\
     149(md.results.LoveSolution.LoveHi[:][1:]),\
     150(md.results.LoveSolution.LoveKi[:][1:]),\
     151(md.results.LoveSolution.LoveLi[:][1:]),\
     152]])
Note: See TracBrowser for help on using the repository browser.