Changeset 4739


Ignore:
Timestamp:
07/22/10 13:50:33 (15 years ago)
Author:
Mathieu Morlighem
Message:

merged prognostic and prognostic2 -> now use TriaRef

Location:
issm/trunk/src
Files:
5 deleted
270 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/EnumDefinitions/EnumAsString.cpp

    r4532 r4739  
    6868                case BalancedvelocitiesSolutionEnum : return "BalancedvelocitiesSolution";
    6969                case BalancedvelocitiesAnalysisEnum : return "BalancedvelocitiesAnalysis";
    70                 case Prognostic2SolutionEnum : return "Prognostic2Solution";
    71                 case Prognostic2AnalysisEnum : return "Prognostic2Analysis";
    7270                case PrognosticSolutionEnum : return "PrognosticSolution";
    7371                case PrognosticAnalysisEnum : return "PrognosticAnalysis";
     
    226224                case P0Enum : return "P0";
    227225                case P1Enum : return "P1";
     226                case P1DGEnum : return "P1DG";
    228227                case DoubleElementResultEnum : return "DoubleElementResult";
    229228                case TriaVertexElementResultEnum : return "TriaVertexElementResult";
  • issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h

    r4532 r4739  
    7272        BalancedvelocitiesSolutionEnum,
    7373        BalancedvelocitiesAnalysisEnum,
    74         Prognostic2SolutionEnum,
    75         Prognostic2AnalysisEnum,
    7674        PrognosticSolutionEnum,
    7775        PrognosticAnalysisEnum,
     
    256254        P0Enum,
    257255        P1Enum,
     256        P1DGEnum,
    258257        /*}}}*/
    259258        /*Results{{{1*/
  • issm/trunk/src/c/EnumDefinitions/StringAsEnum.cpp

    r4532 r4739  
    6666        else if (strcmp(name,"BalancedvelocitiesSolution")==0) return BalancedvelocitiesSolutionEnum;
    6767        else if (strcmp(name,"BalancedvelocitiesAnalysis")==0) return BalancedvelocitiesAnalysisEnum;
    68         else if (strcmp(name,"Prognostic2Solution")==0) return Prognostic2SolutionEnum;
    69         else if (strcmp(name,"Prognostic2Analysis")==0) return Prognostic2AnalysisEnum;
    7068        else if (strcmp(name,"PrognosticSolution")==0) return PrognosticSolutionEnum;
    7169        else if (strcmp(name,"PrognosticAnalysis")==0) return PrognosticAnalysisEnum;
     
    224222        else if (strcmp(name,"P0")==0) return P0Enum;
    225223        else if (strcmp(name,"P1")==0) return P1Enum;
     224        else if (strcmp(name,"P1DG")==0) return P1DGEnum;
    226225        else if (strcmp(name,"DoubleElementResult")==0) return DoubleElementResultEnum;
    227226        else if (strcmp(name,"TriaVertexElementResult")==0) return TriaVertexElementResultEnum;
  • issm/trunk/src/c/Makefile.am

    r4702 r4739  
    390390                                        ./modules/ModelProcessorx/Prognostic/CreateConstraintsPrognostic.cpp\
    391391                                        ./modules/ModelProcessorx/Prognostic/CreateLoadsPrognostic.cpp\
    392                                         ./modules/ModelProcessorx/Prognostic2/UpdateElementsPrognostic2.cpp\
    393                                         ./modules/ModelProcessorx/Prognostic2/CreateNodesPrognostic2.cpp\
    394                                         ./modules/ModelProcessorx/Prognostic2/CreateConstraintsPrognostic2.cpp\
    395                                         ./modules/ModelProcessorx/Prognostic2/CreateLoadsPrognostic2.cpp\
    396392                                        ./modules/ModelProcessorx/Balancedthickness/UpdateElementsBalancedthickness.cpp\
    397393                                        ./modules/ModelProcessorx/Balancedthickness/CreateNodesBalancedthickness.cpp\
     
    950946                                        ./modules/ModelProcessorx/Prognostic/CreateConstraintsPrognostic.cpp\
    951947                                        ./modules/ModelProcessorx/Prognostic/CreateLoadsPrognostic.cpp\
    952                                         ./modules/ModelProcessorx/Prognostic2/UpdateElementsPrognostic2.cpp\
    953                                         ./modules/ModelProcessorx/Prognostic2/CreateNodesPrognostic2.cpp\
    954                                         ./modules/ModelProcessorx/Prognostic2/CreateConstraintsPrognostic2.cpp\
    955                                         ./modules/ModelProcessorx/Prognostic2/CreateLoadsPrognostic2.cpp\
    956948                                        ./modules/ModelProcessorx/Balancedthickness/UpdateElementsBalancedthickness.cpp\
    957949                                        ./modules/ModelProcessorx/Balancedthickness/CreateNodesBalancedthickness.cpp\
     
    11361128                                        ./solutions/adjoint_core.cpp\
    11371129                                        ./solutions/prognostic_core.cpp\
    1138                                         ./solutions/prognostic2_core.cpp\
    11391130                                        ./solutions/balancedthickness_core.cpp\
    11401131                                        ./solutions/balancedthickness2_core.cpp\
  • issm/trunk/src/c/modules/ModelProcessorx/CreateDataSets.cpp

    r4578 r4739  
    9292                        break;
    9393
    94                 case Prognostic2AnalysisEnum:
    95                         CreateNodesPrognostic2(pnodes, iomodel,iomodel_handle);
    96                         CreateConstraintsPrognostic2(pconstraints,iomodel,iomodel_handle);
    97                         CreateLoadsPrognostic2(ploads,iomodel,iomodel_handle);
    98                         UpdateElementsPrognostic2(elements,iomodel,iomodel_handle,analysis_counter,analysis_type);
    99                         break;
    100 
    10194                case BalancedthicknessAnalysisEnum:
    10295                        CreateNodesBalancedthickness(pnodes, iomodel,iomodel_handle);
  • issm/trunk/src/c/modules/ModelProcessorx/ModelProcessorx.h

    r4578 r4739  
    7878void    UpdateElementsPrognostic(Elements* elements,IoModel* iomodel_handle,ConstDataHandle iomodel_handle,int analysis_counter,int analysis_type);
    7979
    80 /*prognostic2:*/
    81 void    CreateNodesPrognostic2(Nodes** pnodes,IoModel* iomodel_handle,ConstDataHandle iomodel_handle);
    82 void    CreateConstraintsPrognostic2(Constraints** pconstraints,IoModel* iomodel,ConstDataHandle iomodel_handle);
    83 void  CreateLoadsPrognostic2(Loads** ploads, IoModel* iomodel, ConstDataHandle iomodel_handle);
    84 void    UpdateElementsPrognostic2(Elements* elements,IoModel* iomodel_handle,ConstDataHandle iomodel_handle,int analysis_counter,int analysis_type);
    85 
    8680/*balancedthickness:*/
    8781void    CreateNodesBalancedthickness(Nodes** pnodes,IoModel* iomodel_handle,ConstDataHandle iomodel_handle);
  • issm/trunk/src/c/modules/ModelProcessorx/Prognostic/CreateConstraintsPrognostic.cpp

    r4236 r4739  
    1818       
    1919        /*Output*/
    20         Constraints* constraints = NULL;
    21         Spc*    spc  = NULL;
     20        Constraints *constraints = NULL;
     21        Spc         *spc         = NULL;
    2222
    2323        /*Recover pointer: */
     
    2727        if(!constraints) constraints = new Constraints(ConstraintsEnum);
    2828
    29         /*Fetch data: */
    30         IoModelFetchData(&iomodel->spcthickness,NULL,NULL,iomodel_handle,"spcthickness");
     29        /*Do not add constraints in DG*/
     30        if(!iomodel->prognostic_DG){
    3131
    32         /*Initialize counter*/
    33         count=0;
     32                /*Fetch data: */
     33                IoModelFetchData(&iomodel->spcthickness,NULL,NULL,iomodel_handle,"spcthickness");
    3434
    35         /*Create spcs from x,y,z, as well as the spc values on those spcs: */
    36         for (i=0;i<iomodel->numberofvertices;i++){
    37                 /*keep only this partition's nodes:*/
    38                 if((iomodel->my_vertices[i]==1)){
     35                /*Initialize counter*/
     36                count=0;
    3937
    40                         if ((int)iomodel->spcthickness[2*i]){
     38                /*Create spcs from x,y,z, as well as the spc values on those spcs: */
     39                for (i=0;i<iomodel->numberofvertices;i++){
     40                        /*keep only this partition's nodes:*/
     41                        if((iomodel->my_vertices[i]==1)){
    4142
    42                                 constraints->AddObject(new Spc(iomodel->constraintcounter+count+1,iomodel->nodecounter+i+1,1,*(iomodel->spcthickness+2*i+1),PrognosticAnalysisEnum));
    43                                 count++;
    44                         }
    45                 } //if((my_vertices[i]==1))
     43                                if ((int)iomodel->spcthickness[2*i]){
     44
     45                                        constraints->AddObject(new Spc(iomodel->constraintcounter+count+1,iomodel->nodecounter+i+1,1,*(iomodel->spcthickness+2*i+1),PrognosticAnalysisEnum));
     46                                        count++;
     47                                }
     48                        } //if((my_vertices[i]==1))
     49                }
     50
     51                /*Free data: */
     52                xfree((void**)&iomodel->spcthickness);
    4653        }
    4754
    48         /*Free data: */
    49         xfree((void**)&iomodel->spcthickness);
    50        
    5155        /*Assign output pointer: */
    5256        *pconstraints=constraints;
  • issm/trunk/src/c/modules/ModelProcessorx/Prognostic/CreateLoadsPrognostic.cpp

    r4236 r4739  
    1313void    CreateLoadsPrognostic(Loads** ploads, IoModel* iomodel,ConstDataHandle iomodel_handle){
    1414
     15        /*Intermediaries*/
     16        int i;
     17        int element;
     18
    1519        /*DataSet*/
    1620        Loads*    loads    = NULL;
     
    2226        if(!loads) loads = new Loads(LoadsEnum);
    2327
     28        /*Loads only in DG*/
     29        if (iomodel->prognostic_DG){
     30
     31                /*Get edges and elements*/
     32                IoModelFetchData(&iomodel->edges,&iomodel->numberofedges,NULL,iomodel_handle,"edges");
     33                IoModelFetchData(&iomodel->elements,NULL,NULL,iomodel_handle,"elements");
     34                IoModelFetchData(&iomodel->thickness,NULL,NULL,iomodel_handle,"thickness");
     35
     36                /*First load data:*/
     37                for (i=0;i<iomodel->numberofedges;i++){
     38
     39                        /*Get left and right elements*/
     40                        element=(int)iomodel->edges[4*i+2]-1; //edges are [node1 node2 elem1 elem2]
     41
     42                        /*Now, if this element is not in the partition, pass: */
     43                        if(!iomodel->my_elements[element]) continue;
     44
     45                        /* Add load */
     46                        loads->AddObject(new Numericalflux(iomodel->loadcounter+i+1,i,iomodel,PrognosticAnalysisEnum));
     47                }
     48
     49                /*Free data: */
     50                xfree((void**)&iomodel->edges);
     51                xfree((void**)&iomodel->elements);
     52                xfree((void**)&iomodel->thickness);
     53
     54        }
     55
    2456        /*Assign output pointer: */
    2557        *ploads=loads;
  • issm/trunk/src/c/modules/ModelProcessorx/Prognostic/CreateNodesPrognostic.cpp

    r4236 r4739  
    1616
    1717        /*Intermediary*/
    18         int i;
     18        int  i,j;
     19        int  node_id;
     20        int  vertex_id;
     21        int  io_index;
    1922        bool continuous_galerkin=true;
    2023
     
    2831        if(!nodes) nodes = new Nodes(NodesEnum);
    2932
    30         /*Continuous Galerkin partition of nodes: */
     33        /*Create partition of nodes: */
     34        if(iomodel->prognostic_DG) continuous_galerkin=false;
    3135        NodesPartitioning(&iomodel->my_nodes,iomodel->my_elements, iomodel->my_vertices, iomodel->my_bordervertices, iomodel, iomodel_handle,continuous_galerkin);
     36
     37        /*Check in 3d*/
     38        if(iomodel->prognostic_DG && iomodel->dim==3) ISSMERROR("DG 3d not implemented yet");
    3239
    3340        /*First fetch data: */
     
    3643                IoModelFetchData(&iomodel->uppernodes,NULL,NULL,iomodel_handle,"uppergrids");
    3744        }
     45        IoModelFetchData(&iomodel->elements,NULL,NULL,iomodel_handle,"elements");
    3846        IoModelFetchData(&iomodel->gridonbed,NULL,NULL,iomodel_handle,"gridonbed");
    3947        IoModelFetchData(&iomodel->gridonsurface,NULL,NULL,iomodel_handle,"gridonsurface");
     
    4149        IoModelFetchData(&iomodel->gridoniceshelf,NULL,NULL,iomodel_handle,"gridoniceshelf");
    4250
    43         for (i=0;i<iomodel->numberofvertices;i++){
     51        if(continuous_galerkin){
    4452
    45                 if(iomodel->my_vertices[i]){
     53                /*Build Nodes dataset (Continuous Galerkin)*/
     54                for (i=0;i<iomodel->numberofvertices;i++){
    4655
    47                         /*Add node to nodes dataset: */
    48                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,i+1,i,iomodel,PrognosticAnalysisEnum));
     56                        if(iomodel->my_vertices[i]){
    4957
     58                                /*Add node to nodes dataset: */
     59                                nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,i+1,i,iomodel,PrognosticAnalysisEnum));
     60
     61                        }
     62                }
     63        }
     64        else{
     65
     66                /*Build Nodes dataset -> 3 for each element (Discontinuous Galerkin)*/
     67                for (i=0;i<iomodel->numberofelements;i++){
     68                        for (j=0;j<3;j++){
     69
     70                                if(iomodel->my_nodes[3*i+j]){
     71
     72                                        //Get index of the vertex on which the current node is located
     73                                        vertex_id=(int)*(iomodel->elements+3*i+j); //(Matlab indexing)
     74                                        io_index=vertex_id-1;                      //(C indexing)
     75                                        ISSMASSERT(vertex_id>0 && vertex_id<=iomodel->numberofvertices);
     76
     77                                        //Compute Node id
     78                                        node_id=iomodel->nodecounter+3*i+j+1;
     79
     80                                        /*Add node to nodes dataset: */
     81                                        nodes->AddObject(new Node(node_id,node_id-1,vertex_id,io_index,iomodel,PrognosticAnalysisEnum));
     82
     83                                }
     84                        }
    5085                }
    5186        }
     
    5893        xfree((void**)&iomodel->gridonicesheet);
    5994        xfree((void**)&iomodel->gridoniceshelf);
     95        xfree((void**)&iomodel->elements);
    6096
    6197        /*Assign output pointer: */
  • issm/trunk/src/c/objects/Elements/Penta.cpp

    r4733 r4739  
    391391        else if (analysis_type==PrognosticAnalysisEnum){
    392392                InputUpdateFromSolutionPrognostic( solution);
    393         }
    394         else if (analysis_type==Prognostic2AnalysisEnum){
    395                 InputUpdateFromSolutionPrognostic2(solution);
    396393        }
    397394        else if (analysis_type==BalancedthicknessAnalysisEnum){
     
    26232620        /*Add thickness as inputs to the tria element: */
    26242621        this->inputs->AddInput(new PentaVertexInput(ThicknessEnum,values));
    2625 }
    2626 /*}}}*/
    2627 /*FUNCTION Penta::InputUpdateFromSolutionPrognostic2 {{{1*/
    2628 void  Penta::InputUpdateFromSolutionPrognostic2(double* solution){
    2629         ISSMERROR(" not supported yet!");
    26302622}
    26312623/*}}}*/
  • issm/trunk/src/c/objects/Elements/Penta.h

    r4733 r4739  
    182182                void    InputUpdateFromSolutionDiagnosticStokes( double* solutiong);
    183183                void    InputUpdateFromSolutionPrognostic( double* solutiong);
    184                 void    InputUpdateFromSolutionPrognostic2(double* solutiong);
    185184                void    InputUpdateFromBedSlopeX( double* solutiong);
    186185                void    InputUpdateFromBedSlopeY( double* solutiong);
  • issm/trunk/src/c/objects/Elements/Tria.cpp

    r4733 r4739  
    3434/*FUNCTION Tria::Tria(int id, int index, IoModel* iomodel,int nummodels){{{1*/
    3535Tria::Tria(int tria_id, int index, IoModel* iomodel,int nummodels)
    36 
    37         :TriaRef(P1Enum) //P1Enum: interpolation type
     36        :TriaRef(nummodels)
    3837        ,TriaHook(nummodels,index+1,iomodel->numberofelements+1) //index+1: matice id, iomodel->numberofelements+1: matpar id
    3938                                                                  { //i is the element index
     
    7372
    7473        //deal with TriaRef mother class
    75         tria->interpolation_type=this->interpolation_type;
     74        tria->element_type_list=(int*)xmalloc(this->numanalyses*sizeof(int));
     75        for(i=0;i<this->numanalyses;i++) tria->element_type_list[i]=this->element_type_list[i];
    7676
    7777        //deal with TriaHook mother class
     
    146146        int i;
    147147        int flaghook;
     148        int type;
    148149
    149150        /*recover marshalled_dataset: */
     
    153154         *object data (thanks to DataSet::Demarshall):*/
    154155        memcpy(&id,marshalled_dataset,sizeof(id));marshalled_dataset+=sizeof(id);
    155         memcpy(&interpolation_type,marshalled_dataset,sizeof(interpolation_type));marshalled_dataset+=sizeof(interpolation_type);
    156156        memcpy(&numanalyses,marshalled_dataset,sizeof(numanalyses));marshalled_dataset+=sizeof(numanalyses);
     157
     158        /*demarshall Ref: */
     159        this->element_type_list=(int*)xmalloc(this->numanalyses*sizeof(int));
     160        for(i=0;i<numanalyses;i++){ memcpy(&element_type_list[i],marshalled_dataset,sizeof(type));marshalled_dataset+=sizeof(type);}
    157161
    158162        /*allocate dynamic memory: */
     
    249253        /*marshall Tria data: */
    250254        memcpy(marshalled_dataset,&id,sizeof(id));marshalled_dataset+=sizeof(id);
    251         memcpy(marshalled_dataset,&interpolation_type,sizeof(interpolation_type));marshalled_dataset+=sizeof(interpolation_type);
    252255        memcpy(marshalled_dataset,&numanalyses,sizeof(numanalyses));marshalled_dataset+=sizeof(numanalyses);
     256
     257        /*Mershall Ref: */
     258        for(i=0;i<numanalyses;i++){
     259                memcpy(marshalled_dataset,&element_type_list[i],sizeof(element_type_list[i]));marshalled_dataset+=sizeof(element_type_list[i]);
     260        }
    253261
    254262 /*Marshall hooks: */
     
    303311        return sizeof(id)
    304312                +hnodes_size
    305                 +sizeof(interpolation_type)
    306313                +sizeof(numanalyses)
     314                +numanalyses*sizeof(int) //element_type_lists
    307315                +hmatice->MarshallSize()
    308316                +hmatpar->MarshallSize()
     
    382390        else if (analysis_type==PrognosticAnalysisEnum){
    383391                InputUpdateFromSolutionPrognostic( solution);
    384         }
    385         else if (analysis_type==Prognostic2AnalysisEnum){
    386                 InputUpdateFromSolutionPrognostic2(solution);
    387392        }
    388393        else if (analysis_type==BalancedthicknessAnalysisEnum){
     
    506511        parametersin->FindParam(&analysis_counter,AnalysisCounterEnum);
    507512
     513        /*Get Element type*/
     514        this->element_type=this->element_type_list[analysis_counter];
     515
    508516        /*Pick up the objects inside the hooks: */
    509517        if(this->hnodes[analysis_counter]) this->nodes=(Node**)this->hnodes[analysis_counter]->deliverp();
     
    524532        /*go into parameters and get the analysis_counter: */
    525533        parametersin->FindParam(&analysis_counter,AnalysisCounterEnum);
     534
     535        /*Get Element type*/
     536        this->element_type=this->element_type_list[analysis_counter];
    526537
    527538        /*Take care of hooking up all objects for this element, ie links the objects in the hooks to their respective
     
    669680        }
    670681        else if (analysis_type==PrognosticAnalysisEnum){
    671                 CreateKMatrixPrognostic( Kgg);
    672         }
    673         else if (analysis_type==Prognostic2AnalysisEnum){
    674                 CreateKMatrixPrognostic2(Kgg);
     682                if (GetElementType()==P1Enum)
     683                 CreateKMatrixPrognostic_CG( Kgg);
     684                else if (GetElementType()==P1DGEnum)
     685                 CreateKMatrixPrognostic_DG( Kgg);
     686                else
     687                 ISSMERROR("Element type %s not supported yet",EnumAsString(GetElementType()));
    675688        }
    676689        else if (analysis_type==BalancedthicknessAnalysisEnum){
     
    716729        }
    717730        else if (analysis_type==PrognosticAnalysisEnum){
    718                 CreatePVectorPrognostic( pg);
    719         }
    720         else if (analysis_type==Prognostic2AnalysisEnum){
    721                 CreatePVectorPrognostic2( pg);
     731                if (GetElementType()==P1Enum)
     732                 CreatePVectorPrognostic_CG( pg);
     733                else if (GetElementType()==P1DGEnum)
     734                CreatePVectorPrognostic_DG( pg);
     735                else
     736                 ISSMERROR("Element type %s not supported yet",EnumAsString(GetElementType()));
    722737        }
    723738        else if (analysis_type==BalancedthicknessAnalysisEnum){
     
    22412256        int    tria_vertex_ids[3];
    22422257        double nodeinputs[3];
     2258        int    tria_type;
    22432259
    22442260        /*Checks if debuging*/
     
    22472263        /*}}}*/
    22482264
     2265        /*Recover element type*/
     2266        if ((analysis_type==PrognosticAnalysisEnum && iomodel->prognostic_DG) || analysis_type==Balancedthickness2AnalysisEnum){
     2267
     2268                /*P1 Discontinuous Galerkin*/
     2269                tria_type=P1DGEnum;
     2270        }
     2271        else{
     2272                /*P1 Continuous Galerkin*/
     2273                tria_type=P1Enum;
     2274        }
     2275
     2276        /*Set TriaRef*/
     2277        this->SetElementType(tria_type,analysis_counter);
     2278
    22492279        /*Recover vertices ids needed to initialize inputs*/
    22502280        for(i=0;i<3;i++){
     
    22532283
    22542284        /*Recover nodes ids needed to initialize the node hook.*/
    2255         if (analysis_type==Prognostic2AnalysisEnum || analysis_type==Balancedthickness2AnalysisEnum){
     2285        if (tria_type==P1DGEnum){
    22562286                /*Discontinuous Galerkin*/
    22572287                tria_node_ids[0]=iomodel->nodecounter+3*index+1;
     
    33593389}
    33603390/*}}}*/
    3361 /*FUNCTION Tria::CreateKMatrixPrognostic {{{1*/
    3362 void  Tria::CreateKMatrixPrognostic(Mat Kgg){
     3391/*FUNCTION Tria::CreateKMatrixPrognostic_CG {{{1*/
     3392void  Tria::CreateKMatrixPrognostic_CG(Mat Kgg){
    33633393
    33643394        /* local declarations */
     
    35293559}
    35303560/*}}}*/
    3531 /*FUNCTION Tria::CreateKMatrixPrognostic2 {{{1*/
    3532 void  Tria::CreateKMatrixPrognostic2(Mat Kgg){
     3561/*FUNCTION Tria::CreateKMatrixPrognostic_DG {{{1*/
     3562void  Tria::CreateKMatrixPrognostic_DG(Mat Kgg){
    35333563
    35343564        /* local declarations */
     
    42864316}
    42874317/*}}}*/
    4288 /*FUNCTION Tria::CreatePVectorPrognostic {{{1*/
    4289 void  Tria::CreatePVectorPrognostic(Vec pg){
     4318/*FUNCTION Tria::CreatePVectorPrognostic_CG {{{1*/
     4319void  Tria::CreatePVectorPrognostic_CG(Vec pg){
    42904320
    42914321
     
    43684398}
    43694399/*}}}*/
    4370 /*FUNCTION Tria::CreatePVectorPrognostic2 {{{1*/
    4371 void  Tria::CreatePVectorPrognostic2(Vec pg){
     4400/*FUNCTION Tria::CreatePVectorPrognostic_DG {{{1*/
     4401void  Tria::CreatePVectorPrognostic_DG(Vec pg){
    43724402
    43734403        /* local declarations */
     
    49494979}
    49504980/*}}}*/
     4981/*FUNCTION Tria::GetElementType {{{1*/
     4982int Tria::GetElementType(){
     4983
     4984        /*return TriaRef field*/
     4985        return this->element_type;
     4986
     4987}
     4988/*}}}*/
    49514989/*FUNCTION Tria::GetDofList {{{1*/
    49524990void  Tria::GetDofList(int* doflist,int* pnumberofdofspernode){
     
    59105948void  Tria::InputUpdateFromSolutionPrognostic(double* solution){
    59115949
    5912         int i;
    5913 
    5914         const int    numvertices=3;
    5915         const int    numdofpervertex=1;
    5916         const int    numdof=numdofpervertex*numvertices;
    5917 
    5918         int          doflist[numdof];
    5919         double       values[numdof];
    5920 
    5921         int          dummy;
    5922 
    5923         /*Get dof list: */
    5924         GetDofList(&doflist[0],&dummy);
    5925 
    5926         /*Use the dof list to index into the solution vector: */
    5927         for(i=0;i<numdof;i++){
    5928                 values[i]=solution[doflist[i]];
    5929         }
    5930 
    5931         /*Add thickness as inputs to the tria element: */
    5932         this->inputs->AddInput(new TriaVertexInput(ThicknessEnum,values));
    5933 }
    5934 /*}}}*/
    5935 /*FUNCTION Tria::InputUpdateFromSolutionPrognostic2 {{{1*/
    5936 void  Tria::InputUpdateFromSolutionPrognostic2(double* solution){
    59375950        int i;
    59385951
  • issm/trunk/src/c/objects/Elements/Tria.h

    r4702 r4739  
    119119                void      CreateKMatrixDiagnosticSurfaceVert(Mat Kgg);
    120120                void      CreateKMatrixMelting(Mat Kgg);
    121                 void      CreateKMatrixPrognostic(Mat Kgg);
    122                 void      CreateKMatrixPrognostic2(Mat Kgg);
     121                void      CreateKMatrixPrognostic_CG(Mat Kgg);
     122                void      CreateKMatrixPrognostic_DG(Mat Kgg);
    123123                void      CreateKMatrixSlope(Mat Kgg);
    124124                void      CreateKMatrixThermal(Mat Kgg);
     
    129129                void      CreatePVectorDiagnosticHoriz(Vec pg);
    130130                void      CreatePVectorDiagnosticHutter(Vec pg);
    131                 void      CreatePVectorPrognostic(Vec pg);
    132                 void      CreatePVectorPrognostic2(Vec pg);
     131                void      CreatePVectorPrognostic_CG(Vec pg);
     132                void      CreatePVectorPrognostic_DG(Vec pg);
    133133                void      CreatePVectorSlope( Vec pg);
    134134                void      CreatePVectorThermalSheet( Vec pg);
    135135                void      CreatePVectorThermalShelf( Vec pg);
    136                 double    GetArea(void);
    137                 double    GetAreaCoordinate(double x, double y, int which_one);
     136                double  GetArea(void);
     137                double  GetAreaCoordinate(double x, double y, int which_one);
    138138                void      GetB(double* B, double* xyz_list, double* gauss_l1l2l3);
    139139                void      GetBPrime(double* Bprime, double* xyz_list, double* gauss_l1l2l3);
    140140                void      GetBPrime_prog(double* Bprime_prog, double* xyz_list, double* gauss_l1l2l3);
    141141                void      GetB_prog(double* B_prog, double* xyz_list, double* gauss_l1l2l3);
     142                int     GetElementType(void);
    142143                void      GetDofList(int* doflist,int* pnumberofdofs);
    143144                void      GetDofList1(int* doflist);
     
    166167                void      InputUpdateFromSolutionSurfaceSlopeY( double* solution);
    167168                void      InputUpdateFromSolutionPrognostic( double* solution);
    168                 void      InputUpdateFromSolutionPrognostic2(double* solution);
    169169                void      InputUpdateFromSolutionBalancedthickness( double* solution);
    170170                void      InputUpdateFromSolutionBalancedthickness2( double* solution);
  • issm/trunk/src/c/objects/Elements/TriaRef.cpp

    r4396 r4739  
    2323/*FUNCTION TriaRef::TriaRef(){{{1*/
    2424TriaRef::TriaRef(){
     25        this->element_type_list=NULL;
    2526}
    2627/*}}}*/
    27 /*FUNCTION TriaRef::TriaRef(int type){{{1*/
    28 TriaRef::TriaRef(int type){
    29         this->interpolation_type=type;
     28/*FUNCTION TriaRef::TriaRef(int* types,int nummodels){{{1*/
     29TriaRef::TriaRef(const int nummodels){
     30
     31        /*Only allocate pointer*/
     32        element_type_list=(int*)xmalloc(nummodels*sizeof(int));
     33
    3034}
    3135/*}}}*/
    3236/*FUNCTION TriaRef::~TriaRef(){{{1*/
    3337TriaRef::~TriaRef(){
     38        xfree((void**)&element_type_list);
    3439}
    3540/*}}}*/
     41
     42/*Management*/
     43/*FUNCTION TriaRef::SetElementType{{{1*/
     44void TriaRef::SetElementType(int type,int type_counter){
     45
     46        ISSMASSERT(type==P1Enum || type==P1DGEnum);
     47
     48        /*initialize element type*/
     49        this->element_type_list[type_counter]=type;
     50}
     51/*}}}*/
  • issm/trunk/src/c/objects/Elements/TriaRef.h

    r4396 r4739  
    1212
    1313        public:
    14                 int interpolation_type; //L1, L2, H1, H2, etc ...
     14                int* element_type_list; //P1CG, P1DG, MINI, P2...
     15                int  element_type;
    1516               
    16                 /*FUNCTION constructors, destructors {{{1*/
    1717                TriaRef();
    18                 TriaRef(int type);
     18                TriaRef(const int nummodels);
    1919                ~TriaRef();
    20                 /*}}}*/
    21                 /*FUNCTION numerics {{{1*/
    22                 /*}}}*/
     20
     21                void SetElementType(int type,int type_counter);
    2322
    2423};
  • issm/trunk/src/c/objects/IoModel.cpp

    r4546 r4739  
    168168        IoModelFetchData(&this->verbose,iomodel_handle,"verbose");
    169169        IoModelFetchData(&this->plot,iomodel_handle,"plot");
    170         IoModelFetchData(&this->artificial_diffusivity,iomodel_handle,"artificial_diffusivity");
     170        IoModelFetchData(&this->prognostic_DG,iomodel_handle,"prognostic_DG");
    171171        IoModelFetchData(&this->nsteps,iomodel_handle,"nsteps");
    172172        IoModelFetchData(&this->eps_cm,iomodel_handle,"eps_cm");
     
    195195        IoModelFetchData(&this->viscosity_overshoot,iomodel_handle,"viscosity_overshoot");
    196196        IoModelFetchData(&this->artdiff,iomodel_handle,"artificial_diffusivity");
     197        IoModelFetchData(&this->prognostic_DG,iomodel_handle,"prognostic_DG");
    197198        IoModelFetchData(&this->stokesreconditioning,iomodel_handle,"stokesreconditioning");
    198199        IoModelFetchData(&this->waitonlock,iomodel_handle,"waitonlock");
     
    308309        this->meanvel=0;
    309310        this->epsvel=0;
    310         this->artificial_diffusivity=0;
    311311        this->nsteps=0;
    312312        this->eps_cm=0;
     
    339339        this->viscosity_overshoot=0;
    340340        this->artdiff=0;
     341        this->prognostic_DG=0;
    341342        this->stokesreconditioning=0;
    342343        this->waitonlock=0;
  • issm/trunk/src/c/objects/IoModel.h

    r4382 r4739  
    109109                double  epsvel;
    110110                int     artdiff;
     111                int     prognostic_DG;
    111112                double  viscosity_overshoot;
    112113                double  stokesreconditioning;
     
    130131                double* weights;
    131132                double* cm_jump;
    132                 int     artificial_diffusivity;
    133133                int     nsteps;
    134134                double  eps_cm;
  • issm/trunk/src/c/objects/Loads/Numericalflux.cpp

    r4575 r4739  
    412412
    413413        /*recover parameters: */
    414         if (analysis_type==Prognostic2AnalysisEnum){
     414        if (analysis_type==PrognosticAnalysisEnum){
    415415                parameters->FindParam(&dt,DtEnum);
    416416        }
     
    529529
    530530        /*recover parameters: */
    531         if (analysis_type==Prognostic2AnalysisEnum){
     531        if (analysis_type==PrognosticAnalysisEnum){
    532532                parameters->FindParam(&dt,DtEnum);
    533533        }
     
    656656
    657657        /*recover parameters: */
    658         if (analysis_type==Prognostic2AnalysisEnum){
     658        if (analysis_type==PrognosticAnalysisEnum){
    659659                parameters->FindParam(&dt,DtEnum);
    660660        }
  • issm/trunk/src/c/shared/Dofs/DistributeNumDofs.cpp

    r4295 r4739  
    3737                numdofs=1;
    3838        }
    39         else if (analysis_type==Prognostic2AnalysisEnum){
    40                 numdofs=1;
    41         }
    4239        else if (analysis_type==BalancedthicknessAnalysisEnum){
    4340                numdofs=1;
  • issm/trunk/src/c/solutions/SolutionConfiguration.cpp

    r4439 r4739  
    6767                        analyses[0]=PrognosticAnalysisEnum;
    6868                        solutioncore=&prognostic_core;
    69                         break;
    70 
    71                 case Prognostic2SolutionEnum:
    72                         numanalyses=1;
    73                         analyses=(int*)xmalloc(numanalyses*sizeof(int));
    74                         analyses[0]=Prognostic2AnalysisEnum;
    75                         solutioncore=&prognostic2_core;
    7669                        break;
    7770
  • issm/trunk/src/c/solutions/solutions.h

    r4445 r4739  
    2222void control_core(FemModel* femmodel);
    2323void prognostic_core(FemModel* femmodel);
    24 void prognostic2_core(FemModel* femmodel);
    2524void balancedthickness_core(FemModel* femmodel);
    2625void balancedthickness2_core(FemModel* femmodel);
  • issm/trunk/src/m/classes/@model/model.m

    r4429 r4739  
    180180        md.ndt=0;
    181181        md.artificial_diffusivity=0;
     182        md.prognostic_DG=0;
    182183        md.timestepping=0;
    183184        md.deltaH=0;
  • issm/trunk/src/m/classes/@model/setdefaultparameters.m

    r4321 r4739  
    131131md.artificial_diffusivity=1;
    132132
     133%Discontinuous Galerkin is used in prognostic to avoid
     134%numerical wiggles of the solution and conserve mass
     135md.prognostic_DG=0;
     136
    133137%do we use an adaptative time stepping?
    134138md.timestepping=0;
  • issm/trunk/src/m/classes/public/marshall.m

    r3994 r4739  
    111111WriteData(fid,md.plot,'Integer','plot');
    112112WriteData(fid,md.artificial_diffusivity,'Integer','artificial_diffusivity');
     113WriteData(fid,md.prognostic_DG,'Integer','prognostic_DG');
    113114WriteData(fid,md.nsteps,'Integer','nsteps');
    114115WriteData(fid,md.eps_cm,'Scalar','eps_cm');
  • issm/trunk/src/m/classes/public/process_solve_options.m

    r4476 r4739  
    2424
    2525%check solution type is supported
    26 if ~ismember(analysis_type,[DiagnosticSolutionEnum,PrognosticSolutionEnum,Prognostic2SolutionEnum,ThermalSolutionEnum,...
     26if ~ismember(analysis_type,[DiagnosticSolutionEnum,PrognosticSolutionEnum,ThermalSolutionEnum,...
    2727                SteadystateSolutionEnum,ParametersSolutionEnum,Transient2DSolutionEnum,Transient3DSolutionEnum...
    2828                BalancedthicknessSolutionEnum,Balancedthickness2SolutionEnum,BalancedvelocitiesSolutionEnum,BedSlopeSolutionEnum,SurfaceSlopeSolutionEnum]),
  • issm/trunk/src/m/enum/AbsoluteEnum.m

    r4521 r4739  
    99%      macro=AbsoluteEnum()
    1010
    11 macro=229;
     11macro=228;
  • issm/trunk/src/m/enum/AccumulationRateEnum.m

    r4475 r4739  
    99%      macro=AccumulationRateEnum()
    1010
    11 macro=123;
     11macro=121;
  • issm/trunk/src/m/enum/AdjointEnum.m

    r4529 r4739  
    99%      macro=AdjointEnum()
    1010
    11 macro=241;
     11macro=240;
  • issm/trunk/src/m/enum/AdjointpEnum.m

    r4475 r4739  
    99%      macro=AdjointpEnum()
    1010
    11 macro=127;
     11macro=125;
  • issm/trunk/src/m/enum/AdjointxEnum.m

    r4475 r4739  
    99%      macro=AdjointxEnum()
    1010
    11 macro=124;
     11macro=122;
  • issm/trunk/src/m/enum/AdjointyEnum.m

    r4475 r4739  
    99%      macro=AdjointyEnum()
    1010
    11 macro=125;
     11macro=123;
  • issm/trunk/src/m/enum/AdjointzEnum.m

    r4475 r4739  
    99%      macro=AdjointzEnum()
    1010
    11 macro=126;
     11macro=124;
  • issm/trunk/src/m/enum/AirEnum.m

    r4475 r4739  
    99%      macro=AirEnum()
    1010
    11 macro=121;
     11macro=119;
  • issm/trunk/src/m/enum/ArtDiffEnum.m

    r4475 r4739  
    99%      macro=ArtDiffEnum()
    1010
    11 macro=128;
     11macro=126;
  • issm/trunk/src/m/enum/BeamEnum.m

    r4475 r4739  
    99%      macro=BeamEnum()
    1010
    11 macro=72;
     11macro=70;
  • issm/trunk/src/m/enum/BeamVertexElementResultEnum.m

    r4521 r4739  
    99%      macro=BeamVertexElementResultEnum()
    1010
    11 macro=214;
     11macro=213;
  • issm/trunk/src/m/enum/BeamVertexInputEnum.m

    r4475 r4739  
    99%      macro=BeamVertexInputEnum()
    1010
    11 macro=95;
     11macro=93;
  • issm/trunk/src/m/enum/BedEnum.m

    r4475 r4739  
    99%      macro=BedEnum()
    1010
    11 macro=129;
     11macro=127;
  • issm/trunk/src/m/enum/BedSlopeXEnum.m

    r4475 r4739  
    99%      macro=BedSlopeXEnum()
    1010
    11 macro=130;
     11macro=128;
  • issm/trunk/src/m/enum/BedSlopeYEnum.m

    r4475 r4739  
    99%      macro=BedSlopeYEnum()
    1010
    11 macro=131;
     11macro=129;
  • issm/trunk/src/m/enum/BetaEnum.m

    r4529 r4739  
    99%      macro=BetaEnum()
    1010
    11 macro=236;
     11macro=235;
  • issm/trunk/src/m/enum/BoolExternalResultEnum.m

    r4521 r4739  
    99%      macro=BoolExternalResultEnum()
    1010
    11 macro=215;
     11macro=214;
  • issm/trunk/src/m/enum/BoolInputEnum.m

    r4475 r4739  
    99%      macro=BoolInputEnum()
    1010
    11 macro=96;
     11macro=94;
  • issm/trunk/src/m/enum/BoolParamEnum.m

    r4475 r4739  
    99%      macro=BoolParamEnum()
    1010
    11 macro=103;
     11macro=101;
  • issm/trunk/src/m/enum/BoundaryEnum.m

    r4475 r4739  
    99%      macro=BoundaryEnum()
    1010
    11 macro=132;
     11macro=130;
  • issm/trunk/src/m/enum/CmGradientEnum.m

    r4529 r4739  
    99%      macro=CmGradientEnum()
    1010
    11 macro=237;
     11macro=236;
  • issm/trunk/src/m/enum/CmJumpEnum.m

    r4529 r4739  
    99%      macro=CmJumpEnum()
    1010
    11 macro=238;
     11macro=237;
  • issm/trunk/src/m/enum/CmMaxDmpSlopeEnum.m

    r4475 r4739  
    99%      macro=CmMaxDmpSlopeEnum()
    1010
    11 macro=133;
     11macro=131;
  • issm/trunk/src/m/enum/CmMaxDmpValueEnum.m

    r4475 r4739  
    99%      macro=CmMaxDmpValueEnum()
    1010
    11 macro=134;
     11macro=132;
  • issm/trunk/src/m/enum/CmMaxEnum.m

    r4529 r4739  
    99%      macro=CmMaxEnum()
    1010
    11 macro=239;
     11macro=238;
  • issm/trunk/src/m/enum/CmMinDmpSlopeEnum.m

    r4475 r4739  
    99%      macro=CmMinDmpSlopeEnum()
    1010
    11 macro=135;
     11macro=133;
  • issm/trunk/src/m/enum/CmMinDmpValueEnum.m

    r4475 r4739  
    99%      macro=CmMinDmpValueEnum()
    1010
    11 macro=136;
     11macro=134;
  • issm/trunk/src/m/enum/CmMinEnum.m

    r4529 r4739  
    99%      macro=CmMinEnum()
    1010
    11 macro=240;
     11macro=239;
  • issm/trunk/src/m/enum/CmNoiseDmpEnum.m

    r4475 r4739  
    99%      macro=CmNoiseDmpEnum()
    1010
    11 macro=137;
     11macro=135;
  • issm/trunk/src/m/enum/CollapseEnum.m

    r4475 r4739  
    99%      macro=CollapseEnum()
    1010
    11 macro=138;
     11macro=136;
  • issm/trunk/src/m/enum/ConnectivityEnum.m

    r4529 r4739  
    99%      macro=ConnectivityEnum()
    1010
    11 macro=244;
     11macro=243;
  • issm/trunk/src/m/enum/ConstantEnum.m

    r4475 r4739  
    99%      macro=ConstantEnum()
    1010
    11 macro=139;
     11macro=137;
  • issm/trunk/src/m/enum/ControlParameterEnum.m

    r4529 r4739  
    99%      macro=ControlParameterEnum()
    1010
    11 macro=245;
     11macro=244;
  • issm/trunk/src/m/enum/ControlSteadyEnum.m

    r4529 r4739  
    99%      macro=ControlSteadyEnum()
    1010
    11 macro=246;
     11macro=245;
  • issm/trunk/src/m/enum/ControlTypeEnum.m

    r4475 r4739  
    99%      macro=ControlTypeEnum()
    1010
    11 macro=140;
     11macro=138;
  • issm/trunk/src/m/enum/ConvergedEnum.m

    r4475 r4739  
    99%      macro=ConvergedEnum()
    1010
    11 macro=141;
     11macro=139;
  • issm/trunk/src/m/enum/DakotaParameterEnum.m

    r4529 r4739  
    99%      macro=DakotaParameterEnum()
    1010
    11 macro=247;
     11macro=246;
  • issm/trunk/src/m/enum/DhDtEnum.m

    r4475 r4739  
    99%      macro=DhDtEnum()
    1010
    11 macro=142;
     11macro=140;
  • issm/trunk/src/m/enum/DimEnum.m

    r4529 r4739  
    99%      macro=DimEnum()
    1010
    11 macro=248;
     11macro=247;
  • issm/trunk/src/m/enum/DofIndexingEnum.m

    r4475 r4739  
    99%      macro=DofIndexingEnum()
    1010
    11 macro=73;
     11macro=71;
  • issm/trunk/src/m/enum/DoubleElementResultEnum.m

    r4521 r4739  
    99%      macro=DoubleElementResultEnum()
    1010
    11 macro=210;
     11macro=209;
  • issm/trunk/src/m/enum/DoubleExternalResultEnum.m

    r4521 r4739  
    99%      macro=DoubleExternalResultEnum()
    1010
    11 macro=216;
     11macro=215;
  • issm/trunk/src/m/enum/DoubleInputEnum.m

    r4475 r4739  
    99%      macro=DoubleInputEnum()
    1010
    11 macro=97;
     11macro=95;
  • issm/trunk/src/m/enum/DoubleMatExternalResultEnum.m

    r4521 r4739  
    99%      macro=DoubleMatExternalResultEnum()
    1010
    11 macro=218;
     11macro=217;
  • issm/trunk/src/m/enum/DoubleMatParamEnum.m

    r4475 r4739  
    99%      macro=DoubleMatParamEnum()
    1010
    11 macro=104;
     11macro=102;
  • issm/trunk/src/m/enum/DoubleParamEnum.m

    r4475 r4739  
    99%      macro=DoubleParamEnum()
    1010
    11 macro=105;
     11macro=103;
  • issm/trunk/src/m/enum/DoubleVecExternalResultEnum.m

    r4521 r4739  
    99%      macro=DoubleVecExternalResultEnum()
    1010
    11 macro=217;
     11macro=216;
  • issm/trunk/src/m/enum/DoubleVecParamEnum.m

    r4475 r4739  
    99%      macro=DoubleVecParamEnum()
    1010
    11 macro=106;
     11macro=104;
  • issm/trunk/src/m/enum/DragCoefficientEnum.m

    r4475 r4739  
    99%      macro=DragCoefficientEnum()
    1010
    11 macro=143;
     11macro=141;
  • issm/trunk/src/m/enum/DragPEnum.m

    r4475 r4739  
    99%      macro=DragPEnum()
    1010
    11 macro=144;
     11macro=142;
  • issm/trunk/src/m/enum/DragQEnum.m

    r4475 r4739  
    99%      macro=DragQEnum()
    1010
    11 macro=145;
     11macro=143;
  • issm/trunk/src/m/enum/DragTypeEnum.m

    r4475 r4739  
    99%      macro=DragTypeEnum()
    1010
    11 macro=146;
     11macro=144;
  • issm/trunk/src/m/enum/DtEnum.m

    r4475 r4739  
    99%      macro=DtEnum()
    1010
    11 macro=147;
     11macro=145;
  • issm/trunk/src/m/enum/ElementEnum.m

    r4475 r4739  
    99%      macro=ElementEnum()
    1010
    11 macro=71;
     11macro=69;
  • issm/trunk/src/m/enum/ElementOnBedEnum.m

    r4475 r4739  
    99%      macro=ElementOnBedEnum()
    1010
    11 macro=148;
     11macro=146;
  • issm/trunk/src/m/enum/ElementOnIceShelfEnum.m

    r4475 r4739  
    99%      macro=ElementOnIceShelfEnum()
    1010
    11 macro=149;
     11macro=147;
  • issm/trunk/src/m/enum/ElementOnSurfaceEnum.m

    r4475 r4739  
    99%      macro=ElementOnSurfaceEnum()
    1010
    11 macro=150;
     11macro=148;
  • issm/trunk/src/m/enum/ElementOnWaterEnum.m

    r4475 r4739  
    99%      macro=ElementOnWaterEnum()
    1010
    11 macro=151;
     11macro=149;
  • issm/trunk/src/m/enum/ElementPropertiesEnum.m

    r4475 r4739  
    99%      macro=ElementPropertiesEnum()
    1010
    11 macro=74;
     11macro=72;
  • issm/trunk/src/m/enum/ElementResultEnum.m

    r4475 r4739  
    99%      macro=ElementResultEnum()
    1010
    11 macro=112;
     11macro=110;
  • issm/trunk/src/m/enum/EnumAsString.m

    r4529 r4739  
    6363        case BalancedvelocitiesSolutionEnum(), string='BalancedvelocitiesSolution'; return
    6464        case BalancedvelocitiesAnalysisEnum(), string='BalancedvelocitiesAnalysis'; return
    65         case Prognostic2SolutionEnum(), string='Prognostic2Solution'; return
    66         case Prognostic2AnalysisEnum(), string='Prognostic2Analysis'; return
    6765        case PrognosticSolutionEnum(), string='PrognosticSolution'; return
    6866        case PrognosticAnalysisEnum(), string='PrognosticAnalysis'; return
     
    221219        case P0Enum(), string='P0'; return
    222220        case P1Enum(), string='P1'; return
     221        case P1DGEnum(), string='P1DG'; return
    223222        case DoubleElementResultEnum(), string='DoubleElementResult'; return
    224223        case TriaVertexElementResultEnum(), string='TriaVertexElementResult'; return
  • issm/trunk/src/m/enum/EpsAbsEnum.m

    r4529 r4739  
    99%      macro=EpsAbsEnum()
    1010
    11 macro=249;
     11macro=248;
  • issm/trunk/src/m/enum/EpsCmEnum.m

    r4529 r4739  
    99%      macro=EpsCmEnum()
    1010
    11 macro=250;
     11macro=249;
  • issm/trunk/src/m/enum/EpsRelEnum.m

    r4529 r4739  
    99%      macro=EpsRelEnum()
    1010
    11 macro=251;
     11macro=250;
  • issm/trunk/src/m/enum/EpsResEnum.m

    r4529 r4739  
    99%      macro=EpsResEnum()
    1010
    11 macro=252;
     11macro=251;
  • issm/trunk/src/m/enum/EpsVelEnum.m

    r4475 r4739  
    99%      macro=EpsVelEnum()
    1010
    11 macro=152;
     11macro=150;
  • issm/trunk/src/m/enum/ExternalResultEnum.m

    r4475 r4739  
    99%      macro=ExternalResultEnum()
    1010
    11 macro=113;
     11macro=111;
  • issm/trunk/src/m/enum/ExtrudeParamEnum.m

    r4529 r4739  
    99%      macro=ExtrudeParamEnum()
    1010
    11 macro=253;
     11macro=252;
  • issm/trunk/src/m/enum/FillEnum.m

    r4475 r4739  
    99%      macro=FillEnum()
    1010
    11 macro=153;
     11macro=151;
  • issm/trunk/src/m/enum/FitEnum.m

    r4475 r4739  
    99%      macro=FitEnum()
    1010
    11 macro=154;
     11macro=152;
  • issm/trunk/src/m/enum/FormulationEnum.m

    r4475 r4739  
    99%      macro=FormulationEnum()
    1010
    11 macro=63;
     11macro=61;
  • issm/trunk/src/m/enum/FractionIncrementEnum.m

    r4475 r4739  
    99%      macro=FractionIncrementEnum()
    1010
    11 macro=155;
     11macro=153;
  • issm/trunk/src/m/enum/FrictionEnum.m

    r4475 r4739  
    99%      macro=FrictionEnum()
    1010
    11 macro=156;
     11macro=154;
  • issm/trunk/src/m/enum/GeographyEnum.m

    r4475 r4739  
    99%      macro=GeographyEnum()
    1010
    11 macro=116;
     11macro=114;
  • issm/trunk/src/m/enum/GeothermalFluxEnum.m

    r4475 r4739  
    99%      macro=GeothermalFluxEnum()
    1010
    11 macro=157;
     11macro=155;
  • issm/trunk/src/m/enum/GradientEnum.m

    r4529 r4739  
    99%      macro=GradientEnum()
    1010
    11 macro=242;
     11macro=241;
  • issm/trunk/src/m/enum/GravityEnum.m

    r4521 r4739  
    99%      macro=GravityEnum()
    1010
    11 macro=232;
     11macro=231;
  • issm/trunk/src/m/enum/HeatCapacityEnum.m

    r4529 r4739  
    99%      macro=HeatCapacityEnum()
    1010
    11 macro=254;
     11macro=253;
  • issm/trunk/src/m/enum/HookEnum.m

    r4475 r4739  
    99%      macro=HookEnum()
    1010
    11 macro=70;
     11macro=68;
  • issm/trunk/src/m/enum/HutterFormulationEnum.m

    r4475 r4739  
    99%      macro=HutterFormulationEnum()
    1010
    11 macro=64;
     11macro=62;
  • issm/trunk/src/m/enum/IceEnum.m

    r4475 r4739  
    99%      macro=IceEnum()
    1010
    11 macro=120;
     11macro=118;
  • issm/trunk/src/m/enum/IceSheetEnum.m

    r4475 r4739  
    99%      macro=IceSheetEnum()
    1010
    11 macro=117;
     11macro=115;
  • issm/trunk/src/m/enum/IceShelfEnum.m

    r4475 r4739  
    99%      macro=IceShelfEnum()
    1010
    11 macro=118;
     11macro=116;
  • issm/trunk/src/m/enum/IcefrontEnum.m

    r4475 r4739  
    99%      macro=IcefrontEnum()
    1010
    11 macro=82;
     11macro=80;
  • issm/trunk/src/m/enum/InputEnum.m

    r4475 r4739  
    99%      macro=InputEnum()
    1010
    11 macro=94;
     11macro=92;
  • issm/trunk/src/m/enum/IntExternalResultEnum.m

    r4521 r4739  
    99%      macro=IntExternalResultEnum()
    1010
    11 macro=219;
     11macro=218;
  • issm/trunk/src/m/enum/IntInputEnum.m

    r4475 r4739  
    99%      macro=IntInputEnum()
    1010
    11 macro=98;
     11macro=96;
  • issm/trunk/src/m/enum/IntParamEnum.m

    r4475 r4739  
    99%      macro=IntParamEnum()
    1010
    11 macro=107;
     11macro=105;
  • issm/trunk/src/m/enum/InternalEnum.m

    r4475 r4739  
    99%      macro=InternalEnum()
    1010
    11 macro=158;
     11macro=156;
  • issm/trunk/src/m/enum/IsHutterEnum.m

    r4529 r4739  
    99%      macro=IsHutterEnum()
    1010
    11 macro=255;
     11macro=254;
  • issm/trunk/src/m/enum/IsMacAyealPattynEnum.m

    r4529 r4739  
    99%      macro=IsMacAyealPattynEnum()
    1010
    11 macro=256;
     11macro=255;
  • issm/trunk/src/m/enum/IsStokesEnum.m

    r4529 r4739  
    99%      macro=IsStokesEnum()
    1010
    11 macro=257;
     11macro=256;
  • issm/trunk/src/m/enum/JEnum.m

    r4521 r4739  
    99%      macro=JEnum()
    1010
    11 macro=222;
     11macro=221;
  • issm/trunk/src/m/enum/KflagEnum.m

    r4475 r4739  
    99%      macro=KflagEnum()
    1010
    11 macro=159;
     11macro=157;
  • issm/trunk/src/m/enum/LatentHeatEnum.m

    r4529 r4739  
    99%      macro=LatentHeatEnum()
    1010
    11 macro=258;
     11macro=257;
  • issm/trunk/src/m/enum/LoadEnum.m

    r4475 r4739  
    99%      macro=LoadEnum()
    1010
    11 macro=81;
     11macro=79;
  • issm/trunk/src/m/enum/LowmemEnum.m

    r4529 r4739  
    99%      macro=LowmemEnum()
    1010
    11 macro=259;
     11macro=258;
  • issm/trunk/src/m/enum/MacAyealFormulationEnum.m

    r4475 r4739  
    99%      macro=MacAyealFormulationEnum()
    1010
    11 macro=65;
     11macro=63;
  • issm/trunk/src/m/enum/MaterialEnum.m

    r4475 r4739  
    99%      macro=MaterialEnum()
    1010
    11 macro=91;
     11macro=89;
  • issm/trunk/src/m/enum/MaticeEnum.m

    r4475 r4739  
    99%      macro=MaticeEnum()
    1010
    11 macro=92;
     11macro=90;
  • issm/trunk/src/m/enum/MatparEnum.m

    r4475 r4739  
    99%      macro=MatparEnum()
    1010
    11 macro=93;
     11macro=91;
  • issm/trunk/src/m/enum/MaxIterEnum.m

    r4529 r4739  
    99%      macro=MaxIterEnum()
    1010
    11 macro=260;
     11macro=259;
  • issm/trunk/src/m/enum/MaxNonlinearIterationsEnum.m

    r4529 r4739  
    99%      macro=MaxNonlinearIterationsEnum()
    1010
    11 macro=261;
     11macro=260;
  • issm/trunk/src/m/enum/MaxPenetrationEnum.m

    r4475 r4739  
    99%      macro=MaxPenetrationEnum()
    1010
    11 macro=160;
     11macro=158;
  • issm/trunk/src/m/enum/MeanVelEnum.m

    r4475 r4739  
    99%      macro=MeanVelEnum()
    1010
    11 macro=161;
     11macro=159;
  • issm/trunk/src/m/enum/MelangeEnum.m

    r4475 r4739  
    99%      macro=MelangeEnum()
    1010
    11 macro=122;
     11macro=120;
  • issm/trunk/src/m/enum/MeltingAnalysisEnum.m

    r4475 r4739  
    99%      macro=MeltingAnalysisEnum()
    1010
    11 macro=56;
     11macro=54;
  • issm/trunk/src/m/enum/MeltingOffsetEnum.m

    r4475 r4739  
    99%      macro=MeltingOffsetEnum()
    1010
    11 macro=162;
     11macro=160;
  • issm/trunk/src/m/enum/MeltingPointEnum.m

    r4529 r4739  
    99%      macro=MeltingPointEnum()
    1010
    11 macro=262;
     11macro=261;
  • issm/trunk/src/m/enum/MeltingRateEnum.m

    r4475 r4739  
    99%      macro=MeltingRateEnum()
    1010
    11 macro=163;
     11macro=161;
  • issm/trunk/src/m/enum/Mesh2gridAnalysisEnum.m

    r4475 r4739  
    99%      macro=Mesh2gridAnalysisEnum()
    1010
    11 macro=57;
     11macro=55;
  • issm/trunk/src/m/enum/MinMechanicalConstraintsEnum.m

    r4529 r4739  
    99%      macro=MinMechanicalConstraintsEnum()
    1010
    11 macro=263;
     11macro=262;
  • issm/trunk/src/m/enum/MinThermalConstraintsEnum.m

    r4529 r4739  
    99%      macro=MinThermalConstraintsEnum()
    1010
    11 macro=264;
     11macro=263;
  • issm/trunk/src/m/enum/MixedLayerCapacityEnum.m

    r4529 r4739  
    99%      macro=MixedLayerCapacityEnum()
    1010
    11 macro=234;
     11macro=233;
  • issm/trunk/src/m/enum/NStepsEnum.m

    r4529 r4739  
    99%      macro=NStepsEnum()
    1010
    11 macro=265;
     11macro=264;
  • issm/trunk/src/m/enum/NdtEnum.m

    r4529 r4739  
    99%      macro=NdtEnum()
    1010
    11 macro=266;
     11macro=265;
  • issm/trunk/src/m/enum/NodeEnum.m

    r4475 r4739  
    99%      macro=NodeEnum()
    1010
    11 macro=79;
     11macro=77;
  • issm/trunk/src/m/enum/NodeOnBedEnum.m

    r4475 r4739  
    99%      macro=NodeOnBedEnum()
    1010
    11 macro=164;
     11macro=162;
  • issm/trunk/src/m/enum/NodeOnIceSheetEnum.m

    r4475 r4739  
    99%      macro=NodeOnIceSheetEnum()
    1010
    11 macro=165;
     11macro=163;
  • issm/trunk/src/m/enum/NodeOnIceShelfEnum.m

    r4475 r4739  
    99%      macro=NodeOnIceShelfEnum()
    1010
    11 macro=166;
     11macro=164;
  • issm/trunk/src/m/enum/NodeOnSurfaceEnum.m

    r4475 r4739  
    99%      macro=NodeOnSurfaceEnum()
    1010
    11 macro=167;
     11macro=165;
  • issm/trunk/src/m/enum/NodePropertiesEnum.m

    r4475 r4739  
    99%      macro=NodePropertiesEnum()
    1010
    11 macro=75;
     11macro=73;
  • issm/trunk/src/m/enum/NoneAnalysisEnum.m

    r4475 r4739  
    99%      macro=NoneAnalysisEnum()
    1010
    11 macro=62;
     11macro=60;
  • issm/trunk/src/m/enum/NoneFormulationEnum.m

    r4475 r4739  
    99%      macro=NoneFormulationEnum()
    1010
    11 macro=66;
     11macro=64;
  • issm/trunk/src/m/enum/NumOutputEnum.m

    r4529 r4739  
    99%      macro=NumOutputEnum()
    1010
    11 macro=267;
     11macro=266;
  • issm/trunk/src/m/enum/NumRiftsEnum.m

    r4529 r4739  
    99%      macro=NumRiftsEnum()
    1010
    11 macro=268;
     11macro=267;
  • issm/trunk/src/m/enum/NumberNodeToElementConnectivityEnum.m

    r4475 r4739  
    99%      macro=NumberNodeToElementConnectivityEnum()
    1010
    11 macro=168;
     11macro=166;
  • issm/trunk/src/m/enum/NumberOfDofsPerNodeEnum.m

    r4529 r4739  
    99%      macro=NumberOfDofsPerNodeEnum()
    1010
    11 macro=269;
     11macro=268;
  • issm/trunk/src/m/enum/NumberOfElementsEnum.m

    r4529 r4739  
    99%      macro=NumberOfElementsEnum()
    1010
    11 macro=270;
     11macro=269;
  • issm/trunk/src/m/enum/NumberOfNodesEnum.m

    r4529 r4739  
    99%      macro=NumberOfNodesEnum()
    1010
    11 macro=271;
     11macro=270;
  • issm/trunk/src/m/enum/NumberOfVerticesEnum.m

    r4529 r4739  
    99%      macro=NumberOfVerticesEnum()
    1010
    11 macro=272;
     11macro=271;
  • issm/trunk/src/m/enum/NumericalfluxEnum.m

    r4475 r4739  
    99%      macro=NumericalfluxEnum()
    1010
    11 macro=83;
     11macro=81;
  • issm/trunk/src/m/enum/ObjectEnum.m

    r4475 r4739  
    99%      macro=ObjectEnum()
    1010
    11 macro=69;
     11macro=67;
  • issm/trunk/src/m/enum/OldGradientEnum.m

    r4529 r4739  
    99%      macro=OldGradientEnum()
    1010
    11 macro=243;
     11macro=242;
  • issm/trunk/src/m/enum/OptScalEnum.m

    r4529 r4739  
    99%      macro=OptScalEnum()
    1010
    11 macro=273;
     11macro=272;
  • issm/trunk/src/m/enum/OutputFileNameEnum.m

    r4529 r4739  
    99%      macro=OutputFileNameEnum()
    1010
    11 macro=274;
     11macro=273;
  • issm/trunk/src/m/enum/P0Enum.m

    r4521 r4739  
    99%      macro=P0Enum()
    1010
    11 macro=208;
     11macro=206;
  • issm/trunk/src/m/enum/P1Enum.m

    r4521 r4739  
    99%      macro=P1Enum()
    1010
    11 macro=209;
     11macro=207;
  • issm/trunk/src/m/enum/ParamEnum.m

    r4475 r4739  
    99%      macro=ParamEnum()
    1010
    11 macro=102;
     11macro=100;
  • issm/trunk/src/m/enum/ParameterOutputEnum.m

    r4529 r4739  
    99%      macro=ParameterOutputEnum()
    1010
    11 macro=275;
     11macro=274;
  • issm/trunk/src/m/enum/ParametersAnalysisEnum.m

    r4475 r4739  
    99%      macro=ParametersAnalysisEnum()
    1010
    11 macro=59;
     11macro=57;
  • issm/trunk/src/m/enum/ParametersSolutionEnum.m

    r4475 r4739  
    99%      macro=ParametersSolutionEnum()
    1010
    11 macro=58;
     11macro=56;
  • issm/trunk/src/m/enum/PatchEnum.m

    r4521 r4739  
    99%      macro=PatchEnum()
    1010
    11 macro=223;
     11macro=222;
  • issm/trunk/src/m/enum/PatchNodesEnum.m

    r4521 r4739  
    99%      macro=PatchNodesEnum()
    1010
    11 macro=225;
     11macro=224;
  • issm/trunk/src/m/enum/PatchVerticesEnum.m

    r4521 r4739  
    99%      macro=PatchVerticesEnum()
    1010
    11 macro=224;
     11macro=223;
  • issm/trunk/src/m/enum/PattynFormulationEnum.m

    r4475 r4739  
    99%      macro=PattynFormulationEnum()
    1010
    11 macro=67;
     11macro=65;
  • issm/trunk/src/m/enum/PenaltyMeltingEnum.m

    r4529 r4739  
    99%      macro=PenaltyMeltingEnum()
    1010
    11 macro=276;
     11macro=275;
  • issm/trunk/src/m/enum/PenaltyOffsetEnum.m

    r4475 r4739  
    99%      macro=PenaltyOffsetEnum()
    1010
    11 macro=169;
     11macro=167;
  • issm/trunk/src/m/enum/PengridEnum.m

    r4475 r4739  
    99%      macro=PengridEnum()
    1010
    11 macro=84;
     11macro=82;
  • issm/trunk/src/m/enum/PenpairEnum.m

    r4475 r4739  
    99%      macro=PenpairEnum()
    1010
    11 macro=85;
     11macro=83;
  • issm/trunk/src/m/enum/PentaEnum.m

    r4475 r4739  
    99%      macro=PentaEnum()
    1010
    11 macro=76;
     11macro=74;
  • issm/trunk/src/m/enum/PentaVertexElementResultEnum.m

    r4521 r4739  
    99%      macro=PentaVertexElementResultEnum()
    1010
    11 macro=212;
     11macro=211;
  • issm/trunk/src/m/enum/PentaVertexInputEnum.m

    r4475 r4739  
    99%      macro=PentaVertexInputEnum()
    1010
    11 macro=99;
     11macro=97;
  • issm/trunk/src/m/enum/PetscMatParamEnum.m

    r4475 r4739  
    99%      macro=PetscMatParamEnum()
    1010
    11 macro=108;
     11macro=106;
  • issm/trunk/src/m/enum/PetscVecExternalResultEnum.m

    r4521 r4739  
    99%      macro=PetscVecExternalResultEnum()
    1010
    11 macro=220;
     11macro=219;
  • issm/trunk/src/m/enum/PetscVecParamEnum.m

    r4475 r4739  
    99%      macro=PetscVecParamEnum()
    1010
    11 macro=109;
     11macro=107;
  • issm/trunk/src/m/enum/PflagEnum.m

    r4475 r4739  
    99%      macro=PflagEnum()
    1010
    11 macro=170;
     11macro=168;
  • issm/trunk/src/m/enum/PressureEnum.m

    r4475 r4739  
    99%      macro=PressureEnum()
    1010
    11 macro=171;
     11macro=169;
  • issm/trunk/src/m/enum/PressureOldEnum.m

    r4475 r4739  
    99%      macro=PressureOldEnum()
    1010
    11 macro=172;
     11macro=170;
  • issm/trunk/src/m/enum/PrognosticAnalysisEnum.m

    r4475 r4739  
    99%      macro=PrognosticAnalysisEnum()
    1010
    11 macro=55;
     11macro=53;
  • issm/trunk/src/m/enum/PrognosticSolutionEnum.m

    r4475 r4739  
    99%      macro=PrognosticSolutionEnum()
    1010
    11 macro=54;
     11macro=52;
  • issm/trunk/src/m/enum/QmuAnalysisEnum.m

    r4529 r4739  
    99%      macro=QmuAnalysisEnum()
    1010
    11 macro=277;
     11macro=276;
  • issm/trunk/src/m/enum/QmuErrNameEnum.m

    r4529 r4739  
    99%      macro=QmuErrNameEnum()
    1010
    11 macro=278;
     11macro=277;
  • issm/trunk/src/m/enum/QmuInNameEnum.m

    r4529 r4739  
    99%      macro=QmuInNameEnum()
    1010
    11 macro=279;
     11macro=278;
  • issm/trunk/src/m/enum/QmuMassFluxNumSegmentsEnum.m

    r4529 r4739  
    99%      macro=QmuMassFluxNumSegmentsEnum()
    1010
    11 macro=280;
     11macro=279;
  • issm/trunk/src/m/enum/QmuMassFluxSegmentsEnum.m

    r4529 r4739  
    99%      macro=QmuMassFluxSegmentsEnum()
    1010
    11 macro=281;
     11macro=280;
  • issm/trunk/src/m/enum/QmuNPartEnum.m

    r4529 r4739  
    99%      macro=QmuNPartEnum()
    1010
    11 macro=282;
     11macro=281;
  • issm/trunk/src/m/enum/QmuOutNameEnum.m

    r4529 r4739  
    99%      macro=QmuOutNameEnum()
    1010
    11 macro=283;
     11macro=282;
  • issm/trunk/src/m/enum/QmuPartEnum.m

    r4529 r4739  
    99%      macro=QmuPartEnum()
    1010
    11 macro=284;
     11macro=283;
  • issm/trunk/src/m/enum/QmuPressureEnum.m

    r4475 r4739  
    99%      macro=QmuPressureEnum()
    1010
    11 macro=173;
     11macro=171;
  • issm/trunk/src/m/enum/QmuVxEnum.m

    r4521 r4739  
    99%      macro=QmuVxEnum()
    1010
    11 macro=196;
     11macro=194;
  • issm/trunk/src/m/enum/QmuVyEnum.m

    r4521 r4739  
    99%      macro=QmuVyEnum()
    1010
    11 macro=201;
     11macro=199;
  • issm/trunk/src/m/enum/QmuVzEnum.m

    r4521 r4739  
    99%      macro=QmuVzEnum()
    1010
    11 macro=206;
     11macro=204;
  • issm/trunk/src/m/enum/QuadIceFrontEnum.m

    r4475 r4739  
    99%      macro=QuadIceFrontEnum()
    1010
    11 macro=86;
     11macro=84;
  • issm/trunk/src/m/enum/QuadRiftFrontEnum.m

    r4475 r4739  
    99%      macro=QuadRiftFrontEnum()
    1010
    11 macro=87;
     11macro=85;
  • issm/trunk/src/m/enum/RelativeEnum.m

    r4521 r4739  
    99%      macro=RelativeEnum()
    1010
    11 macro=227;
     11macro=226;
  • issm/trunk/src/m/enum/ResetPenaltiesEnum.m

    r4521 r4739  
    99%      macro=ResetPenaltiesEnum()
    1010
    11 macro=174;
     11macro=172;
  • issm/trunk/src/m/enum/ResidualEnum.m

    r4521 r4739  
    99%      macro=ResidualEnum()
    1010
    11 macro=228;
     11macro=227;
  • issm/trunk/src/m/enum/ResponseDescriptorsEnum.m

    r4529 r4739  
    99%      macro=ResponseDescriptorsEnum()
    1010
    11 macro=285;
     11macro=284;
  • issm/trunk/src/m/enum/RgbEnum.m

    r4475 r4739  
    99%      macro=RgbEnum()
    1010
    11 macro=114;
     11macro=112;
  • issm/trunk/src/m/enum/RheologyBEnum.m

    r4521 r4739  
    99%      macro=RheologyBEnum()
    1010
    11 macro=175;
     11macro=173;
  • issm/trunk/src/m/enum/RheologyNEnum.m

    r4521 r4739  
    99%      macro=RheologyNEnum()
    1010
    11 macro=176;
     11macro=174;
  • issm/trunk/src/m/enum/RhoIceEnum.m

    r4521 r4739  
    99%      macro=RhoIceEnum()
    1010
    11 macro=230;
     11macro=229;
  • issm/trunk/src/m/enum/RhoWaterEnum.m

    r4521 r4739  
    99%      macro=RhoWaterEnum()
    1010
    11 macro=231;
     11macro=230;
  • issm/trunk/src/m/enum/RiftfrontEnum.m

    r4475 r4739  
    99%      macro=RiftfrontEnum()
    1010
    11 macro=88;
     11macro=86;
  • issm/trunk/src/m/enum/SegmentIcefrontEnum.m

    r4475 r4739  
    99%      macro=SegmentIcefrontEnum()
    1010
    11 macro=89;
     11macro=87;
  • issm/trunk/src/m/enum/SegmentOnIceShelfEnum.m

    r4521 r4739  
    99%      macro=SegmentOnIceShelfEnum()
    1010
    11 macro=177;
     11macro=175;
  • issm/trunk/src/m/enum/SegmentRiftfrontEnum.m

    r4475 r4739  
    99%      macro=SegmentRiftfrontEnum()
    1010
    11 macro=90;
     11macro=88;
  • issm/trunk/src/m/enum/SingEnum.m

    r4475 r4739  
    99%      macro=SingEnum()
    1010
    11 macro=77;
     11macro=75;
  • issm/trunk/src/m/enum/SingVertexElementResultEnum.m

    r4521 r4739  
    99%      macro=SingVertexElementResultEnum()
    1010
    11 macro=213;
     11macro=212;
  • issm/trunk/src/m/enum/SingVertexInputEnum.m

    r4475 r4739  
    99%      macro=SingVertexInputEnum()
    1010
    11 macro=100;
     11macro=98;
  • issm/trunk/src/m/enum/SolverStringEnum.m

    r4529 r4739  
    99%      macro=SolverStringEnum()
    1010
    11 macro=286;
     11macro=285;
  • issm/trunk/src/m/enum/SparsityEnum.m

    r4529 r4739  
    99%      macro=SparsityEnum()
    1010
    11 macro=287;
     11macro=286;
  • issm/trunk/src/m/enum/SpcEnum.m

    r4475 r4739  
    99%      macro=SpcEnum()
    1010
    11 macro=115;
     11macro=113;
  • issm/trunk/src/m/enum/StabilizeConstraintsEnum.m

    r4521 r4739  
    99%      macro=StabilizeConstraintsEnum()
    1010
    11 macro=178;
     11macro=176;
  • issm/trunk/src/m/enum/SteadystateAnalysisEnum.m

    r4475 r4739  
    99%      macro=SteadystateAnalysisEnum()
    1010
    11 macro=61;
     11macro=59;
  • issm/trunk/src/m/enum/SteadystateSolutionEnum.m

    r4475 r4739  
    99%      macro=SteadystateSolutionEnum()
    1010
    11 macro=60;
     11macro=58;
  • issm/trunk/src/m/enum/StokesFormulationEnum.m

    r4475 r4739  
    99%      macro=StokesFormulationEnum()
    1010
    11 macro=68;
     11macro=66;
  • issm/trunk/src/m/enum/StokesReconditioningEnum.m

    r4521 r4739  
    99%      macro=StokesReconditioningEnum()
    1010
    11 macro=179;
     11macro=177;
  • issm/trunk/src/m/enum/StringArrayParamEnum.m

    r4475 r4739  
    99%      macro=StringArrayParamEnum()
    1010
    11 macro=110;
     11macro=108;
  • issm/trunk/src/m/enum/StringAsEnum.m

    r4529 r4739  
    6161elseif (strcmpi(name,'BalancedvelocitiesSolution')), enum=BalancedvelocitiesSolutionEnum(); return
    6262elseif (strcmpi(name,'BalancedvelocitiesAnalysis')), enum=BalancedvelocitiesAnalysisEnum(); return
    63 elseif (strcmpi(name,'Prognostic2Solution')), enum=Prognostic2SolutionEnum(); return
    64 elseif (strcmpi(name,'Prognostic2Analysis')), enum=Prognostic2AnalysisEnum(); return
    6563elseif (strcmpi(name,'PrognosticSolution')), enum=PrognosticSolutionEnum(); return
    6664elseif (strcmpi(name,'PrognosticAnalysis')), enum=PrognosticAnalysisEnum(); return
     
    219217elseif (strcmpi(name,'P0')), enum=P0Enum(); return
    220218elseif (strcmpi(name,'P1')), enum=P1Enum(); return
     219elseif (strcmpi(name,'P1DG')), enum=P1DGEnum(); return
    221220elseif (strcmpi(name,'DoubleElementResult')), enum=DoubleElementResultEnum(); return
    222221elseif (strcmpi(name,'TriaVertexElementResult')), enum=TriaVertexElementResultEnum(); return
  • issm/trunk/src/m/enum/StringExternalResultEnum.m

    r4521 r4739  
    99%      macro=StringExternalResultEnum()
    1010
    11 macro=221;
     11macro=220;
  • issm/trunk/src/m/enum/StringParamEnum.m

    r4475 r4739  
    99%      macro=StringParamEnum()
    1010
    11 macro=111;
     11macro=109;
  • issm/trunk/src/m/enum/SurfaceAreaEnum.m

    r4521 r4739  
    99%      macro=SurfaceAreaEnum()
    1010
    11 macro=180;
     11macro=178;
  • issm/trunk/src/m/enum/SurfaceEnum.m

    r4521 r4739  
    99%      macro=SurfaceEnum()
    1010
    11 macro=181;
     11macro=179;
  • issm/trunk/src/m/enum/SurfaceSlopeXEnum.m

    r4521 r4739  
    99%      macro=SurfaceSlopeXEnum()
    1010
    11 macro=182;
     11macro=180;
  • issm/trunk/src/m/enum/SurfaceSlopeYEnum.m

    r4521 r4739  
    99%      macro=SurfaceSlopeYEnum()
    1010
    11 macro=183;
     11macro=181;
  • issm/trunk/src/m/enum/TemperatureAverageEnum.m

    r4521 r4739  
    99%      macro=TemperatureAverageEnum()
    1010
    11 macro=186;
     11macro=184;
  • issm/trunk/src/m/enum/TemperatureEnum.m

    r4521 r4739  
    99%      macro=TemperatureEnum()
    1010
    11 macro=184;
     11macro=182;
  • issm/trunk/src/m/enum/TemperatureOldEnum.m

    r4521 r4739  
    99%      macro=TemperatureOldEnum()
    1010
    11 macro=185;
     11macro=183;
  • issm/trunk/src/m/enum/ThermalConductivityEnum.m

    r4529 r4739  
    99%      macro=ThermalConductivityEnum()
    1010
    11 macro=233;
     11macro=232;
  • issm/trunk/src/m/enum/ThermalExchangeVelocityEnum.m

    r4529 r4739  
    99%      macro=ThermalExchangeVelocityEnum()
    1010
    11 macro=235;
     11macro=234;
  • issm/trunk/src/m/enum/ThicknessEnum.m

    r4521 r4739  
    99%      macro=ThicknessEnum()
    1010
    11 macro=187;
     11macro=185;
  • issm/trunk/src/m/enum/TimeEnum.m

    r4521 r4739  
    99%      macro=TimeEnum()
    1010
    11 macro=226;
     11macro=225;
  • issm/trunk/src/m/enum/TolXEnum.m

    r4529 r4739  
    99%      macro=TolXEnum()
    1010
    11 macro=288;
     11macro=287;
  • issm/trunk/src/m/enum/TriaEnum.m

    r4475 r4739  
    99%      macro=TriaEnum()
    1010
    11 macro=78;
     11macro=76;
  • issm/trunk/src/m/enum/TriaVertexElementResultEnum.m

    r4521 r4739  
    99%      macro=TriaVertexElementResultEnum()
    1010
    11 macro=211;
     11macro=210;
  • issm/trunk/src/m/enum/TriaVertexInputEnum.m

    r4475 r4739  
    99%      macro=TriaVertexInputEnum()
    1010
    11 macro=101;
     11macro=99;
  • issm/trunk/src/m/enum/TypeEnum.m

    r4521 r4739  
    99%      macro=TypeEnum()
    1010
    11 macro=188;
     11macro=186;
  • issm/trunk/src/m/enum/VariableDescriptorsEnum.m

    r4529 r4739  
    99%      macro=VariableDescriptorsEnum()
    1010
    11 macro=289;
     11macro=288;
  • issm/trunk/src/m/enum/VelEnum.m

    r4521 r4739  
    99%      macro=VelEnum()
    1010
    11 macro=189;
     11macro=187;
  • issm/trunk/src/m/enum/VelObsEnum.m

    r4521 r4739  
    99%      macro=VelObsEnum()
    1010
    11 macro=190;
     11macro=188;
  • issm/trunk/src/m/enum/VerboseEnum.m

    r4529 r4739  
    99%      macro=VerboseEnum()
    1010
    11 macro=290;
     11macro=289;
  • issm/trunk/src/m/enum/VertexEnum.m

    r4475 r4739  
    99%      macro=VertexEnum()
    1010
    11 macro=80;
     11macro=78;
  • issm/trunk/src/m/enum/ViscosityOvershootEnum.m

    r4521 r4739  
    99%      macro=ViscosityOvershootEnum()
    1010
    11 macro=191;
     11macro=189;
  • issm/trunk/src/m/enum/VxAverageEnum.m

    r4521 r4739  
    99%      macro=VxAverageEnum()
    1010
    11 macro=192;
     11macro=190;
  • issm/trunk/src/m/enum/VxEnum.m

    r4521 r4739  
    99%      macro=VxEnum()
    1010
    11 macro=193;
     11macro=191;
  • issm/trunk/src/m/enum/VxObsEnum.m

    r4521 r4739  
    99%      macro=VxObsEnum()
    1010
    11 macro=194;
     11macro=192;
  • issm/trunk/src/m/enum/VxOldEnum.m

    r4521 r4739  
    99%      macro=VxOldEnum()
    1010
    11 macro=195;
     11macro=193;
  • issm/trunk/src/m/enum/VyAverageEnum.m

    r4521 r4739  
    99%      macro=VyAverageEnum()
    1010
    11 macro=197;
     11macro=195;
  • issm/trunk/src/m/enum/VyEnum.m

    r4521 r4739  
    99%      macro=VyEnum()
    1010
    11 macro=198;
     11macro=196;
  • issm/trunk/src/m/enum/VyObsEnum.m

    r4521 r4739  
    99%      macro=VyObsEnum()
    1010
    11 macro=199;
     11macro=197;
  • issm/trunk/src/m/enum/VyOldEnum.m

    r4521 r4739  
    99%      macro=VyOldEnum()
    1010
    11 macro=200;
     11macro=198;
  • issm/trunk/src/m/enum/VzAverageEnum.m

    r4521 r4739  
    99%      macro=VzAverageEnum()
    1010
    11 macro=202;
     11macro=200;
  • issm/trunk/src/m/enum/VzEnum.m

    r4521 r4739  
    99%      macro=VzEnum()
    1010
    11 macro=203;
     11macro=201;
  • issm/trunk/src/m/enum/VzObsEnum.m

    r4521 r4739  
    99%      macro=VzObsEnum()
    1010
    11 macro=204;
     11macro=202;
  • issm/trunk/src/m/enum/VzOldEnum.m

    r4521 r4739  
    99%      macro=VzOldEnum()
    1010
    11 macro=205;
     11macro=203;
  • issm/trunk/src/m/enum/WaitOnLockEnum.m

    r4529 r4739  
    99%      macro=WaitOnLockEnum()
    1010
    11 macro=291;
     11macro=290;
  • issm/trunk/src/m/enum/WaterEnum.m

    r4475 r4739  
    99%      macro=WaterEnum()
    1010
    11 macro=119;
     11macro=117;
  • issm/trunk/src/m/enum/WeightsEnum.m

    r4521 r4739  
    99%      macro=WeightsEnum()
    1010
    11 macro=207;
     11macro=205;
  • issm/trunk/src/m/enum/YtsEnum.m

    r4529 r4739  
    99%      macro=YtsEnum()
    1010
    11 macro=292;
     11macro=291;
  • issm/trunk/src/m/solutions/SolutionConfiguration.m

    r4505 r4739  
    2626                analyses=[PrognosticAnalysisEnum];
    2727                solutioncore='prognostic_core';
    28 
    29         case Prognostic2SolutionEnum,
    30                 numanalyses=1;
    31                 analyses=[Prognostic2AnalysisEnum];
    32                 solutioncore='prognostic2_core';
    3328
    3429        case BalancedthicknessSolutionEnum,
  • issm/trunk/src/m/utils/Nightly/ArchiveName.m

    r4022 r4739  
    1919                        output=[output '_' ];
    2020                        for i=1:length(string),
    21                                 %Keep only capital letters + 3 letters (no numbers! as Prognostic2)
     21                                %Keep only capital letters + 3 letters (no numbers!)
    2222                                if strcmp(string(i),upper(string(i)))
    2323                                        output=[output string(i:min(i+3,length(string)))];
  • issm/trunk/src/m/utils/Nightly/testsgetfields.m

    r4716 r4739  
    6565elseif (analysis_type==PrognosticSolutionEnum),
    6666        field_names={'Thickness'};
    67         thickness=zeros(md.numberofgrids,1);
    68         thickness(md.results.PrognosticSolution.Thickness.index)=md.results.PrognosticSolution.Thickness.value;
    69         field_values={thickness};
    70 
    71 elseif (analysis_type==Prognostic2SolutionEnum),
    72         field_names={'Thickness'};
    7367        thickness   =zeros(md.numberofgrids,1);
    74         connectivity=sparse(md.results.Prognostic2Solution.Thickness.index(:),1,1);
    75         thickness   =sparse(md.results.Prognostic2Solution.Thickness.index(:),1,md.results.Prognostic2Solution.Thickness.value(:));
     68        connectivity=sparse(md.results.PrognosticSolution.Thickness.index(:),1,1);
     69        thickness   =sparse(md.results.PrognosticSolution.Thickness.index(:),1,md.results.PrognosticSolution.Thickness.value(:));
    7670        thickness=full(thickness./connectivity);
    7771        field_values={thickness};
  • issm/trunk/src/m/utils/Nightly/testsgettolerance.m

    r4644 r4739  
    3232
    3333        elseif (analysis_type==PrognosticSolutionEnum),
    34                 tolerance=10^-14;
    35 
    36         elseif (analysis_type==Prognostic2SolutionEnum),
    3734                tolerance=10^-14;
    3835
Note: See TracChangeset for help on using the changeset viewer.