Changeset 19753


Ignore:
Timestamp:
11/18/15 14:53:42 (9 years ago)
Author:
Eric.Larour
Message:

CHG: new modules in javascript, IssmConfig and InterpFromMeshToMesh2d, which uses the
bamg objects.
New classes to be able to finish off the model.js implementation. New paterson, solvers for javacsript.

Location:
issm/trunk-jpl
Files:
22 added
13 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/array/arrayoperations.js

    r19736 r19753  
    7474        return notarray;
    7575} //}}}
     76function ArrayPow(array,coefficient) { //{{{
     77
     78        var powarray=array;
     79        for (i=0;i<array.length;i++)powarray[i]=Math.pow(array[i],coefficient);
     80        return powarray;
     81} //}}}
     82function ArraySqrt(array) { //{{{
     83
     84        var sqrtarray=array;
     85        for (i=0;i<array.length;i++)sqrtarray[i]=Math.sqrt(array[i]);
     86        return sqrtarray;
     87} //}}}
     88function ArrayMag(array1,array2) { //{{{
     89
     90        var arraymag=NewArrayFill(array1.length,0);
     91        for (var i=0;i<array1.length;i++)arraymag[i]=Math.sqrt(Math.pow(array1[i],2)+Math.pow(array2[i],2));
     92        return arraymag;
     93} //}}}
    7694function ArrayAnyNaN(array) { //{{{
    7795
    7896        for(var i=0;i<array.length;i++){
    7997                if (isNaN(array[i])) return 1;
     98        }
     99        return 0;
     100} //}}}
     101function ArrayAnyBelowOrEqual(array,value) { //{{{
     102
     103        for(var i=0;i<array.length;i++){
     104                if (array[i]<=value)return 1;
     105        }
     106        return 0;
     107} //}}}
     108function ArrayAnyBelowStrict(array,value) { //{{{
     109
     110        for(var i=0;i<array.length;i++){
     111                if (array[i]<value)return 1;
    80112        }
    81113        return 0;

  • issm/trunk-jpl/src/m/classes/model.js

    r19719 r19753  
    5252                        this.constants        = new constants();
    5353                        this.geometry         = new geometry();
    54 //                      this.initialization   = initialization();
    55 //                      this.smb              = SMBforcing();
    56 //                      this.basalforcings    = basalforcings();
    57 //                      this.friction         = friction();
    58 //                      this.rifts            = rifts();
    59 //                      this.timestepping     = timestepping();
    60 //                      this.groundingline    = groundingline();
    61 //                      this.materials        = matice();
    62 //                      this.damage           = damage();
    63 //                      this.flowequation     = flowequation();
    64 //                      this.debug            = debug();
    65 //                      this.verbose          = verbose();
    66 //                      this.settings         = settings();
    67 //                      this.toolkits         = toolkits();
    68 //                      this.cluster          = generic();
    69 //                      this.balancethickness = balancethickness();
    70 //                      this.stressbalance    = stressbalance();
    71 //                      this.hydrology        = hydrologyshreve();
    72 //                      this.masstransport    = masstransport();
    73 //                      this.thermal          = thermal();
    74 //                      this.steadystate      = steadystate();
    75 //                      this.transient        = transient();
    76 //                      this.calving          = calving();
    77 //                      this.gia              = gia();
    78 //                      this.autodiff         = autodiff();
    79 //                      this.flaim            = flaim();
    80 //                      this.inversion        = inversion();
    81 //                      this.qmu              = qmu();
    82 //                      this.radaroverlay     = radaroverlay();
    83 //                      this.results          = struct();
    84 //                      this.outputdefinition = outputdefinition();
    85 //                      this.miscellaneous    = miscellaneous();
    86 //                      this.private          = private();
     54                        this.initialization   = new initialization();
     55                        this.smb              = new SMBforcing();
     56                        this.basalforcings    = new basalforcings();
     57                        this.friction         = new friction();
     58                        this.rifts            = new rifts();
     59                        this.timestepping     = new timestepping();
     60                        this.groundingline    = new groundingline();
     61                        this.materials        = new matice();
     62                        this.damage           = new damage();
     63                        this.flowequation     = new flowequation();
     64                        this.debug            = new debug();
     65                        this.verbose          = new verbose();
     66                        this.settings         = new settings();
     67                        this.toolkits         = new toolkits();
     68                        this.toolkits.disp();
     69                        this.cluster          = new generic();
     70                        this.balancethickness = new balancethickness();
     71                        this.stressbalance    = new stressbalance();
     72                        this.hydrology        = new hydrologyshreve();
     73                        this.masstransport    = new masstransport();
     74                        this.thermal          = new thermal();
     75                        this.steadystate      = new steadystate();
     76                        this.transient        = new transient();
     77                        this.calving          = new calving();
     78                        this.gia              = new gia();
     79                        this.autodiff         = new autodiff();
     80                        this.flaim            = new flaim();
     81                        this.inversion        = new inversion();
     82                        this.qmu              = new qmu();
     83                        this.radaroverlay     = new radaroverlay();
     84                        this.results          = [];
     85                        this.outputdefinition = new outputdefinition();
     86                        this.miscellaneous    = new miscellaneous();
     87                        this.private          = new private();
    8788                } //}}}
    8889        //properties

  • issm/trunk-jpl/src/m/classes/pairoptions.js

    r19721 r19753  
    265265
    266266        //properties
    267         list         = [];
     267        this.list         = [];
    268268        this.constructor(args);
    269269}

  • issm/trunk-jpl/src/m/classes/plotoptions.js

    r19721 r19753  
    159159    //properties
    160160         // {{{
    161          numberofplots = 0;
    162          figurenumber  = 1;
    163          list          = [];
     161         this.numberofplots = 0;
     162         this.figurenumber  = 1;
     163         this.list          = [];
    164164         this.constructor(args);
    165165         //}}}

  • issm/trunk-jpl/src/m/miscellaneous/fielddisplay.js

    r19721 r19753  
    2121                }
    2222                else{
    23                         console.log(offset)
    24                         console.log(name)
    25                         console.log(field)
    26                         //displayunit(offset,name,"""" + field """",comment),
     23                        displayunit(offset,name,"'"+field+"'",comment);
    2724                }
    2825        }
    … …  
    8279        else{
    8380                if (typeof(comment) == "string"){
    84                         console.log(sprintf("%s%-23s: %-15s -- %s",offset,name,characterization,comment));
     81                        //console.log(sprintf("%s%-23s: %-15s -- %s",offset,name,characterization,comment));
     82                        console.log(sprintf("%s%s: %-15s -- %s",offset,name,characterization,comment));
    8583                }
    8684                else{

  • issm/trunk-jpl/src/wrappers/InterpFromMeshToMesh2d/InterpFromMeshToMesh2d.h

    r14996 r19753  
    5555#endif
    5656
     57#ifdef _HAVE_JAVASCRIPT_MODULES_
     58/* serial input macros: */
     59#define INDEX          indexin,nelin,3
     60#define X              xin,nodsin
     61#define Y              yin,nodsin
     62#define DATA           datain,nodsin,1
     63#define XINTERP        x_interpin, nods_interpin
     64#define YINTERP        y_interpin, nods_interpin
     65#define ARGUMENTS "default_value",default_value
     66/* serial output macros: */
     67#define DATAINTERP pdata_interp,NULL,NULL
     68#define WRAPPER(modulename) extern "C" { int  InterpFromMeshToMesh2dModule(double** pdata_interp,int* indexin,double* xin,double* yin,double* datain,double* x_interpin,double* y_interpin,int nelin,int nodsin,int nods_interpin,double default_value)
     69#define nrhs  6
     70#endif
     71
    5772/* serial arg counts: */
    5873#undef NLHS

  • issm/trunk-jpl/src/wrappers/IssmConfig/IssmConfig.h

    r15905 r19753  
    3131/* serial output macros: */
    3232#define VALUE (mxArray**)&plhs[0]
     33#define SVALUE (mxArray**)&plhs[0]
    3334#endif
    3435
    … …  
    3839/* serial output macros: */
    3940#define VALUE output,0
     41#define SVALUE output,0
     42#endif
     43
     44
     45#ifdef _HAVE_JAVASCRIPT_MODULES_
     46/* serial input macros: */
     47#define NAME string
     48/* serial output macros: */
     49#define VALUE pvalue
     50#define SVALUE psvalue
     51#define WRAPPER(modulename) extern "C" { int  IssmConfigModule(double* pvalue, char** psvalue, char* string)
    4052#endif
    4153

  • issm/trunk-jpl/src/wrappers/javascript/Makefile.am

    r19739 r19753  
    1010bin_SCRIPTS += ../TriMesh/TriMesh.js  \
    1111                           ../NodeConnectivity/NodeConnectivity.js\
    12                            ../ElementConnectivity/ElementConnectivity.js
     12                           ../ElementConnectivity/ElementConnectivity.js\
     13                           ../InterpFromMeshToMesh2d/InterpFromMeshToMesh2d.js\
     14                           ../IssmConfig/IssmConfig.js
    1315       
    1416#javascript io{{{
    … …  
    7577IssmModule_SOURCES = ../TriMesh/TriMesh.cpp \
    7678                                         ../NodeConnectivity/NodeConnectivity.cpp\
    77                                          ../ElementConnectivity/ElementConnectivity.cpp
     79                                         ../ElementConnectivity/ElementConnectivity.cpp\
     80                                         ../InterpFromMeshToMesh2d/InterpFromMeshToMesh2d.cpp\
     81                                         ../IssmConfig/IssmConfig.cpp
    7882
    79 IssmModule_CXXFLAGS= -fPIC -D_DO_NOT_LOAD_GLOBALS_ --memory-init-file 0 $(AM_CXXFLAGS) $(CXXFLAGS) $(CXXOPTFLAGS) $(COPTFLAGS) -s EXPORTED_FUNCTIONS="['_TriMeshModule','_NodeConnectivityModule','_ElementConnectivityModule']"  -s DISABLE_EXCEPTION_CATCHING=0 -s ALLOW_MEMORY_GROWTH=1
     83IssmModule_CXXFLAGS= -fPIC -D_DO_NOT_LOAD_GLOBALS_ --memory-init-file 0 $(AM_CXXFLAGS) $(CXXFLAGS) $(CXXOPTFLAGS) $(COPTFLAGS) -s EXPORTED_FUNCTIONS="['_TriMeshModule','_NodeConnectivityModule','_ElementConnectivityModule','_InterpFromMeshToMesh2dModule','_IssmConfigModule']"  -s DISABLE_EXCEPTION_CATCHING=0 -s ALLOW_MEMORY_GROWTH=1
    8084IssmModule_LDADD = ${deps} $(TRIANGLELIB)
    8185#}}}

  • issm/trunk-jpl/src/wrappers/javascript/io/FetchJavascriptData.cpp

    r19711 r19753  
    1313
    1414/*Primitive data types*/
     15/*FUNCTION FetchData(char** pstring, char* string){{{*/
     16void FetchData(char** pstring, char* stringin){
     17
     18        char* string=xNew<char>(strlen(stringin)+1);
     19        xMemCpy<char>(string,stringin,strlen(stringin)+1);
     20
     21        *pstring=string;
     22}
     23/*}}}*/
    1524/*FUNCTION FetchData(double* pscalar,double scalar){{{*/
    1625void FetchData(double* pscalar,double scalar){
    … …  
    3544}
    3645/*}}}*/
     46/*FUNCTION FetchData(double **pvector, int* pnods, double* vectorin, int nods){{{*/
     47void FetchData(double** pvector, int* pnods, double* vectorin, int nods){
     48
     49        double* vector=NULL;
     50       
     51        vector=xNew<IssmPDouble>(nods); xMemCpy<IssmPDouble>(vector,vectorin,nods);
     52       
     53        *pvector=vector;
     54        *pnods=nods;
     55}
     56/*}}}*/
    3757/*FUNCTION FetchData(double **pmatrix, int* pM, int* pN, int* matrix, int M, int N){{{*/
    3858void FetchData(double **pmatrix, int* pM, int* pN, int* matrixin, int M, int N){
    … …  
    4464                matrix=xNew<IssmPDouble>(M*N);
    4565                for(int i=0;i<M*N;i++)matrix[i]=(IssmPDouble)matrixin[i];
     66                if (pM)*pM=M;
     67                if (pN)*pN=N;
     68                *pmatrix=matrix;
     69        }
     70}
     71/*}}}*/
     72/*FUNCTION FetchData(double **pmatrix, int* pM, int* pN, double* matrix, int M, int N){{{*/
     73void FetchData(double **pmatrix, int* pM, int* pN, double* matrixin, int M, int N){
     74
     75        double* matrix=NULL;
     76       
     77        if(pmatrix && matrixin){
     78
     79                matrix=xNew<IssmPDouble>(M*N);
     80                for(int i=0;i<M*N;i++)matrix[i]=matrixin[i];
    4681                if (pM)*pM=M;
    4782                if (pN)*pN=N;
    … …  
    89124}
    90125/*}}}*/
     126/*FUNCTION FetchData(Options** poptions,int NRHS, int nrhs, const char* optionname, double optionvalue){{{*/
     127void FetchData(Options** poptions,int NRHS, int nrhs, const char* optionname, double optionvalue){
     128
     129        /*Initialize output*/
     130        Options* options=new Options();
     131       
     132        GenericOption<double> *odouble = NULL;
     133
     134        /*check and parse the name  */
     135        odouble=new GenericOption<double>();
     136        odouble->name =xNew<char>(strlen(optionname)+1);
     137        memcpy(odouble->name,optionname,(strlen(optionname)+1)*sizeof(char));
     138        odouble->value=optionvalue;
     139        odouble->numel=1;
     140        odouble->ndims=1;
     141        odouble->size=NULL;
     142         
     143        options->AddOption((Option*)odouble);
     144
     145        /*Assign output pointers:*/
     146        *poptions=options;
     147}
     148/*}}}*/

  • issm/trunk-jpl/src/wrappers/javascript/io/WriteJavascriptData.cpp

    r19711 r19753  
    4040}
    4141/*}}}*/
     42/*FUNCTION WriteData(IssmPDouble** pmatrix, int* pM, IssmPDouble* pN, , int* matrix, int M,int N){{{*/
     43void WriteData(IssmPDouble** pmatrix,int* pM, int* pN, IssmPDouble* matrix, int M, int N){
     44
     45        if(pmatrix && matrix){
     46
     47                /*Copy matrix: */
     48                IssmPDouble* dmatrix = xNew<IssmPDouble>(M*N);
     49                for (int i=0;i<M*N;i++)dmatrix[i]=matrix[i];
     50                *pmatrix=dmatrix;
     51                if(pM)*pM=M;
     52                if(pN)*pN=N;
     53        }
     54}
     55/*}}}*/
    4256/*FUNCTION WriteData(IssmPDouble** px, int* pnods, double* vector, int M){{{*/
    4357void WriteData(IssmPDouble** px, int* pnods, double* vector, int M){
    … …  
    6478}
    6579/*}}}*/
     80/*FUNCTION WriteData(IssmPDouble* pdouble, IssmPDouble double){{{*/
     81void WriteData(IssmPDouble* pdouble, IssmPDouble doublein){
    6682
    67 /*ISSM objects*/
     83        *pdouble=doublein;
     84}
     85/*}}}*/
     86/*FUNCTION WriteData(char** pstring, char* string){{{*/
     87void WriteData(char** pstring, char* stringin){
    6888
    69 /*Toolkit*/
     89        char* string=xNew<char>(strlen(stringin)+1);
     90        xMemCpy<char>(string,stringin,strlen(stringin)+1);
     91
     92        *pstring=string;
     93}
     94/*}}}*/

  • issm/trunk-jpl/src/wrappers/javascript/io/javascriptio.h

    r19711 r19753  
    2020void WriteData(IssmPDouble** pmatrix,int* pnel, int* matrix, int M,int N);
    2121void WriteData(IssmPDouble** pmatrix,int* pM, int* pN, int* matrix, int M, int N);
     22void WriteData(IssmPDouble** pmatrix,int* pM, int* pN, IssmPDouble* matrix, int M, int N);
    2223void WriteData(IssmPDouble** px, int* pnods, int* vector, int M);
    2324void WriteData(IssmPDouble** px, int* pnods, double* vector, int M);
     25void WriteData(char** pstring, char* stringin);
     26void WriteData(IssmPDouble* pdouble, IssmPDouble doublein);
    2427
     28void FetchData(char** pstring, char* stringin);
    2529void FetchData(double* pscalar,double scalar);
    2630void FetchData(int* pinteger,int integer);
    2731void FetchData(double** pvector, double* vectorin, int nods);
     32void FetchData(double** pvector, int* pnods, double* vectorin, int nods);
    2833void FetchData(double **pmatrix, int* pM, int* pN, int* matrixin, int M, int N);
     34void FetchData(double **pmatrix, int* pM, int* pN, double* matrixin, int M, int N);
    2935void FetchData(int **pmatrix, int* pM, int* pN, int* matrixin, int M, int N);
    3036void FetchData(Contours** pcontours,double* x, double* y, int nods);
     37void FetchData(Options** poptions,int NRHS, int nrhs, const char* optionname, double optionvalue);
    3138
    3239/*Print*/

  • issm/trunk-jpl/test/NightlyRun/test101.html

    r19743 r19753  
    55<!-- Includes {{{-->
    66<script type="text/javascript" src="../Exp/Square.js"></script>
     7<script type="text/javascript" src="../Par/SquareShelfConstrained.js"></script>
     8<script type="text/javascript" src="../Data/SquareShelfConstrained.data.js"></script>
    79<script type="text/javascript" src="../../src/m/miscellaneous/fielddisplay.js"></script>
    810<script type="text/javascript" src="../../src/m/classes/model.js"></script>
    … …  
    1113<script type="text/javascript" src="../../src/m/classes/mesh2d.js"></script>
    1214<script type="text/javascript" src="../../src/m/classes/geometry.js"></script>
     15<script type="text/javascript" src="../../src/m/classes/matice.js"></script>
    1316<script type="text/javascript" src="../../src/m/classes/constants.js"></script>
    1417<script type="text/javascript" src="../../src/m/classes/mask.js"></script>
     18<script type="text/javascript" src="../../src/m/classes/initialization.js"></script>
     19<script type="text/javascript" src="../../src/m/classes/SMBforcing.js"></script>
     20<script type="text/javascript" src="../../src/m/classes/timestepping.js"></script>
     21<script type="text/javascript" src="../../src/m/classes/basalforcings.js"></script>
     22<script type="text/javascript" src="../../src/m/classes/groundingline.js"></script>
     23<script type="text/javascript" src="../../src/m/classes/damage.js"></script>
     24<script type="text/javascript" src="../../src/m/classes/flowequation.js"></script>
     25<script type="text/javascript" src="../../src/m/classes/debug.js"></script>
     26<script type="text/javascript" src="../../src/m/classes/verbose.js"></script>
     27<script type="text/javascript" src="../../src/m/classes/settings.js"></script>
     28<script type="text/javascript" src="../../src/m/classes/toolkits.js"></script>
     29<script type="text/javascript" src="../../src/m/classes/stressbalance.js"></script>
     30<script type="text/javascript" src="../../src/m/classes/rifts.js"></script>
     31<script type="text/javascript" src="../../src/m/classes/friction.js"></script>
    1532<script type="text/javascript" src="../../src/m/classes/plotoptions.js"></script>
    1633<script type="text/javascript" src="../../src/m/classes/pairoptions.js"></script>
     34<script type="text/javascript" src="../../src/m/classes/clusters/generic.js"></script>
    1735<script type="text/javascript" src="../../src/m/parameterization/setmask.js"></script>
    1836<script type="text/javascript" src="../../src/m/miscellaneous/rgbcolor.js"></script>
     37<script type="text/javascript" src="../../src/m/materials/paterson.js"></script>
    1938<script type="text/javascript" src="../../src/m/miscellaneous/colorbars.js"></script>
    2039<script type="text/javascript" src="../../src/m/plot/plotmodel.js"></script>
    … …  
    2645<script type="text/javascript" src="../../src/m/plot/plot_manager.js"></script>
    2746<script type="text/javascript" src="../../src/m/geometry/FlagElements.js"></script>
     47<script type="text/javascript" src="../../src/m/solvers/issmgslsolver.js"></script>
    2848<script type="text/javascript" src="../../src/wrappers/TriMesh/TriMesh.js"></script>
     49<script type="text/javascript" src="../../src/wrappers/IssmConfig/IssmConfig.js"></script>
    2950<script type="text/javascript" src="../../src/wrappers/NodeConnectivity/NodeConnectivity.js"></script>
    3051<script type="text/javascript" src="../../src/wrappers/ElementConnectivity/ElementConnectivity.js"></script>
     52<script type="text/javascript" src="../../src/wrappers/InterpFromMeshToMesh2d/InterpFromMeshToMesh2d.js"></script>
    3153<script type="text/javascript" src="../../build-js/src/wrappers/javascript/IssmModule.js"></script>
    3254<script type="text/javascript" src="../../externalpackages/javascript/src/sprintf.js"></script>
    … …  
    3759<body>
    3860        <canvas id="1" width="640" height="480"> </canvas>
     61        <canvas id="2" width="640" height="480"> </canvas>
     62        <canvas id="3" width="640" height="480"> </canvas>
     63        <canvas id="4" width="640" height="480"> </canvas>
    3964
    4065        <script type="text/javascript" async><!--}}}-->
    4166
    4267        var md = new model();
    43         triangle(md,square,100000);
    44         plotmodel(md,'data',md.mesh.x,'canvasid',1);
     68        triangle(md,square,10000);
     69        parameterize(md);
     70       
     71        plotmodel(md,'data',md.friction.coefficient,'canvasid',1,
     72        'data',md.materials.rheology_B,'canvasid',2,
     73        'data',md.initialization.temperature,'canvasid',3,
     74        'data',md.basalforcings.floatingice_melting_rate,'canvasid',4);
     75
    4576
    4677</script> <!--{{{-->

  • issm/trunk-jpl/test/Par/SquareShelfConstrained.js

    r19718 r19753  
    1818                md.geometry.base[i]=-md.materials.rho_ice/md.materials.rho_water*md.geometry.thickness[i];
    1919                md.geometry.surface[i]=md.geometry.base[i]+md.geometry.thickness[i];
     20                md.geometry.bed[i]=md.geometry.base[i]-10;
    2021        }
    2122       
    22         //Initial velocity 
    23         x     = transpose(ncread('../Data/SquareShelfConstrained.nc','x'));
     23        //Initial velocity: no ncreader for now, so we just load arrays.
     24        /*x     = transpose(ncread('../Data/SquareShelfConstrained.nc','x'));
    2425        y     = transpose(ncread('../Data/SquareShelfConstrained.nc','y'));
    2526        vx    = transpose(ncread('../Data/SquareShelfConstrained.nc','vx'));
    2627        vy    = transpose(ncread('../Data/SquareShelfConstrained.nc','vy'));
    27         index = transpose(ncread('../Data/SquareShelfConstrained.nc','index'));
     28        index = transpose(ncread('../Data/SquareShelfConstrained.nc','index'));*/
     29       
    2830        md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y);
    2931        md.initialization.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y);
    30         clear vx vy x y index;
    31         md.initialization.vz=zeros(md.mesh.numberofvertices,1);
    32         md.initialization.pressure=zeros(md.mesh.numberofvertices,1);
     32        md.initialization.vel=ArrayMag(md.initialization.vx,md.initialization.vy);
     33        md.initialization.vz=NewArrayFill(md.mesh.numberofvertices,0);
     34        md.initialization.pressure=NewArrayFill(md.mesh.numberofvertices,0);
    3335
    3436        //Materials
    35         md.initialization.temperature=(273-20)*ones(md.mesh.numberofvertices,1);
     37        md.initialization.temperature=NewArrayFill(md.mesh.numberofvertices,273-20);
    3638        md.materials.rheology_B=paterson(md.initialization.temperature);
    37         md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
     39        md.materials.rheology_n=NewArrayFill(md.mesh.numberofelements,3);
    3840
    3941        //Surface mass balance and basal melting
    40         md.smb.mass_balance=10*ones(md.mesh.numberofvertices,1);
    41         md.basalforcings.floatingice_melting_rate=5*ones(md.mesh.numberofvertices,1);
    42         md.basalforcings.groundedice_melting_rate=5*ones(md.mesh.numberofvertices,1);
     42        md.smb.mass_balance=NewArrayFill(md.mesh.numberofvertices,10);
     43        md.basalforcings.floatingice_melting_rate=NewArrayFill(md.mesh.numberofvertices,5.0);
     44        md.basalforcings.groundedice_melting_rate=NewArrayFill(md.mesh.numberofvertices,5.0);
    4345
    4446        //Friction
    45         md.friction.coefficient=20*ones(md.mesh.numberofvertices,1);
    46         md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    47         md.friction.p=ones(md.mesh.numberofelements,1);
    48         md.friction.q=ones(md.mesh.numberofelements,1);
     47        md.friction.coefficient=NewArrayFill(md.mesh.numberofvertices,20);
     48        for(var i=0;i<md.mesh.numberofvertices;i++)if(md.mask.groundedice_levelset[i]<0)md.friction.coefficient[i]=0;
     49        md.friction.p=NewArrayFill(md.mesh.numberofelements,1);
     50        md.friction.q=NewArrayFill(md.mesh.numberofelements,1);
    4951
    5052        //Numerical parameters
    5153        md.stressbalance.viscosity_overshoot=0.0;
     54        return;
    5255        md.masstransport.stabilization=1;
    5356        md.thermal.stabilization=1;
    … …  
    6568
    6669        //Change name so that no tests have the same name
    67         A=dbstack;
    68         if (length(A)>2), md.miscellaneous.name=A(3).file(1:end-2); end
     70        //A=dbstack;
     71        //if (length(A)>2), md.miscellaneous.name=A(3).file(1:end-2); end
    6972
    7073}
Note: See TracChangeset for help on using the changeset viewer.