Changeset 3796


Ignore:
Timestamp:
05/18/10 10:05:26 (15 years ago)
Author:
seroussi
Message:

parameter fields names

Location:
issm/trunk/src/m/solutions/jpl
Files:
24 edited

Legend:

Unmodified
Added
Removed

  • TabularUnified issm/trunk/src/m/solutions/jpl/SpawnCore.m

    r2330 r3796  
    77
    88%recover parameters
    9 verbose=models.dh.parameters.verbose;
     9verbose=models.dh.parameters.Verbose;
    1010parameters=models.dh.parameters;
    1111responsedescriptors=models.dh.parameters.responsedescriptors;
    1212npart=models.dh.parameters.qmu_npart;
    1313part=models.dh.parameters.qmu_part+1; %C indexing
    14 numberofnodes=models.dh.parameters.numberofnodes;
     14numberofnodes=models.dh.parameters.NumberOfNodes;
    1515
    1616disp(['   qmu iteration: ',num2str(counter)]);

  • TabularUnified issm/trunk/src/m/solutions/jpl/balancedthickness.m

    r3589 r3796  
    2020        displaystring(md.verbose,'\n%s',['setup inputs...']);
    2121        inputs=inputlist;
    22         inputs=add(inputs,'velocity',models.bt.parameters.u_g,'doublevec',3,models.bt.parameters.numberofnodes);
    23         inputs=add(inputs,'melting',models.bt.parameters.m_g,'doublevec',1,models.bt.parameters.numberofnodes);
    24         inputs=add(inputs,'accumulation',models.bt.parameters.a_g,'doublevec',1,models.bt.parameters.numberofnodes);
     22        inputs=add(inputs,'velocity',models.bt.parameters.u_g,'doublevec',3,models.bt.parameters.NumberOfNodes);
     23        inputs=add(inputs,'melting',models.bt.parameters.m_g,'doublevec',1,models.bt.parameters.NumberOfNodes);
     24        inputs=add(inputs,'accumulation',models.bt.parameters.a_g,'doublevec',1,models.bt.parameters.NumberOfNodes);
    2525
    2626        displaystring(md.verbose,'\n%s',['call computational core:']);

  • TabularUnified issm/trunk/src/m/solutions/jpl/balancedthickness2.m

    r3589 r3796  
    2020        displaystring(md.verbose,'\n%s',['setup inputs...']);
    2121        inputs=inputlist;
    22         inputs=add(inputs,'vx',models.p.parameters.vx_g,'doublevec',1,models.p.parameters.numberofvertices);
    23         inputs=add(inputs,'vy',models.p.parameters.vy_g,'doublevec',1,models.p.parameters.numberofvertices);
    24         inputs=add(inputs,'thickness',models.p.parameters.h_g,'doublevec',1,models.p.parameters.numberofvertices);
    25         inputs=add(inputs,'dhdt',models.p.parameters.dhdt_g,'doublevec',1,models.p.parameters.numberofvertices);
    26         inputs=add(inputs,'melting',models.p.parameters.m_g,'doublevec',1,models.p.parameters.numberofvertices);
    27         inputs=add(inputs,'accumulation',models.p.parameters.a_g,'doublevec',1,models.p.parameters.numberofvertices);
     22        inputs=add(inputs,'vx',models.p.parameters.vx_g,'doublevec',1,models.p.parameters.NumberOfVertices);
     23        inputs=add(inputs,'vy',models.p.parameters.vy_g,'doublevec',1,models.p.parameters.NumberOfVertices);
     24        inputs=add(inputs,'thickness',models.p.parameters.h_g,'doublevec',1,models.p.parameters.NumberOfVertices);
     25        inputs=add(inputs,'dhdt',models.p.parameters.dhdt_g,'doublevec',1,models.p.parameters.NumberOfVertices);
     26        inputs=add(inputs,'melting',models.p.parameters.m_g,'doublevec',1,models.p.parameters.NumberOfVertices);
     27        inputs=add(inputs,'accumulation',models.p.parameters.a_g,'doublevec',1,models.p.parameters.NumberOfVertices);
    2828
    2929        displaystring(md.verbose,'\n%s',['call computational core:']);

  • TabularUnified issm/trunk/src/m/solutions/jpl/balancedthickness2_core.m

    r3589 r3796  
    1010        results.step=1;
    1111
    12         displaystring(m.parameters.verbose,'\n%s',['depth averaging velocity...']);
     12        displaystring(m.parameters.Verbose,'\n%s',['depth averaging velocity...']);
    1313        %Take only the first two dofs of m.parameters.u_g
    1414        vx_g=get(inputs,'vx',1);
    … …  
    1919        %vy_g=FieldDepthAverage(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,vy_g,'vy');
    2020
    21         inputs=add(inputs,'vx_average',vx_g,'doublevec',1,m.parameters.numberofvertices);
    22         inputs=add(inputs,'vy_average',vy_g,'doublevec',1,m.parameters.numberofvertices);
     21        inputs=add(inputs,'vx_average',vx_g,'doublevec',1,m.parameters.NumberOfVertices);
     22        inputs=add(inputs,'vy_average',vy_g,'doublevec',1,m.parameters.NumberOfVertices);
    2323
    24         displaystring(m.parameters.verbose,'\n%s',['call computational core:']);
     24        displaystring(m.parameters.Verbose,'\n%s',['call computational core:']);
    2525        results.h_g=diagnostic_core_linear(m,inputs,analysis_type,sub_analysis_type);
    2626
    27         displaystring(m.parameters.verbose,'\n%s',['averaging over vertices']);
     27        displaystring(m.parameters.Verbose,'\n%s',['averaging over vertices']);
    2828        results.h_g=FieldAverageOntoVertices(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,results.h_g);
    2929
    30         displaystring(m.parameters.verbose,'\n%s',['extrude computed thickness on all layers:']);
     30        displaystring(m.parameters.Verbose,'\n%s',['extrude computed thickness on all layers:']);
    3131        %results.h_g=FieldAverageOntoVertices(m.elements,m.nodes,m.loads,m.materials,m.parameters,results.h_g,'thickness');
    3232

  • TabularUnified issm/trunk/src/m/solutions/jpl/balancedthickness_core.m

    r3589 r3796  
    1010        results.step=1;
    1111
    12         displaystring(m.parameters.verbose,'\n%s',['depth averaging velocity...']);
     12        displaystring(m.parameters.Verbose,'\n%s',['depth averaging velocity...']);
    1313        %Take only the first two dofs of m.parameters.u_g
    1414        u_g=get(inputs,'velocity',[1 1 0 0]);
    1515        u_g=FieldDepthAverage(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,u_g,'velocity');
    16         inputs=add(inputs,'velocity_average',u_g,'doublevec',2,m.parameters.numberofnodes);
     16        inputs=add(inputs,'velocity_average',u_g,'doublevec',2,m.parameters.NumberOfNodes);
    1717
    18         displaystring(m.parameters.verbose,'\n%s',['call computational core:']);
     18        displaystring(m.parameters.Verbose,'\n%s',['call computational core:']);
    1919        results.h_g=diagnostic_core_linear(m,inputs,analysis_type,sub_analysis_type);
    2020
    21         displaystring(m.parameters.verbose,'\n%s',['extrude computed thickness on all layers:']);
     21        displaystring(m.parameters.Verbose,'\n%s',['extrude computed thickness on all layers:']);
    2222        results.h_g=FieldExtrude(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,results.h_g,'thickness',0);
    2323

  • TabularUnified issm/trunk/src/m/solutions/jpl/balancedvelocities.m

    r2722 r3796  
    2020        displaystring(md.verbose,'\n%s',['setup inputs...']);
    2121        inputs=inputlist;
    22         inputs=add(inputs,'velocity',models.p.parameters.u_g,'doublevec',3,models.p.parameters.numberofnodes);
    23         inputs=add(inputs,'thickness',models.p.parameters.h_g,'doublevec',1,models.p.parameters.numberofnodes);
    24         inputs=add(inputs,'melting',models.p.parameters.m_g,'doublevec',1,models.p.parameters.numberofnodes);
    25         inputs=add(inputs,'accumulation',models.p.parameters.a_g,'doublevec',1,models.p.parameters.numberofnodes);
    26         inputs=add(inputs,'dt',models.p.parameters.dt*models.p.parameters.yts,'double');
     22        inputs=add(inputs,'velocity',models.p.parameters.u_g,'doublevec',3,models.p.parameters.NumberOfNodes);
     23        inputs=add(inputs,'thickness',models.p.parameters.h_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     24        inputs=add(inputs,'melting',models.p.parameters.m_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     25        inputs=add(inputs,'accumulation',models.p.parameters.a_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     26        inputs=add(inputs,'dt',models.p.parameters.Dt*models.p.parameters.Yts,'double');
    2727
    2828        displaystring(md.verbose,'\n%s',['call computational core:']);

  • TabularUnified issm/trunk/src/m/solutions/jpl/balancedvelocities_core.m

    r2722 r3796  
    1010        results.step=1;
    1111
    12         displaystring(m.parameters.verbose,'\n%s',['depth averaging velocity...']);
     12        displaystring(m.parameters.Verbose,'\n%s',['depth averaging velocity...']);
    1313        %Take only the first two dofs of m.parameters.u_g
    1414        u_g=get(inputs,'velocity',[1 1 0 0]);
    1515        u_g=FieldDepthAverage(m.elements,m.nodes,m.loads,m.materials,m.parameters,u_g,'velocity');
    16         inputs=add(inputs,'velocity_average',u_g,'doublevec',2,m.parameters.numberofnodes);
     16        inputs=add(inputs,'velocity_average',u_g,'doublevec',2,m.parameters.NumberOfNodes);
    1717
    18         displaystring(m.parameters.verbose,'\n%s',['call computational core:']);
     18        displaystring(m.parameters.Verbose,'\n%s',['call computational core:']);
    1919        results.h_g=diagnostic_core_linear(m,inputs,analysis_type,sub_analysis_type);
    2020
    21         displaystring(m.parameters.verbose,'\n%s',['extrude computed thickness on all layers:']);
     21        displaystring(m.parameters.Verbose,'\n%s',['extrude computed thickness on all layers:']);
    2222        results.v_g=FieldExtrude(m.elements,m.nodes,m.loads,m.materials,m.parameters,results.h_g,'thickness',0);
    2323

  • TabularUnified issm/trunk/src/m/solutions/jpl/control_core.m

    r3524 r3796  
    1212model=models.active;
    1313
     14model.parameters
    1415%recover parameters common to all solutions
    15 verbose=model.parameters.verbose;
    16 dim=model.parameters.dim;
    17 isstokes=model.parameters.isstokes;
     16verbose=model.parameters.Verbose;
     17dim=model.parameters.Dim;
     18isstokes=model.parameters.IsStokes;
    1819
    1920%recover other parameters
    20 eps_cm=model.parameters.eps_cm;
    21 fit=model.parameters.fit;
     21eps_cm=model.parameters.EpsCm;
     22fit=model.parameters.Fit;
    2223
    2324%initialize control parameters
    … …  
    2728options=ControlOptions(model.parameters);
    2829
    29 for n=1:model.parameters.nsteps,
     30for n=1:model.parameters.NSteps,
    3031
    3132        %set options
    32         options=optimset(options,'MaxFunEvals',model.parameters.maxiter(n));
     33        options=optimset(options,'MaxFunEvals',model.parameters.MaxIter(n));
    3334
    34         disp(sprintf('\n%s%s%s%s\n',['   control method step ' num2str(n) '/' num2str(model.parameters.nsteps)]));
     35        disp(sprintf('\n%s%s%s%s\n',['   control method step ' num2str(n) '/' num2str(model.parameters.NSteps)]));
    3536
    3637        %In case we are running a steady state control method, compute new temperature field using new parameter distribution:
    37         if model.parameters.control_steady;
     38        if model.parameters.ControlSteady;
    3839                steadystate_results=steadystate_core(models,inputs); t_g=steadystate_results.t_g;
    39                 inputs=add(inputs,'temperature',t_g,'doublevec',1,model.parameters.numberofnodes);
     40                inputs=add(inputs,'temperature',t_g,'doublevec',1,model.parameters.NumberOfNodes);
    4041        end
    4142
    4243        %update inputs with new fit
    43         inputs=add(inputs,'fit',model.parameters.fit(n),'double');
    44         inputs=add(inputs,model.parameters.control_type,param_g,'doublevec',1,model.parameters.numberofnodes);
     44        inputs=add(inputs,'fit',model.parameters.Fit(n),'double');
     45        inputs=add(inputs,model.parameters.ControlType,param_g,'doublevec',1,model.parameters.NumberOfNodes);
    4546
    4647        %Update inputs in datasets
    … …  
    5253        if dim==3,
    5354                if isstokes,
    54                         inputs=add(inputs,'velocity',u_g,'doublevec',4,model.parameters.numberofnodes);
     55                        inputs=add(inputs,'velocity',u_g,'doublevec',4,model.parameters.NumberOfNodes);
    5556                else
    56                         if model.parameters.control_steady;
    57                                 inputs=add(inputs,'velocity',u_g,'doublevec',3,model.parameters.numberofnodes);
     57                        if model.parameters.ControlSteady;
     58                                inputs=add(inputs,'velocity',u_g,'doublevec',3,model.parameters.NumberOfNodes);
    5859                        else
    59                                 inputs=add(inputs,'velocity',u_g,'doublevec',2,model.parameters.numberofnodes);
     60                                inputs=add(inputs,'velocity',u_g,'doublevec',2,model.parameters.NumberOfNodes);
    6061                        end
    6162                end
    6263        else
    63                 inputs=add(inputs,'velocity',u_g,'doublevec',2,model.parameters.numberofnodes);
     64                inputs=add(inputs,'velocity',u_g,'doublevec',2,model.parameters.NumberOfNodes);
    6465        end
    6566
    … …  
    7273
    7374        %visualize direction.
    74         if model.parameters.plot
     75        if model.parameters.Plot
    7576                plot_direction;
    7677        end
    … …  
    8081
    8182        displaystring(verbose,'\n%s',['      updating parameter using optimized search scalar...']);
    82         param_g=param_g+search_scalar*model.parameters.optscal(n)*c(n).grad_g;
     83        param_g=param_g+search_scalar*model.parameters.Optscal(n)*c(n).grad_g;
    8384
    8485        displaystring(verbose,'\n%s',['      constraining the new distribution...']);
    … …  
    116117
    117118%compute final velocity from diagnostic_core (horiz+vertical)
    118 if model.parameters.control_steady;
    119         inputs=add(inputs,model.parameters.control_type,param_g,'doublevec',1,model.parameters.numberofnodes);
     119if model.parameters.ControlSteady;
     120        inputs=add(inputs,model.parameters.ControlType,param_g,'doublevec',1,model.parameters.NumberOfNodes);
    120121        steadystate_results=steadystate_core(models,inputs); t_g=steadystate_results.t_g;
    121122        u_g=steadystate_results.u_g;
    … …  
    123124        m_g=steadystate_results.m_g;
    124125else
    125         inputs=add(inputs,model.parameters.control_type,param_g,'doublevec',1,model.parameters.numberofnodes);
     126        inputs=add(inputs,model.parameters.ControlType,param_g,'doublevec',1,model.parameters.NumberOfNodes);
    126127        results_diag=diagnostic_core(models,inputs);
    127128        u_g=results_diag.u_g;
    … …  
    140141results.param_g=param_g;
    141142results.u_g=u_g;
    142 if model.parameters.control_steady,
     143if model.parameters.ControlSteady,
    143144        results.t_g=t_g;
    144145        results.m_g=m_g;

  • TabularUnified issm/trunk/src/m/solutions/jpl/diagnostic_core_linear.m

    r3589 r3796  
    1717        %Reduce tangent matrix from g size to f size
    1818        [K_ff, K_fs] = Reducematrixfromgtof( K_gg, m.Gmn, m.nodesets);
    19         displaystring(m.parameters.verbose>1,'%s%g','      condition number of stiffness matrix: ',condest(K_ff));
     19        displaystring(m.parameters.Verbose>1,'%s%g','      condition number of stiffness matrix: ',condest(K_ff));
    2020       
    2121        %Reduce load from g size to f size

  • TabularUnified issm/trunk/src/m/solutions/jpl/diagnostic_core_nonlinear.m

    r3612 r3796  
    1111        m.parameters.kflag=1; m.parameters.pflag=1;
    1212
     13        m.parameters
    1314%   initialize solution vector
    1415        converged=0; count=1;
    15         soln(count).u_g=get(inputs,'velocity',[1 1 (m.parameters.numberofdofspernode>=3) (m.parameters.numberofdofspernode>=3)]); %only keep vz if running with more than 3 dofs per node
     16        soln(count).u_g=get(inputs,'velocity',[1 1 (m.parameters.NumberOfDofsPerNode>=3) (m.parameters.NumberOfDofsPerNode>=3)]); %only keep vz if running with more than 3 dofs per node
    1617
    1718        soln(count).u_f= Reducevectorgtof(soln(count).u_g,m.nodesets);
    … …  
    2021        inputs=add(inputs,'converged',converged,'double');
    2122
    22         displaystring(m.parameters.verbose,'\n%s',['   starting direct shooting method']);
     23        displaystring(m.parameters.Verbose,'\n%s',['   starting direct shooting method']);
    2324        while(~converged),
    2425               
    2526                %add velocity to inputs
    26                 inputs=add(inputs,'velocity',soln(count).u_g,'doublevec',m.parameters.numberofdofspernode,m.parameters.numberofnodes);
     27                inputs=add(inputs,'velocity',soln(count).u_g,'doublevec',m.parameters.NumberOfDofsPerNode,m.parameters.NumberOfNodes);
    2728                if count>1,
    28                         inputs=add(inputs,'old_velocity',soln(count-1).u_g,'doublevec',m.parameters.numberofdofspernode,m.parameters.numberofnodes);
     29                        inputs=add(inputs,'old_velocity',soln(count-1).u_g,'doublevec',m.parameters.NumberOfDofsPerNode,m.parameters.NumberOfNodes);
    2930                else
    30                         inputs=add(inputs,'old_velocity',soln(count).u_g,'doublevec',m.parameters.numberofdofspernode,m.parameters.numberofnodes);
     31                        inputs=add(inputs,'old_velocity',soln(count).u_g,'doublevec',m.parameters.NumberOfDofsPerNode,m.parameters.NumberOfNodes);
    3132                end
    3233               
    … …  
    5960               
    6061                %Deal with penalty loads
    61                 inputs=add(inputs,'velocity',soln(count).u_g,'doublevec',m.parameters.numberofdofspernode,m.parameters.numberofnodes);
     62                inputs=add(inputs,'velocity',soln(count).u_g,'doublevec',m.parameters.NumberOfDofsPerNode,m.parameters.NumberOfNodes);
    6263               
    6364                %penalty constraints
    … …  
    6566                error;
    6667               
    67                 displaystring(m.parameters.verbose,'%s%i','      number of unstable constraints: ',num_unstable_constraints);
     68                displaystring(m.parameters.Verbose,'%s%i','      number of unstable constraints: ',num_unstable_constraints);
    6869               
    6970                %Figure out if convergence have been reached
    … …  
    8788        nout=max(nargout,1)-1;
    8889        if nout==2,
    89                 inputs=add(inputs,'velocity',u_g,'doublevec',m.parameters.numberofdofspernode,m.parameters.numberofnodes);
     90                inputs=add(inputs,'velocity',u_g,'doublevec',m.parameters.NumberOfDofsPerNode,m.parameters.NumberOfNodes);
    9091                m.parameters.kflag=1; m.parameters.pflag=0;
    9192                [K_gg, p_g]=SystemMatrices(m.elements,m.nodes,m.vertices,loads,m.materials,m.parameters,inputs,analysis_type,sub_analysis_type);

  • TabularUnified issm/trunk/src/m/solutions/jpl/gradjcompute_core.m

    r3524 r3796  
    55
    66%recover parameters
    7 verbose=m.parameters.verbose;
    8 dim=m.parameters.dim;
    9 extrude_param=m.parameters.extrude_param;
    10 ishutter=m.parameters.ishutter;
    11 ismacayealpattyn=m.parameters.ismacayealpattyn;
    12 isstokes=m.parameters.isstokes;
    13 analysis_type=m.parameters.analysis_type;
    14 sub_analysis_type=m.parameters.sub_analysis_type;
     7verbose=m.parameters.Verbose;
     8dim=m.parameters.Dim;
     9extrude_param=m.parameters.ExtrudeParam;
     10ishutter=m.parameters.IsHutter;
     11ismacayealpattyn=m.parameters.IsMacAyealPattyn;
     12isstokes=m.parameters.IsStokes;
     13analysis_type=m.parameters.AnalysisType;
     14sub_analysis_type=m.parameters.SubAnalysisType;
    1515
    1616%Recover solution for this stiffness and right hand side:
    1717displaystring(verbose,'%s','         computing velocities...');
    1818[u_g K_ff0 K_fs0 ]=diagnostic_core_nonlinear(m,inputs,analysis_type,sub_analysis_type);
    19 inputs=add(inputs,'velocity',u_g,'doublevec',m.parameters.numberofdofspernode,m.parameters.numberofnodes);
     19inputs=add(inputs,'velocity',u_g,'doublevec',m.parameters.NumberOfDofsPerNode,m.parameters.NumberOfNodes);
    2020
    2121%Buid Du, difference between observed velocity and model velocity.
    … …  
    3232%Merge back to g set
    3333lambda_g= Mergesolutionfromftog( lambda_f, m.Gmn, m.ys0, m.nodesets );
    34 inputs=add(inputs,'adjoint',lambda_g,'doublevec',m.parameters.numberofdofspernode,m.parameters.numberofnodes);
     34inputs=add(inputs,'adjoint',lambda_g,'doublevec',m.parameters.NumberOfDofsPerNode,m.parameters.NumberOfNodes);
    3535
    3636%Compute gradJ
    … …  
    4141end
    4242
    43 if m.parameters.control_steady,
     43if m.parameters.ControlSteady,
    4444        %compute initial velocity from diagnostic_core (horiz+vertical)
    4545        displaystring(verbose,'%s','          compute 3d initial velocity...');

  • TabularUnified issm/trunk/src/m/solutions/jpl/objectivefunctionC.m

    r3502 r3796  
    55
    66%recover some parameters
    7 optscal=m.parameters.optscal(n);
    8 fit=m.parameters.fit(n);
    9 control_type=m.parameters.control_type;
    10 analysis_type=m.parameters.analysis_type;
     7optscal=m.parameters.Optscal(n);
     8fit=m.parameters.Fit(n);
     9control_type=m.parameters.ControlType;
     10analysis_type=m.parameters.AnalysisType;
    1111
    1212%Update along gradient using scalar supplied by fmincon optimization routine
    … …  
    1414
    1515%Plug parameter into inputs
    16 inputs=add(inputs,m.parameters.control_type,parameter,'doublevec',1,m.parameters.numberofnodes);
     16inputs=add(inputs,m.parameters.ControlType,parameter,'doublevec',1,m.parameters.NumberOfNodes);
    1717
    1818%Run diagnostic with updated parameters.
    … …  
    2020
    2121%add velocity to inputs.
    22 inputs=add(inputs,'velocity',u_g,'doublevec',m.parameters.numberofdofspernode,m.parameters.numberofnodes);
     22inputs=add(inputs,'velocity',u_g,'doublevec',m.parameters.NumberOfDofsPerNode,m.parameters.NumberOfNodes);
    2323
    2424%Compute misfit for this velocity field.

  • TabularUnified issm/trunk/src/m/solutions/jpl/processresults.m

    r3552 r3796  
    2121
    2222        %some flags needed
    23         dim=m_dhu.parameters.dim;
    24         ishutter=m_dhu.parameters.ishutter;
    25         ismacayealpattyn=m_dh.parameters.ismacayealpattyn;
    26         isstokes=m_ds.parameters.isstokes;
     23        dim=m_dhu.parameters.Dim;
     24        ishutter=m_dhu.parameters.IsHutter;
     25        ismacayealpattyn=m_dh.parameters.IsMacAyealPattyn;
     26        isstokes=m_ds.parameters.IsStokes;
    2727end
    2828if (analysis_type==ControlAnalysisEnum()),
    … …  
    3030
    3131        %some flags needed
    32         dim=m_dh.parameters.dim;
    33         ishutter=m_dh.parameters.ishutter;
    34         ismacayealpattyn=m_dh.parameters.ismacayealpattyn;
    35         isstokes=m_dh.parameters.isstokes;
     32        dim=m_dh.parameters.Dim;
     33        ishutter=m_dh.parameters.IsHutter;
     34        ismacayealpattyn=m_dh.parameters.IsMacAyealPattyn;
     35        isstokes=m_dh.parameters.IsStokes;
    3636end
    3737if (analysis_type==ThermalAnalysisEnum() | analysis_type==TransientAnalysisEnum() | analysis_type==SteadystateAnalysisEnum()),
    … …  
    5252
    5353                        u_g=results(i).u_g;
    54                         yts=m_dh.parameters.yts;
     54                        yts=m_dh.parameters.Yts;
    5555
    5656                        if dim==2,
    … …  
    9393
    9494                        m_g=results(i).m_g;
    95                         yts=m_m.parameters.yts;
     95                        yts=m_m.parameters.Yts;
    9696                        newresults(i).step=results(i).step;
    9797                        newresults(i).time=results(i).time;

  • TabularUnified issm/trunk/src/m/solutions/jpl/prognostic.m

    r2330 r3796  
    2020        displaystring(md.verbose,'\n%s',['setup inputs...']);
    2121        inputs=inputlist;
    22         inputs=add(inputs,'velocity',models.p.parameters.u_g,'doublevec',3,models.p.parameters.numberofnodes);
    23         inputs=add(inputs,'thickness',models.p.parameters.h_g,'doublevec',1,models.p.parameters.numberofnodes);
    24         inputs=add(inputs,'melting',models.p.parameters.m_g,'doublevec',1,models.p.parameters.numberofnodes);
    25         inputs=add(inputs,'accumulation',models.p.parameters.a_g,'doublevec',1,models.p.parameters.numberofnodes);
    26         inputs=add(inputs,'dt',models.p.parameters.dt*models.p.parameters.yts,'double');
     22        inputs=add(inputs,'velocity',models.p.parameters.u_g,'doublevec',3,models.p.parameters.NumberOfNodes);
     23        inputs=add(inputs,'thickness',models.p.parameters.h_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     24        inputs=add(inputs,'melting',models.p.parameters.m_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     25        inputs=add(inputs,'accumulation',models.p.parameters.a_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     26        inputs=add(inputs,'dt',models.p.parameters.Dt*models.p.parameters.Yts,'double');
    2727
    2828        displaystring(md.verbose,'\n%s',['call computational core:']);

  • TabularUnified issm/trunk/src/m/solutions/jpl/prognostic2_core.m

    r3556 r3796  
    1010        results.step=1;
    1111
    12         displaystring(m.parameters.verbose,'\n%s',['depth averaging velocity...']);
     12        displaystring(m.parameters.Verbose,'\n%s',['depth averaging velocity...']);
    1313        %Take only the first two dofs of m.parameters.u_g
    1414        vx_g=get(inputs,'vx',1);
    … …  
    1919        %vy_g=FieldDepthAverage(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,vy_g,'vy');
    2020
    21         inputs=add(inputs,'vx_average',vx_g,'doublevec',1,m.parameters.numberofvertices);
    22         inputs=add(inputs,'vy_average',vy_g,'doublevec',1,m.parameters.numberofvertices);
     21        inputs=add(inputs,'vx_average',vx_g,'doublevec',1,m.parameters.NumberOfVertices);
     22        inputs=add(inputs,'vy_average',vy_g,'doublevec',1,m.parameters.NumberOfVertices);
    2323
    24         displaystring(m.parameters.verbose,'\n%s',['call computational core:']);
     24        displaystring(m.parameters.Verbose,'\n%s',['call computational core:']);
    2525        results.h_g=diagnostic_core_linear(m,inputs,analysis_type,sub_analysis_type);
    2626
    27         displaystring(m.parameters.verbose,'\n%s',['averaging over vertices']);
     27        displaystring(m.parameters.Verbose,'\n%s',['averaging over vertices']);
    2828        results.h_g=FieldAverageOntoVertices(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,results.h_g);
    2929
    30         displaystring(m.parameters.verbose,'\n%s',['extrude computed thickness on all layers:']);
     30        displaystring(m.parameters.Verbose,'\n%s',['extrude computed thickness on all layers:']);
    3131        %results.h_g=FieldAverageOntoVertices(m.elements,m.nodes,m.loads,m.materials,m.parameters,results.h_g,'thickness');
    3232

  • TabularUnified issm/trunk/src/m/solutions/jpl/prognostic_core.m

    r3485 r3796  
    1010        results.step=1;
    1111
    12         displaystring(m.parameters.verbose,'\n%s',['depth averaging velocity...']);
     12        displaystring(m.parameters.Verbose,'\n%s',['depth averaging velocity...']);
    1313        %Take only the first two dofs of m.parameters.u_g
    1414        u_g=get(inputs,'velocity',[1 1 0 0]);
    1515        u_g=FieldDepthAverage(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,u_g,'velocity');
    16         inputs=add(inputs,'velocity_average',u_g,'doublevec',2,m.parameters.numberofnodes);
     16        inputs=add(inputs,'velocity_average',u_g,'doublevec',2,m.parameters.NumberOfNodes);
    1717
    18         displaystring(m.parameters.verbose,'\n%s',['call computational core:']);
     18        displaystring(m.parameters.Verbose,'\n%s',['call computational core:']);
    1919        results.h_g=diagnostic_core_linear(m,inputs,analysis_type,sub_analysis_type);
    2020
    21         displaystring(m.parameters.verbose,'\n%s',['extrude computed thickness on all layers:']);
     21        displaystring(m.parameters.Verbose,'\n%s',['extrude computed thickness on all layers:']);
    2222        results.h_g=FieldExtrude(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,results.h_g,'thickness',0);
    2323

  • TabularUnified issm/trunk/src/m/solutions/jpl/slopecompute_core.m

    r3485 r3796  
    1010        results.step=1;
    1111
    12         displaystring(m.parameters.verbose,'\n%s',['call computational core:']);
     12        displaystring(m.parameters.Verbose,'\n%s',['call computational core:']);
    1313        results.sx_g=diagnostic_core_linear(m,inputs,analysis_type,SurfaceXAnalysisEnum());
    1414        results.sy_g=diagnostic_core_linear(m,inputs,analysis_type,SurfaceYAnalysisEnum());
    1515
    16         displaystring(m.parameters.verbose,'\n%s',['extrude computed slope on all layers:']);
     16        displaystring(m.parameters.Verbose,'\n%s',['extrude computed slope on all layers:']);
    1717        results.sx_g=FieldExtrude(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,results.sx_g,'slopex',0);
    1818        results.sy_g=FieldExtrude(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,results.sy_g,'slopey',0);

  • TabularUnified issm/trunk/src/m/solutions/jpl/steadystate.m

    r3524 r3796  
    4040        displaystring(md.verbose,'\n%s',['setup inputs...']);
    4141        inputs=inputlist;
    42         inputs=add(inputs,'velocity',models.dh.parameters.u_g,'doublevec',3,models.dh.parameters.numberofnodes);
    43         inputs=add(inputs,'pressure',models.t.parameters.p_g,'doublevec',1,models.t.parameters.numberofnodes);
    44         inputs=add(inputs,'dt',models.t.parameters.dt*models.t.parameters.yts,'double');
     42        inputs=add(inputs,'velocity',models.dh.parameters.u_g,'doublevec',3,models.dh.parameters.NumberOfNodes);
     43        inputs=add(inputs,'pressure',models.t.parameters.p_g,'doublevec',1,models.t.parameters.NumberOfNodes);
     44        inputs=add(inputs,'dt',models.t.parameters.Dt*models.t.parameters.Yts,'double');
    4545        if md.control_analysis,
    46                 inputs=add(inputs,'velocity_obs',models.dh.parameters.u_g_obs,'doublevec',2,models.dh.parameters.numberofnodes);
    47                 inputs=add(inputs,'weights',models.dh.parameters.weights,'doublevec',1,models.dh.parameters.numberofnodes);
     46                inputs=add(inputs,'velocity_obs',models.dh.parameters.u_g_obs,'doublevec',2,models.dh.parameters.NumberOfNodes);
     47                inputs=add(inputs,'weights',models.dh.parameters.Weights,'doublevec',1,models.dh.parameters.NumberOfNodes);
    4848        end
    4949       
    5050        %compute solution
    51         if ~models.dh.parameters.qmu_analysis,
     51        if ~models.dh.parameters.QmuAnalysis,
    5252                if md.control_analysis,
    5353                        %change control_steady to 1
    54                         models.dh.parameters.control_steady=1;
    55                         models.ds.parameters.control_steady=1;
     54                        models.dh.parameters.ControlSteady=1;
     55                        models.ds.parameters.ControlSteady=1;
    5656                        %launch core of control solution.
    5757                        results=control_core(models,inputs);

  • TabularUnified issm/trunk/src/m/solutions/jpl/steadystate_core.m

    r3481 r3796  
    1616
    1717%recover parameters common to all solutions
    18 verbose=m_dhu.parameters.verbose;
    19 dim=m_dhu.parameters.dim;
    20 eps_rel=m_dhu.parameters.eps_rel;
    21 ishutter=m_dhu.parameters.ishutter;
    22 ismacayealpattyn=m_dh.parameters.ismacayealpattyn;
    23 isstokes=m_ds.parameters.isstokes;
     18verbose=m_dhu.parameters.Verbose;
     19dim=m_dhu.parameters.Dim;
     20eps_rel=m_dhu.parameters.EpsRel;
     21ishutter=m_dhu.parameters.IsHutter;
     22ismacayealpattyn=m_dh.parameters.IsMacAyealPattyn;
     23isstokes=m_ds.parameters.IsStokes;
    2424
    2525%convergence
    … …  
    4040        %first compute temperature at steady state.
    4141        if (step>1),
    42                 inputs=add(inputs,'velocity',results_diagnostic.u_g,'doublevec',ndof,m_t.parameters.numberofnodes);
     42                inputs=add(inputs,'velocity',results_diagnostic.u_g,'doublevec',ndof,m_t.parameters.NumberOfNodes);
    4343        end
    4444        results_thermal=thermal_core(models,inputs);
    … …  
    4747        %Compute depth averaged temperature
    4848        temperature_average=FieldDepthAverage(m_t.elements,m_t.nodes,m_t.vertices,m_t.loads,m_t.materials,m_t.parameters,results_thermal.t_g,'temperature');
    49         inputs=add(inputs,'temperature_average',temperature_average,'doublevec',1,m_t.parameters.numberofnodes);
    50         inputs=add(inputs,'temperature',results_thermal.t_g,'doublevec',1,m_t.parameters.numberofnodes);
     49        inputs=add(inputs,'temperature_average',temperature_average,'doublevec',1,m_t.parameters.NumberOfNodes);
     50        inputs=add(inputs,'temperature',results_thermal.t_g,'doublevec',1,m_t.parameters.NumberOfNodes);
    5151       
    5252        %now compute diagnostic velocity using the steady state temperature.

  • TabularUnified issm/trunk/src/m/solutions/jpl/thermal.m

    r2330 r3796  
    2323        displaystring(md.verbose,'\n%s',['setup inputs...']);
    2424        inputs=inputlist;
    25         inputs=add(inputs,'velocity',models.t.parameters.u_g,'doublevec',3,models.t.parameters.numberofnodes);
    26         inputs=add(inputs,'pressure',models.t.parameters.p_g,'doublevec',1,models.t.parameters.numberofnodes);
    27         inputs=add(inputs,'dt',models.t.parameters.dt*models.t.parameters.yts,'double');
     25        inputs=add(inputs,'velocity',models.t.parameters.u_g,'doublevec',3,models.t.parameters.NumberOfNodes);
     26        inputs=add(inputs,'pressure',models.t.parameters.p_g,'doublevec',1,models.t.parameters.NumberOfNodes);
     27        inputs=add(inputs,'dt',models.t.parameters.Dt*models.t.parameters.Yts,'double');
    2828
    2929        %compute solution
    30         if ~models.t.parameters.qmu_analysis,
     30        if ~models.t.parameters.QmuAnalysis,
    3131                %launch core of diagnostic solution.
    3232                results=thermal_core(models,inputs);

  • TabularUnified issm/trunk/src/m/solutions/jpl/thermal_core.m

    r2330 r3796  
    99m_m=models.m;
    1010
    11 if m_t.parameters.dt==0,
     11if m_t.parameters.Dt==0,
    1212
    1313        results.time=0;
    1414        results.step=1;
    1515
    16         displaystring(m_t.parameters.verbose,'\n%s',['computing temperatures...']);
     16        displaystring(m_t.parameters.Verbose,'\n%s',['computing temperatures...']);
    1717        [results.t_g m_t.loads melting_offset]=thermal_core_nonlinear(m_t,inputs,ThermalAnalysisEnum(),NoneAnalysisEnum());
    1818
    19         displaystring(m_t.parameters.verbose,'\n%s',['computing melting...']);
     19        displaystring(m_t.parameters.Verbose,'\n%s',['computing melting...']);
    2020        inputs=add(inputs,'melting_offset',melting_offset,'double');
    21         inputs=add(inputs,'temperature',results.t_g,'doublevec',1,m_t.parameters.numberofnodes);
     21        inputs=add(inputs,'temperature',results.t_g,'doublevec',1,m_t.parameters.NumberOfNodes);
    2222        results.m_g=diagnostic_core_linear(m_m,inputs,MeltingAnalysisEnum(),NoneAnalysisEnum());
    2323
    … …  
    2727        t_g=m_t.parameters.t_g;
    2828        m_g=m_m.parameters.m_g;
    29         nsteps=m_t.parameters.ndt/m_t.parameters.dt;
     29        nsteps=m_t.parameters.Ndt/m_t.parameters.Dt;
    3030
    3131        %initialize temperature and melting
    … …  
    3737        for n=1:nsteps,
    3838
    39                 displaystring(m_t.parameters.verbose,'\n%s%i/%i\n','time step: ',n,nsteps);
     39                displaystring(m_t.parameters.Verbose,'\n%s%i/%i\n','time step: ',n,nsteps);
    4040                results(n+1).step=n+1;
    41                 results(n+1).time=n*m_t.parameters.dt;
     41                results(n+1).time=n*m_t.parameters.Dt;
    4242
    43                 displaystring(m_t.parameters.verbose,'\n%s',['    computing temperatures...']);
    44                 inputs=add(inputs,'temperature',results(n).t_g,'doublevec',1,m_t.parameters.numberofnodes);
     43                displaystring(m_t.parameters.Verbose,'\n%s',['    computing temperatures...']);
     44                inputs=add(inputs,'temperature',results(n).t_g,'doublevec',1,m_t.parameters.NumberOfNodes);
    4545                [results(n+1).t_g m_t.loads melting_offset]=thermal_core_nonlinear(m_t,inputs,ThermalAnalysisEnum(),NoneAnalysisEnum());
    4646
    47                 displaystring(m_t.parameters.verbose,'\n%s',['    computing melting...']);
    48                 inputs=add(inputs,'temperature',results(n+1).t_g,'doublevec',1,m_t.parameters.numberofnodes);
     47                displaystring(m_t.parameters.Verbose,'\n%s',['    computing melting...']);
     48                inputs=add(inputs,'temperature',results(n+1).t_g,'doublevec',1,m_t.parameters.NumberOfNodes);
    4949                inputs=add(inputs,'melting_offset',melting_offset,'double');
    5050                results(n+1).m_g=diagnostic_core_linear(m_m,inputs,MeltingAnalysisEnum(),NoneAnalysisEnum());

  • TabularUnified issm/trunk/src/m/solutions/jpl/thermal_core_nonlinear.m

    r3493 r3796  
    1818        m.parameters.kflag=1; m.parameters.pflag=1;
    1919
    20         displaystring(m.parameters.verbose,'\n%s',['   starting direct shooting method']);
     20        displaystring(m.parameters.Verbose,'\n%s',['   starting direct shooting method']);
    2121       
    2222        while(~converged),
    … …  
    2626
    2727                %system matrices
    28                 if ~m.parameters.lowmem
     28                if ~m.parameters.Lowmem
    2929                        if count==1
    30                                 displaystring(m.parameters.verbose,'%s',['   system matrices']);
     30                                displaystring(m.parameters.Verbose,'%s',['   system matrices']);
    3131                                [K_gg_nopenalty, p_g_nopenalty]=SystemMatrices(m.elements,m.nodes,m.vertices,loads,m.materials,m.parameters,inputs,analysis_type,sub_analysis_type);
    3232                        end
    33                         displaystring(m.parameters.verbose,'%s',['   penalty system matrices']);
     33                        displaystring(m.parameters.Verbose,'%s',['   penalty system matrices']);
    3434                        [K_gg , p_g, melting_offset]=PenaltySystemMatrices(K_gg_nopenalty,p_g_nopenalty,m.elements,m.nodes,m.vertices,loads,m.materials,m.parameters,inputs,analysis_type,sub_analysis_type);
    3535                else
    36                         displaystring(m.parameters.verbose,'%s',['   system matrices']);
     36                        displaystring(m.parameters.Verbose,'%s',['   system matrices']);
    3737                        [K_gg , p_g]=SystemMatrices(m.elements,m.nodes,m.vertices,loads,m.materials,m.parameters,inputs,analysis_type,sub_analysis_type);
    38                         displaystring(m.parameters.verbose,'%s',['   penalty system matrices']);
     38                        displaystring(m.parameters.Verbose,'%s',['   penalty system matrices']);
    3939                        [K_gg , p_g, melting_offset]=PenaltySystemMatrices(K_gg,p_g,m.elements,m.nodes,m.vertices,loads,m.materials,m.parameters,inputs,analysis_type,sub_analysis_type);
    4040                end
    … …  
    4747
    4848                %Solve 
    49                 displaystring(m.parameters.verbose,'%s%g','   condition number of stiffness matrix: ',condest(K_ff));
     49                displaystring(m.parameters.Verbose,'%s%g','   condition number of stiffness matrix: ',condest(K_ff));
    5050                [t_f]=Solver(K_ff,p_f,[],m.parameters);
    5151
    5252                %Merge back to g set
    53                 displaystring(m.parameters.verbose,'%s',['   merging solution back to g set']);
     53                displaystring(m.parameters.Verbose,'%s',['   merging solution back to g set']);
    5454                [t_g]= Mergesolutionfromftog( t_f, m.Gmn, m.ys, m.nodesets );
    5555
    5656                %Update inputs in datasets
    57                 inputs=add(inputs,'temperature',t_g,'doublevec',m.parameters.numberofdofspernode,m.parameters.numberofnodes);
    58                 displaystring(m.parameters.verbose,'%s',['   update inputs']);
     57                inputs=add(inputs,'temperature',t_g,'doublevec',m.parameters.NumberOfDofsPerNode,m.parameters.NumberOfNodes);
     58                displaystring(m.parameters.Verbose,'%s',['   update inputs']);
    5959                [m.elements,m.nodes,m.vertices, loads,m.materials,m.parameters]=UpdateFromInputs(m.elements,m.nodes,m.vertices, loads,m.materials,m.parameters,inputs);
    6060       
    6161                %penalty constraints
    62                 displaystring(m.parameters.verbose,'%s',['   penalty constraints']);
     62                displaystring(m.parameters.Verbose,'%s',['   penalty constraints']);
    6363                [loads,constraints_converged,num_unstable_constraints] =PenaltyConstraints(m.elements,m.nodes,m.vertices,loads, m.materials,m.parameters,inputs,analysis_type,sub_analysis_type);
    6464       
    6565                if ~converged,
    66                         displaystring(m.parameters.verbose,'%s%i','   #unstable constraints ',num_unstable_constraints);
     66                        displaystring(m.parameters.Verbose,'%s%i','   #unstable constraints ',num_unstable_constraints);
    6767                       
    68                         if num_unstable_constraints<=m.parameters.min_thermal_constraints,
     68                        if num_unstable_constraints<=m.parameters.MinThermalConstraints,
    6969                                converged=1;
    7070                        end

  • TabularUnified issm/trunk/src/m/solutions/jpl/transient2d.m

    r3498 r3796  
    4141displaystring(md.verbose,'\n%s',['setup inputs...']);
    4242inputs=inputlist;
    43 inputs=add(inputs,'velocity',solution.u_g,'doublevec',2,models.p.parameters.numberofnodes);
    44 inputs=add(inputs,'melting',models.p.parameters.m_g,'doublevec',1,models.p.parameters.numberofnodes);
    45 inputs=add(inputs,'accumulation',models.p.parameters.a_g,'doublevec',1,models.p.parameters.numberofnodes);
    46 inputs=add(inputs,'dt',models.p.parameters.dt*models.p.parameters.yts,'double'); %change time from yrs to sec
     43inputs=add(inputs,'velocity',solution.u_g,'doublevec',2,models.p.parameters.NumberOfNodes);
     44inputs=add(inputs,'melting',models.p.parameters.m_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     45inputs=add(inputs,'accumulation',models.p.parameters.a_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     46inputs=add(inputs,'dt',models.p.parameters.Dt*models.p.parameters.Yts,'double'); %change time from yrs to sec
    4747
    4848% figure out number of dof: just for information purposes.
    … …  
    5050
    5151%first time step is given by model.
    52 dt=models.p.parameters.dt;
    53 finaltime=models.p.parameters.ndt;
     52dt=models.p.parameters.Dt;
     53finaltime=models.p.parameters.Ndt;
    5454time=dt;
    55 yts=models.p.parameters.yts;
     55yts=models.p.parameters.Yts;
    5656n=1; %counter
    5757
    … …  
    6464
    6565        %update inputs
    66         inputs=add(inputs,'thickness',solution(n).h_g,'doublevec',1,models.p.parameters.numberofnodes);
    67         inputs=add(inputs,'surface',solution(n).s_g,'doublevec',1,models.p.parameters.numberofnodes);
    68         inputs=add(inputs,'bed',solution(n).b_g,'doublevec',1,models.p.parameters.numberofnodes);
     66        inputs=add(inputs,'thickness',solution(n).h_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     67        inputs=add(inputs,'surface',solution(n).s_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     68        inputs=add(inputs,'bed',solution(n).b_g,'doublevec',1,models.p.parameters.NumberOfNodes);
    6969
    7070        %Deal with velocities.
    … …  
    7676        %compute new thickness
    7777        disp(sprintf('%s','   computing new thickness...'));
    78         inputs=add(inputs,'velocity',solution(n+1).u_g,'doublevec',2,models.p.parameters.numberofnodes);
     78        inputs=add(inputs,'velocity',solution(n+1).u_g,'doublevec',2,models.p.parameters.NumberOfNodes);
    7979        rawresults=prognostic_core(models,inputs,PrognosticAnalysisEnum(),NoneAnalysisEnum());
    8080        new_thickness=rawresults.h_g;

  • TabularUnified issm/trunk/src/m/solutions/jpl/transient3d.m

    r3505 r3796  
    5252displaystring(md.verbose,'\n%s',['setup inputs...']);
    5353inputs=inputlist;
    54 inputs=add(inputs,'velocity',results.u_g,'doublevec',3,models.p.parameters.numberofnodes);
    55 inputs=add(inputs,'melting',results.m_g,'doublevec',1,models.p.parameters.numberofnodes);
    56 inputs=add(inputs,'accumulation',results.a_g,'doublevec',1,models.p.parameters.numberofnodes);
    57 inputs=add(inputs,'dt',models.p.parameters.dt*models.p.parameters.yts,'double');
     54inputs=add(inputs,'velocity',results.u_g,'doublevec',3,models.p.parameters.NumberOfNodes);
     55inputs=add(inputs,'melting',results.m_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     56inputs=add(inputs,'accumulation',results.a_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     57inputs=add(inputs,'dt',models.p.parameters.Dt*models.p.parameters.Yts,'double');
    5858
    5959% figure out number of dof: just for information purposes.
    … …  
    6161
    6262%first time step is given by model.
    63 dt=models.p.parameters.dt;
    64 finaltime=models.p.parameters.ndt;
     63dt=models.p.parameters.Dt;
     64finaltime=models.p.parameters.Ndt;
    6565time=dt;
    66 yts=models.p.parameters.yts;
     66yts=models.p.parameters.Yts;
    6767n=1; %counter
    6868
    … …  
    7575
    7676        %update inputs
    77         inputs=add(inputs,'thickness',results(n).h_g,'doublevec',1,models.p.parameters.numberofnodes);
    78         inputs=add(inputs,'surface',results(n).s_g,'doublevec',1,models.p.parameters.numberofnodes);
    79         inputs=add(inputs,'bed',results(n).b_g,'doublevec',1,models.p.parameters.numberofnodes);
    80         inputs=add(inputs,'velocity',results(n).u_g,'doublevec',3,models.p.parameters.numberofnodes);
    81         inputs=add(inputs,'pressure',results(n).p_g,'doublevec',1,models.p.parameters.numberofnodes);
    82         inputs=add(inputs,'temperature',results(n).t_g,'doublevec',1,models.t.parameters.numberofnodes);
     77        inputs=add(inputs,'thickness',results(n).h_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     78        inputs=add(inputs,'surface',results(n).s_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     79        inputs=add(inputs,'bed',results(n).b_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     80        inputs=add(inputs,'velocity',results(n).u_g,'doublevec',3,models.p.parameters.NumberOfNodes);
     81        inputs=add(inputs,'pressure',results(n).p_g,'doublevec',1,models.p.parameters.NumberOfNodes);
     82        inputs=add(inputs,'temperature',results(n).t_g,'doublevec',1,models.t.parameters.NumberOfNodes);
    8383
    8484        %Deal with temperature first
    8585        displaystring(md.verbose,'\n%s',['    computing temperatures...']);
    8686        [results(n+1).t_g models.t.loads melting_offset]=thermal_core_nonlinear(models.t,inputs,ThermalAnalysisEnum(),TransientAnalysisEnum());
    87         inputs=add(inputs,'temperature',results(n+1).t_g,'doublevec',1,models.t.parameters.numberofnodes);
     87        inputs=add(inputs,'temperature',results(n+1).t_g,'doublevec',1,models.t.parameters.NumberOfNodes);
    8888       
    8989        displaystring(md.verbose,'\n%s',['    computing melting...']);
    … …  
    9393        %Compute depth averaged temperature
    9494        temperature_average=FieldDepthAverage(models.t.elements,models.t.nodes,models.t.vertices,models.t.loads,models.t.materials,models.t.parameters,results(n+1).t_g,'temperature');
    95         inputs=add(inputs,'temperature_average',temperature_average,'doublevec',1,models.t.parameters.numberofnodes);
     95        inputs=add(inputs,'temperature_average',temperature_average,'doublevec',1,models.t.parameters.NumberOfNodes);
    9696
    9797        %Deal with velocities.
    … …  
    101101        %compute new thickness
    102102        displaystring(md.verbose,'\n%s',['    computing new thickness...']);
    103         inputs=add(inputs,'velocity',results(n+1).u_g,'doublevec',3,models.p.parameters.numberofnodes);
     103        inputs=add(inputs,'velocity',results(n+1).u_g,'doublevec',3,models.p.parameters.NumberOfNodes);
    104104        rawresults=prognostic_core(models,inputs,PrognosticAnalysisEnum(),NoneAnalysisEnum());
    105105        new_thickness=rawresults.h_g;
Note: See TracChangeset for help on using the changeset viewer.