Ignore:
Timestamp:
10/26/22 06:26:59 (2 years ago)
Author:
vverjans
Message:

NEW: added seasonal capabilities in frontalforcingsrignotarma

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r27318 r27334  
    55md=setflowequation(md,'SSA','all');
    66md.timestepping.start_time = 0;
    7 md.timestepping.time_step  = 1;
    8 md.timestepping.final_time = 10;
     7md.timestepping.time_step  = 0.05;
     8md.timestepping.final_time = 2;
    99
    1010%Basin separation TF
    … …  
    5252trendlin         = [-0.5,-0.2,0.1;0,0,0];
    5353trendquad        = [0,0.0,0;0.1,0.1,0.1];
    54 datebreaks       = [4,7;4,7];
     54datebreaks       = [0.5,1.5;0.5,1.5];
    5555polynomialparams = cat(numparams,intercept,trendlin,trendquad);
     56% Monthly effects params %
     57numbreaksM       = 1;
     58intcpsMp0        = [-0.5,-0.5,0,0,0,0,0.5,0.5,0,0,0,0;
     59                    -1.0,-1.0,0,0,0,0,1.0,1.0,0,0,0,0];
     60intcpsMp1        = [-0.25,-0.25,0,0,0,0,0.,0.,0,0,0,0;
     61                    -0.1,-0.1,0,0,0,0,0.1,0.1,0,0,0,0];
     62intcpsM          = cat(3,intcpsMp0,intcpsMp1);
     63trendsMp0        = [0,0,0,0,0,0,0.,0.0,0,0,0,0;
     64                    0.0,0.0,0,-0.0,0,0,0.0,0.0,0,0,0,0];
     65trendsMp1        = [0,-0.12,0,0,0,0,0.,0.0,0,0.0,0,0;
     66                    0.0,-0.1,0,-0.0,0,0,0.0,0.0,0,0,0,0];
     67trendsM          = cat(3,trendsMp0,trendsMp1);
     68datebreaksM      = [1;1];
    5669
    57 md.frontalforcings.num_basins           = nb_tf;
    58 md.frontalforcings.basin_id             = idb_tf;
    59 md.frontalforcings.num_params           = numparams; %number of parameters in the polynomial
    60 md.frontalforcings.num_breaks           = numbreaks; %number of breakpoints
    61 md.frontalforcings.subglacial_discharge = 0.01*ones(md.mesh.numberofvertices,1);
    62 md.frontalforcings.polynomialparams     = polynomialparams;
    63 md.frontalforcings.datebreaks           = datebreaks;
    64 md.frontalforcings.ar_order             = 4;
    65 md.frontalforcings.ma_order             = 2;
    66 md.frontalforcings.arma_timestep        = 2; %timestep of the ARMA model [yr]
    67 md.frontalforcings.arlag_coefs          = [[0.1,-0.1,0.01,-0.01];[0.2,-0.2,0.1,0.0]]; %autoregressive parameters
    68 md.frontalforcings.malag_coefs          = [[0.1,0.0];[0.0,0.1]]; %moving-average parameters
    69 
     70md.frontalforcings.num_basins              = nb_tf;
     71md.frontalforcings.basin_id                = idb_tf;
     72md.frontalforcings.num_params              = numparams; %number of parameters in the polynomial
     73md.frontalforcings.num_breaks              = numbreaks; %number of breakpoints
     74md.frontalforcings.subglacial_discharge    = 0.01*ones(md.mesh.numberofvertices,1);
     75md.frontalforcings.polynomialparams        = polynomialparams;
     76md.frontalforcings.datebreaks              = datebreaks;
     77md.frontalforcings.ar_order                = 4;
     78md.frontalforcings.ma_order                = 2;
     79md.frontalforcings.arma_timestep           = 2; %timestep of the ARMA model [yr]
     80md.frontalforcings.arlag_coefs             = [[0.1,-0.1,0.01,-0.01];[0.2,-0.2,0.1,0.0]]; %autoregressive parameters
     81md.frontalforcings.malag_coefs             = [[0.1,0.0];[0.0,0.1]]; %moving-average parameters
     82md.frontalforcings.monthlyvals_numbreaks   = numbreaksM;
     83md.frontalforcings.monthlyvals_intercepts  = intcpsM;
     84md.frontalforcings.monthlyvals_trends      = trendsM;
     85md.frontalforcings.monthlyvals_datebreaks  = datebreaksM;
    7086% Floating Ice Melt parameters
    7187md.basalforcings.floatingice_melting_rate = 0.1*ones(md.mesh.numberofvertices,1);
    … …  
    102118   1e-11,2e-11,2e-11,1e-11,1e-9,1e-10,1e-10,1e-10,...
    103119   2e-11,1e-11,1e-11,9e-11,2e-9,1e-10,1e-10,1e-10,...
    104    2e-10,1e-10,1e-10,1e-10,5e-9,1e-10,1e-10,1e-10,...
     120   2e-6,1e-6,1e-6,1e-6,5e-6,1e-6,1e-6,1e-6,...
    105121   };
    106122field_values={...
    … …  
    113129   (md.results.TransientSolution(1).CalvingMeltingrate),...
    114130   (md.results.TransientSolution(1).BasalforcingsFloatingiceMeltingRate),...
    115    (md.results.TransientSolution(5).Vx),...
    116    (md.results.TransientSolution(5).Vy),...
    117    (md.results.TransientSolution(5).Vel),...
    118    (md.results.TransientSolution(5).Thickness),...
    119    (md.results.TransientSolution(5).MaskIceLevelset),...
    120    (md.results.TransientSolution(5).CalvingCalvingrate),...
    121    (md.results.TransientSolution(5).CalvingMeltingrate),...
    122    (md.results.TransientSolution(5).BasalforcingsFloatingiceMeltingRate),...
    123         (md.results.TransientSolution(10).Vx),...
    124         (md.results.TransientSolution(10).Vy),...
    125         (md.results.TransientSolution(10).Vel),...
    126         (md.results.TransientSolution(10).Thickness),...
    127         (md.results.TransientSolution(10).MaskIceLevelset),...
    128         (md.results.TransientSolution(10).CalvingCalvingrate),...
    129         (md.results.TransientSolution(10).CalvingMeltingrate),...
    130         (md.results.TransientSolution(10).BasalforcingsFloatingiceMeltingRate),...
     131   (md.results.TransientSolution(20).Vx),...
     132   (md.results.TransientSolution(20).Vy),...
     133   (md.results.TransientSolution(20).Vel),...
     134   (md.results.TransientSolution(20).Thickness),...
     135   (md.results.TransientSolution(20).MaskIceLevelset),...
     136   (md.results.TransientSolution(20).CalvingCalvingrate),...
     137   (md.results.TransientSolution(20).CalvingMeltingrate),...
     138   (md.results.TransientSolution(20).BasalforcingsFloatingiceMeltingRate),...
     139        (md.results.TransientSolution(40).Vx),...
     140        (md.results.TransientSolution(40).Vy),...
     141        (md.results.TransientSolution(40).Vel),...
     142        (md.results.TransientSolution(40).Thickness),...
     143        (md.results.TransientSolution(40).MaskIceLevelset),...
     144        (md.results.TransientSolution(40).CalvingCalvingrate),...
     145        (md.results.TransientSolution(40).CalvingMeltingrate),...
     146        (md.results.TransientSolution(40).BasalforcingsFloatingiceMeltingRate),...
    131147        };
Note: See TracChangeset for help on using the changeset viewer.