Changeset 3567 for issm/trunk


Ignore:
Timestamp:
04/19/10 09:49:52 (15 years ago)
Author:
Mathieu Morlighem
Message:

all enums are now C Enums, no need to provide a number: the numbering is done automatically

Location:
issm/trunk/src
Files:
178 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/DataSet/DataSet.cpp

    r3556 r3567  
    188188                memcpy(&enum_type,marshalled_dataset,sizeof(int)); marshalled_dataset+=sizeof(int);
    189189
    190                 if(enum_type==NodeEnum()){
     190                if(enum_type==NodeEnum){
    191191                        Node* node=NULL;
    192192                        node=new Node();
     
    194194                        dataset->AddObject(node);
    195195                }
    196                 else if(enum_type==VertexEnum()){
     196                else if(enum_type==VertexEnum){
    197197                        Vertex* vertex=NULL;
    198198                        vertex=new Vertex();
     
    200200                        dataset->AddObject(vertex);
    201201                }
    202                 else if(enum_type==ParamEnum()){
     202                else if(enum_type==ParamEnum){
    203203                        Param* param=NULL;
    204204                        param=new Param();
     
    206206                        dataset->AddObject(param);
    207207                }
    208                 else if(enum_type==NumparEnum()){
     208                else if(enum_type==NumparEnum){
    209209                        Numpar* numpar=NULL;
    210210                        numpar=new Numpar();
     
    212212                        dataset->AddObject(numpar);
    213213                }
    214                 else if(enum_type==TriaEnum()){
     214                else if(enum_type==TriaEnum){
    215215                        Tria* tria=NULL;
    216216                        tria=new Tria();
     
    218218                        dataset->AddObject(tria);
    219219                }
    220                 else if(enum_type==SingEnum()){
     220                else if(enum_type==SingEnum){
    221221                        Sing* sing=NULL;
    222222                        sing=new Sing();
     
    224224                        dataset->AddObject(sing);
    225225                }
    226                 else if(enum_type==BeamEnum()){
     226                else if(enum_type==BeamEnum){
    227227                        Beam* beam=NULL;
    228228                        beam=new Beam();
     
    230230                        dataset->AddObject(beam);
    231231                }
    232                 else if(enum_type==PentaEnum()){
     232                else if(enum_type==PentaEnum){
    233233                        Penta* penta=NULL;
    234234                        penta=new Penta();
     
    236236                        dataset->AddObject(penta);
    237237                }
    238                 else if(enum_type==MaticeEnum()){
     238                else if(enum_type==MaticeEnum){
    239239                        Matice* matice=NULL;
    240240                        matice=new Matice();
     
    242242                        dataset->AddObject(matice);
    243243                }
    244                 else if(enum_type==MatparEnum()){
     244                else if(enum_type==MatparEnum){
    245245                        Matpar* matpar=NULL;
    246246                        matpar=new Matpar();
     
    248248                        dataset->AddObject(matpar);
    249249                }
    250                 else if(enum_type==SpcEnum()){
     250                else if(enum_type==SpcEnum){
    251251                        Spc* spc=NULL;
    252252                        spc=new Spc();
     
    254254                        dataset->AddObject(spc);
    255255                }
    256                 else if(enum_type==PengridEnum()){
     256                else if(enum_type==PengridEnum){
    257257                        Pengrid* pengrid=NULL;
    258258                        pengrid=new Pengrid();
     
    260260                        dataset->AddObject(pengrid);
    261261                }
    262                 else if(enum_type==PenpairEnum()){
     262                else if(enum_type==PenpairEnum){
    263263                        Penpair* penpair=NULL;
    264264                        penpair=new Penpair();
     
    266266                        dataset->AddObject(penpair);
    267267                }
    268                 else if(enum_type==IcefrontEnum()){
     268                else if(enum_type==IcefrontEnum){
    269269                        Icefront* icefront=NULL;
    270270                        icefront=new Icefront();
     
    272272                        dataset->AddObject(icefront);
    273273                }
    274                 else if(enum_type==NumericalfluxEnum()){
     274                else if(enum_type==NumericalfluxEnum){
    275275                        Numericalflux* numericalflux=NULL;
    276276                        numericalflux=new Numericalflux();
     
    278278                        dataset->AddObject(numericalflux);
    279279                }
    280                 else if(enum_type==RgbEnum()){
     280                else if(enum_type==RgbEnum){
    281281                        Rgb* rgb=NULL;
    282282                        rgb=new Rgb();
     
    284284                        dataset->AddObject(rgb);
    285285                }
    286                 else if(enum_type==RiftfrontEnum()){
     286                else if(enum_type==RiftfrontEnum){
    287287                        Riftfront* riftfront=NULL;
    288288                        riftfront=new Riftfront();
     
    385385
    386386                /*Find param type objects: */
    387                 if((*object)->Enum()==ParamEnum()){
     387                if((*object)->Enum()==ParamEnum){
    388388
    389389                        /*Ok, this object is a parameter, recover it and ask which name it has: */
     
    416416
    417417                /*Find param type objects: */
    418                 if((*object)->Enum()==ParamEnum()){
     418                if((*object)->Enum()==ParamEnum){
    419419
    420420                        /*Ok, this object is a parameter, recover it and ask which name it has: */
     
    446446
    447447                /*Find param type objects: */
    448                 if((*object)->Enum()==ParamEnum()){
     448                if((*object)->Enum()==ParamEnum){
    449449
    450450                        /*Ok, this object is a parameter, recover it and ask which name it has: */
     
    477477
    478478                /*Find param type objects: */
    479                 if((*object)->Enum()==ParamEnum()){
     479                if((*object)->Enum()==ParamEnum){
    480480
    481481                        /*Ok, this object is a parameter, recover it and ask which name it has: */
     
    509509
    510510                /*Find param type objects: */
    511                 if((*object)->Enum()==ParamEnum()){
     511                if((*object)->Enum()==ParamEnum){
    512512
    513513                        /*Ok, this object is a parameter, recover it and ask which name it has: */
     
    542542
    543543                /*Find param type objects: */
    544                 if((*object)->Enum()==ParamEnum()){
     544                if((*object)->Enum()==ParamEnum){
    545545
    546546                        /*Ok, this object is a parameter, recover it and ask which name it has: */
     
    573573
    574574                /*Find param type objects: */
    575                 if((*object)->Enum()==ParamEnum()){
     575                if((*object)->Enum()==ParamEnum){
    576576
    577577                        /*Ok, this object is a parameter, recover it and ask which name it has: */
     
    602602
    603603                /*Find param type objects: */
    604                 if((*object)->Enum()==ParamEnum()){
     604                if((*object)->Enum()==ParamEnum){
    605605
    606606                        /*Ok, this object is a parameter, recover it and ask which name it has: */
     
    630630
    631631                /*Find param type objects: */
    632                 if((*object)->Enum()==ResultEnum()){
     632                if((*object)->Enum()==ResultEnum){
    633633
    634634                        /*Ok, this object is a result,recover it and ask which name it has: */
     
    660660
    661661                /*Find param type objects: */
    662                 if((*object)->Enum()==ResultEnum()){
     662                if((*object)->Enum()==ResultEnum){
    663663
    664664                        /*Ok, this object is a result, recover it and ask which name it has: */
     
    678678}
    679679/*}}}*/
    680 /*FUNCTION DataSet::GetEnum(){{{1*/
     680/*FUNCTION DataSet::GetEnum{{{1*/
    681681int  DataSet::GetEnum(){
    682682        return enum_type;
     
    861861                }
    862862
    863                 if((*object)->Enum()==NodeEnum()){
     863                if((*object)->Enum()==NodeEnum){
    864864                        node=(Node*)(*object);
    865865                        node->Configure(nodes,vertices);
    866866                }
    867                 if((*object)->Enum()==NumparEnum()){
     867                if((*object)->Enum()==NumparEnum){
    868868                        numpar=(Numpar*)(*object);
    869869                        numpar->Configure(parameters);
     
    10801080        for ( object=objects.begin() ; object < objects.end(); object++ ){
    10811081
    1082                 if((*object)->Enum()==NodeEnum()){
     1082                if((*object)->Enum()==NodeEnum){
    10831083                        node=(Node*)(*object);
    10841084                        node->FieldAverageOntoVertices(fieldsum,connectivity,field);
     
    10971097        for ( object=objects.begin() ; object < objects.end(); object++ ){
    10981098
    1099                 if((*object)->Enum()==NodeEnum()){
     1099                if((*object)->Enum()==NodeEnum){
    11001100
    11011101                        node=(Node*)(*object);
     
    11151115        for ( object=objects.begin() ; object < objects.end(); object++ ){
    11161116
    1117                 if((*object)->Enum()==PentaEnum()){
     1117                if((*object)->Enum()==PentaEnum){
    11181118
    11191119                        penta=(Penta*)(*object);
     
    11211121
    11221122                }
    1123                 if((*object)->Enum()==NodeEnum()){
     1123                if((*object)->Enum()==NodeEnum){
    11241124
    11251125                        node=(Node*)(*object);
     
    11931193
    11941194                /*Check this is a single point constraint (spc): */
    1195                 if((*object)->Enum()==NodeEnum()){
     1195                if((*object)->Enum()==NodeEnum){
    11961196
    11971197                        node=(Node*)(*object);
     
    12501250        for ( object=objects.begin() ; object < objects.end(); object++ ){
    12511251
    1252                 if((*object)->Enum()==PengridEnum()){
     1252                if((*object)->Enum()==PengridEnum){
    12531253                        found=1;
    12541254                        break;
     
    12821282        for ( object=objects.begin() ; object < objects.end(); object++ ){
    12831283               
    1284                 if((*object)->Enum()==PengridEnum()){
     1284                if((*object)->Enum()==PengridEnum){
    12851285
    12861286                        pengrid=(Pengrid*)(*object);
     
    13431343
    13441344                /*Check this is a node: */
    1345                 if((*object)->Enum()==NodeEnum()){
     1345                if((*object)->Enum()==NodeEnum){
    13461346
    13471347                        node=(Node*)(*object);
     
    13721372
    13731373                /*Check this is a single point constraint (spc): */
    1374                 if((*object)->Enum()==RgbEnum()){ //we assume uniqueness of all Rgbs, no error checking here.
     1374                if((*object)->Enum()==RgbEnum){ //we assume uniqueness of all Rgbs, no error checking here.
    13751375
    13761376                        count++;
     
    13901390        for ( object=objects.begin() ; object < objects.end(); object++ ){
    13911391
    1392                 if((*object)->Enum()==RiftfrontEnum()){
     1392                if((*object)->Enum()==RiftfrontEnum){
    13931393
    13941394                        riftfront=(Riftfront*)(*object);
     
    14481448        for ( object=objects.begin() ; object < objects.end(); object++ ){
    14491449
    1450                 if((*object)->Enum()==PenpairEnum()){
     1450                if((*object)->Enum()==PenpairEnum){
    14511451
    14521452                        penpair=(Penpair*)(*object);
     
    14841484
    14851485                /*Check this is a single point constraint (spc): */
    1486                 if((*object)->Enum()==RgbEnum()){ //we assume uniqueness of all Rgbs, no error checking here.
     1486                if((*object)->Enum()==RgbEnum){ //we assume uniqueness of all Rgbs, no error checking here.
    14871487
    14881488                        /*we found an rgb, increment counter, so that row index for Rmg is up to date: */
     
    15451545
    15461546                /*Check this is a single point constraint (spc): */
    1547                 if((*object)->Enum()==SpcEnum()){
     1547                if((*object)->Enum()==SpcEnum){
    15481548
    15491549                        spc=(Spc*)(*object);
     
    16191619                        material->UpdateFromInputs(inputs);
    16201620                }
    1621                 if((*object)->Enum()==NodeEnum()){
     1621                if((*object)->Enum()==NodeEnum){
    16221622                        node=(Node*)(*object);
    16231623                        node->UpdateFromInputs(inputs);
     
    16351635        for ( object=objects.begin() ; object < objects.end(); object++ ){
    16361636
    1637                 if((*object)->Enum()==VertexEnum()){
     1637                if((*object)->Enum()==VertexEnum){
    16381638
    16391639                        vertex=(Vertex*)(*object);
  • issm/trunk/src/c/EnumDefinitions/EnumDefinitions.cpp

    r3469 r3567  
    66#include "./EnumDefinitions.h"
    77
    8 /*Datasets: */
    9 int DatasetsEnum(void){                     return          100; }
    10 int ElementsEnum(void){                     return          101; }
    11 int NodesEnum(void){                        return          102; }
    12 int ConstraintsEnum(void){                  return          103; }
    13 int LoadsEnum(void){                        return          104; }
    14 int MaterialsEnum(void){                    return          105; }
    15 int ParametersEnum(void){                   return          106; }
    16 int ResultsEnum(void){                      return          107; }
    17 int VerticesEnum(void){                     return          108; }
    18 
    19 /*ANALYSIS TYPES: */
    20 int AnalysisEnum(void){                     return          200; }
    21 //diagnostic
    22 int DiagnosticAnalysisEnum(void){           return          210; }
    23 int HorizAnalysisEnum(void){                return          211; }
    24 int StokesAnalysisEnum(void){               return          212; }
    25 int HutterAnalysisEnum(void){               return          213; }
    26 int VertAnalysisEnum(void){                 return          214; }
    27 //control
    28 int ControlAnalysisEnum(void){              return          220; }
    29 int AdjointAnalysisEnum(void){              return          221; }
    30 int InverseAnalysisEnum(void){              return          222; }
    31 int GradientAnalysisEnum(void){             return          223; }
    32 //thermal
    33 int ThermalAnalysisEnum(void){              return          230; }
    34 //transient
    35 int TransientAnalysisEnum(void){            return          231; }
    36 int SteadyAnalysisEnum(void){               return          232; } //NEEDED by nightly runs
    37 //slope
    38 int SlopecomputeAnalysisEnum(void){         return          240; }
    39 int SurfaceXAnalysisEnum(void){             return          241; }
    40 int SurfaceYAnalysisEnum(void){             return          242; }
    41 int BedXAnalysisEnum(void){                 return          243; }
    42 int BedYAnalysisEnum(void){                 return          244; }
    43 //prognostic
    44 int PrognosticAnalysisEnum(void){           return          250; }
    45 int Prognostic2AnalysisEnum(void){          return          251; }
    46 int BalancedthicknessAnalysisEnum(void){    return          252; }
    47 int BalancedvelocitiesAnalysisEnum(void){   return          253; }
    48 //melting
    49 int MeltingAnalysisEnum(void){              return          260; }
    50 //mesh2grid
    51 int Mesh2gridAnalysisEnum(void){            return          270; }
    52 //parameters
    53 int ParametersAnalysisEnum(void){           return          280; }
    54 //steadystate
    55 int SteadystateAnalysisEnum(void){          return          281; }
    56 //none
    57 int NoneAnalysisEnum(void){                 return          290; }
    58 
    59 
    60 /*Formulations: */
    61 int FormulationEnum(void){                  return          300; }
    62 int NoneFormulationEnum(void){              return          301; }
    63 int HutterFormulationEnum(void){            return          302; }
    64 int MacAyealFormulationEnum(void){          return          303; }
    65 int PattynFormulationEnum(void){            return          304; }
    66 int StokesFormulationEnum(void){            return          305; }
    67 
    68 /*OBJECTS: */
    69 int ObjectEnum(void){                       return          400; }
    70 /*Elements: */
    71 int HookEnum(void){                         return          409; }
    72 int ElementEnum(void){                      return          410; }
    73 int TriaEnum(void){                         return          411; }
    74 int ElementPropertiesEnum(void){            return          412; }
    75 int NodePropertiesEnum(void){               return          413; }
    76 int PentaEnum(void){                        return          414; }
    77 int SingEnum(void){                         return          415; }
    78 int BeamEnum(void){                         return          416; }
    79 int DofIndexingEnum(void){                  return          417; }
    80 
    81 /*Grids: */
    82 int VertexEnum(void){                       return          420; }
    83 int NodeEnum(void){                         return          421; }
    84 /*Loads: */
    85 int LoadEnum(void){                         return          430; }
    86 int IcefrontEnum(void){                     return          431; }
    87 int RiftfrontEnum(void){                    return          432; }
    88 int PenpairEnum(void){                      return          433; }
    89 int PengridEnum(void){                      return          434; }
    90 int NumericalfluxEnum(void){                return          435; }
    91 /*Materials: */
    92 int MaterialEnum(void){                     return          440; }
    93 int MaticeEnum(void){                       return          441; }
    94 int MatparEnum(void){                       return          442; }
    95 int NumparEnum(void){                       return          443; }
    96 /*Inputs: */
    97 int InputEnum(void){                        return          450; }
    98 /*Params: */
    99 int ParamEnum(void){                        return          460; }
    100 /*Results: */
    101 int ResultEnum(void){                       return          470; }
    102 /*Rgb: */
    103 int RgbEnum(void){                          return          480; }
    104 /*Spc: */
    105 int SpcEnum(void){                          return          490; }
    106 /*DofVec: */
    107 int DofVecEnum(void){                       return          495; }
    108 
    109 
    110 /*GEOGRAPHY:*/
    111 int GeographyEnum(void){                    return          500; }
    112 int IceSheetEnum(void){                     return          502; }
    113 int IceShelfEnum(void){                     return          502; }
    114 
    115 /*FILL:*/
    116 int WaterEnum(void){                        return          601; }
    117 int IceEnum(void){                          return          602; }
    118 int AirEnum(void){                          return          603; }
    119 int MelangeEnum(void){                      return          604; }
    120 
    1218/*functions on enums: */
    1229int EnumIsElement(int en){
    123         if ((en>=410) && (en<=419))return 1;
     10        if ((en>=ElementEnum) && (en<=SingEnum))return 1;
    12411        else return 0;
    12512}
    12613
    12714int EnumIsLoad(int en){
    128         if ((en>=430) && (en<=439))return 1;
     15        if ((en>=LoadEnum) && (en<=NumericalfluxEnum))return 1;
    12916        else return 0;
    13017}
    13118int EnumIsMaterial(int en){
    132         if ((en>=440) && (en<=449))return 1;
     19        if ((en>=MaterialEnum) && (en<=NumparEnum))return 1;
    13320        else return 0;
    13421}
  • issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h

    r3417 r3567  
    77#define _ENUM_DEFINITIONS_
    88
    9 /*Datasets: */
    10 int DatasetsEnum(void);
    11 int ElementsEnum(void);
    12 int NodesEnum(void);
    13 int ConstraintsEnum(void);
    14 int LoadsEnum(void);
    15 int MaterialsEnum(void);
    16 int ParametersEnum(void);
    17 int ResultsEnum(void);
    18 int VerticesEnum(void);
     9enum definitions{
    1910
    20 /*ANALYSIS TYPES: */
    21 int AnalysisEnum(void);
    22 //diagnostic
    23 int DiagnosticAnalysisEnum(void);
    24 int HorizAnalysisEnum(void);
    25 int StokesAnalysisEnum(void);
    26 int HutterAnalysisEnum(void);
    27 int VertAnalysisEnum(void);
    28 //control
    29 int ControlAnalysisEnum(void);
    30 int AdjointAnalysisEnum(void);
    31 int InverseAnalysisEnum(void);
    32 int GradientAnalysisEnum(void);
    33 //thermal
    34 int ThermalAnalysisEnum(void);
    35 //transient
    36 int TransientAnalysisEnum(void);
    37 int SteadyAnalysisEnum(void);
    38 //slope
    39 int SlopecomputeAnalysisEnum(void);
    40 int SurfaceXAnalysisEnum(void);
    41 int SurfaceYAnalysisEnum(void);
    42 int BedXAnalysisEnum(void);
    43 int BedYAnalysisEnum(void);
    44 //prognostic
    45 int PrognosticAnalysisEnum(void);
    46 int Prognostic2AnalysisEnum(void);
    47 int BalancedthicknessAnalysisEnum(void);
    48 int BalancedvelocitiesAnalysisEnum(void);
    49 //melting
    50 int MeltingAnalysisEnum(void);
    51 //mesh2grid
    52 int Mesh2gridAnalysisEnum(void);
    53 //parameters
    54 int ParametersAnalysisEnum(void);
    55 //steadystate
    56 int SteadystateAnalysisEnum(void);
    57 //none
    58 int NoneAnalysisEnum(void);
     11        /*Datasets {{{1*/
     12        DatasetsEnum,
     13        ElementsEnum,
     14        NodesEnum,
     15        ConstraintsEnum,
     16        LoadsEnum,
     17        MaterialsEnum,
     18        ParametersEnum,
     19        ResultsEnum,
     20        VerticesEnum,
     21        /*}}}*/
     22        /*Analysis types {{{1 */
     23        AnalysisEnum,
     24        //diagnostic
     25        DiagnosticAnalysisEnum,
     26        HorizAnalysisEnum,
     27        StokesAnalysisEnum,
     28        HutterAnalysisEnum,
     29        VertAnalysisEnum,
     30        //control
     31        ControlAnalysisEnum,
     32        AdjointAnalysisEnum,
     33        InverseAnalysisEnum,
     34        GradientAnalysisEnum,
     35        //thermal
     36        ThermalAnalysisEnum,
     37        //transient
     38        TransientAnalysisEnum,
     39        SteadyAnalysisEnum,
     40        //slope
     41        SlopecomputeAnalysisEnum,
     42        SurfaceXAnalysisEnum,
     43        SurfaceYAnalysisEnum,
     44        BedXAnalysisEnum,
     45        BedYAnalysisEnum,
     46        //prognostic
     47        PrognosticAnalysisEnum,
     48        Prognostic2AnalysisEnum,
     49        BalancedthicknessAnalysisEnum,
     50        BalancedvelocitiesAnalysisEnum,
     51        //melting
     52        MeltingAnalysisEnum,
     53        //mesh2grid
     54        Mesh2gridAnalysisEnum,
     55        //parameters
     56        ParametersAnalysisEnum,
     57        //steadystate
     58        SteadystateAnalysisEnum,
     59        //none
     60        NoneAnalysisEnum,
     61        /*}}}*/
     62        /*Formulations {{{1*/
     63        FormulationEnum,
     64        NoneFormulationEnum,
     65        HutterFormulationEnum,
     66        MacAyealFormulationEnum,
     67        PattynFormulationEnum,
     68        StokesFormulationEnum,
     69        /*}}}*/
     70        /*Objects {{{1*/
     71        ObjectEnum,
     72        HookEnum,
     73        /*Elements: */
     74        ElementEnum,
     75        TriaEnum,
     76        ElementPropertiesEnum,
     77        NodePropertiesEnum,
     78        PentaEnum,
     79        BeamEnum,
     80        SingEnum,
     81        DofIndexingEnum,
     82        /*Grids: */
     83        NodeEnum,
     84        VertexEnum,
     85        /*Loads: */
     86        LoadEnum,
     87        IcefrontEnum,
     88        RiftfrontEnum,
     89        PenpairEnum,
     90        PengridEnum,
     91        NumericalfluxEnum,
     92        /*Materials: */
     93        MaterialEnum,
     94        MaticeEnum,
     95        MatparEnum,
     96        NumparEnum,
     97        /*Inputs: */
     98        InputEnum,
     99        /*Params: */
     100        ParamEnum,
     101        /*Results: */
     102        ResultEnum,
     103        /*Rgb: */
     104        RgbEnum,
     105        /*Spc: */
     106        SpcEnum,
     107        /*DofVec: */
     108        DofVecEnum,
     109        /*}}}*/
     110        /*Geography {{{1*/
     111        GeographyEnum,
     112        IceSheetEnum,
     113        IceShelfEnum,
     114        /*}}}*/
     115        /*Fill {{{1*/
     116        WaterEnum,
     117        IceEnum,
     118        AirEnum,
     119        MelangeEnum,
     120        /*}}}*/
    59121
    60 
    61 /*Formulations: */
    62 int FormulationEnum(void);
    63 int NoneFormulationEnum(void);
    64 int HutterFormulationEnum(void);
    65 int MacAyealFormulationEnum(void);
    66 int PattynFormulationEnum(void);
    67 int StokesFormulationEnum(void);
    68 
    69 /*OBJECTS: */
    70 int ObjectEnum(void);
    71 /*Elements: */
    72 int HookEnum(void);
    73 int ElementEnum(void);
    74 int TriaEnum(void);
    75 int ElementPropertiesEnum(void);
    76 int NodePropertiesEnum(void);
    77 int PentaEnum(void);
    78 int SingEnum(void);
    79 int BeamEnum(void);
    80 int DofIndexingEnum(void);
    81 /*Grids: */
    82 int NodeEnum(void);
    83 int VertexEnum(void);
    84 /*Loads: */
    85 int LoadEnum(void);
    86 int IcefrontEnum(void);
    87 int RiftfrontEnum(void);
    88 int PenpairEnum(void);
    89 int PengridEnum(void);
    90 int NumericalfluxEnum(void);
    91 /*Materials: */
    92 int MaterialEnum(void);
    93 int MaticeEnum(void);
    94 int MatparEnum(void);
    95 int NumparEnum(void);
    96 /*Inputs: */
    97 int InputEnum(void);
    98 /*Params: */
    99 int ParamEnum(void);
    100 /*Results: */
    101 int ResultEnum(void);
    102 /*Rgb: */
    103 int RgbEnum(void);
    104 /*Spc: */
    105 int SpcEnum(void);
    106 /*DofVec: */
    107 int DofVecEnum(void);
    108 
    109 /*GEOGRAPHY:*/
    110 int GeographyEnum(void);
    111 int IceSheetEnum(void);
    112 int IceShelfEnum(void);
    113 
    114 /*FILL: */
    115 int WaterEnum(void);
    116 int IceEnum(void);
    117 int AirEnum(void);
    118 int MelangeEnum(void);
     122};
    119123
    120124/*Functions on enums: */
     
    122126int EnumIsLoad(int en);
    123127int EnumIsMaterial(int en);
     128
    124129#endif
  • issm/trunk/src/c/EnumDefinitions/SynchronizeMatlabEnum.sh

    r1901 r3567  
    77rm $ISSM_DIR/src/m/enum/*.m
    88
    9 #Get the list of interesting lines of EnumDefinitions.cpp and introduce line number
    109echo "Synchronizong with Matlab's Enums..."
    11 cat EnumDefinitions.cpp | grep "return" | grep "int" | awk '{ printf "line%so %s\n", NR, $0 }' | sed -e "s/(void){//"  > temp
     10#Get all lines of EnumDefinitions2.h which hold Enum, | remove all comas | add line number in the first column > put everything in file temp
     11cat EnumDefinitions.h | grep "Enum," |sed -e "s/,//g" | awk '{ printf "%s %s\n", NR, $0 }' > temp
     12
     13#Build File AnalysisTypeFromEnum.m
    1214cat <<END > AnalysisTypeFromEnum.m
    1315function string=AnalysisTypeFromEnum(enum)
     
    2325%initialize output
    2426string='not found';
    25 
    2627END
    2728
    2829#get number of lines in temp
    29 NUMBEROFLINES=`wc -l temp | awk '{printf("%s",$1);}'`;
     30NUMBEROFLINES=$(wc -l temp | awk '{printf("%s",$1);}');
    3031
    3132# go through the lines of temp
     
    3334
    3435        #Get name and enum of the line i
    35         NAME=`cat temp |  awk '/line'$i'o/ {printf("%s\n",$3);}';`;
    36         ENUM=`cat temp |  awk '/line'$i'o/ {printf("%s\n",$5);}';`;
    37         echo writing file $NAME \($i/$NUMBEROFLINES\)
     36        NAME=$(cat temp | grep "^[ ]*$i " | awk '{printf("%s",$2);}');
     37        ENUM=$i;
     38        #offset Enum by one (Enum begins with 0 and not 1!)
     39        let ENUM=$ENUM-1
    3840
    3941        #write corresponding Matlab file
    40         cat <<END > $(echo $NAME.m)
     42        echo "writing file "$NAME".m ($i/$NUMBEROFLINES)"
     43        cat <<END > $(echo $NAME".m")
    4144function macro=$(echo $NAME)()
    4245%$(echo `echo $NAME | sed -e "s/[a-z]/\U&/g"`) - Enum of $(echo `echo $NAME | sed -e "s/Enum//"`)
     
    4750%      macro=$NAME()
    4851
    49 macro=$ENUM
     52macro=$ENUM;
    5053END
    5154
     
    6265
    6366        #move file to m/enum
    64         mv $NAME.m $ISSM_DIR/src/m/enum/
     67        mv $NAME".m" $ISSM_DIR"/src/m/enum/"
    6568
    6669done
  • issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateConstraintsBalancedthickness.cpp

    r3452 r3567  
    1717
    1818        /*Create constraints: */
    19         constraints = new DataSet(ConstraintsEnum());
     19        constraints = new DataSet(ConstraintsEnum);
    2020
    2121        /*Fetch data: */
  • issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateElementsNodesAndMaterialsBalancedthickness.cpp

    r3476 r3567  
    2424
    2525        /*First create the elements, nodes and material properties: */
    26         elements  = new DataSet(ElementsEnum());
    27         nodes     = new DataSet(NodesEnum());
    28         vertices  = new DataSet(VerticesEnum());
    29         materials = new DataSet(MaterialsEnum());
     26        elements  = new DataSet(ElementsEnum);
     27        nodes     = new DataSet(NodesEnum);
     28        vertices  = new DataSet(VerticesEnum);
     29        materials = new DataSet(MaterialsEnum);
    3030
    3131        /*Partition elements and vertices and nodes: */
  • issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateLoadsBalancedthickness.cpp

    r3332 r3567  
    1616
    1717        /*Create loads: */
    18         loads   = new DataSet(LoadsEnum());
     18        loads   = new DataSet(LoadsEnum);
    1919       
    2020       
  • issm/trunk/src/c/ModelProcessorx/Balancedvelocities/CreateConstraintsBalancedvelocities.cpp

    r3456 r3567  
    1919
    2020        /*Create constraints: */
    21         constraints = new DataSet(ConstraintsEnum());
     21        constraints = new DataSet(ConstraintsEnum);
    2222
    2323        /*Fetch data: */
  • issm/trunk/src/c/ModelProcessorx/Balancedvelocities/CreateElementsNodesAndMaterialsBalancedvelocities.cpp

    r3434 r3567  
    2424
    2525        /*First create the elements, nodes and material properties: */
    26         elements  = new DataSet(ElementsEnum());
    27         nodes     = new DataSet(NodesEnum());
    28         vertices  = new DataSet(VerticesEnum());
    29         materials = new DataSet(MaterialsEnum());
     26        elements  = new DataSet(ElementsEnum);
     27        nodes     = new DataSet(NodesEnum);
     28        vertices  = new DataSet(VerticesEnum);
     29        materials = new DataSet(MaterialsEnum);
    3030
    3131        /*Partition elements and vertices and nodes: */
  • issm/trunk/src/c/ModelProcessorx/Balancedvelocities/CreateLoadsBalancedvelocities.cpp

    r3332 r3567  
    1616
    1717        /*Create loads: */
    18         loads   = new DataSet(LoadsEnum());
     18        loads   = new DataSet(LoadsEnum);
    1919       
    2020       
  • issm/trunk/src/c/ModelProcessorx/CreateDataSets.cpp

    r3417 r3567  
    2222
    2323        /*This is just a high level driver: */
    24         if (iomodel->analysis_type==DiagnosticAnalysisEnum()){
     24        if (iomodel->analysis_type==DiagnosticAnalysisEnum){
    2525
    26                 if (iomodel->sub_analysis_type==HorizAnalysisEnum()){
     26                if (iomodel->sub_analysis_type==HorizAnalysisEnum){
    2727
    2828                        CreateElementsNodesAndMaterialsDiagnosticHoriz(pelements,pnodes, pvertices, pmaterials, iomodel,iomodel_handle);
     
    3232
    3333                }
    34                 else if (iomodel->sub_analysis_type==VertAnalysisEnum()){
     34                else if (iomodel->sub_analysis_type==VertAnalysisEnum){
    3535               
    3636                        CreateElementsNodesAndMaterialsDiagnosticVert(pelements,pnodes, pvertices, pmaterials, iomodel,iomodel_handle);
     
    3939                       
    4040                }
    41                 else if (iomodel->sub_analysis_type==StokesAnalysisEnum()){
     41                else if (iomodel->sub_analysis_type==StokesAnalysisEnum){
    4242
    4343                        CreateElementsNodesAndMaterialsDiagnosticStokes(pelements,pnodes, pvertices, pmaterials, iomodel,iomodel_handle);
     
    4646                       
    4747                }
    48                 else if (iomodel->sub_analysis_type==HutterAnalysisEnum()){
     48                else if (iomodel->sub_analysis_type==HutterAnalysisEnum){
    4949
    5050                        CreateElementsNodesAndMaterialsDiagnosticHutter(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle);
     
    5454                }
    5555        }
    56         else if (iomodel->analysis_type==SlopecomputeAnalysisEnum()){
     56        else if (iomodel->analysis_type==SlopecomputeAnalysisEnum){
    5757
    5858                CreateElementsNodesAndMaterialsSlopeCompute(pelements,pnodes, pvertices,pmaterials, iomodel,iomodel_handle);
     
    6161       
    6262        }
    63         else if (iomodel->analysis_type==ThermalAnalysisEnum()){
     63        else if (iomodel->analysis_type==ThermalAnalysisEnum){
    6464
    6565                CreateElementsNodesAndMaterialsThermal(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle);
     
    6969                                       
    7070        }
    71         else if (iomodel->analysis_type==MeltingAnalysisEnum()){
     71        else if (iomodel->analysis_type==MeltingAnalysisEnum){
    7272                       
    7373                CreateElementsNodesAndMaterialsMelting(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle);
     
    7676                CreateParametersMelting(pparameters,iomodel,iomodel_handle);
    7777        }
    78         else if (iomodel->analysis_type==PrognosticAnalysisEnum()){
     78        else if (iomodel->analysis_type==PrognosticAnalysisEnum){
    7979
    8080                CreateElementsNodesAndMaterialsPrognostic(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle);
     
    8484                                       
    8585        }
    86         else if (iomodel->analysis_type==Prognostic2AnalysisEnum()){
     86        else if (iomodel->analysis_type==Prognostic2AnalysisEnum){
    8787
    8888                CreateElementsNodesAndMaterialsPrognostic2(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle);
     
    9292
    9393        }
    94         else if (iomodel->analysis_type==BalancedthicknessAnalysisEnum()){
     94        else if (iomodel->analysis_type==BalancedthicknessAnalysisEnum){
    9595
    9696                CreateElementsNodesAndMaterialsBalancedthickness(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle);
     
    100100
    101101        }
    102         else if (iomodel->analysis_type==BalancedvelocitiesAnalysisEnum()){
     102        else if (iomodel->analysis_type==BalancedvelocitiesAnalysisEnum){
    103103
    104104                CreateElementsNodesAndMaterialsBalancedvelocities(pelements,pnodes,pvertices,pmaterials, iomodel,iomodel_handle);
  • issm/trunk/src/c/ModelProcessorx/CreateParameters.cpp

    r3538 r3567  
    3030
    3131        /*Initialize dataset: */
    32         parameters   = new DataSet(ParametersEnum());
     32        parameters   = new DataSet(ParametersEnum);
    3333
    3434        //solution parameters: always 1'st, to speed-up its lookup by elements..
     
    245245        count++;
    246246        param= new Param(count,"numberofnodes",DOUBLE);
    247         if (iomodel->analysis_type==Prognostic2AnalysisEnum()) param->SetDouble(3*iomodel->numberofelements);
     247        if (iomodel->analysis_type==Prognostic2AnalysisEnum) param->SetDouble(3*iomodel->numberofelements);
    248248        else param->SetDouble(iomodel->numberofvertices);
    249249        parameters->AddObject(param);
  • issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp

    r3507 r3567  
    2323
    2424        /*Create constraints: */
    25         constraints = new DataSet(ConstraintsEnum());
     25        constraints = new DataSet(ConstraintsEnum);
    2626
    2727        /*Now, is the flag macayaealpattyn on? otherwise, do nothing: */
  • issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp

    r3476 r3567  
    2424
    2525        /*First create the elements, nodes and material properties: */
    26         elements  = new DataSet(ElementsEnum());
    27         nodes     = new DataSet(NodesEnum());
    28         vertices  = new DataSet(VerticesEnum());
    29         materials = new DataSet(MaterialsEnum());
     26        elements  = new DataSet(ElementsEnum);
     27        nodes     = new DataSet(NodesEnum);
     28        vertices  = new DataSet(VerticesEnum);
     29        materials = new DataSet(MaterialsEnum);
    3030       
    3131        /*Now, is the flag macayaealpattyn on? otherwise, do nothing: */
     
    5858                        if(iomodel->my_elements[i]){
    5959                               
    60                                 if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum()){ //elements of type 1 are Hutter type Tria. Don't create this elements.
     60                                if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum){ //elements of type 1 are Hutter type Tria. Don't create this elements.
    6161
    6262                                        /*Create and add tria element to elements dataset: */
     
    106106                for (i=0;i<iomodel->numberofelements;i++){
    107107                        if(iomodel->my_elements[i]){
    108                                 if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum() | *(iomodel->elements_type+2*i+0)==PattynFormulationEnum()){ //elements of type 1 are Hutter type Tria. Don't create this elements.
     108                                if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum | *(iomodel->elements_type+2*i+0)==PattynFormulationEnum){ //elements of type 1 are Hutter type Tria. Don't create this elements.
    109109                                        /*Create and add penta element to elements dataset: */
    110110                                        elements->AddObject(new Penta(i,iomodel));
  • issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateLoadsDiagnosticHoriz.cpp

    r3488 r3567  
    2121
    2222        /*Create loads: */
    23         loads   = new DataSet(LoadsEnum());
     23        loads   = new DataSet(LoadsEnum);
    2424
    2525        /*Now, is the flag macayaealpattyn on? otherwise, do nothing: */
     
    4545               
    4646                /*Do not create ice front if Hutter or Stokes elements*/
    47                 if ((int)*(iomodel->elements_type+2*element+0)==(HutterFormulationEnum() || StokesFormulationEnum())) continue;
     47                if ((int)*(iomodel->elements_type+2*element+0)==(HutterFormulationEnum || StokesFormulationEnum)) continue;
    4848
    4949                /*Create and  add load: */
  • issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp

    r3477 r3567  
    1919
    2020        /*Create constraints: */
    21         constraints = new DataSet(ConstraintsEnum());
     21        constraints = new DataSet(ConstraintsEnum);
    2222
    2323        /*Now, is the flag ishutter on? otherwise, do nothing: */
  • issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateElementsNodesAndMaterialsDiagnosticHutter.cpp

    r3508 r3567  
    2424
    2525        /*First create the elements, nodes and material properties: */
    26         elements  = new DataSet(ElementsEnum());
    27         nodes     = new DataSet(NodesEnum());
    28         vertices  = new DataSet(VerticesEnum());
    29         materials = new DataSet(MaterialsEnum());
     26        elements  = new DataSet(ElementsEnum);
     27        nodes     = new DataSet(NodesEnum);
     28        vertices  = new DataSet(VerticesEnum);
     29        materials = new DataSet(MaterialsEnum);
    3030
    3131        /*Now, is the flag ishutter on? otherwise, do nothing: */
  • issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateLoadsDiagnosticHutter.cpp

    r3456 r3567  
    1515
    1616        /*Create loads: */
    17         loads   = new DataSet(LoadsEnum());
     17        loads   = new DataSet(LoadsEnum);
    1818
    1919        /*Now, is the flag ishutter on? otherwise, do nothing: */
  • issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateConstraintsDiagnosticStokes.cpp

    r3456 r3567  
    1818
    1919        /*Create constraints: */
    20         constraints = new DataSet(ConstraintsEnum());
     20        constraints = new DataSet(ConstraintsEnum);
    2121
    2222        /*Now, is the flag ishutter on? otherwise, do nothing: */
  • issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateElementsNodesAndMaterialsDiagnosticStokes.cpp

    r3550 r3567  
    2525
    2626        /*First create the elements, nodes and material properties: */
    27         elements  = new DataSet(ElementsEnum());
    28         nodes     = new DataSet(NodesEnum());
    29         vertices  = new DataSet(VerticesEnum());
    30         materials = new DataSet(MaterialsEnum());
     27        elements  = new DataSet(ElementsEnum);
     28        nodes     = new DataSet(NodesEnum);
     29        vertices  = new DataSet(VerticesEnum);
     30        materials = new DataSet(MaterialsEnum);
    3131
    3232        /*Now, do we have Stokes elements?*/
     
    5959
    6060                if(iomodel->my_elements[i]){
    61                          if (*(iomodel->elements_type+2*i+1)==StokesFormulationEnum()){
     61                         if (*(iomodel->elements_type+2*i+1)==StokesFormulationEnum){
    6262
    6363                                 /*Create and add penta element to elements dataset: */
  • issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateLoadsDiagnosticStokes.cpp

    r3488 r3567  
    2525
    2626        /*Create loads: */
    27         loads   = new DataSet(LoadsEnum());
     27        loads   = new DataSet(LoadsEnum);
    2828
    2929        /*Now, is the flag ishutter on? otherwise, do nothing: */
     
    4949       
    5050                /*Do not create ice front if it is not a Stokes element*/
    51                 if (!(int)*(iomodel->elements_type+2*element+1)==StokesFormulationEnum()){
     51                if (!(int)*(iomodel->elements_type+2*element+1)==StokesFormulationEnum){
    5252                        continue;
    5353                }
  • issm/trunk/src/c/ModelProcessorx/DiagnosticVert/CreateConstraintsDiagnosticVert.cpp

    r3444 r3567  
    1717
    1818        /*Create constraints: */
    19         constraints = new DataSet(ConstraintsEnum());
     19        constraints = new DataSet(ConstraintsEnum);
    2020
    2121        /*return if 2d mesh*/
  • issm/trunk/src/c/ModelProcessorx/DiagnosticVert/CreateElementsNodesAndMaterialsDiagnosticVert.cpp

    r3476 r3567  
    2424
    2525        /*First create the elements, nodes and material properties: */
    26         elements  = new DataSet(ElementsEnum());
    27         nodes     = new DataSet(NodesEnum());
    28         vertices  = new DataSet(VerticesEnum());
    29         materials = new DataSet(MaterialsEnum());
     26        elements  = new DataSet(ElementsEnum);
     27        nodes     = new DataSet(NodesEnum);
     28        vertices  = new DataSet(VerticesEnum);
     29        materials = new DataSet(MaterialsEnum);
    3030
    3131        /*Now, is the flag macayaealpattyn on? otherwise, do nothing: */
  • issm/trunk/src/c/ModelProcessorx/DiagnosticVert/CreateLoadsDiagnosticVert.cpp

    r3443 r3567  
    1515
    1616        /*Create loads: */
    17         loads   = new DataSet(LoadsEnum());
     17        loads   = new DataSet(LoadsEnum);
    1818
    1919        /*Assign output pointer: */
  • issm/trunk/src/c/ModelProcessorx/Melting/CreateConstraintsMelting.cpp

    r3456 r3567  
    1818
    1919        /*Create constraints: */
    20         constraints = new DataSet(ConstraintsEnum());
     20        constraints = new DataSet(ConstraintsEnum);
    2121
    2222        /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these
  • issm/trunk/src/c/ModelProcessorx/Melting/CreateElementsNodesAndMaterialsMelting.cpp

    r3494 r3567  
    2424       
    2525        /*First create the elements, nodes and material properties: */
    26         elements  = new DataSet(ElementsEnum());
    27         nodes     = new DataSet(NodesEnum());
    28         vertices  = new DataSet(VerticesEnum());
    29         materials = new DataSet(MaterialsEnum());
     26        elements  = new DataSet(ElementsEnum);
     27        nodes     = new DataSet(NodesEnum);
     28        vertices  = new DataSet(VerticesEnum);
     29        materials = new DataSet(MaterialsEnum);
    3030
    3131        /*Partition elements and vertices and nodes: */
  • issm/trunk/src/c/ModelProcessorx/Melting/CreateLoadsMelting.cpp

    r3456 r3567  
    2121
    2222        /*Create loads: */
    23         loads   = new DataSet(LoadsEnum());
     23        loads   = new DataSet(LoadsEnum);
    2424
    2525        //create penalties for grids: no grid can have a temperature over the melting point
  • issm/trunk/src/c/ModelProcessorx/Partitioning.cpp

    r3559 r3567  
    2222void  Partitioning(bool** pmy_elements, bool** pmy_vertices, bool** pmy_nodes, bool** pmy_bordervertices, IoModel* iomodel, ConstDataHandle iomodel_handle){
    2323       
    24         if (iomodel->analysis_type==Prognostic2AnalysisEnum())
     24        if (iomodel->analysis_type==Prognostic2AnalysisEnum)
    2525                DiscontinuousGalerkinPartitioning(pmy_elements, pmy_vertices, pmy_nodes, pmy_bordervertices, iomodel, iomodel_handle);
    2626        else
  • issm/trunk/src/c/ModelProcessorx/Prognostic/CreateConstraintsPrognostic.cpp

    r3456 r3567  
    1919
    2020        /*Create constraints: */
    21         constraints = new DataSet(ConstraintsEnum());
     21        constraints = new DataSet(ConstraintsEnum);
    2222
    2323        /*Fetch data: */
  • issm/trunk/src/c/ModelProcessorx/Prognostic/CreateElementsNodesAndMaterialsPrognostic.cpp

    r3434 r3567  
    2424       
    2525        /*First create the elements, nodes and material properties: */
    26         elements  = new DataSet(ElementsEnum());
    27         nodes     = new DataSet(NodesEnum());
    28         vertices  = new DataSet(VerticesEnum());
    29         materials = new DataSet(MaterialsEnum());
     26        elements  = new DataSet(ElementsEnum);
     27        nodes     = new DataSet(NodesEnum);
     28        vertices  = new DataSet(VerticesEnum);
     29        materials = new DataSet(MaterialsEnum);
    3030
    3131        /*Partition elements and vertices and nodes: */
  • issm/trunk/src/c/ModelProcessorx/Prognostic/CreateLoadsPrognostic.cpp

    r3456 r3567  
    1515
    1616        /*Create loads: */
    17         loads   = new DataSet(LoadsEnum());
     17        loads   = new DataSet(LoadsEnum);
    1818       
    1919        /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these
  • issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateConstraintsPrognostic2.cpp

    r3359 r3567  
    1515       
    1616        /*Create constraints: */
    17         constraints = new DataSet(ConstraintsEnum());
     17        constraints = new DataSet(ConstraintsEnum);
    1818       
    1919        /*Assign output pointer: */
  • issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateElementsNodesAndMaterialsPrognostic2.cpp

    r3534 r3567  
    2727
    2828        /*First create the elements, nodes and material properties: */
    29         elements  = new DataSet(ElementsEnum());
    30         nodes     = new DataSet(NodesEnum());
    31         vertices  = new DataSet(VerticesEnum());
    32         materials = new DataSet(MaterialsEnum());
     29        elements  = new DataSet(ElementsEnum);
     30        nodes     = new DataSet(NodesEnum);
     31        vertices  = new DataSet(VerticesEnum);
     32        materials = new DataSet(MaterialsEnum);
    3333
    3434        /*Partition elements and vertices and nodes: */
  • issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateLoadsPrognostic2.cpp

    r3549 r3567  
    2929
    3030        /*Create loads: */
    31         loads   = new DataSet(LoadsEnum());
     31        loads   = new DataSet(LoadsEnum);
    3232
    3333        /*Get edges and elements*/
  • issm/trunk/src/c/ModelProcessorx/SlopeCompute/CreateConstraintsSlopeCompute.cpp

    r3456 r3567  
    1515
    1616        /*Create constraints: */
    17         constraints = new DataSet(ConstraintsEnum());
     17        constraints = new DataSet(ConstraintsEnum);
    1818       
    1919        /*Assign output pointer: */
  • issm/trunk/src/c/ModelProcessorx/SlopeCompute/CreateElementsNodesAndMaterialsSlopeCompute.cpp

    r3425 r3567  
    2424       
    2525        /*First create the elements, nodes and material properties: */
    26         elements  = new DataSet(ElementsEnum());
    27         nodes     = new DataSet(NodesEnum());
    28         vertices  = new DataSet(VerticesEnum());
    29         materials = new DataSet(MaterialsEnum());
     26        elements  = new DataSet(ElementsEnum);
     27        nodes     = new DataSet(NodesEnum);
     28        vertices  = new DataSet(VerticesEnum);
     29        materials = new DataSet(MaterialsEnum);
    3030
    3131        /*Partition elements and vertices and nodes: */
  • issm/trunk/src/c/ModelProcessorx/SlopeCompute/CreateLoadsSlopeCompute.cpp

    r3456 r3567  
    1515
    1616        /*Create loads: */
    17         loads   = new DataSet(LoadsEnum());
     17        loads   = new DataSet(LoadsEnum);
    1818
    1919        /*Now, is the flag isstokes on? otherwise, do nothing: */
  • issm/trunk/src/c/ModelProcessorx/Thermal/CreateConstraintsThermal.cpp

    r3455 r3567  
    1919
    2020        /*Create constraints: */
    21         constraints = new DataSet(ConstraintsEnum());
     21        constraints = new DataSet(ConstraintsEnum);
    2222
    2323        /*return if 2d mesh*/
  • issm/trunk/src/c/ModelProcessorx/Thermal/CreateElementsNodesAndMaterialsThermal.cpp

    r3491 r3567  
    2424       
    2525        /*First create the elements, nodes and material properties: */
    26         elements  = new DataSet(ElementsEnum());
    27         nodes     = new DataSet(NodesEnum());
    28         vertices  = new DataSet(VerticesEnum());
    29         materials = new DataSet(MaterialsEnum());
     26        elements  = new DataSet(ElementsEnum);
     27        nodes     = new DataSet(NodesEnum);
     28        vertices  = new DataSet(VerticesEnum);
     29        materials = new DataSet(MaterialsEnum);
    3030
    3131        /*Partition elements and vertices and nodes: */
  • issm/trunk/src/c/ModelProcessorx/Thermal/CreateLoadsThermal.cpp

    r3456 r3567  
    2222
    2323        /*Create loads: */
    24         loads   = new DataSet(LoadsEnum());
     24        loads   = new DataSet(LoadsEnum);
    2525
    2626        //create penalties for grids: no grid can have a temperature over the melting point
  • issm/trunk/src/c/PenaltyConstraintsx/RiftConstraints.cpp

    r3332 r3567  
    4646        for (i=0;i<loads->Size();i++){
    4747
    48                 if(RiftfrontEnum()==loads->GetEnum(i)){
     48                if(RiftfrontEnum==loads->GetEnum(i)){
    4949
    5050                        riftfront=(Riftfront*)loads->GetObjectByOffset(i);
     
    7777        for (i=0;i<loads->Size();i++){
    7878
    79                 if(RiftfrontEnum()==loads->GetEnum(i)){
     79                if(RiftfrontEnum==loads->GetEnum(i)){
    8080                        found=1;
    8181                        break;
     
    104104        for (i=0;i<loads->Size();i++){
    105105
    106                 if(RiftfrontEnum()==loads->GetEnum(i)){
     106                if(RiftfrontEnum==loads->GetEnum(i)){
    107107
    108108                        riftfront=(Riftfront*)loads->GetObjectByOffset(i);
     
    142142        for (i=0;i<loads->Size();i++){
    143143
    144                 if(RiftfrontEnum()==loads->GetEnum(i)){
     144                if(RiftfrontEnum==loads->GetEnum(i)){
    145145
    146146                        riftfront=(Riftfront*)loads->GetObjectByOffset(i);
     
    164164        for (i=0;i<loads->Size();i++){
    165165
    166                 if (RiftfrontEnum()==loads->GetEnum(i)){
     166                if (RiftfrontEnum==loads->GetEnum(i)){
    167167
    168168                        riftfront=(Riftfront*)loads->GetObjectByOffset(i);
     
    199199        for (i=0;i<loads->Size();i++){
    200200
    201                 if (RiftfrontEnum()==loads->GetEnum(i)){
     201                if (RiftfrontEnum==loads->GetEnum(i)){
    202202
    203203                        riftfront=(Riftfront*)loads->GetObjectByOffset(i);
     
    230230        for (i=0;i<loads->Size();i++){
    231231
    232                 if (RiftfrontEnum()==loads->GetEnum(i)){
     232                if (RiftfrontEnum==loads->GetEnum(i)){
    233233
    234234                        riftfront=(Riftfront*)loads->GetObjectByOffset(i);
     
    253253        for (i=0;i<loads->Size();i++){
    254254
    255                 if (RiftfrontEnum()==loads->GetEnum(i)){
     255                if (RiftfrontEnum==loads->GetEnum(i)){
    256256
    257257                        riftfront=(Riftfront*)loads->GetObjectByOffset(i);
     
    291291        for (i=0;i<loads->Size();i++){
    292292
    293                 if (RiftfrontEnum()==loads->GetEnum(i)){
     293                if (RiftfrontEnum==loads->GetEnum(i)){
    294294
    295295                        riftfront=(Riftfront*)loads->GetObjectByOffset(i);
     
    330330        for (i=0;i<loads->Size();i++){
    331331
    332                 if (RiftfrontEnum()==loads->GetEnum(i)){
     332                if (RiftfrontEnum==loads->GetEnum(i)){
    333333
    334334                        riftfront=(Riftfront*)loads->GetObjectByOffset(i);
  • issm/trunk/src/c/ProcessParamsx/ProcessParamsx.cpp

    r3534 r3567  
    4141                param_type=parameters->GetEnum(i);
    4242
    43                 if(param_type==ParamEnum()){
     43                if(param_type==ParamEnum){
    4444
    4545                        Param* param=(Param*)parameters->GetObjectByOffset(i);
  • issm/trunk/src/c/Qmux/DakotaResponses.cpp

    r3446 r3567  
    238238
    239239                        if(isstokes){
    240                                 femmodel=model->GetFormulation(DiagnosticAnalysisEnum(),StokesAnalysisEnum());
     240                                femmodel=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum);
    241241                        }
    242242                        if(ismacayealpattyn){
    243                                 femmodel=model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     243                                femmodel=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum);
    244244                        }
    245245                        if(ishutter){
    246                                 femmodel=model->GetFormulation(DiagnosticAnalysisEnum(),HutterAnalysisEnum());
     246                                femmodel=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum);
    247247                        }       
    248248
     
    300300
    301301                        if(isstokes){
    302                                 femmodel=model->GetFormulation(DiagnosticAnalysisEnum(),StokesAnalysisEnum());
     302                                femmodel=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum);
    303303                        }
    304304                        if(ismacayealpattyn){
    305                                 femmodel=model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     305                                femmodel=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum);
    306306                        }
    307307                        if(ishutter){
    308                                 femmodel=model->GetFormulation(DiagnosticAnalysisEnum(),HutterAnalysisEnum());
     308                                femmodel=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum);
    309309                        }
    310310
  • issm/trunk/src/c/Qmux/SpawnCoreParallel.cpp

    r3332 r3567  
    6060               
    6161        /*First off, recover the response descriptors for the response functions: */
    62         param=(Param*)model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum())->parameters->FindParamObject("responsedescriptors");
     62        param=(Param*)model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum)->parameters->FindParamObject("responsedescriptors");
    6363        if(!param)ISSMERROR(" could not find response descriptors!");
    6464
     
    123123
    124124        //initialize results:
    125         results=new DataSet(ResultsEnum());
     125        results=new DataSet(ResultsEnum);
    126126
    127127        /*Modify core inputs to reflect the dakota variables inputs: */
    128         inputs->UpdateFromDakota(variables,variables_descriptors,numvariables,model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum())->parameters,qmu_part,qmu_npart); //diagnostic horiz model is the one holding the parameters for Dakota.
     128        inputs->UpdateFromDakota(variables,variables_descriptors,numvariables,model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum)->parameters,qmu_part,qmu_npart); //diagnostic horiz model is the one holding the parameters for Dakota.
    129129
    130130        /*Run the analysis core solution sequence, with the updated inputs: */
    131         if(analysis_type==DiagnosticAnalysisEnum()){
     131        if(analysis_type==DiagnosticAnalysisEnum){
    132132                       
    133133                if(verbose)_printf_("Starting diagnostic core\n");
     
    136136
    137137        }
    138         else if(analysis_type==ThermalAnalysisEnum()){
     138        else if(analysis_type==ThermalAnalysisEnum){
    139139               
    140140                if(verbose)_printf_("Starting thermal core\n");
     
    142142
    143143        }
    144         else if(analysis_type==PrognosticAnalysisEnum()){
     144        else if(analysis_type==PrognosticAnalysisEnum){
    145145
    146146                if(verbose)_printf_("Starting prognostic core\n");
     
    148148
    149149        }
    150         else if(analysis_type==TransientAnalysisEnum()){
     150        else if(analysis_type==TransientAnalysisEnum){
    151151
    152152                if(verbose)_printf_("Starting transient core\n");
  • issm/trunk/src/c/UpdateGeometryx/UpdateGeometryx.cpp

    r3500 r3567  
    5757        for(i=0;i<materials->Size();i++){
    5858                Object* object=materials->GetObjectByOffset(i);
    59                 if  (object->Enum()==MatparEnum()){
     59                if  (object->Enum()==MatparEnum){
    6060                        matpar=(Matpar*)object;
    6161                        break;
  • issm/trunk/src/c/objects/Beam.cpp

    r3543 r3567  
    205205
    206206        /*get enum type of Beam: */
    207         enum_type=BeamEnum();
     207        enum_type=BeamEnum;
    208208       
    209209        /*marshall enum: */
     
    377377
    378378        /*Just branch to the correct element stiffness matrix generator, according to the type of analysis we are carrying out: */
    379         if (analysis_type==DiagnosticAnalysisEnum()) {
    380        
    381                 if (sub_analysis_type==HutterAnalysisEnum()) {
     379        if (analysis_type==DiagnosticAnalysisEnum) {
     380       
     381                if (sub_analysis_type==HutterAnalysisEnum) {
    382382
    383383                        CreateKMatrixDiagnosticHutter( Kgg,inputs,analysis_type,sub_analysis_type);
     
    430430       
    431431        /*Just branch to the correct load generator, according to the type of analysis we are carrying out: */
    432         if (analysis_type==DiagnosticAnalysisEnum()) {
    433                 if (sub_analysis_type==HutterAnalysisEnum()) {
     432        if (analysis_type==DiagnosticAnalysisEnum) {
     433                if (sub_analysis_type==HutterAnalysisEnum) {
    434434                        CreatePVectorDiagnosticHutter( pg,inputs,analysis_type,sub_analysis_type);
    435435                }
     
    577577int Beam::Enum(void){
    578578
    579         return BeamEnum();
     579        return BeamEnum;
    580580
    581581}
  • issm/trunk/src/c/objects/DofIndexing.cpp

    r3483 r3567  
    150150
    151151        /*get enum type of DofIndexing: */
    152         enum_type=DofIndexingEnum();
     152        enum_type=DofIndexingEnum;
    153153       
    154154        /*marshall enum: */
  • issm/trunk/src/c/objects/DofVec.cpp

    r3470 r3567  
    191191
    192192        /*get enum type of DofVec: */
    193         enum_type=DofVecEnum();
     193        enum_type=DofVecEnum;
    194194       
    195195        /*marshall enum: */
     
    227227int DofVec::Enum(void){
    228228
    229         return DofVecEnum();
     229        return DofVecEnum;
    230230
    231231}
  • issm/trunk/src/c/objects/Element.cpp

    r3470 r3567  
    99int Element::Enum(void){
    1010
    11         return ElementEnum();
     11        return ElementEnum;
    1212}
  • issm/trunk/src/c/objects/ElementProperties.cpp

    r3483 r3567  
    257257
    258258        /*get enum type of ElementProperties: */
    259         enum_type=ElementPropertiesEnum();
     259        enum_type=ElementPropertiesEnum;
    260260       
    261261        /*marshall enum: */
  • issm/trunk/src/c/objects/Hook.cpp

    r3483 r3567  
    117117
    118118        /*get enum type of Hook: */
    119         enum_type=HookEnum();
     119        enum_type=HookEnum;
    120120       
    121121        /*marshall enum: */
  • issm/trunk/src/c/objects/Icefront.cpp

    r3511 r3567  
    8585        if (strcmp(iomodel->meshtype,"2d")==0){
    8686                segment_width=4;
    87                 element_type=TriaEnum();
     87                element_type=TriaEnum;
    8888        }
    8989        else{
    9090                segment_width=6;
    91                 element_type=PentaEnum();
     91                element_type=PentaEnum;
    9292        }
    9393        element=(int)(*(iomodel->pressureload+segment_width*i+segment_width-2)-1); //element is in the penultimate column (grid1 grid2 ... elem fill)
     
    112112        icefront_b[1]=iomodel->bed[i2-1];
    113113
    114         if (iomodel->sub_analysis_type==HorizAnalysisEnum()){
    115                 if ((int)*(iomodel->elements_type+2*element+0)==MacAyealFormulationEnum()){ //this is a collapsed 3d element, icefront will be 2d
     114        if (iomodel->sub_analysis_type==HorizAnalysisEnum){
     115                if ((int)*(iomodel->elements_type+2*element+0)==MacAyealFormulationEnum){ //this is a collapsed 3d element, icefront will be 2d
    116116                        strcpy(icefront_type,"segment");
    117117                }
    118                 else if ((int)*(iomodel->elements_type+2*element+0)==PattynFormulationEnum()){ //this is a real 3d element, icefront will be 3d.
     118                else if ((int)*(iomodel->elements_type+2*element+0)==PattynFormulationEnum){ //this is a real 3d element, icefront will be 3d.
    119119                        strcpy(icefront_type,"quad");
    120120                        i3=(int)*(iomodel->pressureload+segment_width*i+2);
     
    133133                }
    134134        }
    135         else if (iomodel->sub_analysis_type==StokesAnalysisEnum()){
     135        else if (iomodel->sub_analysis_type==StokesAnalysisEnum){
    136136                //We have a Stokes element, so we need a 3d Icefront
    137137                strcpy(icefront_type,"quad");
     
    209209
    210210        /*get enum type of Icefront: */
    211         enum_type=IcefrontEnum();
     211        enum_type=IcefrontEnum;
    212212       
    213213        /*marshall enum: */
     
    306306
    307307        /*Just branch to the correct element icefront vector generator, according to the type of analysis we are carrying out: */
    308         if (analysis_type==ControlAnalysisEnum()){
     308        if (analysis_type==ControlAnalysisEnum){
    309309               
    310310                CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type);
    311311
    312312        }
    313         else if (analysis_type==DiagnosticAnalysisEnum()){
    314        
    315                 if (sub_analysis_type==HorizAnalysisEnum()){
     313        else if (analysis_type==DiagnosticAnalysisEnum){
     314       
     315                if (sub_analysis_type==HorizAnalysisEnum){
    316316               
    317317                        CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type);
    318318
    319319                }
    320                 else if (sub_analysis_type==StokesAnalysisEnum()){
     320                else if (sub_analysis_type==StokesAnalysisEnum){
    321321                       
    322322                        CreatePVectorDiagnosticStokes( pg,inputs,analysis_type,sub_analysis_type);
     
    380380
    381381        //check that the element is onbed (collapsed formulation) otherwise:pe=0
    382         if(element_type==PentaEnum()){
     382        if(element_type==PentaEnum){
    383383                if  (!element->GetOnBed()){
    384384                        return;
     
    389389               
    390390        /*Identify which grids are comprised in the segment. */
    391         if(element_type==TriaEnum())element_nodes=(Node**)xmalloc(3*sizeof(Node*));
    392         if(element_type==PentaEnum())element_nodes=(Node**)xmalloc(6*sizeof(Node*));
     391        if(element_type==TriaEnum)element_nodes=(Node**)xmalloc(3*sizeof(Node*));
     392        if(element_type==PentaEnum)element_nodes=(Node**)xmalloc(6*sizeof(Node*));
    393393        element->GetNodes((void**)element_nodes);
    394394
     
    496496
    497497        /*check icefront is associated to a pentaelem: */
    498         if(element_type!=PentaEnum()){
     498        if(element_type!=PentaEnum){
    499499                ISSMERROR(" quad icefront is associated to a TriaElem!");
    500500        }
     
    510510       
    511511        //Identify which grids are comprised in the quad:
    512         if(element_type==PentaEnum())element_nodes=(Node**)xmalloc(6*sizeof(Node*));
     512        if(element_type==PentaEnum)element_nodes=(Node**)xmalloc(6*sizeof(Node*));
    513513        element->GetNodes((void**)element_nodes);
    514514
     
    653653
    654654        /*check icefront is associated to a pentaelem: */
    655         if(element_type!=PentaEnum()){
     655        if(element_type!=PentaEnum){
    656656                ISSMERROR(" quad icefront is associated to a TriaElem!");
    657657        }
     
    667667       
    668668        //Identify which grids are comprised in the quad:
    669         if(element_type==PentaEnum())element_nodes=(Node**)xmalloc(6*sizeof(Node*));
     669        if(element_type==PentaEnum)element_nodes=(Node**)xmalloc(6*sizeof(Node*));
    670670        element->GetNodes((void**)element_nodes);
    671671
     
    855855int Icefront::Enum(void){
    856856
    857         return IcefrontEnum();
     857        return IcefrontEnum;
    858858
    859859}
     
    11281128
    11291129                        //Now deal with water pressure:
    1130                         if(fill==WaterEnum()){ //icefront ends in water
     1130                        if(fill==WaterEnum){ //icefront ends in water
    11311131                                water_level_above_g_tria=min(0,z_g[i]);//0 if the gaussian point is above water level
    11321132                                water_pressure_tria=rho_water*gravity*water_level_above_g_tria;
    11331133                        }
    1134                         else if(fill==AirEnum()){
     1134                        else if(fill==AirEnum){
    11351135                                water_pressure_tria=0;
    11361136                        }
     
    13821382
    13831383                        //Now deal with water pressure:
    1384                         if(fill==WaterEnum()){ //icefront ends in water
     1384                        if(fill==WaterEnum){ //icefront ends in water
    13851385                                water_level_above_g_tria=min(0,z_g[i]);//0 if the gaussian point is above water level
    13861386                                water_pressure_tria=rho_water*gravity*water_level_above_g_tria;
    13871387                        }
    1388                         else if(fill==AirEnum()){
     1388                        else if(fill==AirEnum){
    13891389                                water_pressure_tria=0;
    13901390                        }
     
    14631463                bed=b1*(1+segment_gauss_coord[ig])/2+b2*(1-segment_gauss_coord[ig])/2;
    14641464
    1465                 if (fill==WaterEnum()){
     1465                if (fill==WaterEnum){
    14661466                        //icefront ends in water:
    14671467                        ice_pressure=1.0/2.0*gravity*rho_ice*pow(thickness,2);
     
    14731473                        water_pressure=1.0/2.0*gravity*rho_water*(pow(surface_under_water,2) - pow(base_under_water,2));
    14741474                }
    1475                 else if (fill==AirEnum()){
     1475                else if (fill==AirEnum){
    14761476                        ice_pressure=1.0/2.0*gravity*rho_ice*pow(thickness,2);
    14771477                        air_pressure=0;
  • issm/trunk/src/c/objects/Input.cpp

    r3463 r3567  
    153153int Input::Enum(void){
    154154
    155         return InputEnum();
     155        return InputEnum;
    156156
    157157}
  • issm/trunk/src/c/objects/Load.cpp

    r246 r3567  
    1010int Load::Enum(void){
    1111
    12         return LoadEnum();
     12        return LoadEnum;
    1313}
    1414               
  • issm/trunk/src/c/objects/Material.cpp

    r246 r3567  
    1010int Material::Enum(void){
    1111
    12         return MaterialEnum();
     12        return MaterialEnum;
    1313}
    1414               
  • issm/trunk/src/c/objects/Matice.cpp

    r3541 r3567  
    102102
    103103        /*get enum type of Matice: */
    104         enum_type=MaticeEnum();
     104        enum_type=MaticeEnum;
    105105       
    106106        /*marshall enum: */
     
    175175int Matice::Enum(void){
    176176
    177         return MaticeEnum();
     177        return MaticeEnum;
    178178
    179179}
  • issm/trunk/src/c/objects/Matpar.cpp

    r3509 r3567  
    4646        double  matpar_g;
    4747
    48         if (iomodel->analysis_type==DiagnosticAnalysisEnum() && iomodel->sub_analysis_type==HutterAnalysisEnum()){
     48        if (iomodel->analysis_type==DiagnosticAnalysisEnum && iomodel->sub_analysis_type==HutterAnalysisEnum){
    4949                if (strcmp(iomodel->meshtype,"2d")==0){
    5050                        matpar_mid=iomodel->numberofvertices+1; //put it at the end of the materials
     
    133133
    134134        /*get enum type of Matpar: */
    135         enum_type=MatparEnum();
     135        enum_type=MatparEnum;
    136136       
    137137        /*marshall enum: */
     
    221221int Matpar::Enum(void){
    222222
    223         return MatparEnum();
     223        return MatparEnum;
    224224
    225225}
  • issm/trunk/src/c/objects/Model.cpp

    r3417 r3567  
    143143        }
    144144        /*sub_analysis_type should default to none, as it was not specified: */
    145         sub_analysis_type=NoneAnalysisEnum();
     145        sub_analysis_type=NoneAnalysisEnum;
    146146       
    147147        if (sub_analysis_type!=0){
  • issm/trunk/src/c/objects/Node.cpp

    r3557 r3567  
    9595
    9696        /*Diagnostic Horiz*/
    97         if (iomodel->analysis_type==DiagnosticAnalysisEnum() && iomodel->sub_analysis_type==HorizAnalysisEnum()){
     97        if (iomodel->analysis_type==DiagnosticAnalysisEnum && iomodel->sub_analysis_type==HorizAnalysisEnum){
    9898                if (strcmp(iomodel->meshtype,"3d")==0){
    9999                        /*We have a  3d mesh, we may have collapsed elements, hence dead grids. Freeze them out: */
     
    115115
    116116        /*Diagnostic Stokes*/
    117         if (iomodel->analysis_type==DiagnosticAnalysisEnum() && iomodel->sub_analysis_type==StokesAnalysisEnum()){
     117        if (iomodel->analysis_type==DiagnosticAnalysisEnum && iomodel->sub_analysis_type==StokesAnalysisEnum){
    118118                /*On a 3d mesh, in stokes formualtions, only stokes grids are free, the others are frozen: */
    119119                if (!iomodel->borderstokes) ISSMERROR("iomodel->borderstokes is NULL");
     
    132132
    133133        /*Diagnostic Hutter*/
    134         if (iomodel->analysis_type==DiagnosticAnalysisEnum() && iomodel->sub_analysis_type==HutterAnalysisEnum()){
     134        if (iomodel->analysis_type==DiagnosticAnalysisEnum && iomodel->sub_analysis_type==HutterAnalysisEnum){
    135135                /*Spc all nodes that are not Hutter*/
    136136                if (!iomodel->gridonhutter) ISSMERROR("iomodel->gridonhutter is NULL");
     
    144144        /*Prognostic/ Melting/ Slopecompute/ Balancedvelocities/ Balancedthickness*/
    145145        if (
    146                                 iomodel->analysis_type==PrognosticAnalysisEnum() ||
    147                                 iomodel->analysis_type==MeltingAnalysisEnum() ||
    148                                 iomodel->analysis_type==SlopecomputeAnalysisEnum() ||
    149                                 iomodel->analysis_type==BalancedvelocitiesAnalysisEnum() ||
    150                                 iomodel->analysis_type==BalancedthicknessAnalysisEnum()
     146                                iomodel->analysis_type==PrognosticAnalysisEnum ||
     147                                iomodel->analysis_type==MeltingAnalysisEnum ||
     148                                iomodel->analysis_type==SlopecomputeAnalysisEnum ||
     149                                iomodel->analysis_type==BalancedvelocitiesAnalysisEnum ||
     150                                iomodel->analysis_type==BalancedthicknessAnalysisEnum
    151151                                ){
    152152                if (strcmp(iomodel->meshtype,"3d")==0){
     
    295295int Node::Enum(void){
    296296
    297         return NodeEnum();
     297        return NodeEnum;
    298298
    299299}
     
    354354
    355355        /*get enum type of Node: */
    356         enum_type=NodeEnum();
     356        enum_type=NodeEnum;
    357357       
    358358        /*marshall enum: */
  • issm/trunk/src/c/objects/NodeProperties.cpp

    r3483 r3567  
    8080
    8181        /*get enum type of NodeProperties: */
    82         enum_type=NodePropertiesEnum();
     82        enum_type=NodePropertiesEnum;
    8383       
    8484        /*marshall enum: */
  • issm/trunk/src/c/objects/Numericalflux.cpp

    r3565 r3567  
    9696
    9797        /*get enum type of Numericalflux: */
    98         enum_type=NumericalfluxEnum();
     98        enum_type=NumericalfluxEnum;
    9999       
    100100        /*marshall enum: */
     
    571571int Numericalflux::Enum(void){
    572572
    573         return NumericalfluxEnum();
     573        return NumericalfluxEnum;
    574574
    575575}
  • issm/trunk/src/c/objects/Numpar.cpp

    r3463 r3567  
    6060
    6161        /*get enum type of Numpar: */
    62         enum_type=NumparEnum();
     62        enum_type=NumparEnum;
    6363       
    6464        /*marshall enum: */
     
    204204int Numpar::Enum(void){
    205205
    206         return NumparEnum();
     206        return NumparEnum;
    207207
    208208}
  • issm/trunk/src/c/objects/Param.cpp

    r3463 r3567  
    233233
    234234        /*get enum type of Param: */
    235         enum_type=ParamEnum();
     235        enum_type=ParamEnum;
    236236       
    237237        /*marshall enum: */
     
    437437int Param::Enum(void){
    438438
    439         return ParamEnum();
     439        return ParamEnum;
    440440}
    441441/*}}}*/
  • issm/trunk/src/c/objects/Pengrid.cpp

    r3463 r3567  
    6464
    6565        /*get enum type of Pengrid: */
    66         enum_type=PengridEnum();
     66        enum_type=PengridEnum;
    6767       
    6868        /*marshall enum: */
     
    219219int Pengrid::Enum(void){
    220220
    221         return PengridEnum();
     221        return PengridEnum;
    222222}
    223223/*}}}1*/
     
    255255void  Pengrid::PenaltyConstrain(int* punstable,void* vinputs,int analysis_type,int sub_analysis_type){
    256256
    257         if ((analysis_type==DiagnosticAnalysisEnum()) && ((sub_analysis_type==StokesAnalysisEnum()))){
     257        if ((analysis_type==DiagnosticAnalysisEnum) && ((sub_analysis_type==StokesAnalysisEnum))){
    258258
    259259                /*No penalty to check*/
     
    261261
    262262        }
    263         else if (analysis_type==ThermalAnalysisEnum()){
     263        else if (analysis_type==ThermalAnalysisEnum){
    264264               
    265265                PenaltyConstrainThermal(punstable,vinputs,analysis_type,sub_analysis_type);
    266266               
    267267        }
    268         else if (analysis_type==MeltingAnalysisEnum()){
     268        else if (analysis_type==MeltingAnalysisEnum){
    269269                       
    270270                /*No penalty to check*/
     
    362362void  Pengrid::PenaltyCreateKMatrix(Mat Kgg,void* inputs,double kmax,int analysis_type,int sub_analysis_type){
    363363
    364         if ((analysis_type==DiagnosticAnalysisEnum()) && ((sub_analysis_type==StokesAnalysisEnum()))){
     364        if ((analysis_type==DiagnosticAnalysisEnum) && ((sub_analysis_type==StokesAnalysisEnum))){
    365365
    366366                PenaltyCreateKMatrixDiagnosticStokes( Kgg,inputs,kmax,analysis_type,sub_analysis_type);
    367367        }
    368         else if (analysis_type==ThermalAnalysisEnum()){
     368        else if (analysis_type==ThermalAnalysisEnum){
    369369               
    370370                PenaltyCreateKMatrixThermal( Kgg,inputs,kmax,analysis_type,sub_analysis_type);
    371371               
    372372        }
    373         else if (analysis_type==MeltingAnalysisEnum()){
     373        else if (analysis_type==MeltingAnalysisEnum){
    374374                       
    375375                PenaltyCreateKMatrixMelting( Kgg,inputs,kmax,analysis_type,sub_analysis_type);
     
    500500void  Pengrid::PenaltyCreatePVector(Vec pg,void* inputs,double kmax,int analysis_type,int sub_analysis_type){
    501501
    502         if (analysis_type==ThermalAnalysisEnum()){
     502        if (analysis_type==ThermalAnalysisEnum){
    503503               
    504504                PenaltyCreatePVectorThermal( pg,inputs,kmax,analysis_type,sub_analysis_type);
    505505               
    506506        }
    507         else if (analysis_type==MeltingAnalysisEnum()){
     507        else if (analysis_type==MeltingAnalysisEnum){
    508508                       
    509509                PenaltyCreatePVectorMelting( pg,inputs,kmax,analysis_type,sub_analysis_type);
    510510
    511511        }
    512         else if (analysis_type==DiagnosticAnalysisEnum()){
     512        else if (analysis_type==DiagnosticAnalysisEnum){
    513513
    514514                /*No loads applied, do nothing: */
  • issm/trunk/src/c/objects/Penpair.cpp

    r3463 r3567  
    6161
    6262        /*get enum type of Penpair: */
    63         enum_type=PenpairEnum();
     63        enum_type=PenpairEnum;
    6464       
    6565        /*marshall enum: */
     
    198198int Penpair::Enum(void){
    199199
    200         return PenpairEnum();
     200        return PenpairEnum;
    201201}
    202202/*}}}1*/
  • issm/trunk/src/c/objects/Penta.cpp

    r3529 r3567  
    8181        //elements of type 3 are macayeal type penta. we collapse the formulation on their base.
    8282        if(iomodel->elements_type){
    83                 if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum()){
     83                if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum){
    8484                        this->properties.collapse=1;
    8585                }
     
    196196
    197197        /*get enum type of Penta: */
    198         enum_type=PentaEnum();
     198        enum_type=PentaEnum;
    199199
    200200        /*marshall enum: */
     
    387387
    388388        /*Check analysis_types*/
    389         if (analysis_type!=DiagnosticAnalysisEnum() || sub_analysis_type!=StokesAnalysisEnum()) ISSMERROR("Not supported yet!");
     389        if (analysis_type!=DiagnosticAnalysisEnum || sub_analysis_type!=StokesAnalysisEnum) ISSMERROR("Not supported yet!");
    390390
    391391        /*recover pointers: */
     
    555555
    556556        /*Just branch to the correct element stiffness matrix generator, according to the type of analysis we are carrying out: */
    557         if (analysis_type==ControlAnalysisEnum()){
     557        if (analysis_type==ControlAnalysisEnum){
    558558
    559559                CreateKMatrixDiagnosticHoriz( Kgg,inputs,analysis_type,sub_analysis_type);
    560560
    561561        }
    562         else if (analysis_type==DiagnosticAnalysisEnum()){
    563 
    564                 if (sub_analysis_type==HorizAnalysisEnum()){
     562        else if (analysis_type==DiagnosticAnalysisEnum){
     563
     564                if (sub_analysis_type==HorizAnalysisEnum){
    565565
    566566                        CreateKMatrixDiagnosticHoriz( Kgg,inputs,analysis_type,sub_analysis_type);
    567567                }
    568                 else if (sub_analysis_type==VertAnalysisEnum()){
     568                else if (sub_analysis_type==VertAnalysisEnum){
    569569
    570570                        CreateKMatrixDiagnosticVert( Kgg,inputs,analysis_type,sub_analysis_type);
    571571                }
    572                 else if (sub_analysis_type==StokesAnalysisEnum()){
     572                else if (sub_analysis_type==StokesAnalysisEnum){
    573573
    574574                        CreateKMatrixDiagnosticStokes( Kgg,inputs,analysis_type,sub_analysis_type);
     
    577577                else ISSMERROR(exprintf("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet"));
    578578        }
    579         else if (analysis_type==SlopecomputeAnalysisEnum()){
     579        else if (analysis_type==SlopecomputeAnalysisEnum){
    580580
    581581                CreateKMatrixSlopeCompute( Kgg,inputs,analysis_type,sub_analysis_type);
    582582        }
    583         else if (analysis_type==PrognosticAnalysisEnum()){
     583        else if (analysis_type==PrognosticAnalysisEnum){
    584584
    585585                CreateKMatrixPrognostic( Kgg,inputs,analysis_type,sub_analysis_type);
    586586        }
    587         else if (analysis_type==BalancedthicknessAnalysisEnum()){
     587        else if (analysis_type==BalancedthicknessAnalysisEnum){
    588588
    589589                CreateKMatrixBalancedthickness( Kgg,inputs,analysis_type,sub_analysis_type);
    590590        }
    591         else if (analysis_type==BalancedvelocitiesAnalysisEnum()){
     591        else if (analysis_type==BalancedvelocitiesAnalysisEnum){
    592592
    593593                CreateKMatrixBalancedvelocities( Kgg,inputs,analysis_type,sub_analysis_type);
    594594        }
    595         else if (analysis_type==ThermalAnalysisEnum()){
     595        else if (analysis_type==ThermalAnalysisEnum){
    596596
    597597                CreateKMatrixThermal( Kgg,inputs,analysis_type,sub_analysis_type);
    598598        }
    599         else if (analysis_type==MeltingAnalysisEnum()){
     599        else if (analysis_type==MeltingAnalysisEnum){
    600600
    601601                CreateKMatrixMelting( Kgg,inputs,analysis_type,sub_analysis_type);
     
    16831683
    16841684        /*Just branch to the correct element stiffness matrix generator, according to the type of analysis we are carrying out: */
    1685         if (analysis_type==ControlAnalysisEnum()){
     1685        if (analysis_type==ControlAnalysisEnum){
    16861686
    16871687                CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type);
    16881688        }
    1689         else if (analysis_type==DiagnosticAnalysisEnum()){
    1690 
    1691                 if (sub_analysis_type==HorizAnalysisEnum()){
     1689        else if (analysis_type==DiagnosticAnalysisEnum){
     1690
     1691                if (sub_analysis_type==HorizAnalysisEnum){
    16921692
    16931693                        CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type);
    16941694                }
    1695                 else if (sub_analysis_type==VertAnalysisEnum()){
     1695                else if (sub_analysis_type==VertAnalysisEnum){
    16961696
    16971697                        CreatePVectorDiagnosticVert( pg,inputs,analysis_type,sub_analysis_type);
    16981698                }
    1699                 else if (sub_analysis_type==StokesAnalysisEnum()){
     1699                else if (sub_analysis_type==StokesAnalysisEnum){
    17001700
    17011701                        CreatePVectorDiagnosticStokes( pg,inputs,analysis_type,sub_analysis_type);
     
    17031703                else ISSMERROR(exprintf("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet"));
    17041704        }
    1705         else if (analysis_type==SlopecomputeAnalysisEnum()){
     1705        else if (analysis_type==SlopecomputeAnalysisEnum){
    17061706
    17071707                CreatePVectorSlopeCompute( pg,inputs,analysis_type,sub_analysis_type);
    17081708        }
    1709         else if (analysis_type==PrognosticAnalysisEnum()){
     1709        else if (analysis_type==PrognosticAnalysisEnum){
    17101710
    17111711                CreatePVectorPrognostic( pg,inputs,analysis_type,sub_analysis_type);
    17121712        }
    1713         else if (analysis_type==BalancedthicknessAnalysisEnum()){
     1713        else if (analysis_type==BalancedthicknessAnalysisEnum){
    17141714
    17151715                CreatePVectorBalancedthickness( pg,inputs,analysis_type,sub_analysis_type);
    17161716        }
    1717         else if (analysis_type==BalancedvelocitiesAnalysisEnum()){
     1717        else if (analysis_type==BalancedvelocitiesAnalysisEnum){
    17181718
    17191719                CreatePVectorBalancedvelocities( pg,inputs,analysis_type,sub_analysis_type);
    17201720        }
    1721         else if (analysis_type==ThermalAnalysisEnum()){
     1721        else if (analysis_type==ThermalAnalysisEnum){
    17221722
    17231723                CreatePVectorThermal( pg,inputs,analysis_type,sub_analysis_type);
    17241724        }
    1725         else if (analysis_type==MeltingAnalysisEnum()){
     1725        else if (analysis_type==MeltingAnalysisEnum){
    17261726
    17271727                CreatePVectorMelting( pg,inputs,analysis_type,sub_analysis_type);
     
    26202620int Penta::Enum(void){
    26212621
    2622         return PentaEnum();
     2622        return PentaEnum;
    26232623
    26242624}
     
    40594059        if (!this->properties.onbed) return;
    40604060
    4061         if (sub_analysis_type==HorizAnalysisEnum()){
     4061        if (sub_analysis_type==HorizAnalysisEnum){
    40624062
    40634063                /*MacAyeal or Pattyn*/
     
    40674067                return;
    40684068        }
    4069         else if (sub_analysis_type==StokesAnalysisEnum()){
     4069        else if (sub_analysis_type==StokesAnalysisEnum){
    40704070
    40714071                /*Stokes*/
  • issm/trunk/src/c/objects/Result.cpp

    r3463 r3567  
    202202int Result::Enum(void){
    203203
    204         return ResultEnum();
     204        return ResultEnum;
    205205
    206206}
  • issm/trunk/src/c/objects/Rgb.cpp

    r3463 r3567  
    5252
    5353        /*get enum type of Rgb: */
    54         enum_type=RgbEnum();
     54        enum_type=RgbEnum;
    5555       
    5656        /*marshall enum: */
     
    130130int Rgb::Enum(void){
    131131
    132         return RgbEnum();
     132        return RgbEnum;
    133133
    134134}
  • issm/trunk/src/c/objects/Riftfront.cpp

    r3463 r3567  
    182182
    183183        /*get enum type of Riftfront: */
    184         enum_type=RiftfrontEnum();
     184        enum_type=RiftfrontEnum;
    185185       
    186186        /*marshall enum: */
     
    453453int Riftfront::Enum(void){
    454454
    455         return RiftfrontEnum();
     455        return RiftfrontEnum;
    456456
    457457}
     
    729729
    730730                /*Ok, now compute the pressure (in norm) that is being applied to the flanks, depending on the type of fill: */
    731                 if(fill==WaterEnum()){
     731                if(fill==WaterEnum){
    732732                        if(shelf){
    733733                                /*We are on an ice shelf, hydrostatic equilibrium is used to determine the pressure for water fill: */
     
    739739                        }
    740740                }
    741                 else if(fill==AirEnum()){
     741                else if(fill==AirEnum){
    742742                        pressure=rho_ice*gravity*pow(thickness,(double)2)/(double)2;   //icefront on an ice sheet, pressure imbalance ice vs air.
    743743                }
    744                 else if(fill==IceEnum()){ //icefront finding itself against another icefront (pressure imbalance is fully compensated, ice vs ice)
     744                else if(fill==IceEnum){ //icefront finding itself against another icefront (pressure imbalance is fully compensated, ice vs ice)
    745745                        pressure=0;
    746746                }
    747                 else if(fill==MelangeEnum()){ //icefront finding itself against another icefront (pressure imbalance is fully compensated, ice vs ice)
     747                else if(fill==MelangeEnum){ //icefront finding itself against another icefront (pressure imbalance is fully compensated, ice vs ice)
    748748                       
    749749                        if(!shelf) ISSMERROR(exprintf("%s%i%s","fill type ",fill," not supported on ice sheets yet."));
  • issm/trunk/src/c/objects/Sing.cpp

    r3529 r3567  
    193193
    194194        /*get enum type of Sing: */
    195         enum_type=SingEnum();
     195        enum_type=SingEnum;
    196196       
    197197        /*marshall enum: */
     
    348348
    349349        /*Just branch to the correct element stiffness matrix generator, according to the type of analysis we are carrying out: */
    350         if ((analysis_type==DiagnosticAnalysisEnum()) && (sub_analysis_type==HutterAnalysisEnum())){
     350        if ((analysis_type==DiagnosticAnalysisEnum) && (sub_analysis_type==HutterAnalysisEnum)){
    351351
    352352                CreateKMatrixDiagnosticHutter( Kgg,inputs,analysis_type,sub_analysis_type);
     
    385385       
    386386        /*Just branch to the correct load generator, according to the type of analysis we are carrying out: */
    387         if ((analysis_type==DiagnosticAnalysisEnum()) && (sub_analysis_type==HutterAnalysisEnum())){
     387        if ((analysis_type==DiagnosticAnalysisEnum) && (sub_analysis_type==HutterAnalysisEnum)){
    388388       
    389389                        CreatePVectorDiagnosticHutter( pg,inputs,analysis_type,sub_analysis_type);
     
    472472int Sing::Enum(void){
    473473
    474         return SingEnum();
     474        return SingEnum;
    475475
    476476}
  • issm/trunk/src/c/objects/Spc.cpp

    r3463 r3567  
    5353
    5454        /*get enum type of Spc: */
    55         enum_type=SpcEnum();
     55        enum_type=SpcEnum;
    5656       
    5757        /*marshall enum: */
     
    130130int Spc::Enum(void){
    131131
    132         return SpcEnum();
     132        return SpcEnum;
    133133
    134134}
  • issm/trunk/src/c/objects/Tria.cpp

    r3556 r3567  
    7777        /*hooks: */
    7878        //go recover node ids, needed to initialize the node hook.
    79         if (iomodel->analysis_type==Prognostic2AnalysisEnum()){
     79        if (iomodel->analysis_type==Prognostic2AnalysisEnum){
    8080                /*Discontinuous Galerkin*/
    8181                tria_node_ids[0]=3*i+1;
     
    209209
    210210        /*get enum type of Tria: */
    211         enum_type=TriaEnum();
     211        enum_type=TriaEnum;
    212212
    213213        /*marshall enum: */
     
    569569
    570570        /*Just branch to the correct element stiffness matrix generator, according to the type of analysis we are carrying out: */
    571         if (analysis_type==ControlAnalysisEnum()){
     571        if (analysis_type==ControlAnalysisEnum){
    572572               
    573573                CreateKMatrixDiagnosticHoriz( Kgg,inputs,analysis_type,sub_analysis_type);
    574574        }
    575         else if (analysis_type==DiagnosticAnalysisEnum()){
     575        else if (analysis_type==DiagnosticAnalysisEnum){
    576576       
    577                 if (sub_analysis_type==HorizAnalysisEnum()){
     577                if (sub_analysis_type==HorizAnalysisEnum){
    578578
    579579                        CreateKMatrixDiagnosticHoriz( Kgg,inputs,analysis_type,sub_analysis_type);
     
    582582
    583583        }
    584         else if (analysis_type==SlopecomputeAnalysisEnum()){
     584        else if (analysis_type==SlopecomputeAnalysisEnum){
    585585
    586586                CreateKMatrixSlopeCompute( Kgg,inputs,analysis_type,sub_analysis_type);
    587587        }
    588         else if (analysis_type==PrognosticAnalysisEnum()){
     588        else if (analysis_type==PrognosticAnalysisEnum){
    589589
    590590                CreateKMatrixPrognostic( Kgg,inputs,analysis_type,sub_analysis_type);
    591591        }
    592         else if (analysis_type==Prognostic2AnalysisEnum()){
     592        else if (analysis_type==Prognostic2AnalysisEnum){
    593593
    594594                CreateKMatrixPrognostic2(Kgg,inputs,analysis_type,sub_analysis_type);
    595595        }
    596         else if (analysis_type==BalancedthicknessAnalysisEnum()){
     596        else if (analysis_type==BalancedthicknessAnalysisEnum){
    597597
    598598                CreateKMatrixBalancedthickness( Kgg,inputs,analysis_type,sub_analysis_type);
    599599        }
    600         else if (analysis_type==BalancedvelocitiesAnalysisEnum()){
     600        else if (analysis_type==BalancedvelocitiesAnalysisEnum){
    601601
    602602                CreateKMatrixBalancedvelocities( Kgg,inputs,analysis_type,sub_analysis_type);
     
    21772177       
    21782178        /*Just branch to the correct load generator, according to the type of analysis we are carrying out: */
    2179         if (analysis_type==ControlAnalysisEnum()){
     2179        if (analysis_type==ControlAnalysisEnum){
    21802180               
    21812181                CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type);
    21822182       
    21832183        }
    2184         else if (analysis_type==DiagnosticAnalysisEnum()){
    2185                 if (sub_analysis_type==HorizAnalysisEnum()){
     2184        else if (analysis_type==DiagnosticAnalysisEnum){
     2185                if (sub_analysis_type==HorizAnalysisEnum){
    21862186               
    21872187                        CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type);
     
    21892189                else ISSMERROR(exprintf("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet"));
    21902190        }
    2191         else if (analysis_type==SlopecomputeAnalysisEnum()){
     2191        else if (analysis_type==SlopecomputeAnalysisEnum){
    21922192               
    21932193                CreatePVectorSlopeCompute( pg,inputs,analysis_type,sub_analysis_type);
    21942194        }
    2195         else if (analysis_type==PrognosticAnalysisEnum()){
     2195        else if (analysis_type==PrognosticAnalysisEnum){
    21962196
    21972197                CreatePVectorPrognostic( pg,inputs,analysis_type,sub_analysis_type);
    21982198        }
    2199         else if (analysis_type==Prognostic2AnalysisEnum()){
     2199        else if (analysis_type==Prognostic2AnalysisEnum){
    22002200
    22012201                CreatePVectorPrognostic2( pg,inputs,analysis_type,sub_analysis_type);
    22022202        }
    2203         else if (analysis_type==BalancedthicknessAnalysisEnum()){
     2203        else if (analysis_type==BalancedthicknessAnalysisEnum){
    22042204
    22052205                CreatePVectorBalancedthickness( pg,inputs,analysis_type,sub_analysis_type);
    22062206        }
    2207         else if (analysis_type==BalancedvelocitiesAnalysisEnum()){
     2207        else if (analysis_type==BalancedvelocitiesAnalysisEnum){
    22082208
    22092209                CreatePVectorBalancedvelocities( pg,inputs,analysis_type,sub_analysis_type);
     
    30043004        for(i=0;i<numdof;i++) pe_g[i]=0.0;
    30053005
    3006         if ( (sub_analysis_type==SurfaceXAnalysisEnum()) || (sub_analysis_type==SurfaceYAnalysisEnum())){
     3006        if ( (sub_analysis_type==SurfaceXAnalysisEnum) || (sub_analysis_type==SurfaceYAnalysisEnum)){
    30073007                for(i=0;i<numdof;i++) param[i]=this->properties.s[i];
    30083008        }
    3009         if ( (sub_analysis_type==BedXAnalysisEnum()) || (sub_analysis_type==BedYAnalysisEnum())){
     3009        if ( (sub_analysis_type==BedXAnalysisEnum) || (sub_analysis_type==BedYAnalysisEnum)){
    30103010                for(i=0;i<numdof;i++) param[i]=this->properties.b[i];
    30113011        }
     
    30323032
    30333033                /*Build pe_g_gaussian vector: */
    3034                 if ( (sub_analysis_type==SurfaceXAnalysisEnum()) || (sub_analysis_type==BedXAnalysisEnum())){
     3034                if ( (sub_analysis_type==SurfaceXAnalysisEnum) || (sub_analysis_type==BedXAnalysisEnum)){
    30353035                        for(i=0;i<numdof;i++) pe_g_gaussian[i]=Jdet*gauss_weight*slope[0]*l1l2l3[i];
    30363036                }
    3037                 if ( (sub_analysis_type==SurfaceYAnalysisEnum()) || (sub_analysis_type==BedYAnalysisEnum())){
     3037                if ( (sub_analysis_type==SurfaceYAnalysisEnum) || (sub_analysis_type==BedYAnalysisEnum)){
    30383038                        for(i=0;i<numdof;i++) pe_g_gaussian[i]=Jdet*gauss_weight*slope[1]*l1l2l3[i];
    30393039                }
     
    36133613int Tria::Enum(void){
    36143614
    3615         return TriaEnum();
     3615        return TriaEnum;
    36163616
    36173617}
  • issm/trunk/src/c/objects/Vertex.cpp

    r3479 r3567  
    116116int Vertex::Enum(void){
    117117
    118         return VertexEnum();
     118        return VertexEnum;
    119119
    120120}
     
    138138
    139139        /*get enum type of Vertex: */
    140         enum_type=VertexEnum();
     140        enum_type=VertexEnum;
    141141
    142142        /*marshall enum: */
  • issm/trunk/src/c/parallel/ControlInitialization.cpp

    r3529 r3567  
    6060
    6161        /*recover fem models: */
    62         fem_dh=model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum());
    63         fem_dv=model->GetFormulation(DiagnosticAnalysisEnum(),VertAnalysisEnum());
    64         fem_ds=model->GetFormulation(DiagnosticAnalysisEnum(),StokesAnalysisEnum());
    65         fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum(),HutterAnalysisEnum());
    66         fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum());
     62        fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum);
     63        fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum);
     64        fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum);
     65        fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum);
     66        fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum);
    6767
    6868        //specific parameters for specific models
     
    8181        //compute slopes
    8282        if(verbose)_printf_("%s\n","computing bed slope (x and y derivatives)...");
    83         diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum(),BedXAnalysisEnum());
    84         diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum(),BedYAnalysisEnum());
     83        diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum,BedXAnalysisEnum);
     84        diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum,BedYAnalysisEnum);
    8585        FieldExtrudex( slopex, fem_sl->elements,fem_sl->nodes,fem_sl->vertices,fem_sl->loads,fem_sl->materials,fem_sl->parameters,"slopex",0);
    8686        FieldExtrudex( slopey, fem_sl->elements,fem_sl->nodes,fem_sl->vertices,fem_sl->loads,fem_sl->materials,fem_sl->parameters,"slopey",0);
     
    9595        //horizontal velocity
    9696        if(verbose)_printf_("%s\n"," computing horizontal velocities...");
    97         diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_dh,inputs,DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     97        diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_dh,inputs,DiagnosticAnalysisEnum,HorizAnalysisEnum);
    9898        if(verbose)_printf_("%s\n"," extruding horizontal velocities...");
    9999        VecDuplicatePatch(&ug_horiz,ug); FieldExtrudex( ug_horiz,fem_dh->elements,fem_dh->nodes, fem_dh->vertices,fem_dh->loads,fem_dh-> materials,fem_dh->parameters,"velocity",1);
     
    102102        if(verbose)_printf_("%s\n"," computing vertical velocities...");
    103103        inputs->Add("velocity",ug_horiz,numberofdofspernode_dh,numberofnodes);
    104         diagnostic_core_linear(&ug_vert,fem_dv,inputs,DiagnosticAnalysisEnum(),VertAnalysisEnum());
     104        diagnostic_core_linear(&ug_vert,fem_dv,inputs,DiagnosticAnalysisEnum,VertAnalysisEnum);
    105105
    106106        //Create 3d u_g
     
    113113        //Create 4d u_g
    114114        if(verbose)_printf_("%s\n"," computing pressure according to Pattyn...");
    115         ComputePressurex( &pg,fem_dh->elements, fem_dh->nodes, fem_dh->vertices,fem_dh->loads,  fem_dh->materials, fem_dh->parameters,inputs,DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     115        ComputePressurex( &pg,fem_dh->elements, fem_dh->nodes, fem_dh->vertices,fem_dh->loads,  fem_dh->materials, fem_dh->parameters,inputs,DiagnosticAnalysisEnum,HorizAnalysisEnum);
    116116        VecScale(pg,1.0/stokesreconditioning);
    117117        ug_stokes=NewVec(fem_ds->nodesets->GetGSize());
     
    133133        if(verbose)_printf_("%s\n"," computing stokes velocities and pressure ...");
    134134        VecFree(&ug);
    135         diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_ds,inputs,DiagnosticAnalysisEnum(),StokesAnalysisEnum());
     135        diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_ds,inputs,DiagnosticAnalysisEnum,StokesAnalysisEnum);
    136136
    137137        //Add in inputs
  • issm/trunk/src/c/parallel/ControlRestart.cpp

    r3332 r3567  
    3434        for(i=0;i<numberofnodes;i++) param_g_copy[i]=param_g[i];
    3535
    36         temporary_results=new DataSet(ResultsEnum());
     36        temporary_results=new DataSet(ResultsEnum);
    3737        result=new Result(temporary_results->Size()+1,0,1,"param_g",param_g_copy,numberofnodes);
    3838        temporary_results->AddObject(result);
     
    4242
    4343        //process results
    44         ProcessResults(&results,temporary_results,model,ControlAnalysisEnum());
     44        ProcessResults(&results,temporary_results,model,ControlAnalysisEnum);
    4545
    4646        //Write results on disk
  • issm/trunk/src/c/parallel/ProcessResults.cpp

    r3556 r3567  
    109109
    110110        /*Initialize new results: */
    111         newresults=new DataSet(ResultsEnum());
     111        newresults=new DataSet(ResultsEnum);
    112112       
    113113        /*some flags needed: */
     
    118118
    119119        /*Recover femmodels first: */
    120         fem_dh=model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     120        fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum);
    121121        fem_c=fem_dh;
    122         fem_dv=model->GetFormulation(DiagnosticAnalysisEnum(),VertAnalysisEnum());
    123         fem_ds=model->GetFormulation(DiagnosticAnalysisEnum(),StokesAnalysisEnum());
    124         fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum(),HutterAnalysisEnum());
    125         fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum());
    126         if(analysis_type==PrognosticAnalysisEnum()){
    127                 fem_p=model->GetFormulation(PrognosticAnalysisEnum());
     122        fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum);
     123        fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum);
     124        fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum);
     125        fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum);
     126        if(analysis_type==PrognosticAnalysisEnum){
     127                fem_p=model->GetFormulation(PrognosticAnalysisEnum);
    128128        }
    129         if(analysis_type==Prognostic2AnalysisEnum()){
    130                 fem_p=model->GetFormulation(Prognostic2AnalysisEnum());
     129        if(analysis_type==Prognostic2AnalysisEnum){
     130                fem_p=model->GetFormulation(Prognostic2AnalysisEnum);
    131131        }
    132         if(analysis_type==TransientAnalysisEnum()){
    133                 fem_p=model->GetFormulation(PrognosticAnalysisEnum());
     132        if(analysis_type==TransientAnalysisEnum){
     133                fem_p=model->GetFormulation(PrognosticAnalysisEnum);
    134134        }
    135         if(analysis_type==BalancedthicknessAnalysisEnum()){
    136                 fem_p=model->GetFormulation(BalancedthicknessAnalysisEnum());
     135        if(analysis_type==BalancedthicknessAnalysisEnum){
     136                fem_p=model->GetFormulation(BalancedthicknessAnalysisEnum);
    137137        }
    138         if(analysis_type==BalancedvelocitiesAnalysisEnum()){
    139                 fem_p=model->GetFormulation(BalancedvelocitiesAnalysisEnum());
     138        if(analysis_type==BalancedvelocitiesAnalysisEnum){
     139                fem_p=model->GetFormulation(BalancedvelocitiesAnalysisEnum);
    140140        }
    141         fem_t=model->GetFormulation(ThermalAnalysisEnum());
     141        fem_t=model->GetFormulation(ThermalAnalysisEnum);
    142142
    143143        for(n=0;n<results->Size();n++){
  • issm/trunk/src/c/parallel/balancedthickness.cpp

    r3332 r3567  
    7676
    7777        _printf_("read and create finite element model:\n");
    78         model->AddFormulation(fid,BalancedthicknessAnalysisEnum());
     78        model->AddFormulation(fid,BalancedthicknessAnalysisEnum);
    7979
    8080        /*recover parameters: */
     
    8484        _printf_("initialize inputs:\n");
    8585       
    86         model->FindParam(&u_g_serial,NULL,NULL,"u_g",BalancedthicknessAnalysisEnum());
    87         model->FindParam(&melting_g,NULL,NULL,"m_g",BalancedthicknessAnalysisEnum());
    88         model->FindParam(&accumulation_g,NULL,NULL,"a_g",BalancedthicknessAnalysisEnum());
     86        model->FindParam(&u_g_serial,NULL,NULL,"u_g",BalancedthicknessAnalysisEnum);
     87        model->FindParam(&melting_g,NULL,NULL,"m_g",BalancedthicknessAnalysisEnum);
     88        model->FindParam(&accumulation_g,NULL,NULL,"a_g",BalancedthicknessAnalysisEnum);
    8989        model->FindParam(&numberofnodes,"numberofnodes");
    9090       
     
    9595
    9696        _printf_("initialize results:\n");
    97         results=new DataSet(ResultsEnum());
     97        results=new DataSet(ResultsEnum);
    9898        MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime();
    9999
     
    115115                #ifdef _HAVE_DAKOTA_
    116116                MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
    117                 Qmux(model,inputs,BalancedthicknessAnalysisEnum(),NoneAnalysisEnum());
     117                Qmux(model,inputs,BalancedthicknessAnalysisEnum,NoneAnalysisEnum);
    118118                MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
    119119                #else
     
    127127       
    128128        _printf_("process results:\n");
    129         ProcessResults(&processedresults,results,model,BalancedthicknessAnalysisEnum());
     129        ProcessResults(&processedresults,results,model,BalancedthicknessAnalysisEnum);
    130130       
    131131        _printf_("write results to disk:\n");
  • issm/trunk/src/c/parallel/balancedthickness_core.cpp

    r3446 r3567  
    3535
    3636        /*recover fem model: */
    37         fem_p=model->GetFormulation(BalancedthicknessAnalysisEnum());
     37        fem_p=model->GetFormulation(BalancedthicknessAnalysisEnum);
    3838
    3939        //first recover parameters common to all solutions
     
    4848       
    4949        _printf_("call computational core:\n");
    50         diagnostic_core_linear(&h_g,fem_p,inputs,BalancedthicknessAnalysisEnum(),NoneAnalysisEnum());
     50        diagnostic_core_linear(&h_g,fem_p,inputs,BalancedthicknessAnalysisEnum,NoneAnalysisEnum);
    5151
    5252        _printf_("extrude computed thickness on all layers:\n");
  • issm/trunk/src/c/parallel/balancedvelocities.cpp

    r3332 r3567  
    7373
    7474        _printf_("read and create finite element model:\n");
    75         model->AddFormulation(fid,BalancedvelocitiesAnalysisEnum());
     75        model->AddFormulation(fid,BalancedvelocitiesAnalysisEnum);
    7676
    7777        /*recover parameters: */
     
    8181        _printf_("initialize inputs:\n");
    8282       
    83         model->FindParam(&u_g_serial,NULL,NULL,"u_g",BalancedvelocitiesAnalysisEnum());
    84         model->FindParam(&melting_g,NULL,NULL,"m_g",BalancedvelocitiesAnalysisEnum());
    85         model->FindParam(&accumulation_g,NULL,NULL,"a_g",BalancedvelocitiesAnalysisEnum());
     83        model->FindParam(&u_g_serial,NULL,NULL,"u_g",BalancedvelocitiesAnalysisEnum);
     84        model->FindParam(&melting_g,NULL,NULL,"m_g",BalancedvelocitiesAnalysisEnum);
     85        model->FindParam(&accumulation_g,NULL,NULL,"a_g",BalancedvelocitiesAnalysisEnum);
    8686        model->FindParam(&numberofnodes,"numberofnodes");
    8787       
     
    9292
    9393        _printf_("initialize results:\n");
    94         results=new DataSet(ResultsEnum());
     94        results=new DataSet(ResultsEnum);
    9595        MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime();
    9696
     
    112112                #ifdef _HAVE_DAKOTA_
    113113                MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
    114                 Qmux(model,inputs,BalancedvelocitiesAnalysisEnum(),NoneAnalysisEnum());
     114                Qmux(model,inputs,BalancedvelocitiesAnalysisEnum,NoneAnalysisEnum);
    115115                MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
    116116                #else
     
    124124       
    125125        _printf_("process results:\n");
    126         ProcessResults(&processedresults,results,model,BalancedvelocitiesAnalysisEnum());
     126        ProcessResults(&processedresults,results,model,BalancedvelocitiesAnalysisEnum);
    127127       
    128128        _printf_("write results to disk:\n");
  • issm/trunk/src/c/parallel/balancedvelocities_core.cpp

    r3446 r3567  
    3434
    3535        /*recover fem model: */
    36         fem_p=model->GetFormulation(BalancedvelocitiesAnalysisEnum());
     36        fem_p=model->GetFormulation(BalancedvelocitiesAnalysisEnum);
    3737
    3838        //first recover parameters common to all solutions
     
    4747       
    4848        _printf_("call computational core:\n");
    49         diagnostic_core_linear(&v_g,fem_p,inputs,BalancedvelocitiesAnalysisEnum(),NoneAnalysisEnum());
     49        diagnostic_core_linear(&v_g,fem_p,inputs,BalancedvelocitiesAnalysisEnum,NoneAnalysisEnum);
    5050
    5151        _printf_("extrude computed thickness on all layers:\n");
  • issm/trunk/src/c/parallel/control_core.cpp

    r3446 r3567  
    9797                 * distribution: */
    9898                if (control_steady){
    99                         steadystate_results=new DataSet(ResultsEnum());
     99                        steadystate_results=new DataSet(ResultsEnum);
    100100                        steadystate_core(steadystate_results,model,inputs);
    101101                        VecFree(&t_g); steadystate_results->FindResult(&t_g,"t_g");
     
    108108
    109109                _printf_("%s\n","      computing gradJ...");
    110                 gradjcompute_results=new DataSet(ResultsEnum());
     110                gradjcompute_results=new DataSet(ResultsEnum);
    111111                gradjcompute_core(gradjcompute_results,model, inputs);
    112112                gradjcompute_results->FindResult(&grad_g,"grad_g");
     
    205205        if (control_steady){
    206206                inputs->Add(control_type,param_g,1,numberofnodes);
    207                 steadystate_results=new DataSet(ResultsEnum());
     207                steadystate_results=new DataSet(ResultsEnum);
    208208                steadystate_core(steadystate_results,model,inputs);
    209209
     
    216216        else{
    217217                inputs->Add(control_type,param_g,1,numberofnodes);
    218                 diagnostic_results=new DataSet(ResultsEnum());
     218                diagnostic_results=new DataSet(ResultsEnum);
    219219                diagnostic_core(diagnostic_results,model, inputs);
    220220
  • issm/trunk/src/c/parallel/diagnostic.cpp

    r3332 r3567  
    7676        _printf_("read and create finite element model:\n");
    7777        _printf_("\n   reading diagnostic horiz model data:\n");
    78         model->AddFormulation(fid,DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     78        model->AddFormulation(fid,DiagnosticAnalysisEnum,HorizAnalysisEnum);
    7979
    8080        _printf_("\n   reading diagnostic vert model data:\n");
    81         model->AddFormulation(fid,DiagnosticAnalysisEnum(),VertAnalysisEnum());
     81        model->AddFormulation(fid,DiagnosticAnalysisEnum,VertAnalysisEnum);
    8282       
    8383        _printf_("\n   reading diagnostic stokes model data:\n");
    84         model->AddFormulation(fid,DiagnosticAnalysisEnum(),StokesAnalysisEnum());
     84        model->AddFormulation(fid,DiagnosticAnalysisEnum,StokesAnalysisEnum);
    8585       
    8686        _printf_("\n   reading diagnostic hutter model data:\n");
    87         model->AddFormulation(fid,DiagnosticAnalysisEnum(),HutterAnalysisEnum());
     87        model->AddFormulation(fid,DiagnosticAnalysisEnum,HutterAnalysisEnum);
    8888       
    8989        _printf_("\n   reading surface and bed slope computation model data:\n");
    90         model->AddFormulation(fid,SlopecomputeAnalysisEnum());
     90        model->AddFormulation(fid,SlopecomputeAnalysisEnum);
    9191
    9292        /*get parameters: */
     
    9696
    9797        _printf_("initialize inputs:\n");
    98         model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     98        model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum,HorizAnalysisEnum);
    9999        model->FindParam(&numberofnodes,"numberofnodes");
    100100
     
    103103
    104104        if(control_analysis){
    105                 model->FindParam(&u_g_obs,NULL,NULL,"u_g_obs",DiagnosticAnalysisEnum(),HorizAnalysisEnum());
    106                 model->FindParam(&weights,NULL,NULL,"weights",DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     105                model->FindParam(&u_g_obs,NULL,NULL,"u_g_obs",DiagnosticAnalysisEnum,HorizAnalysisEnum);
     106                model->FindParam(&weights,NULL,NULL,"weights",DiagnosticAnalysisEnum,HorizAnalysisEnum);
    107107                inputs->Add("velocity_obs",u_g_obs,2,numberofnodes);
    108108                inputs->Add("weights",weights,1,numberofnodes);
     
    110110       
    111111        _printf_("initialize results:\n");
    112         results=new DataSet(ResultsEnum());
     112        results=new DataSet(ResultsEnum);
    113113        MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime();
    114114
     
    127127                       
    128128                        _printf_("process results:\n");
    129                         ProcessResults(&processed_results,results,model,DiagnosticAnalysisEnum());
     129                        ProcessResults(&processed_results,results,model,DiagnosticAnalysisEnum);
    130130                }
    131131                else{
     
    145145
    146146                        _printf_("process results:\n");
    147                         ProcessResults(&processed_results,results,model,ControlAnalysisEnum());
     147                        ProcessResults(&processed_results,results,model,ControlAnalysisEnum);
    148148                }
    149149
     
    157157                #ifdef _HAVE_DAKOTA_
    158158                MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
    159                 Qmux(model,inputs,DiagnosticAnalysisEnum(),NoneAnalysisEnum());
     159                Qmux(model,inputs,DiagnosticAnalysisEnum,NoneAnalysisEnum);
    160160                MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
    161161                #else
  • issm/trunk/src/c/parallel/diagnostic_core.cpp

    r3529 r3567  
    7070
    7171        /*recover fem models: */
    72         fem_dh=model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum());
    73         fem_dv=model->GetFormulation(DiagnosticAnalysisEnum(),VertAnalysisEnum());
    74         fem_ds=model->GetFormulation(DiagnosticAnalysisEnum(),StokesAnalysisEnum());
    75         fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum(),HutterAnalysisEnum());
    76         fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum());
     72        fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum);
     73        fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum);
     74        fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum);
     75        fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum);
     76        fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum);
    7777
    7878        //specific parameters for specific models
     
    8383        //for qmu analysis, be sure the velocity input we are starting from  is the one in the parameters: */
    8484        if(qmu_analysis){
    85                 model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     85                model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum,HorizAnalysisEnum);
    8686                inputs->Add("velocity",u_g_initial,3,numberofnodes);
    8787        }
     
    9090                       
    9191                if(verbose)_printf_("%s\n","computing surface slope (x and y derivatives)...");
    92                 diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum(),SurfaceXAnalysisEnum());
    93                 diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum(),SurfaceYAnalysisEnum());
     92                diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum,SurfaceXAnalysisEnum);
     93                diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum,SurfaceYAnalysisEnum);
    9494
    9595                if (dim==3){
     
    106106
    107107                if(verbose)_printf_("%s\n"," computing hutter velocities...");
    108                 diagnostic_core_linear(&ug,fem_dhu,inputs,DiagnosticAnalysisEnum(),HutterAnalysisEnum());
     108                diagnostic_core_linear(&ug,fem_dhu,inputs,DiagnosticAnalysisEnum,HutterAnalysisEnum);
    109109
    110110                if(verbose)_printf_("%s\n"," computing pressure according to MacAyeal...");
    111                 ComputePressurex(&pg,fem_dhu->elements,fem_dhu->nodes, fem_dhu->vertices,fem_dhu->loads,fem_dhu->materials,fem_dhu->parameters,inputs,DiagnosticAnalysisEnum(),HutterAnalysisEnum());
     111                ComputePressurex(&pg,fem_dhu->elements,fem_dhu->nodes, fem_dhu->vertices,fem_dhu->loads,fem_dhu->materials,fem_dhu->parameters,inputs,DiagnosticAnalysisEnum,HutterAnalysisEnum);
    112112
    113113                if(verbose)_printf_("%s\n"," update boundary conditions for macyeal pattyn using hutter results...");
     
    123123               
    124124                if(verbose)_printf_("%s\n"," computing horizontal velocities...");
    125                 diagnostic_core_nonlinear(&ug,NULL,NULL,fem_dh->loads,fem_dh,inputs,DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     125                diagnostic_core_nonlinear(&ug,NULL,NULL,fem_dh->loads,fem_dh,inputs,DiagnosticAnalysisEnum,HorizAnalysisEnum);
    126126
    127127                if(dim==2){
    128128                        if(verbose)_printf_("%s\n"," computing pressure according to MacAyeal...");
    129                         ComputePressurex(&pg,fem_dh->elements,fem_dh->nodes, fem_dh->vertices,fem_dh->loads,fem_dh->materials,fem_dh->parameters,inputs,DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     129                        ComputePressurex(&pg,fem_dh->elements,fem_dh->nodes, fem_dh->vertices,fem_dh->loads,fem_dh->materials,fem_dh->parameters,inputs,DiagnosticAnalysisEnum,HorizAnalysisEnum);
    130130                }
    131131
     
    140140                if(verbose)_printf_("%s\n"," computing vertical velocities...");
    141141                inputs->Add("velocity",ug_horiz,numberofdofspernode_dh,numberofnodes);
    142                 diagnostic_core_linear(&ug_vert,fem_dv,inputs,DiagnosticAnalysisEnum(),VertAnalysisEnum());
     142                diagnostic_core_linear(&ug_vert,fem_dv,inputs,DiagnosticAnalysisEnum,VertAnalysisEnum);
    143143
    144144                if(verbose)_printf_("%s\n"," combining horizontal and vertical velocities...");
     
    150150
    151151                if(verbose)_printf_("%s\n"," computing pressure according to Pattyn...");
    152                 ComputePressurex(&pg,fem_dh->elements, fem_dh->nodes, fem_dh->vertices,fem_dh->loads,  fem_dh->materials,fem_dh->parameters,inputs,DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     152                ComputePressurex(&pg,fem_dh->elements, fem_dh->nodes, fem_dh->vertices,fem_dh->loads,  fem_dh->materials,fem_dh->parameters,inputs,DiagnosticAnalysisEnum,HorizAnalysisEnum);
    153153               
    154154                if (isstokes){
     
    158158
    159159                        if(verbose)_printf_("%s\n","computing bed slope (x and y derivatives)...");
    160                         diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum(),BedXAnalysisEnum());
    161                         diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum(),BedYAnalysisEnum());
     160                        diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum,BedXAnalysisEnum);
     161                        diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum,BedYAnalysisEnum);
    162162                        FieldExtrudex( slopex, fem_sl->elements,fem_sl->nodes,fem_sl->vertices,fem_sl->loads,fem_sl->materials,fem_sl->parameters,"slopex",0);
    163163                        FieldExtrudex( slopey, fem_sl->elements,fem_sl->nodes,fem_sl->vertices,fem_sl->loads,fem_sl->materials,fem_sl->parameters,"slopey",0);
     
    182182                        if(verbose)_printf_("%s\n"," computing stokes velocities and pressure ...");
    183183                        VecFree(&ug);
    184                         diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_ds,inputs,DiagnosticAnalysisEnum(),StokesAnalysisEnum());
     184                        diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_ds,inputs,DiagnosticAnalysisEnum,StokesAnalysisEnum);
    185185               
    186186                        //decondition" pressure
  • issm/trunk/src/c/parallel/gradjcompute_core.cpp

    r3446 r3567  
    6262
    6363        _printf_("%s\n","      recover solution for this stiffness and right hand side:");
    64         diagnostic_core_nonlinear(&u_g,&K_ff0,&K_fs0,NULL, femmodel,inputs,DiagnosticAnalysisEnum(),sub_analysis_type);
     64        diagnostic_core_nonlinear(&u_g,&K_ff0,&K_fs0,NULL, femmodel,inputs,DiagnosticAnalysisEnum,sub_analysis_type);
    6565        VecToMPISerial(&u_g_double,u_g); VecFree(&u_g);
    6666        inputs->Add("velocity",u_g_double,numberofdofspernode,numberofnodes);
     
    9595
    9696        if(control_steady){
    97                 diagnostic_results=new DataSet(ResultsEnum());
     97                diagnostic_results=new DataSet(ResultsEnum);
    9898                diagnostic_core(diagnostic_results,model, inputs);
    9999
  • issm/trunk/src/c/parallel/objectivefunctionC.cpp

    r3563 r3567  
    9090        /*Run diagnostic with updated parameters.*/
    9191        if(!control_steady){
    92                 diagnostic_core_nonlinear(&u_g,NULL,NULL,NULL,femmodel,inputs,DiagnosticAnalysisEnum(),sub_analysis_type);
     92                diagnostic_core_nonlinear(&u_g,NULL,NULL,NULL,femmodel,inputs,DiagnosticAnalysisEnum,sub_analysis_type);
    9393                VecToMPISerial(&u_g_double,u_g); VecFree(&u_g);
    9494                inputs->Add("velocity",u_g_double,numberofdofspernode,numberofnodes);
     
    9696        else{
    9797                //We need a 3D velocity!! (vz is required for the next thermal run)
    98                 diagnostic_results=new DataSet(ResultsEnum());
     98                diagnostic_results=new DataSet(ResultsEnum);
    9999                diagnostic_core(diagnostic_results,model, inputs);
    100100
  • issm/trunk/src/c/parallel/prognostic.cpp

    r3381 r3567  
    7474
    7575        _printf_("read and create finite element model:\n");
    76         model->AddFormulation(fid,PrognosticAnalysisEnum());
     76        model->AddFormulation(fid,PrognosticAnalysisEnum);
    7777
    7878        /*recover parameters: */
     
    8282        _printf_("initialize inputs:\n");
    8383       
    84         model->FindParam(&u_g_serial,NULL,NULL,"u_g",PrognosticAnalysisEnum());
    85         model->FindParam(&h_g_initial,NULL,NULL,"h_g",PrognosticAnalysisEnum());
    86         model->FindParam(&melting_g,NULL,NULL,"m_g",PrognosticAnalysisEnum());
    87         model->FindParam(&accumulation_g,NULL,NULL,"a_g",PrognosticAnalysisEnum());
     84        model->FindParam(&u_g_serial,NULL,NULL,"u_g",PrognosticAnalysisEnum);
     85        model->FindParam(&h_g_initial,NULL,NULL,"h_g",PrognosticAnalysisEnum);
     86        model->FindParam(&melting_g,NULL,NULL,"m_g",PrognosticAnalysisEnum);
     87        model->FindParam(&accumulation_g,NULL,NULL,"a_g",PrognosticAnalysisEnum);
    8888        model->FindParam(&dt,"dt");
    8989        model->FindParam(&yts,"yts");
     
    9898
    9999        _printf_("initialize results:\n");
    100         results=new DataSet(ResultsEnum());
     100        results=new DataSet(ResultsEnum);
    101101        MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime();
    102102
     
    118118                #ifdef _HAVE_DAKOTA_
    119119                MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
    120                 Qmux(model,inputs,PrognosticAnalysisEnum(),NoneAnalysisEnum());
     120                Qmux(model,inputs,PrognosticAnalysisEnum,NoneAnalysisEnum);
    121121                MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
    122122                #else
     
    130130       
    131131        _printf_("process results:\n");
    132         ProcessResults(&processedresults,results,model,PrognosticAnalysisEnum());
     132        ProcessResults(&processedresults,results,model,PrognosticAnalysisEnum);
    133133       
    134134        _printf_("write results to disk:\n");
  • issm/trunk/src/c/parallel/prognostic2.cpp

    r3381 r3567  
    7474
    7575        _printf_("read and create finite element model:\n");
    76         model->AddFormulation(fid,Prognostic2AnalysisEnum());
     76        model->AddFormulation(fid,Prognostic2AnalysisEnum);
    7777
    7878        /*recover parameters: */
     
    8282        _printf_("initialize inputs:\n");
    8383       
    84         model->FindParam(&vx_g,NULL,NULL,"vx_g",Prognostic2AnalysisEnum());
    85         model->FindParam(&vy_g,NULL,NULL,"vy_g",Prognostic2AnalysisEnum());
    86         model->FindParam(&h_g_initial,NULL,NULL,"h_g",Prognostic2AnalysisEnum());
    87         model->FindParam(&m_g,NULL,NULL,"m_g",Prognostic2AnalysisEnum());
    88         model->FindParam(&a_g,NULL,NULL,"a_g",Prognostic2AnalysisEnum());
     84        model->FindParam(&vx_g,NULL,NULL,"vx_g",Prognostic2AnalysisEnum);
     85        model->FindParam(&vy_g,NULL,NULL,"vy_g",Prognostic2AnalysisEnum);
     86        model->FindParam(&h_g_initial,NULL,NULL,"h_g",Prognostic2AnalysisEnum);
     87        model->FindParam(&m_g,NULL,NULL,"m_g",Prognostic2AnalysisEnum);
     88        model->FindParam(&a_g,NULL,NULL,"a_g",Prognostic2AnalysisEnum);
    8989        model->FindParam(&dt,"dt");
    9090        model->FindParam(&yts,"yts");
     
    100100
    101101        _printf_("initialize results:\n");
    102         results=new DataSet(ResultsEnum());
     102        results=new DataSet(ResultsEnum);
    103103        MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime();
    104104
     
    120120                #ifdef _HAVE_DAKOTA_
    121121                MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
    122                 Qmux(model,inputs,Prognostic2AnalysisEnum(),NoneAnalysisEnum());
     122                Qmux(model,inputs,Prognostic2AnalysisEnum,NoneAnalysisEnum);
    123123                MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
    124124                #else
     
    132132       
    133133        _printf_("process results:\n");
    134         ProcessResults(&processedresults,results,model,Prognostic2AnalysisEnum());
     134        ProcessResults(&processedresults,results,model,Prognostic2AnalysisEnum);
    135135       
    136136        _printf_("write results to disk:\n");
  • issm/trunk/src/c/parallel/prognostic2_core.cpp

    r3556 r3567  
    3535
    3636        /*recover fem model: */
    37         fem_p=model->GetFormulation(Prognostic2AnalysisEnum());
     37        fem_p=model->GetFormulation(Prognostic2AnalysisEnum);
    3838
    3939        //first recover parameters common to all solutions
     
    5454       
    5555        _printf_("call computational core:\n");
    56         diagnostic_core_linear(&h_g,fem_p,inputs,Prognostic2AnalysisEnum(),NoneAnalysisEnum());
     56        diagnostic_core_linear(&h_g,fem_p,inputs,Prognostic2AnalysisEnum,NoneAnalysisEnum);
    5757
    5858        _printf_("Averaging over vertices:\n");
  • issm/trunk/src/c/parallel/prognostic_core.cpp

    r3446 r3567  
    3434
    3535        /*recover fem model: */
    36         fem_p=model->GetFormulation(PrognosticAnalysisEnum());
     36        fem_p=model->GetFormulation(PrognosticAnalysisEnum);
    3737
    3838        //first recover parameters common to all solutions
     
    4747       
    4848        _printf_("call computational core:\n");
    49         diagnostic_core_linear(&h_g,fem_p,inputs,PrognosticAnalysisEnum(),NoneAnalysisEnum());
     49        diagnostic_core_linear(&h_g,fem_p,inputs,PrognosticAnalysisEnum,NoneAnalysisEnum);
    5050
    5151        _printf_("extrude computed thickness on all layers:\n");
  • issm/trunk/src/c/parallel/slopecompute.cpp

    r3332 r3567  
    6868
    6969        _printf_("read and create finite element model:\n");
    70         model->AddFormulation(fid,SlopecomputeAnalysisEnum());
     70        model->AddFormulation(fid,SlopecomputeAnalysisEnum);
    7171
    7272        /*recover parameters: */
     
    8080
    8181        _printf_("initialize results:\n");
    82         results=new DataSet(ResultsEnum());
     82        results=new DataSet(ResultsEnum);
    8383        MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime();
    8484
     
    9494       
    9595        _printf_("process results:\n");
    96         ProcessResults(&processedresults,results,model,SlopecomputeAnalysisEnum());
     96        ProcessResults(&processedresults,results,model,SlopecomputeAnalysisEnum);
    9797       
    9898        _printf_("write results to disk:\n");
  • issm/trunk/src/c/parallel/slopecompute_core.cpp

    r3446 r3567  
    3131
    3232        /*recover fem model: */
    33         fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum());
     33        fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum);
    3434
    3535        //first recover parameters common to all solutions
     
    3939
    4040        _printf_("call computational core:\n");
    41         diagnostic_core_linear(&sx_g,fem_sl,inputs,SlopecomputeAnalysisEnum(),SurfaceXAnalysisEnum());
    42         diagnostic_core_linear(&sy_g,fem_sl,inputs,SlopecomputeAnalysisEnum(),SurfaceYAnalysisEnum());
     41        diagnostic_core_linear(&sx_g,fem_sl,inputs,SlopecomputeAnalysisEnum,SurfaceXAnalysisEnum);
     42        diagnostic_core_linear(&sy_g,fem_sl,inputs,SlopecomputeAnalysisEnum,SurfaceYAnalysisEnum);
    4343
    4444        _printf_("extrude computed slope on all layers:\n");
  • issm/trunk/src/c/parallel/steadystate.cpp

    r3332 r3567  
    7979       
    8080        _printf_("\n   reading diagnostic horiz model data:\n");
    81         model->AddFormulation(fid,DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     81        model->AddFormulation(fid,DiagnosticAnalysisEnum,HorizAnalysisEnum);
    8282
    8383        _printf_("\n   reading diagnostic vert model data:\n");
    84         model->AddFormulation(fid,DiagnosticAnalysisEnum(),VertAnalysisEnum());
     84        model->AddFormulation(fid,DiagnosticAnalysisEnum,VertAnalysisEnum);
    8585       
    8686        _printf_("\n   reading diagnostic stokes model data:\n");
    87         model->AddFormulation(fid,DiagnosticAnalysisEnum(),StokesAnalysisEnum());
     87        model->AddFormulation(fid,DiagnosticAnalysisEnum,StokesAnalysisEnum);
    8888       
    8989        _printf_("\n   reading diagnostic hutter model data:\n");
    90         model->AddFormulation(fid,DiagnosticAnalysisEnum(),HutterAnalysisEnum());
     90        model->AddFormulation(fid,DiagnosticAnalysisEnum,HutterAnalysisEnum);
    9191       
    9292        _printf_("\n   reading surface and bed slope computation model data:\n");
    93         model->AddFormulation(fid,SlopecomputeAnalysisEnum());
     93        model->AddFormulation(fid,SlopecomputeAnalysisEnum);
    9494
    9595        _printf_("\n   read and create thermal finite element model:\n");
    96         model->AddFormulation(fid,ThermalAnalysisEnum());
     96        model->AddFormulation(fid,ThermalAnalysisEnum);
    9797        _printf_("\n   read and create melting finite element model:\n");
    98         model->AddFormulation(fid,MeltingAnalysisEnum());
     98        model->AddFormulation(fid,MeltingAnalysisEnum);
    9999
    100100        /*recover parameters: */
     
    104104
    105105        _printf_("initialize inputs:\n");
    106         model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum(),HorizAnalysisEnum());
    107         model->FindParam(&p_g_initial,NULL,NULL,"p_g",ThermalAnalysisEnum());
     106        model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum,HorizAnalysisEnum);
     107        model->FindParam(&p_g_initial,NULL,NULL,"p_g",ThermalAnalysisEnum);
    108108        model->FindParam(&dt,"dt");
    109109        model->FindParam(&numberofnodes,"numberofnodes");
     
    115115
    116116        if(control_analysis){
    117                 model->FindParam(&u_g_obs,NULL,NULL,"u_g_obs",DiagnosticAnalysisEnum(),HorizAnalysisEnum());
    118                 model->FindParam(&weights,NULL,NULL,"weights",DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     117                model->FindParam(&u_g_obs,NULL,NULL,"u_g_obs",DiagnosticAnalysisEnum,HorizAnalysisEnum);
     118                model->FindParam(&weights,NULL,NULL,"weights",DiagnosticAnalysisEnum,HorizAnalysisEnum);
    119119                inputs->Add("velocity_obs",u_g_obs,2,numberofnodes);
    120120                inputs->Add("weights",weights,1,numberofnodes);
     
    122122       
    123123        _printf_("initialize results:\n");
    124         results=new DataSet(ResultsEnum());
     124        results=new DataSet(ResultsEnum);
    125125        MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime();
    126126
     
    140140
    141141                        _printf_("process results:\n");
    142                         ProcessResults(&processed_results,results,model,SteadystateAnalysisEnum());
     142                        ProcessResults(&processed_results,results,model,SteadystateAnalysisEnum);
    143143                }
    144144                else{
    145145                        /*change control_steady to 1 to know we are doing steadystate*/
    146                         fem_dh=model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum());
    147                         fem_ds=model->GetFormulation(DiagnosticAnalysisEnum(),StokesAnalysisEnum());
     146                        fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum);
     147                        fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum);
    148148                        param=(Param*)fem_dh->parameters->FindParamObject("control_steady");
    149149                        param->SetDouble(1);
     
    165165
    166166                        _printf_("process results:\n");
    167                         ProcessResults(&processed_results,results,model,ControlAnalysisEnum());
     167                        ProcessResults(&processed_results,results,model,ControlAnalysisEnum);
    168168                }
    169169
     
    177177                #ifdef _HAVE_DAKOTA_
    178178                MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
    179                 Qmux(model,inputs,SteadystateAnalysisEnum(),NoneAnalysisEnum());
     179                Qmux(model,inputs,SteadystateAnalysisEnum,NoneAnalysisEnum);
    180180                MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
    181181                #else
  • issm/trunk/src/c/parallel/steadystate_core.cpp

    r3446 r3567  
    5151
    5252        /*recover fem models: */
    53         fem_dh=model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum());
    54         fem_dv=model->GetFormulation(DiagnosticAnalysisEnum(),VertAnalysisEnum());
    55         fem_ds=model->GetFormulation(DiagnosticAnalysisEnum(),StokesAnalysisEnum());
    56         fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum(),HutterAnalysisEnum());
    57         fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum());
    58         fem_t=model->GetFormulation(ThermalAnalysisEnum());
    59         fem_m=model->GetFormulation(MeltingAnalysisEnum());
     53        fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum);
     54        fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum);
     55        fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum);
     56        fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum);
     57        fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum);
     58        fem_t=model->GetFormulation(ThermalAnalysisEnum);
     59        fem_m=model->GetFormulation(MeltingAnalysisEnum);
    6060
    6161
     
    8181                        inputs->Add("velocity",u_g,ndof,numberofnodes);
    8282                }
    83                 results_thermal=new DataSet(ResultsEnum());
     83                results_thermal=new DataSet(ResultsEnum);
    8484                thermal_core(results_thermal,model,inputs);
    8585       
     
    9898
    9999                //now compute diagnostic velocity using the steady state temperature.
    100                 results_diagnostic=new DataSet(ResultsEnum());
     100                results_diagnostic=new DataSet(ResultsEnum);
    101101                diagnostic_core(results_diagnostic,model, inputs);
    102102
  • issm/trunk/src/c/parallel/thermal.cpp

    r3332 r3567  
    7373
    7474        _printf_("read and create thermal finite element model:\n");
    75         model->AddFormulation(fid,ThermalAnalysisEnum());
     75        model->AddFormulation(fid,ThermalAnalysisEnum);
    7676        _printf_("read and create melting finite element model:\n");
    77         model->AddFormulation(fid,MeltingAnalysisEnum());
     77        model->AddFormulation(fid,MeltingAnalysisEnum);
    7878
    7979        /*recover parameters: */
     
    8282       
    8383        _printf_("initialize inputs:\n");
    84         model->FindParam(&u_g,NULL,NULL,"u_g",ThermalAnalysisEnum());
    85         model->FindParam(&p_g,NULL,NULL,"p_g",ThermalAnalysisEnum());
     84        model->FindParam(&u_g,NULL,NULL,"u_g",ThermalAnalysisEnum);
     85        model->FindParam(&p_g,NULL,NULL,"p_g",ThermalAnalysisEnum);
    8686        model->FindParam(&numberofnodes,"numberofnodes");
    8787        model->FindParam(&dt,"dt");
     
    9595
    9696        _printf_("initialize results:\n");
    97         results=new DataSet(ResultsEnum());
     97        results=new DataSet(ResultsEnum);
    9898        MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime();
    9999
     
    112112                       
    113113                _printf_("process results:\n");
    114                 ProcessResults(&processed_results,results,model,ThermalAnalysisEnum());
     114                ProcessResults(&processed_results,results,model,ThermalAnalysisEnum);
    115115               
    116116                _printf_("write results to disk:\n");
     
    123123                #ifdef _HAVE_DAKOTA_
    124124                MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
    125                 Qmux(model,inputs,ThermalAnalysisEnum(),NoneAnalysisEnum());
     125                Qmux(model,inputs,ThermalAnalysisEnum,NoneAnalysisEnum);
    126126                MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
    127127                #else
  • issm/trunk/src/c/parallel/thermal_core.cpp

    r3332 r3567  
    4545
    4646        /*recover fem models: */
    47         fem_t=model->GetFormulation(ThermalAnalysisEnum());
    48         fem_m=model->GetFormulation(MeltingAnalysisEnum());
     47        fem_t=model->GetFormulation(ThermalAnalysisEnum);
     48        fem_m=model->GetFormulation(MeltingAnalysisEnum);
    4949
    5050        //first recover parameters common to all solutions
     
    6565
    6666                if(verbose)_printf_("computing temperatures:\n");
    67                 thermal_core_nonlinear(&t_g[0],&melting_offset,fem_t,inputs,ThermalAnalysisEnum(),NoneAnalysisEnum());
     67                thermal_core_nonlinear(&t_g[0],&melting_offset,fem_t,inputs,ThermalAnalysisEnum,NoneAnalysisEnum);
    6868                inputs->Add("temperature",t_g[0],1,numberofnodes);
    6969                inputs->Add("melting_offset",melting_offset);
    7070               
    7171                if(verbose)_printf_("computing melting:\n");
    72                 diagnostic_core_linear(&m_g[0],fem_m,inputs,MeltingAnalysisEnum(),NoneAnalysisEnum());
     72                diagnostic_core_linear(&m_g[0],fem_m,inputs,MeltingAnalysisEnum,NoneAnalysisEnum);
    7373        }
    7474        else{
     
    9696                        if(verbose)_printf_("computing temperatures:\n");
    9797                        inputs->Add("temperature",t_g[i],1,numberofnodes);
    98                         thermal_core_nonlinear(&t_g[i+1],&melting_offset,fem_t,inputs,ThermalAnalysisEnum(),NoneAnalysisEnum());
     98                        thermal_core_nonlinear(&t_g[i+1],&melting_offset,fem_t,inputs,ThermalAnalysisEnum,NoneAnalysisEnum);
    9999                       
    100100                        if(verbose)_printf_("computing melting:\n");
    101101                        inputs->Add("temperature",t_g[i+1],1,numberofnodes);
    102102                        inputs->Add("melting_offset",melting_offset);
    103                         diagnostic_core_linear(&m_g[i+1],fem_m,inputs,MeltingAnalysisEnum(),NoneAnalysisEnum());
     103                        diagnostic_core_linear(&m_g[i+1],fem_m,inputs,MeltingAnalysisEnum,NoneAnalysisEnum);
    104104                }
    105105        }
  • issm/trunk/src/c/parallel/transient.cpp

    r3332 r3567  
    7676        _printf_("read and create finite element model:\n");
    7777        _printf_("\n   reading diagnostic horiz model data:\n");
    78         model->AddFormulation(fid,DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     78        model->AddFormulation(fid,DiagnosticAnalysisEnum,HorizAnalysisEnum);
    7979
    8080        _printf_("\n   reading diagnostic vert model data:\n");
    81         model->AddFormulation(fid,DiagnosticAnalysisEnum(),VertAnalysisEnum());
     81        model->AddFormulation(fid,DiagnosticAnalysisEnum,VertAnalysisEnum);
    8282       
    8383        _printf_("\n   reading diagnostic stokes model data:\n");
    84         model->AddFormulation(fid,DiagnosticAnalysisEnum(),StokesAnalysisEnum());
     84        model->AddFormulation(fid,DiagnosticAnalysisEnum,StokesAnalysisEnum);
    8585       
    8686        _printf_("\n   reading diagnostic hutter model data:\n");
    87         model->AddFormulation(fid,DiagnosticAnalysisEnum(),HutterAnalysisEnum());
     87        model->AddFormulation(fid,DiagnosticAnalysisEnum,HutterAnalysisEnum);
    8888       
    8989        _printf_("\n   reading surface and bed slope computation model data:\n");
    90         model->AddFormulation(fid,SlopecomputeAnalysisEnum());
     90        model->AddFormulation(fid,SlopecomputeAnalysisEnum);
    9191
    9292        _printf_("\n   reading prognositc model data:\n");
    93         model->AddFormulation(fid,PrognosticAnalysisEnum());
     93        model->AddFormulation(fid,PrognosticAnalysisEnum);
    9494       
    9595        /*Do we run in 3d?, in which case we need thermal and melting also:*/
     
    9797        if(dim==3){
    9898                _printf_("read and create thermal finite element model:\n");
    99                 model->AddFormulation(fid,ThermalAnalysisEnum(),TransientAnalysisEnum());
     99                model->AddFormulation(fid,ThermalAnalysisEnum,TransientAnalysisEnum);
    100100                _printf_("read and create melting finite element model:\n");
    101                 model->AddFormulation(fid,MeltingAnalysisEnum(),TransientAnalysisEnum());
     101                model->AddFormulation(fid,MeltingAnalysisEnum,TransientAnalysisEnum);
    102102        }
    103103
     
    107107
    108108        _printf_("initialize inputs:\n");
    109         model->FindParam(&u_g,NULL,NULL,"u_g",PrognosticAnalysisEnum());
    110         model->FindParam(&m_g,NULL,NULL,"m_g",PrognosticAnalysisEnum());
    111         model->FindParam(&a_g,NULL,NULL,"a_g",PrognosticAnalysisEnum());
     109        model->FindParam(&u_g,NULL,NULL,"u_g",PrognosticAnalysisEnum);
     110        model->FindParam(&m_g,NULL,NULL,"m_g",PrognosticAnalysisEnum);
     111        model->FindParam(&a_g,NULL,NULL,"a_g",PrognosticAnalysisEnum);
    112112        model->FindParam(&numberofnodes,"numberofnodes");
    113113        model->FindParam(&dt,"dt");
     
    122122       
    123123        _printf_("initialize results:\n");
    124         results=new DataSet(ResultsEnum());
     124        results=new DataSet(ResultsEnum);
    125125        MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime();
    126126
     
    139139
    140140                _printf_("process results:\n");
    141                 ProcessResults(&processed_results,results,model,TransientAnalysisEnum());
     141                ProcessResults(&processed_results,results,model,TransientAnalysisEnum);
    142142               
    143143                _printf_("write results to disk:\n");
     
    150150                #ifdef _HAVE_DAKOTA_
    151151                MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
    152                 Qmux(model,inputs,TransientAnalysisEnum(),NoneAnalysisEnum());
     152                Qmux(model,inputs,TransientAnalysisEnum,NoneAnalysisEnum);
    153153                MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
    154154                #else
  • issm/trunk/src/c/parallel/transient_core_2d.cpp

    r3446 r3567  
    6363
    6464        /*recover fem models: */
    65         fem_dh=model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum());
    66         fem_dv=model->GetFormulation(DiagnosticAnalysisEnum(),VertAnalysisEnum());
    67         fem_ds=model->GetFormulation(DiagnosticAnalysisEnum(),StokesAnalysisEnum());
    68         fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum(),HutterAnalysisEnum());
    69         fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum());
    70         fem_p=model->GetFormulation(PrognosticAnalysisEnum());
     65        fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum);
     66        fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum);
     67        fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum);
     68        fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum);
     69        fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum);
     70        fem_p=model->GetFormulation(PrognosticAnalysisEnum);
    7171
    7272
     
    119119                //Deal with velocities.
    120120                _printf_("%s\n","computing new velocity");
    121                 diagnostic_results=new DataSet(ResultsEnum());
     121                diagnostic_results=new DataSet(ResultsEnum);
    122122                diagnostic_core(diagnostic_results,model, inputs);
    123123
     
    131131               
    132132                inputs->Add("velocity",u_g,2,numberofnodes);
    133                 prognostic_results=new DataSet(ResultsEnum());
     133                prognostic_results=new DataSet(ResultsEnum);
    134134                prognostic_core(prognostic_results,model, inputs);
    135135       
  • issm/trunk/src/c/parallel/transient_core_3d.cpp

    r3544 r3567  
    6767
    6868        /*recover fem models: */
    69         fem_dh=model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum());
    70         fem_dv=model->GetFormulation(DiagnosticAnalysisEnum(),VertAnalysisEnum());
    71         fem_ds=model->GetFormulation(DiagnosticAnalysisEnum(),StokesAnalysisEnum());
    72         fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum(),HutterAnalysisEnum());
    73         fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum());
    74         fem_p=model->GetFormulation(PrognosticAnalysisEnum());
    75         fem_t=model->GetFormulation(ThermalAnalysisEnum());
    76         fem_m=model->GetFormulation(MeltingAnalysisEnum());
     69        fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum);
     70        fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum);
     71        fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum);
     72        fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum);
     73        fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum);
     74        fem_p=model->GetFormulation(PrognosticAnalysisEnum);
     75        fem_t=model->GetFormulation(ThermalAnalysisEnum);
     76        fem_m=model->GetFormulation(MeltingAnalysisEnum);
    7777
    7878
     
    134134                //Deal with temperature first
    135135                if(verbose)_printf_("%s\n","computing temperature");
    136                 thermal_core_nonlinear(&t_g,&melting_offset,fem_t,inputs,ThermalAnalysisEnum(),TransientAnalysisEnum());
     136                thermal_core_nonlinear(&t_g,&melting_offset,fem_t,inputs,ThermalAnalysisEnum,TransientAnalysisEnum);
    137137                if(verbose)_printf_("%s\n","computing melting");
    138138                inputs->Add("temperature",t_g,1,numberofnodes);
    139139                inputs->Add("melting_offset",melting_offset);
    140                 diagnostic_core_linear(&m_g,fem_m,inputs,MeltingAnalysisEnum(),TransientAnalysisEnum());
     140                diagnostic_core_linear(&m_g,fem_m,inputs,MeltingAnalysisEnum,TransientAnalysisEnum);
    141141
    142142                //Compute depth averaged temperature and add to inputs
     
    148148
    149149                //Deal with velocities.
    150                 diagnostic_results=new DataSet(ResultsEnum());
     150                diagnostic_results=new DataSet(ResultsEnum);
    151151                diagnostic_core(diagnostic_results,model, inputs);
    152152
     
    160160               
    161161                inputs->Add("velocity",u_g,3,numberofnodes);
    162                 prognostic_results=new DataSet(ResultsEnum()); prognostic_core(prognostic_results,model, inputs);
     162                prognostic_results=new DataSet(ResultsEnum); prognostic_core(prognostic_results,model, inputs);
    163163       
    164164                //extract h_g prognostic_results, and erase prognostic_results;
     
    198198
    199199                        model->FindParam(&outputfilename,"outputfilename");
    200                         ProcessResults(&processed_results,results,model,TransientAnalysisEnum());
     200                        ProcessResults(&processed_results,results,model,TransientAnalysisEnum);
    201201                        OutputResults(processed_results,outputfilename);
    202202
  • issm/trunk/src/c/shared/Dofs/DistributeNumDofs.cpp

    r3354 r3567  
    1313
    1414        /*ok, according to analysis type: */
    15         if (analysis_type==ControlAnalysisEnum()){
    16                 if (sub_analysis_type==HorizAnalysisEnum()){
     15        if (analysis_type==ControlAnalysisEnum){
     16                if (sub_analysis_type==HorizAnalysisEnum){
    1717                        numdofs=2;
    1818                }
    19                 else if (sub_analysis_type==VertAnalysisEnum()){
     19                else if (sub_analysis_type==VertAnalysisEnum){
    2020                        numdofs=1;
    2121                }
    22                 else if (sub_analysis_type==StokesAnalysisEnum()){
     22                else if (sub_analysis_type==StokesAnalysisEnum){
    2323                        numdofs=4;
    2424                }
    25                 else if (sub_analysis_type==HutterAnalysisEnum()){
     25                else if (sub_analysis_type==HutterAnalysisEnum){
    2626                        numdofs=2;
    2727                }
    2828        }
    29         else if (analysis_type==DiagnosticAnalysisEnum()){
    30                 if (sub_analysis_type==HorizAnalysisEnum()){
     29        else if (analysis_type==DiagnosticAnalysisEnum){
     30                if (sub_analysis_type==HorizAnalysisEnum){
    3131                        numdofs=2;
    3232                }
    33                 else if (sub_analysis_type==VertAnalysisEnum()){
     33                else if (sub_analysis_type==VertAnalysisEnum){
    3434                        numdofs=1;
    3535                }
    36                 else if (sub_analysis_type==StokesAnalysisEnum()){
     36                else if (sub_analysis_type==StokesAnalysisEnum){
    3737                        numdofs=4;
    3838                }
    39                 else if (sub_analysis_type==HutterAnalysisEnum()){
     39                else if (sub_analysis_type==HutterAnalysisEnum){
    4040                        numdofs=2;
    4141                }
    4242        }
    43         else if (analysis_type==SlopecomputeAnalysisEnum()){
     43        else if (analysis_type==SlopecomputeAnalysisEnum){
    4444                numdofs=1;
    4545        }
    46         else if (analysis_type==ThermalAnalysisEnum()){
     46        else if (analysis_type==ThermalAnalysisEnum){
    4747                numdofs=1;
    4848        }
    49         else if (analysis_type==MeltingAnalysisEnum()){
     49        else if (analysis_type==MeltingAnalysisEnum){
    5050                numdofs=1;
    5151        }
    52         else if (analysis_type==PrognosticAnalysisEnum()){
     52        else if (analysis_type==PrognosticAnalysisEnum){
    5353                numdofs=1;
    5454        }
    55         else if (analysis_type==Prognostic2AnalysisEnum()){
     55        else if (analysis_type==Prognostic2AnalysisEnum){
    5656                numdofs=1;
    5757        }
    58         else if (analysis_type==BalancedthicknessAnalysisEnum()){
     58        else if (analysis_type==BalancedthicknessAnalysisEnum){
    5959                numdofs=1;
    6060        }
    61         else if (analysis_type==BalancedvelocitiesAnalysisEnum()){
     61        else if (analysis_type==BalancedvelocitiesAnalysisEnum){
    6262                numdofs=1;
    6363        }
  • issm/trunk/src/m/enum/AdjointAnalysisEnum.m

    r1714 r3567  
    77%      macro=AdjointAnalysisEnum()
    88
    9 macro=221;
     9macro=16;
  • issm/trunk/src/m/enum/AirEnum.m

    r1749 r3567  
    77%      macro=AirEnum()
    88
    9 macro=603;
     9macro=75;
  • issm/trunk/src/m/enum/AnalysisEnum.m

    r1714 r3567  
    77%      macro=AnalysisEnum()
    88
    9 macro=200;
     9macro=9;
  • issm/trunk/src/m/enum/AnalysisTypeFromEnum.m

    r3362 r3567  
    1111%initialize output
    1212string='not found';
    13 
    1413
    1514if enum==AnalysisEnum(),
  • issm/trunk/src/m/enum/BalancedthicknessAnalysisEnum.m

    r3354 r3567  
    77%      macro=BalancedthicknessAnalysisEnum()
    88
    9 macro=252;
     9macro=29;
  • issm/trunk/src/m/enum/BalancedvelocitiesAnalysisEnum.m

    r3354 r3567  
    77%      macro=BalancedvelocitiesAnalysisEnum()
    88
    9 macro=253;
     9macro=30;
  • issm/trunk/src/m/enum/BeamEnum.m

    r3417 r3567  
    77%      macro=BeamEnum()
    88
    9 macro=416;
     9macro=49;
  • issm/trunk/src/m/enum/BedXAnalysisEnum.m

    r1714 r3567  
    77%      macro=BedXAnalysisEnum()
    88
    9 macro=243;
     9macro=25;
  • issm/trunk/src/m/enum/BedYAnalysisEnum.m

    r1714 r3567  
    77%      macro=BedYAnalysisEnum()
    88
    9 macro=244;
     9macro=26;
  • issm/trunk/src/m/enum/ConstraintsEnum.m

    r1714 r3567  
    77%      macro=ConstraintsEnum()
    88
    9 macro=103;
     9macro=3;
  • issm/trunk/src/m/enum/ControlAnalysisEnum.m

    r1714 r3567  
    77%      macro=ControlAnalysisEnum()
    88
    9 macro=220;
     9macro=15;
  • issm/trunk/src/m/enum/DatasetsEnum.m

    r1714 r3567  
    77%      macro=DatasetsEnum()
    88
    9 macro=100;
     9macro=0;
  • issm/trunk/src/m/enum/DiagnosticAnalysisEnum.m

    r1714 r3567  
    77%      macro=DiagnosticAnalysisEnum()
    88
    9 macro=210;
     9macro=10;
  • issm/trunk/src/m/enum/DofIndexingEnum.m

    r3417 r3567  
    77%      macro=DofIndexingEnum()
    88
    9 macro=417;
     9macro=51;
  • issm/trunk/src/m/enum/DofVecEnum.m

    r2311 r3567  
    77%      macro=DofVecEnum()
    88
    9 macro=495;
     9macro=69;
  • issm/trunk/src/m/enum/ElementEnum.m

    r1714 r3567  
    77%      macro=ElementEnum()
    88
    9 macro=410;
     9macro=44;
  • issm/trunk/src/m/enum/ElementPropertiesEnum.m

    r3383 r3567  
    77%      macro=ElementPropertiesEnum()
    88
    9 macro=412;
     9macro=46;
  • issm/trunk/src/m/enum/ElementsEnum.m

    r1714 r3567  
    77%      macro=ElementsEnum()
    88
    9 macro=101;
     9macro=1;
  • issm/trunk/src/m/enum/FormulationEnum.m

    r1714 r3567  
    77%      macro=FormulationEnum()
    88
    9 macro=300;
     9macro=36;
  • issm/trunk/src/m/enum/GeographyEnum.m

    r1729 r3567  
    77%      macro=GeographyEnum()
    88
    9 macro=500;
     9macro=70;
  • issm/trunk/src/m/enum/GradientAnalysisEnum.m

    r1714 r3567  
    77%      macro=GradientAnalysisEnum()
    88
    9 macro=223;
     9macro=18;
  • issm/trunk/src/m/enum/HookEnum.m

    r3383 r3567  
    77%      macro=HookEnum()
    88
    9 macro=409;
     9macro=43;
  • issm/trunk/src/m/enum/HorizAnalysisEnum.m

    r1714 r3567  
    77%      macro=HorizAnalysisEnum()
    88
    9 macro=211;
     9macro=11;
  • issm/trunk/src/m/enum/HutterAnalysisEnum.m

    r1714 r3567  
    77%      macro=HutterAnalysisEnum()
    88
    9 macro=213;
     9macro=13;
  • issm/trunk/src/m/enum/HutterFormulationEnum.m

    r1714 r3567  
    77%      macro=HutterFormulationEnum()
    88
    9 macro=302;
     9macro=38;
  • issm/trunk/src/m/enum/IceEnum.m

    r1749 r3567  
    77%      macro=IceEnum()
    88
    9 macro=602;
     9macro=74;
  • issm/trunk/src/m/enum/IceSheetEnum.m

    r1729 r3567  
    77%      macro=IceSheetEnum()
    88
    9 macro=502;
     9macro=71;
  • issm/trunk/src/m/enum/IceShelfEnum.m

    r1729 r3567  
    77%      macro=IceShelfEnum()
    88
    9 macro=502;
     9macro=72;
  • issm/trunk/src/m/enum/IcefrontEnum.m

    r1714 r3567  
    77%      macro=IcefrontEnum()
    88
    9 macro=431;
     9macro=55;
  • issm/trunk/src/m/enum/InputEnum.m

    r1714 r3567  
    77%      macro=InputEnum()
    88
    9 macro=450;
     9macro=64;
  • issm/trunk/src/m/enum/InverseAnalysisEnum.m

    r1714 r3567  
    77%      macro=InverseAnalysisEnum()
    88
    9 macro=222;
     9macro=17;
  • issm/trunk/src/m/enum/LoadEnum.m

    r1714 r3567  
    77%      macro=LoadEnum()
    88
    9 macro=430;
     9macro=54;
  • issm/trunk/src/m/enum/LoadsEnum.m

    r1714 r3567  
    77%      macro=LoadsEnum()
    88
    9 macro=104;
     9macro=4;
  • issm/trunk/src/m/enum/MacAyealFormulationEnum.m

    r1714 r3567  
    77%      macro=MacAyealFormulationEnum()
    88
    9 macro=303;
     9macro=39;
  • issm/trunk/src/m/enum/MaterialEnum.m

    r1714 r3567  
    77%      macro=MaterialEnum()
    88
    9 macro=440;
     9macro=60;
  • issm/trunk/src/m/enum/MaterialsEnum.m

    r1714 r3567  
    77%      macro=MaterialsEnum()
    88
    9 macro=105;
     9macro=5;
  • issm/trunk/src/m/enum/MaticeEnum.m

    r1714 r3567  
    77%      macro=MaticeEnum()
    88
    9 macro=441;
     9macro=61;
  • issm/trunk/src/m/enum/MatparEnum.m

    r1714 r3567  
    77%      macro=MatparEnum()
    88
    9 macro=442;
     9macro=62;
  • issm/trunk/src/m/enum/MelangeEnum.m

    r1749 r3567  
    77%      macro=MelangeEnum()
    88
    9 macro=604;
     9macro=76;
  • issm/trunk/src/m/enum/MeltingAnalysisEnum.m

    r1714 r3567  
    77%      macro=MeltingAnalysisEnum()
    88
    9 macro=260;
     9macro=31;
  • issm/trunk/src/m/enum/Mesh2gridAnalysisEnum.m

    r1714 r3567  
    77%      macro=Mesh2gridAnalysisEnum()
    88
    9 macro=270;
     9macro=32;
  • issm/trunk/src/m/enum/NodeEnum.m

    r3417 r3567  
    77%      macro=NodeEnum()
    88
    9 macro=421;
     9macro=52;
  • issm/trunk/src/m/enum/NodePropertiesEnum.m

    r3417 r3567  
    77%      macro=NodePropertiesEnum()
    88
    9 macro=413;
     9macro=47;
  • issm/trunk/src/m/enum/NodesEnum.m

    r1714 r3567  
    77%      macro=NodesEnum()
    88
    9 macro=102;
     9macro=2;
  • issm/trunk/src/m/enum/NoneAnalysisEnum.m

    r1714 r3567  
    77%      macro=NoneAnalysisEnum()
    88
    9 macro=290;
     9macro=35;
  • issm/trunk/src/m/enum/NoneFormulationEnum.m

    r1714 r3567  
    77%      macro=NoneFormulationEnum()
    88
    9 macro=301;
     9macro=37;
  • issm/trunk/src/m/enum/NumericalfluxEnum.m

    r3359 r3567  
    77%      macro=NumericalfluxEnum()
    88
    9 macro=435;
     9macro=59;
  • issm/trunk/src/m/enum/NumparEnum.m

    r2333 r3567  
    77%      macro=NumparEnum()
    88
    9 macro=443;
     9macro=63;
  • issm/trunk/src/m/enum/ObjectEnum.m

    r1714 r3567  
    77%      macro=ObjectEnum()
    88
    9 macro=400;
     9macro=42;
  • issm/trunk/src/m/enum/ParamEnum.m

    r1714 r3567  
    77%      macro=ParamEnum()
    88
    9 macro=460;
     9macro=65;
  • issm/trunk/src/m/enum/ParametersAnalysisEnum.m

    r1714 r3567  
    77%      macro=ParametersAnalysisEnum()
    88
    9 macro=280;
     9macro=33;
  • issm/trunk/src/m/enum/ParametersEnum.m

    r1714 r3567  
    77%      macro=ParametersEnum()
    88
    9 macro=106;
     9macro=6;
  • issm/trunk/src/m/enum/PattynFormulationEnum.m

    r1714 r3567  
    77%      macro=PattynFormulationEnum()
    88
    9 macro=304;
     9macro=40;
  • issm/trunk/src/m/enum/PengridEnum.m

    r1714 r3567  
    77%      macro=PengridEnum()
    88
    9 macro=434;
     9macro=58;
  • issm/trunk/src/m/enum/PenpairEnum.m

    r1714 r3567  
    77%      macro=PenpairEnum()
    88
    9 macro=433;
     9macro=57;
  • issm/trunk/src/m/enum/PentaEnum.m

    r3417 r3567  
    77%      macro=PentaEnum()
    88
    9 macro=414;
     9macro=48;
  • issm/trunk/src/m/enum/Prognostic2AnalysisEnum.m

    r3354 r3567  
    77%      macro=Prognostic2AnalysisEnum()
    88
    9 macro=251;
     9macro=28;
  • issm/trunk/src/m/enum/PrognosticAnalysisEnum.m

    r1714 r3567  
    77%      macro=PrognosticAnalysisEnum()
    88
    9 macro=250;
     9macro=27;
  • issm/trunk/src/m/enum/ResultEnum.m

    r1714 r3567  
    77%      macro=ResultEnum()
    88
    9 macro=470;
     9macro=66;
  • issm/trunk/src/m/enum/ResultsEnum.m

    r1714 r3567  
    77%      macro=ResultsEnum()
    88
    9 macro=107;
     9macro=7;
  • issm/trunk/src/m/enum/RgbEnum.m

    r1714 r3567  
    77%      macro=RgbEnum()
    88
    9 macro=480;
     9macro=67;
  • issm/trunk/src/m/enum/RiftfrontEnum.m

    r1714 r3567  
    77%      macro=RiftfrontEnum()
    88
    9 macro=432;
     9macro=56;
  • issm/trunk/src/m/enum/SingEnum.m

    r3417 r3567  
    77%      macro=SingEnum()
    88
    9 macro=415;
     9macro=50;
  • issm/trunk/src/m/enum/SlopecomputeAnalysisEnum.m

    r3087 r3567  
    77%      macro=SlopecomputeAnalysisEnum()
    88
    9 macro=240;
     9macro=22;
  • issm/trunk/src/m/enum/SpcEnum.m

    r1714 r3567  
    77%      macro=SpcEnum()
    88
    9 macro=490;
     9macro=68;
  • issm/trunk/src/m/enum/SteadyAnalysisEnum.m

    r3362 r3567  
    77%      macro=SteadyAnalysisEnum()
    88
    9 macro=232;
     9macro=21;
  • issm/trunk/src/m/enum/SteadystateAnalysisEnum.m

    r1865 r3567  
    77%      macro=SteadystateAnalysisEnum()
    88
    9 macro=281;
     9macro=34;
  • issm/trunk/src/m/enum/StokesAnalysisEnum.m

    r1714 r3567  
    77%      macro=StokesAnalysisEnum()
    88
    9 macro=212;
     9macro=12;
  • issm/trunk/src/m/enum/StokesFormulationEnum.m

    r1714 r3567  
    77%      macro=StokesFormulationEnum()
    88
    9 macro=305;
     9macro=41;
  • issm/trunk/src/m/enum/SurfaceXAnalysisEnum.m

    r1714 r3567  
    77%      macro=SurfaceXAnalysisEnum()
    88
    9 macro=241;
     9macro=23;
  • issm/trunk/src/m/enum/SurfaceYAnalysisEnum.m

    r1714 r3567  
    77%      macro=SurfaceYAnalysisEnum()
    88
    9 macro=242;
     9macro=24;
  • issm/trunk/src/m/enum/ThermalAnalysisEnum.m

    r1714 r3567  
    77%      macro=ThermalAnalysisEnum()
    88
    9 macro=230;
     9macro=19;
  • issm/trunk/src/m/enum/TransientAnalysisEnum.m

    r1901 r3567  
    77%      macro=TransientAnalysisEnum()
    88
    9 macro=231;
     9macro=20;
  • issm/trunk/src/m/enum/TriaEnum.m

    r1714 r3567  
    77%      macro=TriaEnum()
    88
    9 macro=411;
     9macro=45;
  • issm/trunk/src/m/enum/VertAnalysisEnum.m

    r1714 r3567  
    77%      macro=VertAnalysisEnum()
    88
    9 macro=214;
     9macro=14;
  • issm/trunk/src/m/enum/VertexEnum.m

    r3417 r3567  
    77%      macro=VertexEnum()
    88
    9 macro=420;
     9macro=53;
  • issm/trunk/src/m/enum/VerticesEnum.m

    r3417 r3567  
    77%      macro=VerticesEnum()
    88
    9 macro=108;
     9macro=8;
  • issm/trunk/src/m/enum/WaterEnum.m

    r1750 r3567  
    77%      macro=WaterEnum()
    88
    9 macro=601;
     9macro=73;
  • issm/trunk/src/mex/TriMeshProcessRifts/TriMeshProcessRifts.cpp

    r1904 r3567  
    287287                        /*Friction fraction, fractionincrement  and fill: */
    288288                        mxSetField(pmxa_array,i,"friction",mxCreateDoubleScalar(0));
    289                         mxSetField(pmxa_array,i,"fill",mxCreateDoubleScalar(IceEnum())); //default is ice
     289                        mxSetField(pmxa_array,i,"fill",mxCreateDoubleScalar(IceEnum)); //default is ice
    290290                        mxSetField(pmxa_array,i,"fraction",mxCreateDoubleScalar(0)); //default is ice
    291291                        mxSetField(pmxa_array,i,"fractionincrement",mxCreateDoubleScalar(0.1));
     
    295295                /*output NaN :*/
    296296                pNaN=(double*)xmalloc(sizeof(double));
    297                 *pNaN=0.0/0.0;
     297                *pNaN=NAN;
    298298                pmxa_array= mxCreateDoubleMatrix(0,0,mxREAL);
    299299                mxSetM(pmxa_array,1);
Note: See TracChangeset for help on using the changeset viewer.