Changeset 24959


Ignore:
Timestamp:
06/03/20 14:59:24 (5 years ago)
Author:
Eric.Larour
Message:

CHG: fixed large issues in computation of centroid. Also introduced planetradius instead
of area.

Location:
issm/trunk-jpl/src/c
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/analyses/EsaAnalysis.cpp

    r24950 r24959  
    4747        int         M,m,lower_row,upper_row;
    4848        IssmDouble  degacc=.01;
     49        IssmDouble  planetradius=0;
     50        IssmDouble  planetarea=0;
    4951
    5052        int     numoutputs;
    … …  
    5961        /*some constant parameters: */
    6062        parameters->AddObject(iomodel->CopyConstantObject("md.esa.hemisphere",EsaHemisphereEnum));
    61         parameters->AddObject(iomodel->CopyConstantObject("md.slr.eartharea",SealevelEarthAreaEnum));
    6263        parameters->AddObject(iomodel->CopyConstantObject("md.slr.rigid",SealevelriseRigidEnum));
    6364        parameters->AddObject(iomodel->CopyConstantObject("md.slr.horiz",SealevelriseHorizEnum));
    6465        parameters->AddObject(iomodel->CopyConstantObject("md.slr.elastic",SealevelriseElasticEnum));
    6566        parameters->AddObject(iomodel->CopyConstantObject("md.slr.rotation",SealevelriseRotationEnum));
    66        
     67
     68        /*deal with planet radius and area: */
     69        parameters->AddObject(iomodel->CopyConstantObject("md.slr.planetradius",SealevelPlanetRadiusEnum));
     70        iomodel->FetchData(&planetradius,"md.slr.planetradius");
     71        planetarea=4*PI*planetradius*planetradius;
     72        parameters->AddObject(new DoubleParam(SealevelPlanetAreaEnum,planetarea));
     73
    6774        /*love numbers: */
    6875        iomodel->FetchData(&love_h,&nl,NULL,"md.esa.love_h");

  • issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp

    r24947 r24959  
    11#include "./SealevelriseAnalysis.h"
     2#include <math.h>
    23#include "../toolkits/toolkits.h"
    34#include "../classes/classes.h"
    … …  
    292293        int         M,m,lower_row,upper_row;
    293294        IssmDouble  degacc=.01;
     295        IssmDouble  planetradius=0;
     296        IssmDouble  planetarea=0;
    294297
    295298        int     numoutputs;
    … …  
    320323        parameters->AddObject(iomodel->CopyConstantObject("md.slr.ocean_area_scaling",SealevelriseOceanAreaScalingEnum));
    321324        parameters->AddObject(iomodel->CopyConstantObject("md.slr.geodetic",SealevelriseGeodeticEnum));
    322         parameters->AddObject(iomodel->CopyConstantObject("md.slr.eartharea",SealevelEarthAreaEnum));
     325        parameters->AddObject(iomodel->CopyConstantObject("md.slr.planetradius",SealevelPlanetRadiusEnum));
     326
     327        /*compute planet area and plug into parameters:*/
     328        iomodel->FetchData(&planetradius,"md.slr.planetradius");
     329        planetarea=4*PI*planetradius*planetradius;
     330        parameters->AddObject(new DoubleParam(SealevelPlanetAreaEnum,planetarea));
    323331
    324332        /*Deal with dsl multi-model ensembles: {{{*/

  • issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

    r24951 r24959  
    52975297        IssmDouble llr_list[NUMVERTICES][3];
    52985298        IssmDouble xyz_list[NUMVERTICES][3];
    5299         IssmDouble area,eartharea;
     5299        IssmDouble area,planetarea;
    53005300        IssmDouble I;           //ice/water loading
    53015301        IssmDouble late,longe,re;
    … …  
    53315331       
    53325332        /*recover earth area: */
    5333         this->parameters->FindParam(&eartharea,SealevelEarthAreaEnum);
     5333        this->parameters->FindParam(&planetarea,SealevelPlanetAreaEnum);
    53345334
    53355335        /*how many dofs are we working with here? */
    … …  
    53565356        }
    53575357
    5358         // correction at the north pole
     5358        // correction at the north pole: given longitude of the North pole a definition
     5359        // closer to the other two vertices.
    53595360        if(llr_list[0][0]==0)llr_list[0][1]=(llr_list[1][1]+llr_list[2][1])/2.0;
    53605361        if(llr_list[1][0]==0)llr_list[1][1]=(llr_list[0][1]+llr_list[2][1])/2.0;
    53615362        if(llr_list[2][0]==0)llr_list[2][1]=(llr_list[0][1]+llr_list[1][1])/2.0;
    53625363
    5363         //correction at the south pole
     5364        // correction at the north pole: given longitude of the North pole a definition
     5365        // closer to the other two vertices.
    53645366        if(llr_list[0][0]==180)llr_list[0][1]=(llr_list[1][1]+llr_list[2][1])/2.0;
    53655367        if(llr_list[1][0]==180)llr_list[1][1]=(llr_list[0][1]+llr_list[2][1])/2.0;
    … …  
    53705372
    53715373        late=90-late;
    5372         if(longe>180)longe=(longe-180)-180;
     5374        if(longe>180)longe=longe-360;
    53735375
    53745376        late=late/180*PI;
    … …  
    54325434
    54335435                /*Add all components to the pUp solution vectors:*/
    5434                 U_values[i]+=3*rho_ice/rho_earth*area/eartharea*I*U_elastic[i];
    5435                 N_values[i]+=3*rho_ice/rho_earth*area/eartharea*I*N_elastic[i];
    5436                 E_values[i]+=3*rho_ice/rho_earth*area/eartharea*I*E_elastic[i];
     5436                U_values[i]+=3*rho_ice/rho_earth*area/planetarea*I*U_elastic[i];
     5437                N_values[i]+=3*rho_ice/rho_earth*area/planetarea*I*N_elastic[i];
     5438                E_values[i]+=3*rho_ice/rho_earth*area/planetarea*I*E_elastic[i];
    54375439        }
    54385440        pUp->SetValues(gsize,indices,U_values,ADD_VAL);
    … …  
    54795481
    54805482        /*Compute area of element:*/
    5481         IssmDouble area,eartharea;
     5483        IssmDouble area,planetarea;
    54825484        this->GetInput2Value(&area,AreaEnum);
    54835485
    54845486        /*recover earth area: */
    5485         this->parameters->FindParam(&eartharea,SealevelEarthAreaEnum);
     5487        this->parameters->FindParam(&planetarea,SealevelPlanetAreaEnum);
    54865488
    54875489        /*Compute lat,long,radius of elemental centroid: */
    … …  
    55415543                 * ALL in geographic coordinates
    55425544                 * */
    5543                 dI_list[0] = -4*PI*(rho_water*S*area)*pow(re,4)*(sin(late)*cos(late)*cos(longe))/eartharea;
    5544                 dI_list[1] = -4*PI*(rho_water*S*area)*pow(re,4)*(sin(late)*cos(late)*sin(longe))/eartharea;
    5545                 dI_list[2] = +4*PI*(rho_water*S*area)*pow(re,4)*(1-pow(sin(late),2))/eartharea;
     5545                dI_list[0] = -4*PI*(rho_water*S*area)*pow(re,4)*(sin(late)*cos(late)*cos(longe))/planetarea;
     5546                dI_list[1] = -4*PI*(rho_water*S*area)*pow(re,4)*(sin(late)*cos(late)*sin(longe))/planetarea;
     5547                dI_list[2] = +4*PI*(rho_water*S*area)*pow(re,4)*(1-pow(sin(late),2))/planetarea;
    55465548        }
    55475549        else if(masks->isiceonly[this->lid]){
    … …  
    55565558                deltathickness_input->GetInputAverage(&I);
    55575559
    5558                 dI_list[0] = -4*PI*(rho_ice*I*area)*pow(re,4)*(sin(late)*cos(late)*cos(longe))/eartharea;
    5559                 dI_list[1] = -4*PI*(rho_ice*I*area)*pow(re,4)*(sin(late)*cos(late)*sin(longe))/eartharea;
    5560                 dI_list[2] = +4*PI*(rho_ice*I*area)*pow(re,4)*(1-pow(sin(late),2))/eartharea;
     5560                dI_list[0] = -4*PI*(rho_ice*I*area)*pow(re,4)*(sin(late)*cos(late)*cos(longe))/planetarea;
     5561                dI_list[1] = -4*PI*(rho_ice*I*area)*pow(re,4)*(sin(late)*cos(late)*sin(longe))/planetarea;
     5562                dI_list[2] = +4*PI*(rho_ice*I*area)*pow(re,4)*(1-pow(sin(late),2))/planetarea;
    55615563        }
    55625564
    … …  
    55855587        IssmDouble llr_list[NUMVERTICES][3];
    55865588        IssmDouble xyz_list[NUMVERTICES][3];
    5587         IssmDouble area,eartharea;
     5589        IssmDouble area,planetarea,planetradius;
    55885590        IssmDouble I;  //change in ice thickness or water level(Farrel and Clarke, Equ. 4)
    55895591        IssmDouble rho_earth;
    … …  
    56175619        this->parameters->FindParam(&computeelastic,SealevelriseElasticEnum);
    56185620        this->parameters->FindParam(&gsize,MeshNumberofverticesEnum);
    5619         this->parameters->FindParam(&eartharea,SealevelEarthAreaEnum);
     5621        this->parameters->FindParam(&planetarea,SealevelPlanetAreaEnum);
     5622        this->parameters->FindParam(&planetradius,SealevelPlanetRadiusEnum);
    56205623        this->parameters->FindParam(&horiz,SealevelriseHorizEnum);
    56215624
    … …  
    56455648        area=GetAreaSpherical();
    56465649
     5650
     5651        /* Where is the centroid of this element:*/
     5652
     5653        /*retrieve coordinates: */
     5654        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5655
    56475656        /*figure out gravity center of our element (Cartesian): */
    5648         ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    56495657        x_element=(xyz_list[0][0]+xyz_list[1][0]+xyz_list[2][0])/3.0;
    56505658        y_element=(xyz_list[0][1]+xyz_list[1][1]+xyz_list[2][1])/3.0;
    56515659        z_element=(xyz_list[0][2]+xyz_list[1][2]+xyz_list[2][2])/3.0;
    56525660
    5653         /* Where is the centroid of this element?:{{{*/
    5654 
    5655         /*retrieve coordinates: */
    5656         ::GetVerticesCoordinates(&llr_list[0][0],this->vertices,NUMVERTICES,spherical);
    5657 
    5658         IssmDouble minlong=400;
    5659         IssmDouble maxlong=-20;
    5660         for (int i=0;i<NUMVERTICES;i++){
    5661                 llr_list[i][0]=(90-llr_list[i][0]);
    5662                 if(llr_list[i][1]<0)llr_list[i][1]=180+(180+llr_list[i][1]);
    5663                 if(llr_list[i][1]>maxlong)maxlong=llr_list[i][1];
    5664                 if(llr_list[i][1]<minlong)minlong=llr_list[i][1];
    5665         }
    5666         if(minlong==0 && maxlong>180){
    5667                 if (llr_list[0][1]==0)llr_list[0][1]=360;
    5668                 if (llr_list[1][1]==0)llr_list[1][1]=360;
    5669                 if (llr_list[2][1]==0)llr_list[2][1]=360;
    5670         }
    5671 
    5672         // correction at the north pole
    5673         if(llr_list[0][0]==0)llr_list[0][1]=(llr_list[1][1]+llr_list[2][1])/2.0;
    5674         if(llr_list[1][0]==0)llr_list[1][1]=(llr_list[0][1]+llr_list[2][1])/2.0;
    5675         if(llr_list[2][0]==0)llr_list[2][1]=(llr_list[0][1]+llr_list[1][1])/2.0;
    5676 
    5677         //correction at the south pole
    5678         if(llr_list[0][0]==180)llr_list[0][1]=(llr_list[1][1]+llr_list[2][1])/2.0;
    5679         if(llr_list[1][0]==180)llr_list[1][1]=(llr_list[0][1]+llr_list[2][1])/2.0;
    5680         if(llr_list[2][0]==180)llr_list[2][1]=(llr_list[0][1]+llr_list[1][1])/2.0;
    5681 
    5682         late=(llr_list[0][0]+llr_list[1][0]+llr_list[2][0])/3.0;
    5683         longe=(llr_list[0][1]+llr_list[1][1]+llr_list[2][1])/3.0;
    5684 
    5685         late=90-late;
    5686         if(longe>180)longe=(longe-180)-180;
    5687 
    5688         late=late/180*PI;
    5689         longe=longe/180*PI;
    5690         /*}}}*/
    5691                
    5692         constant=3/rho_earth*area/eartharea;
     5661        /*compute gravity center in lat,long: */
     5662        late= asin(z_element/planetradius);
     5663        longe = atan2(y_element,x_element);
     5664
     5665        constant=3/rho_earth*area/planetarea;
    56935666       
    56945667        for(int i=0;i<gsize;i++){
    … …  
    56985671                /*Compute alpha angle between centroid and current vertex and index into precomputed tables: */
    56995672                lati=latitude[i]/180*PI; longi=longitude[i]/180*PI;
    5700                 delPhi=fabs(lati-late); delLambda=fabs(longi-longe);
     5673                delPhi=fabs(lati-late); delLambda=fabs(longi-longe); if (delLambda>PI)delLambda=2*PI-delLambda;
    57015674                alpha=2.*asin(sqrt(pow(sin(delPhi/2),2)+cos(lati)*cos(late)*pow(sin(delLambda/2),2)));
    57025675                indices[i]=alpha/PI*reCast<IssmDouble,int>(M-1);
    … …  
    59045877        bool spherical=true;
    59055878        IssmDouble llr_list[NUMVERTICES][3];
    5906         IssmDouble area,eartharea;
     5879        IssmDouble area,planetarea;
    59075880        IssmDouble I;  //change in ice thickness or water level(Farrel and Clarke, Equ. 4)
    59085881        IssmDouble rho;
    … …  
    59395912
    59405913        /*recover earth area: */
    5941         this->parameters->FindParam(&eartharea,SealevelEarthAreaEnum);
     5914        this->parameters->FindParam(&planetarea,SealevelPlanetAreaEnum);
    59425915
    59435916        /*recover love numbers and computational flags: */
    … …  
    60386011
    60396012                        /*Add all components to the pSgi or pSgo solution vectors:*/
    6040                         Sgi[i]+=3*rho_water/rho_earth*area/eartharea*I*(G_rigid+G_elastic);
     6013                        Sgi[i]+=3*rho_water/rho_earth*area/planetarea*I*(G_rigid+G_elastic);
    60416014                }
    60426015        }

  • issm/trunk-jpl/src/c/shared/Enum/Enum.vim

    r24947 r24959  
    302302syn keyword cConstant RootPathEnum
    303303syn keyword cConstant SaveResultsEnum
    304 syn keyword cConstant SealevelEarthAreaEnum
     304syn keyword cConstant SealevelPlanetRadiusEnum
     305syn keyword cConstant SealevelPlanetAreaEnum
    305306syn keyword cConstant SealevelEustaticEnum
    306307syn keyword cConstant SealevelriseAbstolEnum

  • issm/trunk-jpl/src/c/shared/Enum/EnumDefinitions.h

    r24947 r24959  
    296296        RootPathEnum,
    297297        SaveResultsEnum,
    298         SealevelEarthAreaEnum,
     298        SealevelPlanetRadiusEnum,
     299        SealevelPlanetAreaEnum,
    299300        SealevelEustaticEnum,
    300301        SealevelriseAbstolEnum,

  • issm/trunk-jpl/src/c/shared/Enum/EnumToStringx.cpp

    r24947 r24959  
    304304                case RootPathEnum : return "RootPath";
    305305                case SaveResultsEnum : return "SaveResults";
    306                 case SealevelEarthAreaEnum : return "SealevelEarthArea";
     306                case SealevelPlanetRadiusEnum : return "SealevelPlanetRadius";
     307                case SealevelPlanetAreaEnum : return "SealevelPlanetArea";
    307308                case SealevelEustaticEnum : return "SealevelEustatic";
    308309                case SealevelriseAbstolEnum : return "SealevelriseAbstol";

  • issm/trunk-jpl/src/c/shared/Enum/StringToEnumx.cpp

    r24947 r24959  
    310310              else if (strcmp(name,"RootPath")==0) return RootPathEnum;
    311311              else if (strcmp(name,"SaveResults")==0) return SaveResultsEnum;
    312               else if (strcmp(name,"SealevelEarthArea")==0) return SealevelEarthAreaEnum;
     312              else if (strcmp(name,"SealevelPlanetRadius")==0) return SealevelPlanetRadiusEnum;
     313              else if (strcmp(name,"SealevelPlanetArea")==0) return SealevelPlanetAreaEnum;
    313314              else if (strcmp(name,"SealevelEustatic")==0) return SealevelEustaticEnum;
    314315              else if (strcmp(name,"SealevelriseAbstol")==0) return SealevelriseAbstolEnum;
    … …  
    382383              else if (strcmp(name,"SmbNumRequestedOutputs")==0) return SmbNumRequestedOutputsEnum;
    383384              else if (strcmp(name,"SmbPfac")==0) return SmbPfacEnum;
    384               else if (strcmp(name,"SmbRdl")==0) return SmbRdlEnum;
    385385         else stage=4;
    386386   }
    387387   if(stage==4){
    388               if (strcmp(name,"SmbRequestedOutputs")==0) return SmbRequestedOutputsEnum;
     388              if (strcmp(name,"SmbRdl")==0) return SmbRdlEnum;
     389              else if (strcmp(name,"SmbRequestedOutputs")==0) return SmbRequestedOutputsEnum;
    389390              else if (strcmp(name,"SmbRlaps")==0) return SmbRlapsEnum;
    390391              else if (strcmp(name,"SmbRlapslgm")==0) return SmbRlapslgmEnum;
    … …  
    505506              else if (strcmp(name,"BedSlopeY")==0) return BedSlopeYEnum;
    506507              else if (strcmp(name,"CalvingCalvingrate")==0) return CalvingCalvingrateEnum;
    507               else if (strcmp(name,"CalvingHabFraction")==0) return CalvingHabFractionEnum;
    508508         else stage=5;
    509509   }
    510510   if(stage==5){
    511               if (strcmp(name,"CalvingMeltingrate")==0) return CalvingMeltingrateEnum;
     511              if (strcmp(name,"CalvingHabFraction")==0) return CalvingHabFractionEnum;
     512              else if (strcmp(name,"CalvingMeltingrate")==0) return CalvingMeltingrateEnum;
    512513              else if (strcmp(name,"CalvingStressThresholdFloatingice")==0) return CalvingStressThresholdFloatingiceEnum;
    513514              else if (strcmp(name,"CalvingStressThresholdGroundedice")==0) return CalvingStressThresholdGroundediceEnum;
    … …  
    628629              else if (strcmp(name,"Ice")==0) return IceEnum;
    629630              else if (strcmp(name,"IceMaskNodeActivation")==0) return IceMaskNodeActivationEnum;
    630               else if (strcmp(name,"Input")==0) return InputEnum;
    631631         else stage=6;
    632632   }
    633633   if(stage==6){
    634               if (strcmp(name,"InversionCostFunctionsCoefficients")==0) return InversionCostFunctionsCoefficientsEnum;
     634              if (strcmp(name,"Input")==0) return InputEnum;
     635              else if (strcmp(name,"InversionCostFunctionsCoefficients")==0) return InversionCostFunctionsCoefficientsEnum;
    635636              else if (strcmp(name,"InversionSurfaceObs")==0) return InversionSurfaceObsEnum;
    636637              else if (strcmp(name,"InversionThicknessObs")==0) return InversionThicknessObsEnum;
    … …  
    751752              else if (strcmp(name,"SmbGsp")==0) return SmbGspEnum;
    752753              else if (strcmp(name,"SmbGspini")==0) return SmbGspiniEnum;
    753               else if (strcmp(name,"SmbHref")==0) return SmbHrefEnum;
    754754         else stage=7;
    755755   }
    756756   if(stage==7){
    757               if (strcmp(name,"SmbIsInitialized")==0) return SmbIsInitializedEnum;
     757              if (strcmp(name,"SmbHref")==0) return SmbHrefEnum;
     758              else if (strcmp(name,"SmbIsInitialized")==0) return SmbIsInitializedEnum;
    758759              else if (strcmp(name,"SmbMAdd")==0) return SmbMAddEnum;
    759760              else if (strcmp(name,"SmbMassBalanceClimate")==0) return SmbMassBalanceClimateEnum;
    … …  
    874875              else if (strcmp(name,"Waterheight")==0) return WaterheightEnum;
    875876              else if (strcmp(name,"WeightsSurfaceObservation")==0) return WeightsSurfaceObservationEnum;
    876               else if (strcmp(name,"Outputdefinition1")==0) return Outputdefinition1Enum;
    877877         else stage=8;
    878878   }
    879879   if(stage==8){
    880               if (strcmp(name,"Outputdefinition10")==0) return Outputdefinition10Enum;
     880              if (strcmp(name,"Outputdefinition1")==0) return Outputdefinition1Enum;
     881              else if (strcmp(name,"Outputdefinition10")==0) return Outputdefinition10Enum;
    881882              else if (strcmp(name,"Outputdefinition11")==0) return Outputdefinition11Enum;
    882883              else if (strcmp(name,"Outputdefinition12")==0) return Outputdefinition12Enum;
    … …  
    997998              else if (strcmp(name,"BalancethicknessSolution")==0) return BalancethicknessSolutionEnum;
    998999              else if (strcmp(name,"BalancethicknessSpcthickness")==0) return BalancethicknessSpcthicknessEnum;
    999               else if (strcmp(name,"BalancevelocityAnalysis")==0) return BalancevelocityAnalysisEnum;
    10001000         else stage=9;
    10011001   }
    10021002   if(stage==9){
    1003               if (strcmp(name,"BalancevelocitySolution")==0) return BalancevelocitySolutionEnum;
     1003              if (strcmp(name,"BalancevelocityAnalysis")==0) return BalancevelocityAnalysisEnum;
     1004              else if (strcmp(name,"BalancevelocitySolution")==0) return BalancevelocitySolutionEnum;
    10041005              else if (strcmp(name,"BasalforcingsIsmip6")==0) return BasalforcingsIsmip6Enum;
    10051006              else if (strcmp(name,"BasalforcingsPico")==0) return BasalforcingsPicoEnum;
    … …  
    11201121              else if (strcmp(name,"HydrologySolution")==0) return HydrologySolutionEnum;
    11211122              else if (strcmp(name,"Hydrologydc")==0) return HydrologydcEnum;
    1122               else if (strcmp(name,"Hydrologypism")==0) return HydrologypismEnum;
    11231123         else stage=10;
    11241124   }
    11251125   if(stage==10){
    1126               if (strcmp(name,"Hydrologyshakti")==0) return HydrologyshaktiEnum;
     1126              if (strcmp(name,"Hydrologypism")==0) return HydrologypismEnum;
     1127              else if (strcmp(name,"Hydrologyshakti")==0) return HydrologyshaktiEnum;
    11271128              else if (strcmp(name,"Hydrologyshreve")==0) return HydrologyshreveEnum;
    11281129              else if (strcmp(name,"IceMass")==0) return IceMassEnum;
    … …  
    12431244              else if (strcmp(name,"P2bubblecondensed")==0) return P2bubblecondensedEnum;
    12441245              else if (strcmp(name,"P2xP1")==0) return P2xP1Enum;
    1245               else if (strcmp(name,"P2xP4")==0) return P2xP4Enum;
    12461246         else stage=11;
    12471247   }
    12481248   if(stage==11){
    1249               if (strcmp(name,"Paterson")==0) return PatersonEnum;
     1249              if (strcmp(name,"P2xP4")==0) return P2xP4Enum;
     1250              else if (strcmp(name,"Paterson")==0) return PatersonEnum;
    12501251              else if (strcmp(name,"Pengrid")==0) return PengridEnum;
    12511252              else if (strcmp(name,"Penpair")==0) return PenpairEnum;
    … …  
    13661367              else if (strcmp(name,"StressbalanceViscosityOvershoot")==0) return StressbalanceViscosityOvershootEnum;
    13671368              else if (strcmp(name,"SubelementMigration4")==0) return SubelementMigration4Enum;
    1368               else if (strcmp(name,"TimesteppingTimeAdapt")==0) return TimesteppingTimeAdaptEnum;
    13691369         else stage=12;
    13701370   }
    13711371   if(stage==12){
    1372               if (strcmp(name,"TriangleInterp")==0) return TriangleInterpEnum;
     1372              if (strcmp(name,"TimesteppingTimeAdapt")==0) return TimesteppingTimeAdaptEnum;
     1373              else if (strcmp(name,"TriangleInterp")==0) return TriangleInterpEnum;
    13731374              else if (strcmp(name,"MaximumNumberOfDefinitions")==0) return MaximumNumberOfDefinitionsEnum;
    13741375         else stage=13;
Note: See TracChangeset for help on using the changeset viewer.