Changeset 25290


Ignore:
Timestamp:
07/18/20 23:19:58 (5 years ago)
Author:
adhikari
Message:

CHG: GIA ivins nightly merged, 2051,52,53 => 2051; 2071,72,73=> 2052; 2081,2082,2083 => 2053

Location:
issm/trunk-jpl/test
Files:
10 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/test/NightlyRun/test2051.m

    r25134 r25290  
    1 %Test Name: GiaIvinsBenchmarksAB2dA1
     1%Test Name: GiaIvinsBenchmarksAB2dA
    22% Benchmark experiments (Figure A2a Ivins and James, 1999, Geophys. J. Int.)
    33md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
     
    88md.gia.cross_section_shape=1;    % for square-edged x-section
    99
    10 %% define loading history
     10% evaluation time (termed start_time)
    1111md.timestepping.start_time=2002100; % after 2 kyr of deglaciation
    1212md.timestepping.final_time=2500000; % 2,500 kyr
     13
     14%% define loading history {{{
    1315md.geometry.thickness=[...
    1416        [md.geometry.thickness*0.0; 0.0],...
     
    1820        [md.geometry.thickness*0.0; md.timestepping.start_time],...
    1921        ];
     22% }}}
     23
     24% find out elements that have zero loads throughout the loading history.
     25pos = find(sum(abs(md.geometry.thickness(1:end-1,:)),2)==0);
     26md.mask.ice_levelset(pos)=1; % no-ice.
     27
     28md.cluster=generic('name',oshostname(),'np',3);
     29md.verbose=verbose('1111111');
    2030
    2131%% solve for GIA deflection
    22 md.cluster=generic('name',oshostname(),'np',3);
    23 md.verbose=verbose('1111111');
    2432md=solve(md,'Gia');
    2533
     34%Test Name: GiaIvinsBenchmarksAB2dA1
     35U_AB2dA1 = md.results.GiaSolution.UGia;
     36URate_AB2dA1 = md.results.GiaSolution.UGiaRate;
     37
     38%Test Name: GiaIvinsBenchmarksAB2dA2
     39%% different evaluation time. {{{
     40md.timestepping.start_time=2005100; % after 5 kyr of deglaciation
     41md.geometry.thickness(end,end) = md.timestepping.start_time;
     42
     43md=solve(md,'Gia');
     44
     45U_AB2dA2 = md.results.GiaSolution.UGia;
     46URate_AB2dA2 = md.results.GiaSolution.UGiaRate;
     47% }}}
     48
     49%Test Name: GiaIvinsBenchmarksAB2dA3
     50%% different evaluation time. {{{
     51md.timestepping.start_time=2010100; % after 10 kyr of deglaciation
     52md.geometry.thickness(end,end) = md.timestepping.start_time;
     53
     54md=solve(md,'Gia');
     55
     56U_AB2dA3 = md.results.GiaSolution.UGia;
     57URate_AB2dA3 = md.results.GiaSolution.UGiaRate;
     58% }}}
     59
    2660%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         };
     61field_names     ={'U_AB2dA1','URate_AB2dA1','U_AB2dA2','URate_AB2dA2','U_AB2dA3','URate_AB2dA3'};
     62field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     63field_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: GiaIvinsBenchmarksAB2dA2
     1%Test Name: GiaIvinsBenchmarksAB2dC
    22% Benchmark experiments (Figure A2a Ivins and James, 1999, Geophys. J. Int.)
    33md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
    44md=setmask(md,'','');
    5 md=parameterize(md,'../Par/GiaIvinsBenchmarksAB.par');
     5md=parameterize(md,'../Par/GiaIvinsBenchmarksCD.par');
    66
    77%% indicate what you want to compute
    88md.gia.cross_section_shape=1;    % for square-edged x-section
    99
    10 %% define loading history
    11 md.timestepping.start_time=2005100; % after 5 kyr of deglaciation
     10% evaluation time (termed start_time)
     11md.timestepping.start_time=0.3;     % for t \approx 0 kyr : to get eleastic response! 
    1212md.timestepping.final_time=2500000; % 2,500 kyr
     13
     14%% define loading history {{{
    1315md.geometry.thickness=[...
    1416        [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],...
    1920        ];
     21% }}}
     22
     23% find out elements that have zero loads throughout the loading history.
     24pos = find(sum(abs(md.geometry.thickness(1:end-1,:)),2)==0);
     25md.mask.ice_levelset(pos)=1; % no-ice.
     26
     27md.cluster=generic('name',oshostname(),'np',3);
     28md.verbose=verbose('1111111');
    2029
    2130%% solve for GIA deflection
    22 md.cluster=generic('name',oshostname(),'np',3);
    23 md.verbose=verbose('1111111');
    2431md=solve(md,'Gia');
    2532
     33%Test Name: GiaIvinsBenchmarksAB2dC1
     34U_AB2dC1 = md.results.GiaSolution.UGia;
     35URate_AB2dC1 = md.results.GiaSolution.UGiaRate;
     36
     37%Test Name: GiaIvinsBenchmarksAB2dC2
     38%% different evaluation time. {{{
     39md.timestepping.start_time=1000.3;  % for t \approx 1 kyr
     40md.geometry.thickness(end,end) = md.timestepping.start_time;
     41
     42md=solve(md,'Gia');
     43
     44U_AB2dC2 = md.results.GiaSolution.UGia;
     45URate_AB2dC2 = md.results.GiaSolution.UGiaRate;
     46% }}}
     47
     48%Test Name: GiaIvinsBenchmarksAB2dC3
     49%% different evaluation time. {{{
     50md.timestepping.start_time=2400000; % for t \approx \infty
     51md.geometry.thickness(end,end) = md.timestepping.start_time;
     52
     53md=solve(md,'Gia');
     54
     55U_AB2dC3 = md.results.GiaSolution.UGia;
     56URate_AB2dC3 = md.results.GiaSolution.UGiaRate;
     57% }}}
     58
    2659%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         };
     60field_names     ={'U_AB2dC1','URate_AB2dC1','U_AB2dC2','URate_AB2dC2','U_AB2dC3','URate_AB2dC3'};
     61field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     62field_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: GiaIvinsBenchmarksAB2dA2
     1#Test Name: GiaIvinsBenchmarksAB2dC1
    22from socket import gethostname
    33
     
    1111
    1212
    13 #Benchmark experiments (Figure A2a Ivins and James, 1999, Geophys. J. Int.)
     13#Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.)
    1414md = triangle(model(), '../Exp/RoundFrontEISMINT.exp', 200000)
    1515md = setmask(md, '', '')
    16 md = parameterize(md, '../Par/GiaIvinsBenchmarksAB.py')
     16md = parameterize(md, '../Par/GiaIvinsBenchmarksCD.py')
    1717
    1818#indicate what you want to compute
    19 md.gia.cross_section_shape = 1  # for square-edged x - section
     19md.gia.cross_section_shape = 1  # for square-edged x-section
    2020
    2121#define loading history
    22 md.timestepping.start_time = 2005100  # after 5 kyr of deglaciation
    23 md.timestepping.final_time = 2500000  # 2, 500 kyr
     22md.timestepping.start_time = 0.3  # for t \approx 0 kyr : to get eleastic response!
     23md.timestepping.final_time = 2500000  # 2,500 kyr
    2424md.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
    2928
    3029#solve for GIA deflection
     
    3433
    3534#Fields and tolerances to track changes
    36 field_names = ['Ugia', 'UGiaRate']
     35field_names = ['UGia', 'UGiaRate']
    3736field_tolerances = [1e-13, 1e-13]
    3837field_values = [md.results.GiaSolution.UGia, md.results.GiaSolution.UGiaRate]
  • issm/trunk-jpl/test/NightlyRun/test2053.m

    r25131 r25290  
    1 %Test Name: GiaIvinsBenchmarksAB2dA3
     1%Test Name: GiaIvinsBenchmarksAB2dD
    22% Benchmark experiments (Figure A2a Ivins and James, 1999, Geophys. J. Int.)
    33md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
    44md=setmask(md,'','');
    5 md=parameterize(md,'../Par/GiaIvinsBenchmarksAB.par');
     5md=parameterize(md,'../Par/GiaIvinsBenchmarksCD.par');
    66
    77%% indicate what you want to compute
    8 md.gia.cross_section_shape=1;    % for square-edged x-section
     8md.gia.cross_section_shape=2;    % for elliptical edge
    99
    10 %% define loading history
    11 md.timestepping.start_time=2010100; % after 10 kyr of deglaciation
     10% evaluation time (termed start_time)
     11md.timestepping.start_time=0.3;     % for t \approx 0 kyr : to get eleastic response! 
    1212md.timestepping.final_time=2500000; % 2,500 kyr
     13
     14%% define loading history {{{
    1315md.geometry.thickness=[...
    1416        [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],...
    1920        ];
     21% }}}
     22
     23% find out elements that have zero loads throughout the loading history.
     24pos = find(sum(abs(md.geometry.thickness(1:end-1,:)),2)==0);
     25md.mask.ice_levelset(pos)=1; % no-ice.
     26
     27md.cluster=generic('name',oshostname(),'np',3);
     28md.verbose=verbose('1111111');
    2029
    2130%% solve for GIA deflection
    22 md.cluster=generic('name',oshostname(),'np',3);
    23 md.verbose=verbose('1111111');
    2431md=solve(md,'Gia');
    2532
     33%Test Name: GiaIvinsBenchmarksAB2dD1
     34U_AB2dD1 = md.results.GiaSolution.UGia;
     35URate_AB2dD1 = md.results.GiaSolution.UGiaRate;
     36
     37%Test Name: GiaIvinsBenchmarksAB2dD2
     38%% different evaluation time. {{{
     39md.timestepping.start_time=1000.3;  % for t \approx 1 kyr
     40md.geometry.thickness(end,end) = md.timestepping.start_time;
     41
     42md=solve(md,'Gia');
     43
     44U_AB2dD2 = md.results.GiaSolution.UGia;
     45URate_AB2dD2 = md.results.GiaSolution.UGiaRate;
     46% }}}
     47
     48%Test Name: GiaIvinsBenchmarksAB2dD3
     49%% different evaluation time. {{{
     50md.timestepping.start_time=2400000; % for t \approx \infty
     51md.geometry.thickness(end,end) = md.timestepping.start_time;
     52
     53md=solve(md,'Gia');
     54
     55U_AB2dD3 = md.results.GiaSolution.UGia;
     56URate_AB2dD3 = md.results.GiaSolution.UGiaRate;
     57% }}}
     58
    2659%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         };
     60field_names     ={'U_AB2dD1','URate_AB2dD1','U_AB2dD2','URate_AB2dD2','U_AB2dD3','URate_AB2dD3'};
     61field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     62field_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: GiaIvinsBenchmarksAB2dA3
     1#Test Name: GiaIvinsBenchmarksAB2dD1
    22from socket import gethostname
    33
     
    1111
    1212
    13 #Benchmark experiments (Figure A2a Ivins and James, 1999, Geophys. J. Int.)
     13#Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.)
    1414md = triangle(model(), '../Exp/RoundFrontEISMINT.exp', 200000)
    1515md = setmask(md, '', '')
    16 md = parameterize(md, '../Par/GiaIvinsBenchmarksAB.py')
     16md = parameterize(md, '../Par/GiaIvinsBenchmarksCD.py')
    1717
    1818#indicate what you want to compute
    19 md.gia.cross_section_shape = 1  # for square-edged x - section
     19md.gia.cross_section_shape = 2  # for square-edged x - section
    2020
    2121#define loading history
    22 md.timestepping.start_time = 2010100  # after 10 kyr of deglaciation
     22md.timestepping.start_time = 0.3  # for t \approx 0 kyr : to get eleastic response!
    2323md.timestepping.final_time = 2500000  # 2, 500 kyr
    2424md.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
    2928
    3029#solve for GIA deflection
     
    3433
    3534#Fields and tolerances to track changes
    36 field_names = ['Ugia', 'UgiaRate']
     35field_names = ['UGia', 'UGiaRate']
    3736field_tolerances = [1e-13, 1e-13]
    3837field_values = [md.results.GiaSolution.UGia, md.results.GiaSolution.UGiaRate]
  • issm/trunk-jpl/test/Par/GiaIvinsBenchmarksAB.par

    r25075 r25290  
    77                md.geometry.thickness(i)=2000.0;
    88        else
    9                 md.geometry.thickness(i)=1.0; % non-zero thickness
     9                md.geometry.thickness(i)=0;
    1010        end
    1111end
  • issm/trunk-jpl/test/Par/GiaIvinsBenchmarksCD.par

    r25070 r25290  
    77                md.geometry.thickness(i)=3000.0;
    88        else
    9                 md.geometry.thickness(i)=1.0; % non-zero thickness
     9                md.geometry.thickness(i)=0;
    1010        end
    1111end
Note: See TracChangeset for help on using the changeset viewer.