Changeset 25290
- Timestamp:
- 07/18/20 23:19:58 (5 years ago)
- Location:
- issm/trunk-jpl/test
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test2051.m
r25134 r25290 1 %Test Name: GiaIvinsBenchmarksAB2dA 11 %Test Name: GiaIvinsBenchmarksAB2dA 2 2 % Benchmark experiments (Figure A2a Ivins and James, 1999, Geophys. J. Int.) 3 3 md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.); … … 8 8 md.gia.cross_section_shape=1; % for square-edged x-section 9 9 10 % % define loading history10 % evaluation time (termed start_time) 11 11 md.timestepping.start_time=2002100; % after 2 kyr of deglaciation 12 12 md.timestepping.final_time=2500000; % 2,500 kyr 13 14 %% define loading history {{{ 13 15 md.geometry.thickness=[... 14 16 [md.geometry.thickness*0.0; 0.0],... … … 18 20 [md.geometry.thickness*0.0; md.timestepping.start_time],... 19 21 ]; 22 % }}} 23 24 % find out elements that have zero loads throughout the loading history. 25 pos = find(sum(abs(md.geometry.thickness(1:end-1,:)),2)==0); 26 md.mask.ice_levelset(pos)=1; % no-ice. 27 28 md.cluster=generic('name',oshostname(),'np',3); 29 md.verbose=verbose('1111111'); 20 30 21 31 %% solve for GIA deflection 22 md.cluster=generic('name',oshostname(),'np',3);23 md.verbose=verbose('1111111');24 32 md=solve(md,'Gia'); 25 33 34 %Test Name: GiaIvinsBenchmarksAB2dA1 35 U_AB2dA1 = md.results.GiaSolution.UGia; 36 URate_AB2dA1 = md.results.GiaSolution.UGiaRate; 37 38 %Test Name: GiaIvinsBenchmarksAB2dA2 39 %% different evaluation time. {{{ 40 md.timestepping.start_time=2005100; % after 5 kyr of deglaciation 41 md.geometry.thickness(end,end) = md.timestepping.start_time; 42 43 md=solve(md,'Gia'); 44 45 U_AB2dA2 = md.results.GiaSolution.UGia; 46 URate_AB2dA2 = md.results.GiaSolution.UGiaRate; 47 % }}} 48 49 %Test Name: GiaIvinsBenchmarksAB2dA3 50 %% different evaluation time. {{{ 51 md.timestepping.start_time=2010100; % after 10 kyr of deglaciation 52 md.geometry.thickness(end,end) = md.timestepping.start_time; 53 54 md=solve(md,'Gia'); 55 56 U_AB2dA3 = md.results.GiaSolution.UGia; 57 URate_AB2dA3 = md.results.GiaSolution.UGiaRate; 58 % }}} 59 26 60 %Fields and tolerances to track changes 27 field_names ={'UGia','UGiaRate'}; 28 field_tolerances={1e-13,1e-13}; 29 field_values={... 30 (md.results.GiaSolution.UGia),... 31 (md.results.GiaSolution.UGiaRate),... 32 }; 61 field_names ={'U_AB2dA1','URate_AB2dA1','U_AB2dA2','URate_AB2dA2','U_AB2dA3','URate_AB2dA3'}; 62 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 63 field_values={U_AB2dA1,URate_AB2dA1,U_AB2dA2,URate_AB2dA2,U_AB2dA3,URate_AB2dA3}; 64 -
issm/trunk-jpl/test/NightlyRun/test2052.m
r25131 r25290 1 %Test Name: GiaIvinsBenchmarksAB2d A21 %Test Name: GiaIvinsBenchmarksAB2dC 2 2 % Benchmark experiments (Figure A2a Ivins and James, 1999, Geophys. J. Int.) 3 3 md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.); 4 4 md=setmask(md,'',''); 5 md=parameterize(md,'../Par/GiaIvinsBenchmarks AB.par');5 md=parameterize(md,'../Par/GiaIvinsBenchmarksCD.par'); 6 6 7 7 %% indicate what you want to compute 8 8 md.gia.cross_section_shape=1; % for square-edged x-section 9 9 10 % % define loading history11 md.timestepping.start_time= 2005100; % after 5 kyr of deglaciation10 % evaluation time (termed start_time) 11 md.timestepping.start_time=0.3; % for t \approx 0 kyr : to get eleastic response! 12 12 md.timestepping.final_time=2500000; % 2,500 kyr 13 14 %% define loading history {{{ 13 15 md.geometry.thickness=[... 14 16 [md.geometry.thickness*0.0; 0.0],... 15 [md.geometry.thickness; 1000],... 16 [md.geometry.thickness; 2000000],... 17 [md.geometry.thickness*0.0; 2000100],... 18 [md.geometry.thickness*0.0; md.timestepping.start_time],... 17 [md.geometry.thickness/2.0; 0.1],... 18 [md.geometry.thickness; 0.2],... 19 [md.geometry.thickness; md.timestepping.start_time],... 19 20 ]; 21 % }}} 22 23 % find out elements that have zero loads throughout the loading history. 24 pos = find(sum(abs(md.geometry.thickness(1:end-1,:)),2)==0); 25 md.mask.ice_levelset(pos)=1; % no-ice. 26 27 md.cluster=generic('name',oshostname(),'np',3); 28 md.verbose=verbose('1111111'); 20 29 21 30 %% solve for GIA deflection 22 md.cluster=generic('name',oshostname(),'np',3);23 md.verbose=verbose('1111111');24 31 md=solve(md,'Gia'); 25 32 33 %Test Name: GiaIvinsBenchmarksAB2dC1 34 U_AB2dC1 = md.results.GiaSolution.UGia; 35 URate_AB2dC1 = md.results.GiaSolution.UGiaRate; 36 37 %Test Name: GiaIvinsBenchmarksAB2dC2 38 %% different evaluation time. {{{ 39 md.timestepping.start_time=1000.3; % for t \approx 1 kyr 40 md.geometry.thickness(end,end) = md.timestepping.start_time; 41 42 md=solve(md,'Gia'); 43 44 U_AB2dC2 = md.results.GiaSolution.UGia; 45 URate_AB2dC2 = md.results.GiaSolution.UGiaRate; 46 % }}} 47 48 %Test Name: GiaIvinsBenchmarksAB2dC3 49 %% different evaluation time. {{{ 50 md.timestepping.start_time=2400000; % for t \approx \infty 51 md.geometry.thickness(end,end) = md.timestepping.start_time; 52 53 md=solve(md,'Gia'); 54 55 U_AB2dC3 = md.results.GiaSolution.UGia; 56 URate_AB2dC3 = md.results.GiaSolution.UGiaRate; 57 % }}} 58 26 59 %Fields and tolerances to track changes 27 field_names ={'UGia','UGiaRate'}; 28 field_tolerances={1e-13,1e-13}; 29 field_values={... 30 (md.results.GiaSolution.UGia),... 31 (md.results.GiaSolution.UGiaRate),... 32 }; 60 field_names ={'U_AB2dC1','URate_AB2dC1','U_AB2dC2','URate_AB2dC2','U_AB2dC3','URate_AB2dC3'}; 61 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 62 field_values={U_AB2dC1,URate_AB2dC1,U_AB2dC2,URate_AB2dC2,U_AB2dC3,URate_AB2dC3}; 63 -
issm/trunk-jpl/test/NightlyRun/test2052.py
r25131 r25290 1 #Test Name: GiaIvinsBenchmarksAB2d A21 #Test Name: GiaIvinsBenchmarksAB2dC1 2 2 from socket import gethostname 3 3 … … 11 11 12 12 13 #Benchmark experiments (Figure A2 aIvins and James, 1999, Geophys. J. Int.)13 #Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.) 14 14 md = triangle(model(), '../Exp/RoundFrontEISMINT.exp', 200000) 15 15 md = setmask(md, '', '') 16 md = parameterize(md, '../Par/GiaIvinsBenchmarks AB.py')16 md = parameterize(md, '../Par/GiaIvinsBenchmarksCD.py') 17 17 18 18 #indicate what you want to compute 19 md.gia.cross_section_shape = 1 # for square-edged x -section19 md.gia.cross_section_shape = 1 # for square-edged x-section 20 20 21 21 #define loading history 22 md.timestepping.start_time = 2005100 # after 5 kyr of deglaciation23 md.timestepping.final_time = 2500000 # 2, 22 md.timestepping.start_time = 0.3 # for t \approx 0 kyr : to get eleastic response! 23 md.timestepping.final_time = 2500000 # 2,500 kyr 24 24 md.geometry.thickness = np.array([np.append(md.geometry.thickness * 0.0, 0.0), 25 np.append(md.geometry.thickness, 1000), 26 np.append(md.geometry.thickness, 2000000), 27 np.append(md.geometry.thickness * 0.0, 2000100), 28 np.append(md.geometry.thickness * 0.0, md.timestepping.start_time)]).T 25 np.append(md.geometry.thickness / 2.0, 0.1), 26 np.append(md.geometry.thickness, 0.2), 27 np.append(md.geometry.thickness, md.timestepping.start_time)]).T 29 28 30 29 #solve for GIA deflection … … 34 33 35 34 #Fields and tolerances to track changes 36 field_names = ['U gia', 'UGiaRate']35 field_names = ['UGia', 'UGiaRate'] 37 36 field_tolerances = [1e-13, 1e-13] 38 37 field_values = [md.results.GiaSolution.UGia, md.results.GiaSolution.UGiaRate] -
issm/trunk-jpl/test/NightlyRun/test2053.m
r25131 r25290 1 %Test Name: GiaIvinsBenchmarksAB2d A31 %Test Name: GiaIvinsBenchmarksAB2dD 2 2 % Benchmark experiments (Figure A2a Ivins and James, 1999, Geophys. J. Int.) 3 3 md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.); 4 4 md=setmask(md,'',''); 5 md=parameterize(md,'../Par/GiaIvinsBenchmarks AB.par');5 md=parameterize(md,'../Par/GiaIvinsBenchmarksCD.par'); 6 6 7 7 %% indicate what you want to compute 8 md.gia.cross_section_shape= 1; % for square-edged x-section8 md.gia.cross_section_shape=2; % for elliptical edge 9 9 10 % % define loading history11 md.timestepping.start_time= 2010100; % after 10 kyr of deglaciation10 % evaluation time (termed start_time) 11 md.timestepping.start_time=0.3; % for t \approx 0 kyr : to get eleastic response! 12 12 md.timestepping.final_time=2500000; % 2,500 kyr 13 14 %% define loading history {{{ 13 15 md.geometry.thickness=[... 14 16 [md.geometry.thickness*0.0; 0.0],... 15 [md.geometry.thickness; 1000],... 16 [md.geometry.thickness; 2000000],... 17 [md.geometry.thickness*0.0; 2000100],... 18 [md.geometry.thickness*0.0; md.timestepping.start_time],... 17 [md.geometry.thickness/2.0; 0.1],... 18 [md.geometry.thickness; 0.2],... 19 [md.geometry.thickness; md.timestepping.start_time],... 19 20 ]; 21 % }}} 22 23 % find out elements that have zero loads throughout the loading history. 24 pos = find(sum(abs(md.geometry.thickness(1:end-1,:)),2)==0); 25 md.mask.ice_levelset(pos)=1; % no-ice. 26 27 md.cluster=generic('name',oshostname(),'np',3); 28 md.verbose=verbose('1111111'); 20 29 21 30 %% solve for GIA deflection 22 md.cluster=generic('name',oshostname(),'np',3);23 md.verbose=verbose('1111111');24 31 md=solve(md,'Gia'); 25 32 33 %Test Name: GiaIvinsBenchmarksAB2dD1 34 U_AB2dD1 = md.results.GiaSolution.UGia; 35 URate_AB2dD1 = md.results.GiaSolution.UGiaRate; 36 37 %Test Name: GiaIvinsBenchmarksAB2dD2 38 %% different evaluation time. {{{ 39 md.timestepping.start_time=1000.3; % for t \approx 1 kyr 40 md.geometry.thickness(end,end) = md.timestepping.start_time; 41 42 md=solve(md,'Gia'); 43 44 U_AB2dD2 = md.results.GiaSolution.UGia; 45 URate_AB2dD2 = md.results.GiaSolution.UGiaRate; 46 % }}} 47 48 %Test Name: GiaIvinsBenchmarksAB2dD3 49 %% different evaluation time. {{{ 50 md.timestepping.start_time=2400000; % for t \approx \infty 51 md.geometry.thickness(end,end) = md.timestepping.start_time; 52 53 md=solve(md,'Gia'); 54 55 U_AB2dD3 = md.results.GiaSolution.UGia; 56 URate_AB2dD3 = md.results.GiaSolution.UGiaRate; 57 % }}} 58 26 59 %Fields and tolerances to track changes 27 field_names ={'UGia','UGiaRate'}; 28 field_tolerances={1e-13,1e-13}; 29 field_values={... 30 (md.results.GiaSolution.UGia),... 31 (md.results.GiaSolution.UGiaRate),... 32 }; 60 field_names ={'U_AB2dD1','URate_AB2dD1','U_AB2dD2','URate_AB2dD2','U_AB2dD3','URate_AB2dD3'}; 61 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 62 field_values={U_AB2dD1,URate_AB2dD1,U_AB2dD2,URate_AB2dD2,U_AB2dD3,URate_AB2dD3}; 63 -
issm/trunk-jpl/test/NightlyRun/test2053.py
r25131 r25290 1 #Test Name: GiaIvinsBenchmarksAB2d A31 #Test Name: GiaIvinsBenchmarksAB2dD1 2 2 from socket import gethostname 3 3 … … 11 11 12 12 13 #Benchmark experiments (Figure A2 aIvins and James, 1999, Geophys. J. Int.)13 #Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.) 14 14 md = triangle(model(), '../Exp/RoundFrontEISMINT.exp', 200000) 15 15 md = setmask(md, '', '') 16 md = parameterize(md, '../Par/GiaIvinsBenchmarks AB.py')16 md = parameterize(md, '../Par/GiaIvinsBenchmarksCD.py') 17 17 18 18 #indicate what you want to compute 19 md.gia.cross_section_shape = 1# for square-edged x - section19 md.gia.cross_section_shape = 2 # for square-edged x - section 20 20 21 21 #define loading history 22 md.timestepping.start_time = 2010100 # after 10 kyr of deglaciation22 md.timestepping.start_time = 0.3 # for t \approx 0 kyr : to get eleastic response! 23 23 md.timestepping.final_time = 2500000 # 2, 500 kyr 24 24 md.geometry.thickness = np.array([np.append(md.geometry.thickness * 0.0, 0.0), 25 np.append(md.geometry.thickness, 1000), 26 np.append(md.geometry.thickness, 2000000), 27 np.append(md.geometry.thickness * 0.0, 2000100), 28 np.append(md.geometry.thickness * 0.0, md.timestepping.start_time)]).T 25 np.append(md.geometry.thickness / 2.0, 0.1), 26 np.append(md.geometry.thickness, 0.2), 27 np.append(md.geometry.thickness, md.timestepping.start_time)]).T 29 28 30 29 #solve for GIA deflection … … 34 33 35 34 #Fields and tolerances to track changes 36 field_names = ['U gia', 'UgiaRate']35 field_names = ['UGia', 'UGiaRate'] 37 36 field_tolerances = [1e-13, 1e-13] 38 37 field_values = [md.results.GiaSolution.UGia, md.results.GiaSolution.UGiaRate] -
issm/trunk-jpl/test/Par/GiaIvinsBenchmarksAB.par
r25075 r25290 7 7 md.geometry.thickness(i)=2000.0; 8 8 else 9 md.geometry.thickness(i)= 1.0; % non-zero thickness9 md.geometry.thickness(i)=0; 10 10 end 11 11 end -
issm/trunk-jpl/test/Par/GiaIvinsBenchmarksCD.par
r25070 r25290 7 7 md.geometry.thickness(i)=3000.0; 8 8 else 9 md.geometry.thickness(i)= 1.0; % non-zero thickness9 md.geometry.thickness(i)=0; 10 10 end 11 11 end
Note:
See TracChangeset
for help on using the changeset viewer.