Changeset 5747


Ignore:
Timestamp:
09/10/10 11:45:33 (15 years ago)
Author:
Mathieu Morlighem
Message:

Cleaned elements by taking out Inputs declaration and added assertion for all inputs if they are used

Location:
issm/trunk/src/c/objects/Elements
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/objects/Elements/Penta.cpp

    r5746 r5747  
    556556
    557557        /*Retrieve all inputs we will be needing: */
    558         Input* pressure_input=inputs->GetInput(PressureEnum);
    559         Input* vx_input=inputs->GetInput(VxEnum);
    560         Input* vy_input=inputs->GetInput(VyEnum);
    561         Input* vz_input=inputs->GetInput(VzEnum);
     558        Input* pressure_input=inputs->GetInput(PressureEnum); ISSMASSERT(pressure_input);
     559        Input* vx_input=inputs->GetInput(VxEnum);             ISSMASSERT(vx_input);
     560        Input* vy_input=inputs->GetInput(VyEnum);             ISSMASSERT(vy_input);
     561        Input* vz_input=inputs->GetInput(VzEnum);             ISSMASSERT(vz_input);
    562562
    563563        /* Start  looping on the number of gaussian points: */
    … …  
    20852085        Penta* pentabase=NULL;
    20862086
    2087         /*inputs: */
    2088         Input* vx_input=NULL;
    2089         Input* vy_input=NULL;
    2090         Input* vxold_input=NULL;
    2091         Input* vyold_input=NULL;
    2092 
    20932087        /*retrieve some parameters: */
    20942088        this->parameters->FindParam(&viscosity_overshoot,ViscosityOvershootEnum);
    … …  
    21072101
    21082102        /*Retrieve all inputs we will be needing: */
    2109         vx_input=inputs->GetInput(VxEnum);
    2110         vy_input=inputs->GetInput(VyEnum);
    2111         vxold_input=inputs->GetInput(VxOldEnum);
    2112         vyold_input=inputs->GetInput(VyOldEnum);
     2103        Input* vx_input=inputs->GetInput(VxEnum);       ISSMASSERT(vx_input);
     2104        Input* vy_input=inputs->GetInput(VyEnum);       ISSMASSERT(vy_input);
     2105        Input* vxold_input=inputs->GetInput(VxOldEnum); ISSMASSERT(vxold_input);
     2106        Input* vyold_input=inputs->GetInput(VyOldEnum); ISSMASSERT(vyold_input);
    21132107
    21142108        /* Start  looping on the number of gaussian points: */
    … …  
    23012295        Tria*  tria=NULL;
    23022296        Penta* pentabase=NULL;
    2303 
    2304         /*inputs: */
    23052297        int  approximation;
    2306         Input* vx_input=NULL;
    2307         Input* vy_input=NULL;
    2308         Input* vxold_input=NULL;
    2309         Input* vyold_input=NULL;
    23102298
    23112299        inputs->GetParameterValue(&approximation,ApproximationEnum);
    … …  
    23592347
    23602348                /*Retrieve all inputs we will be needing: */
    2361                 vx_input=inputs->GetInput(VxEnum);
    2362                 vy_input=inputs->GetInput(VyEnum);
    2363                 vxold_input=inputs->GetInput(VxOldEnum);
    2364                 vyold_input=inputs->GetInput(VyOldEnum);
     2349                Input* vx_input=inputs->GetInput(VxEnum);       ISSMASSERT(vx_input);
     2350                Input* vy_input=inputs->GetInput(VyEnum);       ISSMASSERT(vy_input);
     2351                Input* vxold_input=inputs->GetInput(VxOldEnum); ISSMASSERT(vxold_input);
     2352                Input* vyold_input=inputs->GetInput(VyOldEnum); ISSMASSERT(vyold_input);
    23652353
    23662354                /* Start  looping on the number of gaussian points: */
    … …  
    24792467        Tria*  tria=NULL;
    24802468
    2481         /*inputs: */
    2482         Input* vx_input=NULL;
    2483         Input* vy_input=NULL;
    2484         Input* vxold_input=NULL;
    2485         Input* vyold_input=NULL;
    2486 
    24872469        /*retrieve some parameters: */
    24882470        this->parameters->FindParam(&viscosity_overshoot,ViscosityOvershootEnum);
    … …  
    24982480
    24992481        /*Retrieve all inputs we will be needing: */
    2500         vx_input=inputs->GetInput(VxEnum);
    2501         vy_input=inputs->GetInput(VyEnum);
    2502         vxold_input=inputs->GetInput(VxOldEnum);
    2503         vyold_input=inputs->GetInput(VyOldEnum);
     2482        Input* vx_input=inputs->GetInput(VxEnum);       ISSMASSERT(vx_input);
     2483        Input* vy_input=inputs->GetInput(VyEnum);       ISSMASSERT(vy_input);
     2484        Input* vxold_input=inputs->GetInput(VxOldEnum); ISSMASSERT(vxold_input);
     2485        Input* vyold_input=inputs->GetInput(VyOldEnum); ISSMASSERT(vyold_input);
    25042486
    25052487        /* Start  looping on the number of gaussian points: */
    … …  
    26182600
    26192601        /*Retrieve all inputs we will be needing: */
    2620         Input* vx_input=inputs->GetInput(VxEnum);
    2621         Input* vy_input=inputs->GetInput(VyEnum);
    2622         Input* vz_input=inputs->GetInput(VzEnum);
     2602        Input* vx_input=inputs->GetInput(VxEnum); ISSMASSERT(vx_input);
     2603        Input* vy_input=inputs->GetInput(VyEnum); ISSMASSERT(vy_input);
     2604        Input* vz_input=inputs->GetInput(VzEnum); ISSMASSERT(vz_input);
    26232605
    26242606        /* Start  looping on the number of gaussian points: */
    … …  
    29342916        Tria*  tria=NULL;
    29352917
    2936         /*inputs: */
    2937         Input* vx_input=NULL;
    2938         Input* vy_input=NULL;
    2939         Input* vz_input=NULL;
    2940 
    29412918        /*If on water, skip: */
    29422919        if(IsOnWater())return;
    … …  
    29582935        this->parameters->FindParam(&epsvel,EpsVelEnum);
    29592936
    2960         vx_input=inputs->GetInput(VxEnum);
    2961         vy_input=inputs->GetInput(VyEnum);
    2962         vz_input=inputs->GetInput(VzEnum);
     2937        Input* vx_input=inputs->GetInput(VxEnum); ISSMASSERT(vx_input);
     2938        Input* vy_input=inputs->GetInput(VyEnum); ISSMASSERT(vy_input);
     2939        Input* vz_input=inputs->GetInput(VzEnum); ISSMASSERT(vz_input);
    29632940
    29642941        /* Start  looping on the number of gaussian points: */
    … …  
    32113188
    32123189        //recover extra inputs
    3213         Input* thickness_input=inputs->GetInput(ThicknessEnum);
    3214         Input* surface_input=inputs->GetInput(SurfaceEnum);
    3215         Input* slopex_input=inputs->GetInput(SurfaceSlopeXEnum);
    3216         Input* slopey_input=inputs->GetInput(SurfaceSlopeYEnum);
     3190        Input* thickness_input=inputs->GetInput(ThicknessEnum);  ISSMASSERT(thickness_input);
     3191        Input* surface_input=inputs->GetInput(SurfaceEnum);      ISSMASSERT(surface_input);
     3192        Input* slopex_input=inputs->GetInput(SurfaceSlopeXEnum); ISSMASSERT(slopex_input);
     3193        Input* slopey_input=inputs->GetInput(SurfaceSlopeYEnum); ISSMASSERT(slopey_input);
    32173194
    32183195        GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
    … …  
    33343311        double  pe_g_gaussian[numdof];
    33353312
    3336         /*inputs: */
    3337         Input* surface_input=NULL;
    3338         Input* thickness_input=NULL;
    3339 
    33403313        /*If on water, skip load: */
    33413314        if(IsOnWater())return;
    … …  
    33483321
    33493322        /*Retrieve all inputs we will be needing: */
    3350         thickness_input=inputs->GetInput(ThicknessEnum);
    3351         surface_input=inputs->GetInput(SurfaceEnum);
     3323        Input* thickness_input=inputs->GetInput(ThicknessEnum); ISSMASSERT(thickness_input);
     3324        Input* surface_input=inputs->GetInput(SurfaceEnum);     ISSMASSERT(surface_input);
    33523325
    33533326        /* Start  looping on the number of gaussian points: */
    … …  
    34403413        /*inputs: */
    34413414        int  approximation;
    3442         Input* vx_input=NULL;
    3443         Input* vy_input=NULL;
    3444         Input* vz_input=NULL;
    3445         Input* bed_input=NULL;
    34463415
    34473416        /*retrieve inputs :*/
    … …  
    34643433
    34653434        /*Retrieve all inputs we will be needing: */
    3466         vx_input=inputs->GetInput(VxEnum);
    3467         vy_input=inputs->GetInput(VyEnum);
    3468         vz_input=inputs->GetInput(VzEnum);
    3469         bed_input=inputs->GetInput(BedEnum);
     3435        Input* vx_input=inputs->GetInput(VxEnum);   ISSMASSERT(vx_input);
     3436        Input* vy_input=inputs->GetInput(VyEnum);   ISSMASSERT(vy_input);
     3437        Input* vz_input=inputs->GetInput(VzEnum);   ISSMASSERT(vz_input);
     3438        Input* bed_input=inputs->GetInput(BedEnum); ISSMASSERT(bed_input);
    34703439
    34713440        /* Start  looping on the number of gaussian points: */
    … …  
    36103579        double dv[3];
    36113580        double dudx,dvdy;
    3612 
    3613         /*inputs: */
    3614         Input* vx_input=NULL;
    3615         Input* vy_input=NULL;
    36163581
    36173582        /*If on water, skip: */
    … …  
    36323597
    36333598        /*Retrieve all inputs we will be needing: */
    3634         vx_input=inputs->GetInput(VxEnum);
    3635         vy_input=inputs->GetInput(VyEnum);
     3599        Input* vx_input=inputs->GetInput(VxEnum); ISSMASSERT(vx_input);
     3600        Input* vy_input=inputs->GetInput(VyEnum); ISSMASSERT(vy_input);
    36363601
    36373602        /* Start  looping on the number of gaussian points: */
    … …  
    37703735        /*Collapsed formulation: */
    37713736        Tria*  tria=NULL;
    3772 
    3773         /*parameters: */
    37743737        double dt;
    3775 
    3776         /*inputs: */
    3777         Input* vx_input=NULL;
    3778         Input* vy_input=NULL;
    3779         Input* vz_input=NULL;
    3780         Input* temperature_input=NULL;
    37813738
    37823739        /*retrieve some parameters: */
    … …  
    37993756
    38003757        /*Retrieve all inputs we will be needing: */
    3801         vx_input=inputs->GetInput(VxEnum);
    3802         vy_input=inputs->GetInput(VyEnum);
    3803         vz_input=inputs->GetInput(VzEnum);
    3804         if (dt) temperature_input=inputs->GetInput(TemperatureEnum);
     3758        Input* vx_input=inputs->GetInput(VxEnum); ISSMASSERT(vx_input);
     3759        Input* vy_input=inputs->GetInput(VyEnum); ISSMASSERT(vy_input);
     3760        Input* vz_input=inputs->GetInput(VzEnum); ISSMASSERT(vz_input);
     3761        Input* temperature_input=NULL;
     3762        if (dt){
     3763                temperature_input=inputs->GetInput(TemperatureEnum); ISSMASSERT(inputs);
     3764        }
    38053765
    38063766        /* Start  looping on the number of gaussian points: */
    … …  
    45034463        double       xyz_list[NUMVERTICES][3];
    45044464
    4505         Input  *vz_input        = NULL;
    45064465        double *vz_ptr          = NULL;
    45074466        Penta  *penta          = NULL;
    … …  
    45374496
    45384497                /*Now Compute vel*/
    4539                 vz_input=inputs->GetInput(VzEnum);
     4498                Input* vz_input=inputs->GetInput(VzEnum);
    45404499                if (vz_input){
    45414500                        if (vz_input->Enum()!=PentaVertexInputEnum) ISSMERROR("Cannot compute Vel as Vz is of type %s",EnumToString(vz_input->Enum()));
    … …  
    45954554        double       xyz_list[NUMVERTICES][3];
    45964555
    4597         Input  *vz_input       = NULL;
    45984556        double *vz_ptr         = NULL;
    45994557        Penta  *penta          = NULL;
    … …  
    46274585
    46284586        /*Get Vz*/
    4629         vz_input=inputs->GetInput(VzEnum);
     4587        Input* vz_input=inputs->GetInput(VzEnum);
    46304588        if (vz_input){
    46314589                if (vz_input->Enum()!=PentaVertexInputEnum){
    … …  
    46864644        double       rho_ice,g;
    46874645        double       xyz_list[NUMVERTICES][3];
    4688        
    4689         Input  *vz_input        = NULL;
    46904646        double *vz_ptr          = NULL;
    46914647
    … …  
    47084664
    47094665        /*Get Vz*/
    4710         vz_input=inputs->GetInput(VzEnum);
     4666        Input* vz_input=inputs->GetInput(VzEnum);
    47114667        if (vz_input){
    47124668                if (vz_input->Enum()!=PentaVertexInputEnum){
    … …  
    47654721        double       rho_ice,g;
    47664722        double       xyz_list[NUMVERTICES][3];
    4767        
    4768         Input*       vz_input=NULL;
    47694723        double*      vz_ptr=NULL;
    47704724
    … …  
    47874741
    47884742        /*Get Vz*/
    4789         vz_input=inputs->GetInput(VzEnum);
     4743        Input* vz_input=inputs->GetInput(VzEnum);
    47904744        if (vz_input){
    47914745                if (vz_input->Enum()!=PentaVertexInputEnum){
    … …  
    48454799        double       xyz_list[NUMVERTICES][3];
    48464800
    4847         Input*       vx_input=NULL;
    48484801        double*      vx_ptr=NULL;
    4849         Input*       vy_input=NULL;
    48504802        double*      vy_ptr=NULL;
    48514803
    … …  
    48674819
    48684820        /*Get Vx and Vy*/
    4869         vx_input=inputs->GetInput(VxEnum);
     4821        Input* vx_input=inputs->GetInput(VxEnum);
    48704822        if (vx_input){
    48714823                if (vx_input->Enum()!=PentaVertexInputEnum) ISSMERROR("Cannot compute Vel as Vx is of type %s",EnumToString(vx_input->Enum()));
    … …  
    48754827        else for(i=0;i<NUMVERTICES;i++) vx[i]=0.0;
    48764828
    4877         vy_input=inputs->GetInput(VyEnum);
     4829        Input* vy_input=inputs->GetInput(VyEnum);
    48784830        if (vy_input){
    48794831                if (vy_input->Enum()!=PentaVertexInputEnum) ISSMERROR("Cannot compute Vel as Vy is of type %s",EnumToString(vy_input->Enum()));

  • issm/trunk/src/c/objects/Elements/Tria.cpp

    r5746 r5747  
    647647        this->parameters->FindParam(&control_type,ControlTypeEnum);
    648648        this->parameters->FindParam(&cm_noisedmp,CmNoiseDmpEnum);
    649         Input* drag_input=inputs->GetInput(DragCoefficientEnum);
    650         Input* dhdt_input=inputs->GetInput(DhDtEnum);
    651         Input* Bbar_input=matice->inputs->GetInput(RheologyBbarEnum);
     649        Input* drag_input=inputs->GetInput(DragCoefficientEnum);      ISSMASSERT(drag_input);
     650        Input* dhdt_input=inputs->GetInput(DhDtEnum);                 ISSMASSERT(dhdt_input);
     651        Input* Bbar_input=matice->inputs->GetInput(RheologyBbarEnum); ISSMASSERT(Bbar_input);
    652652
    653653        /*If on water, return 0: */
    … …  
    10051005        double epsilon[3]; /* epsilon=[exx,eyy,exy];*/
    10061006
    1007         /*inputs: */
    1008         Input* vx_input=NULL;
    1009         Input* vy_input=NULL;
    1010         Input* adjointx_input=NULL;
    1011         Input* adjointy_input=NULL;
    1012         Input* dragcoefficient_input=NULL;
    1013 
    10141007        /*parameters: */
    10151008        double  cm_noisedmp;
    1016 
    10171009        int analysis_type;
    10181010
    … …  
    10351027
    10361028        /*Retrieve all inputs we will be needing: */
    1037         adjointx_input=inputs->GetInput(AdjointxEnum);
    1038         adjointy_input=inputs->GetInput(AdjointyEnum);
    1039         vx_input=inputs->GetInput(VxEnum);
    1040         vy_input=inputs->GetInput(VyEnum);
    1041         dragcoefficient_input=inputs->GetInput(DragCoefficientEnum);
     1029        Input* adjointx_input=inputs->GetInput(AdjointxEnum);               ISSMASSERT(adjointx_input);
     1030        Input* adjointy_input=inputs->GetInput(AdjointyEnum);               ISSMASSERT(adjointy_input);
     1031        Input* vx_input=inputs->GetInput(VxEnum);                           ISSMASSERT(vx_input);
     1032        Input* vy_input=inputs->GetInput(VyEnum);                           ISSMASSERT(vy_input);
     1033        Input* dragcoefficient_input=inputs->GetInput(DragCoefficientEnum); ISSMASSERT(dragcoefficient_input);
    10421034
    10431035        /* Start  looping on the number of gaussian points: */
    … …  
    15441536        Input* thickness_input   =inputs->GetInput(ThicknessEnum);   ISSMASSERT(thickness_input);
    15451537        Input* thicknessobs_input=inputs->GetInput(ThicknessObsEnum);ISSMASSERT(thicknessobs_input);
    1546         Input* weights_input     =inputs->GetInput(WeightsEnum);ISSMASSERT(weights_input);
     1538        Input* weights_input     =inputs->GetInput(WeightsEnum);     ISSMASSERT(weights_input);
    15471539
    15481540        /* Start  looping on the number of gaussian points: */
    … …  
    24552447        /*Retrieve all inputs we will be needed*/
    24562448        if(dim==2){
    2457                 vxaverage_input=inputs->GetInput(VxEnum);
    2458                 vyaverage_input=inputs->GetInput(VyEnum);
     2449                vxaverage_input=inputs->GetInput(VxEnum); ISSMASSERT(vxaverage_input);
     2450                vyaverage_input=inputs->GetInput(VyEnum); ISSMASSERT(vyaverage_input);
    24592451        }
    24602452        else{
    2461                 vxaverage_input=inputs->GetInput(VxAverageEnum);
    2462                 vyaverage_input=inputs->GetInput(VyAverageEnum);
     2453                vxaverage_input=inputs->GetInput(VxAverageEnum); ISSMASSERT(vxaverage_input);
     2454                vyaverage_input=inputs->GetInput(VyAverageEnum); ISSMASSERT(vyaverage_input);
    24632455        }
    24642456
    … …  
    25832575        int     dim;
    25842576
    2585         /*inputs: */
    2586         Input* vx_input=NULL;
    2587         Input* vy_input=NULL;
    2588 
    25892577        /* Get node coordinates and dof list: */
    25902578        GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
    … …  
    25932581
    25942582        /*Retrieve all inputs we will be needed*/
    2595         if(dim==2){
    2596                 vx_input=inputs->GetInput(VxEnum);
    2597                 vy_input=inputs->GetInput(VyEnum);
    2598         }
     2583        Input* vx_input=inputs->GetInput(VxEnum); ISSMASSERT(vx_input);
     2584        Input* vy_input=inputs->GetInput(VyEnum); ISSMASSERT(vy_input);
    25992585
    26002586        /* Start  looping on the number of gaussian points: */
    … …  
    26852671        Input* vxaverage_input=NULL;
    26862672        Input* vyaverage_input=NULL;
    2687         Input* surface_input=NULL;
    26882673        bool artdiff;
    26892674
    … …  
    26972682
    26982683        /*Retrieve all inputs we will be needed*/
    2699         surface_input=inputs->GetInput(SurfaceEnum);
     2684        Input* surface_input=inputs->GetInput(SurfaceEnum); ISSMASSERT(surface_input);
    27002685        if(dim==2){
    2701                 vxaverage_input=inputs->GetInput(VxEnum);
    2702                 vyaverage_input=inputs->GetInput(VyEnum);
     2686                vxaverage_input=inputs->GetInput(VxEnum); ISSMASSERT(vxaverage_input);
     2687                vyaverage_input=inputs->GetInput(VyEnum); ISSMASSERT(vyaverage_input);
    27032688        }
    27042689        else{
    2705                 vxaverage_input=inputs->GetInput(VxAverageEnum);
    2706                 vyaverage_input=inputs->GetInput(VyAverageEnum);
     2690                vxaverage_input=inputs->GetInput(VxAverageEnum); ISSMASSERT(vxaverage_input);
     2691                vyaverage_input=inputs->GetInput(VyAverageEnum); ISSMASSERT(vyaverage_input);
    27072692        }
    27082693
    … …  
    28492834        double  thickness;
    28502835
    2851         /*inputs: */
    2852         Input* thickness_input=NULL;
    2853         Input* vx_input=NULL;
    2854         Input* vy_input=NULL;
    2855         Input* vxold_input=NULL;
    2856         Input* vyold_input=NULL;
    2857 
    28582836        /*retrieve some parameters: */
    28592837        this->parameters->FindParam(&viscosity_overshoot,ViscosityOvershootEnum);
    … …  
    28672845
    28682846        /*Retrieve all inputs we will be needing: */
    2869         thickness_input=inputs->GetInput(ThicknessEnum);
    2870         vx_input=inputs->GetInput(VxEnum);
    2871         vy_input=inputs->GetInput(VyEnum);
    2872         vxold_input=inputs->GetInput(VxOldEnum);
    2873         vyold_input=inputs->GetInput(VyOldEnum);
     2847        Input* thickness_input=inputs->GetInput(ThicknessEnum); ISSMASSERT(thickness_input);
     2848        Input* vx_input=inputs->GetInput(VxEnum);               ISSMASSERT(vx_input);
     2849        Input* vy_input=inputs->GetInput(VyEnum);               ISSMASSERT(vy_input);
     2850        Input* vxold_input=inputs->GetInput(VxOldEnum);         ISSMASSERT(vxold_input);
     2851        Input* vyold_input=inputs->GetInput(VyOldEnum);         ISSMASSERT(vyold_input);
    28742852
    28752853        /* Start  looping on the number of gaussian points: */
    … …  
    29732951        /*inputs: */
    29742952        int  drag_type;
    2975         Input* surface_input=NULL;
    2976         Input* vx_input=NULL;
    2977         Input* vy_input=NULL;
    2978         Input* vz_input=NULL;
    29792953
    29802954        /*retrive parameters: */
    … …  
    29832957        /*retrieve inputs :*/
    29842958        inputs->GetParameterValue(&drag_type,DragTypeEnum);
    2985         surface_input=inputs->GetInput(SurfaceEnum);
    2986         vx_input=inputs->GetInput(VxEnum);
    2987         vy_input=inputs->GetInput(VyEnum);
    2988         vz_input=inputs->GetInput(VzEnum);
     2959        Input* surface_input=inputs->GetInput(SurfaceEnum); ISSMASSERT(surface_input);
     2960        Input* vx_input=inputs->GetInput(VxEnum);           ISSMASSERT(vx_input);
     2961        Input* vy_input=inputs->GetInput(VyEnum);           ISSMASSERT(vy_input);
     2962        Input* vz_input=inputs->GetInput(VzEnum);           ISSMASSERT(vz_input);
    29892963       
    29902964        /* Get node coordinates and dof list: */
    … …  
    30943068        /*inputs: */
    30953069        int  drag_type;
    3096         Input* surface_input=NULL;
    3097         Input* vx_input=NULL;
    3098         Input* vy_input=NULL;
    3099         Input* vz_input=NULL;
    31003070
    31013071        /*retrive parameters: */
    … …  
    31043074        /*retrieve inputs :*/
    31053075        inputs->GetParameterValue(&drag_type,DragTypeEnum);
    3106         surface_input=inputs->GetInput(SurfaceEnum);
    3107         vx_input=inputs->GetInput(VxEnum);
    3108         vy_input=inputs->GetInput(VyEnum);
    3109         vz_input=inputs->GetInput(VzEnum);
     3076        Input* surface_input=inputs->GetInput(SurfaceEnum); ISSMASSERT(surface_input);
     3077        Input* vx_input=inputs->GetInput(VxEnum);           ISSMASSERT(vx_input);
     3078        Input* vy_input=inputs->GetInput(VyEnum);           ISSMASSERT(vy_input);
     3079        Input* vz_input=inputs->GetInput(VzEnum);           ISSMASSERT(vz_input);
    31103080       
    31113081        /* Get node coordinates and dof list: */
    … …  
    32123182        /*inputs: */
    32133183        int  drag_type;
    3214         Input* surface_input=NULL;
    3215         Input* vx_input=NULL;
    3216         Input* vy_input=NULL;
    3217         Input* vz_input=NULL;
    32183184
    32193185        /*retrive parameters: */
    … …  
    32223188        /*retrieve inputs :*/
    32233189        inputs->GetParameterValue(&drag_type,DragTypeEnum);
    3224         surface_input=inputs->GetInput(SurfaceEnum);
    3225         vx_input=inputs->GetInput(VxEnum);
    3226         vy_input=inputs->GetInput(VyEnum);
    3227         vz_input=inputs->GetInput(VzEnum);
     3190        Input* surface_input=inputs->GetInput(SurfaceEnum); ISSMASSERT(surface_input);
     3191        Input* vx_input=inputs->GetInput(VxEnum);           ISSMASSERT(vx_input);
     3192        Input* vy_input=inputs->GetInput(VyEnum);           ISSMASSERT(vy_input);
     3193        Input* vz_input=inputs->GetInput(VzEnum);           ISSMASSERT(vz_input);
    32283194       
    32293195        /* Get node coordinates and dof list: */
    … …  
    35443510        /*Retrieve all inputs we will be needing: */
    35453511        if(dim==2){
    3546                 vxaverage_input=inputs->GetInput(VxEnum);
    3547                 vyaverage_input=inputs->GetInput(VyEnum);
     3512                vxaverage_input=inputs->GetInput(VxEnum); ISSMASSERT(vxaverage_input);
     3513                vyaverage_input=inputs->GetInput(VyEnum); ISSMASSERT(vyaverage_input);
    35483514        }
    35493515        else{
    3550                 vxaverage_input=inputs->GetInput(VxAverageEnum);
    3551                 vyaverage_input=inputs->GetInput(VyAverageEnum);
     3516                vxaverage_input=inputs->GetInput(VxAverageEnum); ISSMASSERT(vxaverage_input);
     3517                vyaverage_input=inputs->GetInput(VyAverageEnum); ISSMASSERT(vyaverage_input);
    35523518        }
    35533519
    … …  
    37033669        /*Retrieve all inputs we will be needed*/
    37043670        if(dim==2){
    3705                 vxaverage_input=inputs->GetInput(VxEnum);
    3706                 vyaverage_input=inputs->GetInput(VyEnum);
     3671                vxaverage_input=inputs->GetInput(VxEnum); ISSMASSERT(vxaverage_input);
     3672                vyaverage_input=inputs->GetInput(VyEnum); ISSMASSERT(vyaverage_input);
    37073673        }
    37083674        else{
    3709                 vxaverage_input=inputs->GetInput(VxAverageEnum);
    3710                 vyaverage_input=inputs->GetInput(VyAverageEnum);
     3675                vxaverage_input=inputs->GetInput(VxAverageEnum); ISSMASSERT(vxaverage_input);
     3676                vyaverage_input=inputs->GetInput(VyAverageEnum); ISSMASSERT(vyaverage_input);
    37113677        }
    37123678
    … …  
    39203886        double  dhdt_g;
    39213887
    3922         /*inputs: */
    3923         Input* accumulation_input=NULL;
    3924         Input* melting_input=NULL;
    3925         Input* dhdt_input=NULL;
    3926 
    39273888        /* Get node coordinates and dof list: */
    39283889        GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
    … …  
    39303891
    39313892        /*retrieve inputs :*/
    3932         accumulation_input=inputs->GetInput(AccumulationRateEnum); ISSMASSERT(accumulation_input);
    3933         melting_input=inputs->GetInput(MeltingRateEnum);           ISSMASSERT(melting_input);
    3934         dhdt_input=inputs->GetInput(DhDtEnum);                     ISSMASSERT(dhdt_input);
     3893        Input* accumulation_input=inputs->GetInput(AccumulationRateEnum); ISSMASSERT(accumulation_input);
     3894        Input* melting_input=inputs->GetInput(MeltingRateEnum);           ISSMASSERT(melting_input);
     3895        Input* dhdt_input=inputs->GetInput(DhDtEnum);                     ISSMASSERT(dhdt_input);
    39353896       
    39363897        /* Start  looping on the number of gaussian points: */
    … …  
    39903951        double  dhdt_g;
    39913952
    3992         /*inputs: */
    3993         Input* accumulation_input=NULL;
    3994         Input* melting_input=NULL;
    3995         Input* dhdt_input=NULL;
    3996 
    39973953        /* Get node coordinates and dof list: */
    39983954        GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
    … …  
    40003956
    40013957        /*Retrieve all inputs we will be needing: */
    4002         accumulation_input=inputs->GetInput(AccumulationRateEnum);
    4003         melting_input=inputs->GetInput(MeltingRateEnum);
    4004         dhdt_input=inputs->GetInput(DhDtEnum);
     3958        Input* accumulation_input=inputs->GetInput(AccumulationRateEnum); ISSMASSERT(accumulation_input);
     3959        Input* melting_input=inputs->GetInput(MeltingRateEnum);           ISSMASSERT(melting_input);
     3960        Input* dhdt_input=inputs->GetInput(DhDtEnum);                     ISSMASSERT(dhdt_input);
    40053961
    40063962        /* Start  looping on the number of gaussian points: */
    … …  
    40594015        double  melting_g;
    40604016
    4061         /*inputs: */
    4062         Input* accumulation_input=NULL;
    4063         Input* melting_input=NULL;
    4064 
    40654017        /* Get node coordinates and dof list: */
    40664018        GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
    … …  
    40684020
    40694021        /*Retrieve all inputs we will be needing: */
    4070         accumulation_input=inputs->GetInput(AccumulationRateEnum);
    4071         melting_input=inputs->GetInput(MeltingRateEnum);
     4022        Input* accumulation_input=inputs->GetInput(AccumulationRateEnum); ISSMASSERT(accumulation_input);
     4023        Input* melting_input=inputs->GetInput(MeltingRateEnum);           ISSMASSERT(melting_input);
    40724024
    40734025        /* Start  looping on the number of gaussian points: */
    … …  
    41334085        double  dbdx,dbdy;
    41344086
    4135         /*inputs: */
    4136         Input* melting_input=NULL;
    4137         Input* vx_input=NULL;
    4138         Input* vy_input=NULL;
    4139         Input* bed_input=NULL;
    4140 
    41414087        /* Get node coordinates and dof list: */
    41424088        GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
    … …  
    41444090
    41454091        /*Retrieve all inputs we will be needing: */
    4146         melting_input=inputs->GetInput(MeltingRateEnum);
    4147         vx_input=inputs->GetInput(VxEnum);
    4148         vy_input=inputs->GetInput(VyEnum);
    4149         bed_input=inputs->GetInput(BedEnum);
     4092        Input* melting_input=inputs->GetInput(MeltingRateEnum); ISSMASSERT(melting_input);
     4093        Input* vx_input=inputs->GetInput(VxEnum);               ISSMASSERT(vx_input);
     4094        Input* vy_input=inputs->GetInput(VyEnum);               ISSMASSERT(vy_input);
     4095        Input* bed_input=inputs->GetInput(BedEnum);             ISSMASSERT(bed_input);
    41504096
    41514097        /*For icesheets: */
    … …  
    43084254        double l1l2l3[3];
    43094255
    4310         /*inputs: */
    4311         Input* thickness_input=NULL;
    4312         Input* thicknessobs_input=NULL;
    4313         Input* weights_input=NULL;
    4314 
    43154256        /* Get node coordinates and dof list: */
    43164257        GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
    … …  
    43184259
    43194260        /*Retrieve all inputs we will be needing: */
    4320         thickness_input   =inputs->GetInput(ThicknessEnum);   ISSMASSERT(thickness_input);
    4321         thicknessobs_input=inputs->GetInput(ThicknessObsEnum);ISSMASSERT(thicknessobs_input);
    4322         weights_input     =inputs->GetInput(WeightsEnum);ISSMASSERT(weights_input);
     4261        Input* thickness_input   =inputs->GetInput(ThicknessEnum);   ISSMASSERT(thickness_input);
     4262        Input* thicknessobs_input=inputs->GetInput(ThicknessObsEnum);ISSMASSERT(thicknessobs_input);
     4263        Input* weights_input     =inputs->GetInput(WeightsEnum);     ISSMASSERT(weights_input);
    43234264
    43244265        /* Start  looping on the number of gaussian points: */
    … …  
    48084749        double    slope2;
    48094750        int       connectivity;
    4810 
    4811         /*flags: */
    4812         Input* slopex_input=NULL;
    4813         Input* slopey_input=NULL;
    4814         Input* thickness_input=NULL;
    48154751        GaussTria* gauss=NULL;
    48164752
    … …  
    48274763
    48284764        /* Get slopes and thickness */
    4829         slopex_input=inputs->GetInput(SurfaceSlopeXEnum);
    4830         slopey_input=inputs->GetInput(SurfaceSlopeYEnum);
    4831         thickness_input=inputs->GetInput(ThicknessEnum);
     4765        Input* slopex_input=inputs->GetInput(SurfaceSlopeXEnum); ISSMASSERT(slopex_input);
     4766        Input* slopey_input=inputs->GetInput(SurfaceSlopeYEnum); ISSMASSERT(slopey_input);
     4767        Input* thickness_input=inputs->GetInput(ThicknessEnum);  ISSMASSERT(thickness_input);
    48324768
    48334769        /*Spawn 3 sing elements: */
    … …  
    48884824        double  melting_g;
    48894825        double  thickness_g;
    4890 
    4891         /*inputs: */
    48924826        double  dt;
    4893         Input* accumulation_input=NULL;
    4894         Input* melting_input=NULL;
    4895         Input* thickness_input=NULL;
    48964827
    48974828        /*retrieve some parameters: */
    … …  
    49034834
    49044835        /*Retrieve all inputs we will be needing: */
    4905         accumulation_input=inputs->GetInput(AccumulationRateEnum);
    4906         melting_input=inputs->GetInput(MeltingRateEnum);
    4907         thickness_input=inputs->GetInput(ThicknessEnum);
     4836        Input* accumulation_input=inputs->GetInput(AccumulationRateEnum); ISSMASSERT(accumulation_input);
     4837        Input* melting_input=inputs->GetInput(MeltingRateEnum);           ISSMASSERT(melting_input);
     4838        Input* thickness_input=inputs->GetInput(ThicknessEnum);           ISSMASSERT(thickness_input);
    49084839
    49094840        /* Start  looping on the number of gaussian points: */
    … …  
    49644895        double  thickness_g;
    49654896        double  dt;
    4966         Input*  accumulation_input=NULL;
    4967         Input*  melting_input=NULL;
    4968         Input*  thickness_input=NULL;
    49694897
    49704898        /*retrieve some parameters: */
    … …  
    49764904
    49774905        /*Retrieve all inputs we will be needing: */
    4978         accumulation_input=inputs->GetInput(AccumulationRateEnum);
    4979         melting_input=inputs->GetInput(MeltingRateEnum);
    4980         thickness_input=inputs->GetInput(ThicknessEnum);
     4906        Input* accumulation_input=inputs->GetInput(AccumulationRateEnum); ISSMASSERT(accumulation_input);
     4907        Input* melting_input=inputs->GetInput(MeltingRateEnum);           ISSMASSERT(melting_input);
     4908        Input* thickness_input=inputs->GetInput(ThicknessEnum);           ISSMASSERT(thickness_input);
    49814909
    49824910        /* Start  looping on the number of gaussian points: */
    … …  
    50484976        /*Retrieve all inputs we will be needing: */
    50494977        if ( (analysis_type==SurfaceSlopeXAnalysisEnum) || (analysis_type==SurfaceSlopeYAnalysisEnum)){
    5050                 slope_input=inputs->GetInput(SurfaceEnum);
     4978                slope_input=inputs->GetInput(SurfaceEnum); ISSMASSERT(slope_input);
    50514979        }
    50524980        if ( (analysis_type==BedSlopeXAnalysisEnum) || (analysis_type==BedSlopeYAnalysisEnum)){
    5053                 slope_input=inputs->GetInput(BedEnum);
     4981                slope_input=inputs->GetInput(BedEnum);     ISSMASSERT(slope_input);
    50544982        }
    50554983               
    … …  
    51185046        double  P_terms[numdof]={0.0};
    51195047        double  l1l2l3[NUMVERTICES];
    5120 
    51215048        double  t_pmp;
    51225049        double  scalar_ocean;
    5123 
    5124         /*inputs: */
    5125         Input* pressure_input=NULL;
    51265050
    51275051        /* Get node coordinates and dof list: */
    … …  
    51445068
    51455069        /*Retrieve all inputs we will be needing: */
    5146         pressure_input=inputs->GetInput(PressureEnum);
     5070        Input* pressure_input=inputs->GetInput(PressureEnum); ISSMASSERT(pressure_input);
    51475071
    51485072        /* Start looping on the number of gauss 2d (nodes on the bedrock) */
    … …  
    52185142        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    52195143
    5220         /*inputs: */
    5221         Input* vx_input=NULL;
    5222         Input* vy_input=NULL;
    5223         Input* vz_input=NULL;
    5224         Input* geothermalflux_input=NULL;
    5225        
    52265144        /* Get node coordinates and dof list: */
    52275145        GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
    … …  
    52415159
    52425160        /*Retrieve all inputs we will be needing: */
    5243         vx_input=inputs->GetInput(VxEnum);
    5244         vy_input=inputs->GetInput(VyEnum);
    5245         vz_input=inputs->GetInput(VzEnum);
    5246         geothermalflux_input=inputs->GetInput(GeothermalFluxEnum);
     5161        Input* vx_input=inputs->GetInput(VxEnum);                         ISSMASSERT(vx_input);
     5162        Input* vy_input=inputs->GetInput(VyEnum);                         ISSMASSERT(vy_input);
     5163        Input* vz_input=inputs->GetInput(VzEnum);                         ISSMASSERT(vz_input);
     5164        Input* geothermalflux_input=inputs->GetInput(GeothermalFluxEnum); ISSMASSERT(geothermalflux_input);
    52475165
    52485166        /* Ice/ocean heat exchange flux on ice shelf base */
Note: See TracChangeset for help on using the changeset viewer.