Changeset 22401


Ignore:
Timestamp:
02/07/18 21:07:09 (7 years ago)
Author:
Mathieu Morlighem
Message:

CHG: cosmetics

Location:
issm/trunk-jpl/test/NightlyRun
Files:
2 edited

Legend:

Unmodified
Added
Removed

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

    r22081 r22401  
    11
    22%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
     3%Forward Love number solution for a viscoelastic earth, model M3-L70-V01 from
     4%Spada, G., Barletta, V. R., Klemann, V., Riva, R. E. M., Martinec, Z.,
     5%Gasperini, P., Lund, B., Wolf, D., Vermeersen, L. L. A. and King, M. A.
     6%(2011), A benchmark study for glacial isostatic adjustment codes. Geophysical
     7%Journal International, 185: 106--132. doi:10.1111/j.1365-246X.2011.04952.x
    48
    59md=model();
    … …  
    913md.miscellaneous.name='FourierLoveTest';
    1014
    11 md.verbose=verbose('111111101');
     15md.verbose=verbose('all');
    1216cst=365.25*24*3600*1000;
    1317
    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);
     18md.materials.numlayers=6;
     19md.materials.radius =  [10 1222.5 3.4800e+03   5.7010e+03   5.9510e+03   6.3010e+03   6.3710e+03]'*1e3;
     20md.materials.density=  [1.0750e4 1.0750e+04   4.9780e+03   3.8710e+03   3.4380e+03   3.0370e+03]';
     21md.materials.lame_mu=  [1e-5         0   2.2834e+00   1.0549e+00   7.0363e-01   5.0605e-01]'*1e11;
     22md.materials.viscosity=[0            0   2.0000e+00   1.0000e+00   1.0000e+00   1.0000e+25]'*1e21;
     23md.materials.lame_lambda=md.materials.lame_mu*0+5e14;
     24md.materials.issolid=[1 0 1 1 1 1]';
     25md.materials.isburgers=zeros(md.materials.numlayers,1);
    2226
    2327md.love.allow_layer_deletion=1;
    … …  
    3034
    3135md=solve(md,'lv');
    32 
    33 
    3436
    3537%Fields and tolerances to track changes
    … …  
    5153
    5254md=solve(md,'lv');
     55
    5356%Fields and tolerances to track changes
    54 %loading love numbers
    5557field_names     ={field_names{:},'LoveH_loading_realpart','LoveK_loading_realpart','LoveL_loading_realpart','LoveH_loading_imagpart','LoveK_loading_imagpart','LoveL_loading_imagpart'};
    5658field_tolerances={field_tolerances{:},1e-10,1e-10,1e-10,1e-10,1e-10,1e-10};
    … …  
    6567
    6668md.love.forcing_type=9;
    67 md.love.sh_nmin=2
     69md.love.sh_nmin=2;
    6870md.love.frequencies=([0 1e-3 1e-2 1e-1 1 -1e-3 -1e-2 -1e-1 -1]*2*pi)'/cst;
    6971md.love.nfreq=length(md.love.frequencies);
    7072md=solve(md,'lv');
     73
    7174%tidal love numbers
    7275field_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'};
    … …  
    8487        };
    8588
    86 
    87 
    88  %Many layers PREM-based model
     89%Many layers PREM-based model
    8990%data=load('../Data/PREM_500layers');
    9091%md.love.sh_nmin=1;
    … …  
    126127%       };
    127128
    128 
    129129%Model VSS96 from Vermeersen, L.L.A., Sabadini, R. & Spada, G., 1996a. Analytical visco-elastic relaxation models, Geophys. Res. Lett., 23, 697–700.
    130130md.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;
    … …  
    145145md.love.nfreq=length(md.love.frequencies);
    146146
    147 
    148147field_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'};
    149148field_tolerances={field_tolerances{:},1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10};

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

    r22096 r22401  
    1919md.miscellaneous.name='FourierLoveTest'
    2020
    21 md.verbose=verbose('111111101')
     21md.verbose=verbose('all')
    2222cst=365.25*24*3600*1000
    2323
    … …  
    4545md=solve(md,'lv')
    4646
    47 
    48 
    4947#Fields and tolerances to track changes
    5048#loading love numbers
    … …  
    6563
    6664md=solve(md,'lv')
    67 
    68 
    6965
    7066#Fields and tolerances to track changes
    … …  
    8985md=solve(md,'lv')
    9086
    91 
    92 
    9387#tidal love numbers, check
    9488field_names+=['LoveHtr','LoveKtr','LoveLtr','LoveHti','LoveKti','LoveLti']
    … …  
    10599np.array(md.results.LoveSolution.LoveLi)[:,1:]
    106100]
    107 
    108 
    109101
    110102#Many layers PREM-based model
    … …  
    189181md=solve(md,'lv')
    190182
    191 
    192 
    193183md.love.frequencies=(np.array([0,1e-3,1e-2,1,-1e-3,-1e-2,
    194184                                -1])*2*pi).reshape(-1,1)/cst
Note: See TracChangeset for help on using the changeset viewer.