Changeset 22079
- Timestamp:
- 09/13/17 15:31:03 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test2084.m
r22068 r22079 36 36 %loading love numbers 37 37 field_names ={'LoveH_loading_elastic','LoveK_loading_elastic','LoveL_loading_elastic'}; 38 field_tolerances={1e-1 3,1e-13,1e-13};38 field_tolerances={1e-10,1e-10,1e-10}; 39 39 field_values={... 40 40 (md.results.LoveSolution.LoveHr(:,1)),... … … 54 54 %loading love numbers 55 55 field_names ={field_names{:},'LoveH_loading_realpart','LoveK_loading_realpart','LoveL_loading_realpart','LoveH_loading_imagpart','LoveK_loading_imagpart','LoveL_loading_imagpart'}; 56 field_tolerances={field_tolerances{:},1e-1 3,1e-13,1e-13,1e-13,1e-13,1e-13};56 field_tolerances={field_tolerances{:},1e-10,1e-10,1e-10,1e-10,1e-10,1e-10}; 57 57 field_values={field_values{:},... 58 58 (md.results.LoveSolution.LoveHr(:,:)),... … … 66 66 md.love.forcing_type=9; 67 67 md.love.sh_nmin=2 68 md.love.frequencies=([0 1e-4 1e-3 1e-2 1e-1 1 -1e-4 -1e-3 -1e-2 -1e-1 -1]*2*pi)'/cst;68 md.love.frequencies=([0 1e-4 1e-3 1e-2 1e-1 1 1e-4 -1e-3 -1e-2 -1e-1 -1]*2*pi)'/cst; 69 69 md.love.nfreq=length(md.love.frequencies); 70 70 md=solve(md,'lv'); 71 71 %tidal love numbers 72 72 field_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'}; 73 field_tolerances={field_tolerances{:},1e-1 3,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};73 field_tolerances={field_tolerances{:},1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10}; 74 74 field_values={field_values{:},... 75 75 (md.results.LoveSolution.LoveHr(:,1)),... … … 113 113 % 114 114 %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'}; 115 %field_tolerances={field_tolerances{:},1e-1 3,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};115 %field_tolerances={field_tolerances{:},1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10}; 116 116 %field_values={field_values{:},... 117 117 % (md.results.LoveSolution.LoveHr(:,1)),... … … 147 147 148 148 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'}; 149 field_tolerances={field_tolerances{:},1e-1 3,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};149 field_tolerances={field_tolerances{:},1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10}; 150 150 field_values={field_values{:},... 151 151 (md.results.LoveSolution.LoveHr(:,1)),...
Note:
See TracChangeset
for help on using the changeset viewer.