Changeset 1104


Ignore:
Timestamp:
06/25/09 17:42:16 (16 years ago)
Author:
Eric.Larour
Message:

Transformed some Tria and Penta fields in booleans

Location:
issm/trunk/src/c
Files:
17 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp

    r1025 r1104  
    7171        double tria_viscosity_overshoot;
    7272        int    tria_artdiff;
     73        bool   tria_onwater;
    7374       
    7475        /*matice constructor input: */
    … …  
    231232                ModelFetchData((void**)&model->q,NULL,NULL,model_handle,"q","Matrix","Mat");
    232233                ModelFetchData((void**)&model->elementoniceshelf,NULL,NULL,model_handle,"elementoniceshelf","Matrix","Mat");
     234                ModelFetchData((void**)&model->elementonwater,NULL,NULL,model_handle,"elementonwater","Matrix","Mat");
    233235                ModelFetchData((void**)&model->elements_type,NULL,NULL,model_handle,"elements_type","Matrix","Mat");
    234236                ModelFetchData((void**)&model->B,NULL,NULL,model_handle,"B","Matrix","Mat");
    … …  
    288290                                tria_accumulation[2]=*(model->accumulation+        ((int)*(model->elements+elements_width*i+2)-1));
    289291
    290                                 /*element on iceshelf?:*/
     292                                /*element on iceshelf, water?:*/
    291293                                tria_shelf=(int)*(model->elementoniceshelf+i);
     294                                tria_onwater=(bool)*(model->elementonwater+i);
    292295
    293296                                tria_meanvel=model->meanvel;
    … …  
    298301
    299302                                /*Create tria element using its constructor:*/
    300                                 tria=new Tria(tria_id, tria_mid, tria_mparid, tria_g, tria_h, tria_s, tria_b, tria_k, tria_melting,tria_accumulation,tria_geothermalflux,tria_friction_type, tria_p, tria_q, tria_shelf, tria_meanvel, tria_epsvel, tria_viscosity_overshoot,tria_artdiff);
     303                                tria=new Tria(tria_id, tria_mid, tria_mparid, tria_g, tria_h, tria_s, tria_b, tria_k, tria_melting,tria_accumulation,tria_geothermalflux,tria_friction_type, tria_p, tria_q, tria_shelf, tria_meanvel, tria_epsvel, tria_viscosity_overshoot,tria_artdiff,tria_onwater);
    301304
    302305                                /*Add tria element to elements dataset: */
    … …  
    349352                xfree((void**)&model->q);
    350353                xfree((void**)&model->elementoniceshelf);
     354                xfree((void**)&model->elementonwater);
    351355                xfree((void**)&model->B);
    352356                xfree((void**)&model->n);
    … …  
    425429                                                penta_p,penta_q,penta_shelf,penta_onbed,penta_onsurface,penta_meanvel,penta_epsvel,
    426430                                                penta_collapse,penta_melting,penta_accumulation,penta_geothermalflux,penta_artdiff,
    427                                                 penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning);
     431                                                penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning,0);
    428432
    429433                                /*Add penta element to elements dataset: */

  • issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateElementsNodesAndMaterialsDiagnosticHutter.cpp

    r1025 r1104  
    7373        double beam_b[2];
    7474        double beam_k[2];
    75         int    beam_onbed;
     75        bool    beam_onbed;
    7676                                       
    7777        /*matpar constructor input: */
    … …  
    291291                                        beam_k[1]=model->drag[(int)(model->uppernodes[i]-1)];
    292292
    293                                         beam_onbed=(int)model->gridonbed[i];
     293                                        beam_onbed=(bool)model->gridonbed[i];
    294294
    295295                                        /*Create beam element ubeam its constructor:*/

  • issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateElementsNodesAndMaterialsDiagnosticStokes.cpp

    r1025 r1104  
    270270                                                penta_p,penta_q,penta_shelf,penta_onbed,penta_onsurface,penta_meanvel,penta_epsvel,
    271271                                                penta_collapse,penta_melting,penta_accumulation,penta_geothermalflux,penta_artdiff,
    272                                                 penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning);
     272                                                penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning,0);
    273273
    274274                                /*Add penta element to elements dataset: */

  • issm/trunk/src/c/ModelProcessorx/DiagnosticVert/CreateElementsNodesAndMaterialsDiagnosticVert.cpp

    r1025 r1104  
    238238                                penta_p,penta_q,penta_shelf,penta_onbed,penta_onsurface,penta_meanvel,penta_epsvel,
    239239                                penta_collapse,penta_melting,penta_accumulation,penta_geothermalflux,penta_artdiff,
    240                                 penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning);
     240                                penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning,0);
    241241
    242242                /*Add penta element to elements dataset: */

  • issm/trunk/src/c/ModelProcessorx/Melting/CreateElementsNodesAndMaterialsMelting.cpp

    r1025 r1104  
    242242                                penta_p,penta_q,penta_shelf,penta_onbed,penta_onsurface,penta_meanvel,penta_epsvel,
    243243                                penta_collapse,penta_melting,penta_accumulation,penta_geothermalflux,penta_artdiff,
    244                                 penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning);
     244                                penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning,0);
    245245
    246246                /*Add penta element to elements dataset: */

  • issm/trunk/src/c/ModelProcessorx/Model.cpp

    r962 r1104  
    7676        model->bed=NULL;
    7777        model->elementoniceshelf=NULL;
     78        model->elementonwater=NULL;
    7879        model->gridonicesheet=NULL;
    7980        model->gridoniceshelf=NULL;
    … …  
    237238        xfree((void**)&model->q);
    238239        xfree((void**)&model->elementoniceshelf);
     240        xfree((void**)&model->elementonwater);
    239241        xfree((void**)&model->gridonicesheet);
    240242        xfree((void**)&model->gridoniceshelf);

  • issm/trunk/src/c/ModelProcessorx/Model.h

    r962 r1104  
    7373        double* bed;
    7474        double* elementoniceshelf;
     75        double* elementonwater;
    7576        double* gridonicesheet;
    7677        double* gridoniceshelf;

  • issm/trunk/src/c/ModelProcessorx/Prognostic/CreateElementsNodesAndMaterialsPrognostic.cpp

    r1025 r1104  
    7171        double tria_viscosity_overshoot;
    7272        int    tria_artdiff;
     73        bool   tria_onwater;
    7374
    7475        /*matice constructor input: */
    … …  
    225226                ModelFetchData((void**)&model->bed,NULL,NULL,model_handle,"bed","Matrix","Mat");
    226227                ModelFetchData((void**)&model->elementoniceshelf,NULL,NULL,model_handle,"elementoniceshelf","Matrix","Mat");
     228                ModelFetchData((void**)&model->elementonwater,NULL,NULL,model_handle,"elementonwater","Matrix","Mat");
    227229               
    228230                for (i=0;i<model->numberofelements;i++){
    … …  
    259261                        /*element on iceshelf?:*/
    260262                        tria_shelf=(int)*(model->elementoniceshelf+i);
     263                        tria_onwater=(bool)*(model->elementonwater+i);
    261264                        tria_artdiff=model->artificial_diffusivity;
    262265
    263266                        /*Create tria element using its constructor:*/
    264                         tria=new Tria(tria_id, tria_mid, tria_mparid, tria_g, tria_h, tria_s, tria_b, tria_k, tria_melting,tria_accumulation,tria_geothermalflux,tria_friction_type, tria_p, tria_q, tria_shelf, tria_meanvel, tria_epsvel, tria_viscosity_overshoot,tria_artdiff);
     267                        tria=new Tria(tria_id, tria_mid, tria_mparid, tria_g, tria_h, tria_s, tria_b, tria_k, tria_melting,tria_accumulation,tria_geothermalflux,tria_friction_type, tria_p, tria_q, tria_shelf, tria_meanvel, tria_epsvel, tria_viscosity_overshoot,tria_artdiff,tria_onwater);
    265268
    266269                        /*Add tria element to elements dataset: */
    … …  
    290293                xfree((void**)&model->bed);
    291294                xfree((void**)&model->elementoniceshelf);
     295                xfree((void**)&model->elementonwater);
    292296
    293297        }
    … …  
    335339                                        penta_p,penta_q,penta_shelf,penta_onbed,penta_onsurface,penta_meanvel,penta_epsvel,
    336340                                        penta_collapse,penta_melting,penta_accumulation,penta_geothermalflux,penta_artdiff,
    337                                         penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning);
     341                                        penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning,0);
    338342
    339343                        /*Add penta element to elements dataset: */

  • issm/trunk/src/c/ModelProcessorx/Qmu/CreateParametersQmu.cpp

    r962 r1104  
    156156                #endif
    157157
    158                 /*partition grids in model->qmu_npart parts: */
    159 
    160                 if(strcmp(model->meshtype,"2d")==0){
    161                         ModelFetchData((void**)&model->elements,NULL,NULL,model_handle,"elements","Matrix","Mat");
    162                         elements_width=3; //tria elements
    163                 }
    164                 else{
    165                         ModelFetchData((void**)&model->elements2d,NULL,NULL,model_handle,"elements2d","Matrix","Mat");
    166                         elements_width=6; //penta elements
    167                 }
    168 
    169                 MeshPartitionx(&epart, &part,model->numberofelements,model->numberofnodes,model->elements, model->numberofelements2d,model->numberofnodes2d,model->elements2d,model->numlayers,elements_width, model->meshtype,model->qmu_npart);
    170 
    171                 dpart=(double*)xmalloc(model->numberofnodes*sizeof(double));
    172                 for(i=0;i<model->numberofnodes;i++)dpart[i]=part[i];
     158                /*partition grids in model->qmu_npart parts, unless a partition is already present: */
     159                ModelFetchData((void**)&dpart,NULL,NULL,model_handle,"part","Matrix","Mat");
     160
     161                if(!dpart){
     162
     163                        if(strcmp(model->meshtype,"2d")==0){
     164                                ModelFetchData((void**)&model->elements,NULL,NULL,model_handle,"elements","Matrix","Mat");
     165                                elements_width=3; //tria elements
     166                        }
     167                        else{
     168                                ModelFetchData((void**)&model->elements2d,NULL,NULL,model_handle,"elements2d","Matrix","Mat");
     169                                elements_width=6; //penta elements
     170                        }
     171
     172                        MeshPartitionx(&epart, &part,model->numberofelements,model->numberofnodes,model->elements, model->numberofelements2d,model->numberofnodes2d,model->elements2d,model->numlayers,elements_width, model->meshtype,model->qmu_npart);
     173
     174                        dpart=(double*)xmalloc(model->numberofnodes*sizeof(double));
     175                        for(i=0;i<model->numberofnodes;i++)dpart[i]=part[i];
     176                }
    173177
    174178                count++;

  • issm/trunk/src/c/ModelProcessorx/SlopeCompute/CreateElementsNodesAndMaterialsSlopeCompute.cpp

    r1036 r1104  
    6767        double tria_viscosity_overshoot;
    6868        int    tria_artdiff;
     69        bool   tria_onwater;
    6970
    7071        /*penta constructor input: */
    … …  
    201202                ModelFetchData((void**)&model->elements,NULL,NULL,model_handle,"elements","Matrix","Mat");
    202203                ModelFetchData((void**)&model->surface,NULL,NULL,model_handle,"surface","Matrix","Mat");
     204                ModelFetchData((void**)&model->elementonwater,NULL,NULL,model_handle,"elementonwater","Matrix","Mat");
    203205                ModelFetchData((void**)&model->bed,NULL,NULL,model_handle,"bed","Matrix","Mat");
    204206               
    … …  
    229231                        tria_b[1]=*(model->bed+        ((int)*(model->elements+elements_width*i+1)-1));
    230232                        tria_b[2]=*(model->bed+        ((int)*(model->elements+elements_width*i+2)-1));
     233                       
     234                        /*element on water? : */
     235                        tria_onwater=(bool)*(model->elementonwater+i);
    231236
    232237                        /*Create tria element using its constructor:*/
    233                         tria=new Tria(tria_id, tria_mid, tria_mparid, tria_g, tria_h, tria_s, tria_b, tria_k, tria_melting,tria_accumulation,tria_geothermalflux,tria_friction_type, tria_p, tria_q, tria_shelf, tria_meanvel, tria_epsvel, tria_viscosity_overshoot,tria_artdiff);
     238                        tria=new Tria(tria_id, tria_mid, tria_mparid, tria_g, tria_h, tria_s, tria_b, tria_k, tria_melting,tria_accumulation,tria_geothermalflux,tria_friction_type, tria_p, tria_q, tria_shelf, tria_meanvel, tria_epsvel, tria_viscosity_overshoot,tria_artdiff,tria_onwater);
    234239
    235240                        /*Add tria element to elements dataset: */
    … …  
    294299                                        penta_p,penta_q,penta_shelf,penta_onbed,penta_onsurface,penta_meanvel,penta_epsvel,
    295300                                        penta_collapse,penta_melting,penta_accumulation,penta_geothermalflux,penta_artdiff,
    296                                         penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning);
     301                                        penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning,0);
    297302
    298303                        /*Add penta element to elements dataset: */

  • issm/trunk/src/c/ModelProcessorx/Thermal/CreateElementsNodesAndMaterialsThermal.cpp

    r1025 r1104  
    246246                                penta_p,penta_q,penta_shelf,penta_onbed,penta_onsurface,penta_meanvel,penta_epsvel,
    247247                                penta_collapse,penta_melting,penta_accumulation,penta_geothermalflux,penta_artdiff,
    248                                 penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning);
     248                                penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning,0);
    249249
    250250                /*Add penta element to elements dataset: */

  • issm/trunk/src/c/objects/Beam.cpp

    r803 r1104  
    2626}
    2727               
    28 Beam::Beam(int beam_id, int beam_mid, int beam_mparid, int beam_g[2], double beam_h[2], double beam_s[2],double beam_b[2],double beam_k[2],int beam_onbed){
     28Beam::Beam(int beam_id, int beam_mid, int beam_mparid, int beam_g[2], double beam_h[2], double beam_s[2],double beam_b[2],double beam_k[2],bool beam_onbed){
    2929
    3030        int i;

  • issm/trunk/src/c/objects/Beam.h

    r803 r1104  
    3939                double k[2];
    4040
    41                 int onbed;
     41                bool onbed;
    4242
    4343        public:
    4444
    4545                Beam();
    46                 Beam(int beam_id, int beam_mid, int beam_mparid, int beam_g[2], double beam_h[2], double beam_s[2],double beam_b[2],double beam_k[2],int beam_onbed);
     46                Beam(int beam_id, int beam_mid, int beam_mparid, int beam_g[2], double beam_h[2], double beam_s[2],double beam_b[2],double beam_k[2],bool beam_onbed);
    4747                ~Beam();
    4848

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

    r995 r1104  
    2222                                double penta_p, double penta_q, int penta_shelf, int penta_onbed, int penta_onsurface, double penta_meanvel,double penta_epsvel,
    2323                                int penta_collapse, double penta_melting[6], double penta_accumulation[6], double penta_geothermalflux[6],
    24                                 int penta_artdiff, int penta_thermal_steadystate,double penta_viscosity_overshoot,double penta_stokesreconditioning){
     24                                int penta_artdiff, int penta_thermal_steadystate,double penta_viscosity_overshoot,double penta_stokesreconditioning,bool penta_onwater){
    2525       
    2626        int i;
    … …  
    5252        shelf = penta_shelf;
    5353        onbed = penta_onbed;
     54        onwater = penta_onwater;
    5455        onsurface = penta_onsurface;
    5556        meanvel = penta_meanvel;
    … …  
    8990        printf("   shelf: %i\n",shelf);
    9091        printf("   onbed: %i\n",onbed);
     92        printf("   onwater: %i\n",onwater);
    9193        printf("   onsurface: %i\n",onsurface);
    9294        printf("   meanvel: %g\n",meanvel);
    … …  
    125127        printf("   shelf: %i\n",shelf);
    126128        printf("   onbed: %i\n",onbed);
     129        printf("   onwater: %i\n",onwater);
    127130        printf("   onsurface: %i\n",onsurface);
    128131        printf("   meanvel: %g\n",meanvel);
    … …  
    173176        memcpy(marshalled_dataset,&shelf,sizeof(shelf));marshalled_dataset+=sizeof(shelf);
    174177        memcpy(marshalled_dataset,&onbed,sizeof(onbed));marshalled_dataset+=sizeof(onbed);
     178        memcpy(marshalled_dataset,&onwater,sizeof(onwater));marshalled_dataset+=sizeof(onwater);
    175179        memcpy(marshalled_dataset,&onsurface,sizeof(onsurface));marshalled_dataset+=sizeof(onsurface);
    176180        memcpy(marshalled_dataset,&meanvel,sizeof(meanvel));marshalled_dataset+=sizeof(meanvel);
    … …  
    210214                sizeof(shelf)+
    211215                sizeof(onbed)+
     216                sizeof(onwater)+
    212217                sizeof(onsurface)+
    213218                sizeof(meanvel)+
    … …  
    258263        memcpy(&shelf,marshalled_dataset,sizeof(shelf));marshalled_dataset+=sizeof(shelf);
    259264        memcpy(&onbed,marshalled_dataset,sizeof(onbed));marshalled_dataset+=sizeof(onbed);
     265        memcpy(&onwater,marshalled_dataset,sizeof(onwater));marshalled_dataset+=sizeof(onwater);
    260266        memcpy(&onsurface,marshalled_dataset,sizeof(onsurface));marshalled_dataset+=sizeof(onsurface);
    261267        memcpy(&meanvel,marshalled_dataset,sizeof(meanvel));marshalled_dataset+=sizeof(meanvel);
    … …  
    452458        /*Collapsed formulation: */
    453459        Tria*  tria=NULL;
     460
     461        /*If on water, skip stiffness: */
     462        if(onwater)return;
     463
    454464
    455465        /*recover pointers: */
    … …  
    13251335        tria_node_offsets[2]=node_offsets[g2];
    13261336
    1327         tria= new Tria(id,mid,mparid,tria_node_ids,tria_h,tria_s,tria_b,tria_k, tria_melting, tria_accumulation, tria_geothermalflux,friction_type,p,q,shelf,meanvel,epsvel,viscosity_overshoot,artdiff);
     1337        tria= new Tria(id,mid,mparid,tria_node_ids,tria_h,tria_s,tria_b,tria_k, tria_melting, tria_accumulation, tria_geothermalflux,friction_type,p,q,shelf,meanvel,epsvel,viscosity_overshoot,artdiff,onwater);
    13281338
    13291339        tria->NodeConfiguration(tria_node_ids,tria_nodes,tria_node_offsets);
    … …  
    17681778        /*Spawning: */
    17691779        Tria* tria=NULL;
     1780
     1781        /*If on water, skip load: */
     1782        if(onwater)return;
    17701783
    17711784        /*recover pointers: */

  • issm/trunk/src/c/objects/Penta.h

    r848 r1104  
    4646                int    onbed;
    4747                int    onsurface;
     48                bool    onwater;
    4849                double meanvel;/*!scaling ratio for velocities*/
    4950                double epsvel; /*!minimum velocity to avoid infinite velocity ratios*/
    … …  
    6162                                double p, double q, int    shelf, int    onbed, int    onsurface, double meanvel,double epsvel,
    6263                                int    collapse, double melting[6], double accumulation[6], double geothermalflux[6],
    63                                 int    artdiff, int    thermal_steadystate,double viscosity_overshoot,double stokesreconditioning);
     64                                int    artdiff, int    thermal_steadystate,double viscosity_overshoot,double stokesreconditioning,bool onwater);
    6465                ~Penta();
    6566

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

    r1051 r1104  
    3434Tria::Tria(int tria_id,int tria_mid,int tria_mparid,int tria_node_ids[3],double tria_h[3],double tria_s[3],double tria_b[3],double tria_k[3],double tria_melting[3],
    3535                                double tria_accumulation[3],double tria_geothermalflux[3],int tria_friction_type,double tria_p,double tria_q,int tria_shelf,double tria_meanvel,double tria_epsvel,
    36                                 double tria_viscosity_overshoot,int tria_artdiff){
     36                                double tria_viscosity_overshoot,int tria_artdiff,bool tria_onwater){
    3737       
    3838        int i;
    … …  
    6464        epsvel=tria_epsvel;
    6565        onbed=1;
     66        onwater=tria_onwater;
    6667        viscosity_overshoot=tria_viscosity_overshoot;
    6768        artdiff=tria_artdiff;
    … …  
    9798        printf("   epsvel: %g\n",epsvel);
    9899        printf("   onbed: %i\n",onbed);
     100        printf("   onwater: %i\n",onwater);
    99101        printf("   viscosity_overshoot=%g\n",viscosity_overshoot);
    100102        printf("   artdiff=%g\n",artdiff);
    … …  
    135137        printf("   epsvel: %g\n",epsvel);
    136138        printf("   onbed: %i\n",onbed);
     139        printf("   onwater: %i\n",onwater);
    137140        printf("   viscosity_overshoot=%g\n",viscosity_overshoot);
    138141        printf("   artdiff=%g\n",artdiff);
    … …  
    180183        memcpy(marshalled_dataset,&friction_type,sizeof(friction_type));marshalled_dataset+=sizeof(friction_type);
    181184        memcpy(marshalled_dataset,&onbed,sizeof(onbed));marshalled_dataset+=sizeof(onbed);
     185        memcpy(marshalled_dataset,&onwater,sizeof(onwater));marshalled_dataset+=sizeof(onwater);
    182186        memcpy(marshalled_dataset,&p,sizeof(p));marshalled_dataset+=sizeof(p);
    183187        memcpy(marshalled_dataset,&q,sizeof(q));marshalled_dataset+=sizeof(q);
    … …  
    212216                +sizeof(friction_type)
    213217                +sizeof(onbed)
     218                +sizeof(onwater)
    214219                +sizeof(p)
    215220                +sizeof(q)
    … …  
    256261        memcpy(&friction_type,marshalled_dataset,sizeof(friction_type));marshalled_dataset+=sizeof(friction_type);
    257262        memcpy(&onbed,marshalled_dataset,sizeof(onbed));marshalled_dataset+=sizeof(onbed);
     263        memcpy(&onwater,marshalled_dataset,sizeof(onwater));marshalled_dataset+=sizeof(onwater);
    258264        memcpy(&p,marshalled_dataset,sizeof(p));marshalled_dataset+=sizeof(p);
    259265        memcpy(&q,marshalled_dataset,sizeof(q));marshalled_dataset+=sizeof(q);
    … …  
    399405
    400406        ParameterInputs* inputs=NULL;
     407
     408        /*First, if we are on water, return empty matrix: */
     409        if(onwater)return;
    401410
    402411        /*recover pointers: */
    … …  
    12041213        ParameterInputs* inputs=NULL;
    12051214
     1215        /*First, if we are on water, return empty vector: */
     1216        if(onwater)return;
     1217
    12061218        /*recover pointers: */
    12071219        inputs=(ParameterInputs*)vinputs;

  • issm/trunk/src/c/objects/Tria.h

    r803 r1104  
    4848                double epsvel; /*!minimum velocity to avoid infinite velocity ratios*/
    4949                int    onbed;
     50                bool   onwater;
    5051                double viscosity_overshoot;
    5152                int    artdiff;
    … …  
    5556                Tria();
    5657                Tria(int id,int mid,int mparid,int node_ids[3],double h[3],double s[3],double b[3],double k[3],double melting[3],double accumulation[3],double geothermalflux[3],
    57                                 int friction_type,double p,double q,int shelf,double meanvel,double epsvel,double viscosity_overshoot,int artdiff);
     58                                int friction_type,double p,double q,int shelf,double meanvel,double epsvel,double viscosity_overshoot,int artdiff,bool onwater);
    5859                ~Tria();
    5960
Note: See TracChangeset for help on using the changeset viewer.