Changeset 15988


Ignore:
Timestamp:
08/28/13 08:05:45 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: updated parameter files with new mask that only has levelsets

Location:
issm/trunk-jpl/test/Par
Files:
34 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/test/Par/79North.par

    r15771 r15988  
    2222
    2323%Friction
    24 pos=find(md.mask.elementonfloatingice);
    2524md.friction.coefficient=50.*ones(md.mesh.numberofvertices,1);
    26 md.friction.coefficient(md.mesh.elements(pos,:))=0.;
     25md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    2726md.friction.p=ones(md.mesh.numberofelements,1);
    2827md.friction.q=ones(md.mesh.numberofelements,1);
    2928
    3029%Ice shelf melting and surface mass balance
    31 pos=find(md.mask.elementonfloatingice);
    3230md.basalforcings.melting_rate=zeros(md.mesh.numberofvertices,1);
    33 md.basalforcings.melting_rate(md.mesh.elements(pos,:))=10;
     31md.basalforcings.melting_rate(find(md.mask.groundedice_levelset<0.))=0.;
    3432md.surfaceforcings.mass_balance=15.*ones(md.mesh.numberofvertices,1);
    3533

  • issm/trunk-jpl/test/Par/79North.py

    r15771 r15988  
    3434
    3535#Friction
    36 pos=numpy.nonzero(md.mask.elementonfloatingice)
    3736md.friction.coefficient=50.*ones((md.mesh.numberofvertices,1))
    38 md.friction.coefficient[md.mesh.elements[pos].astype(int)-1]=0.
     37md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    3938md.friction.p=ones((md.mesh.numberofelements,1))
    4039md.friction.q=ones((md.mesh.numberofelements,1))
    4140
    4241#Ice shelf melting and surface mass balance
    43 pos=numpy.nonzero(md.mask.elementonfloatingice)
    4442md.basalforcings.melting_rate=zeros((md.mesh.numberofvertices,1))
    45 md.basalforcings.melting_rate[md.mesh.elements[pos].astype(int)-1]=10.
     43md.basalforcings.melting_rate[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    4644md.surfaceforcings.mass_balance=15*numpy.ones((md.mesh.numberofvertices,1))
    4745

  • issm/trunk-jpl/test/Par/GiaBenchmarksAB.par

    r15771 r15988  
    4545
    4646%Friction
    47 pos=find(md.mask.elementonfloatingice);
    4847md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
    49 md.friction.coefficient(md.mesh.elements(pos,:))=0;
     48md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    5049md.friction.p=ones(md.mesh.numberofelements,1);
    5150md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/GiaBenchmarksCD.par

    r15771 r15988  
    4444
    4545%Friction
    46 pos=find(md.mask.elementonfloatingice);
    4746md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
    48 md.friction.coefficient(md.mesh.elements(pos,:))=0;
     47md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    4948md.friction.p=ones(md.mesh.numberofelements,1);
    5049md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/ISMIPA.par

    r15771 r15988  
    77
    88disp('      creating drag');
    9 md.friction.coefficient=200.*ones(md.mesh.numberofvertices,1); %q=1.
    10 %Take care of iceshelves: no basal drag
    11 pos=find(md.mask.elementonfloatingice);
    12 md.friction.coefficient(md.mesh.elements(pos,:))=0.;
     9md.friction.coefficient=200.*ones(md.mesh.numberofvertices,1);
     10md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    1311md.friction.p=ones(md.mesh.numberofelements,1);
    1412md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/ISMIPA.py

    r15771 r15988  
    1010
    1111print "      creating drag"
    12 md.friction.coefficient=200.*numpy.ones((md.mesh.numberofvertices,1))    #q=1.
    13 #Take care of iceshelves: no basal drag
    14 pos=numpy.nonzero(md.mask.elementonfloatingice)[0]
    15 md.friction.coefficient[md.mesh.elements[pos,:]-1]=0.
     12md.friction.coefficient=200.*numpy.ones((md.mesh.numberofvertices,1))
     13md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    1614md.friction.p=numpy.ones((md.mesh.numberofelements,1))
    1715md.friction.q=numpy.ones((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/ISMIPB.par

    r15771 r15988  
    77
    88disp('      creating drag');
    9 md.friction.coefficient=200.*ones(md.mesh.numberofvertices,1); %q=1.
    10 %Take care of iceshelves: no basal drag
    11 pos=find(md.mask.elementonfloatingice);
    12 md.friction.coefficient(md.mesh.elements(pos,:))=0.;
     9md.friction.coefficient=200.*ones(md.mesh.numberofvertices,1);
     10md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    1311md.friction.p=ones(md.mesh.numberofelements,1);
    1412md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/ISMIPB.py

    r15771 r15988  
    1010
    1111print "      creating drag"
    12 md.friction.coefficient=200.*numpy.ones((md.mesh.numberofvertices,1))    #q=1.
    13 #Take care of iceshelves: no basal drag
    14 pos=numpy.nonzero(md.mask.elementonfloatingice)[0]
    15 md.friction.coefficient[md.mesh.elements[pos,:]-1]=0.
     12md.friction.coefficient=200.*numpy.ones((md.mesh.numberofvertices,1))
     13md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    1614md.friction.p=numpy.ones((md.mesh.numberofelements,1))
    1715md.friction.q=numpy.ones((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/ISMIPC.par

    r15771 r15988  
    99%md.friction.coefficient=sqrt(md.constants.yts.*(1000.+1000.*sin(md.mesh.x*2.*pi/max(md.mesh.x/2.)).*sin(md.mesh.y*2.*pi/max(md.mesh.x/2.)))./(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.bed)));
    1010md.friction.coefficient=sqrt(md.constants.yts.*(1000.+1000.*sin(md.mesh.x*2.*pi/max(md.mesh.x)).*sin(md.mesh.y*2.*pi/max(md.mesh.x))));
    11 %Take care of iceshelves: no basal drag
    12 pos=find(md.mask.elementonfloatingice);
    13 md.friction.coefficient(md.mesh.elements(pos,:))=0.;
     11md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    1412md.friction.p=ones(md.mesh.numberofelements,1);
    1513md.friction.q=zeros(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/ISMIPC.py

    r15771 r15988  
    1212#md.friction.coefficient=sqrt(md.constants.yts.*(1000.+1000.*sin(md.mesh.x*2.*pi/max(md.mesh.x/2.)).*sin(md.mesh.y*2.*pi/max(md.mesh.x/2.)))./(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.bed)));
    1313md.friction.coefficient=numpy.sqrt(md.constants.yts*(1000.+1000.*numpy.sin(md.mesh.x.reshape(-1,1)*2.*numpy.pi/numpy.max(md.mesh.x))*numpy.sin(md.mesh.y.reshape(-1,1)*2.*numpy.pi/numpy.max(md.mesh.x))))
    14 #Take care of iceshelves: no basal drag
    15 pos=numpy.nonzero(md.mask.elementonfloatingice)[0]
    16 md.friction.coefficient[md.mesh.elements[pos,:]-1]=0.
     14md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    1715md.friction.p=numpy.ones((md.mesh.numberofelements,1))
    1816md.friction.q=numpy.zeros((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/ISMIPD.par

    r15771 r15988  
    88disp('      creating drag');
    99md.friction.coefficient=sqrt(md.constants.yts.*(1000.+1000.*sin(md.mesh.x*2.*pi/max(md.mesh.x))));
    10 %Take care of iceshelves: no basal drag
    11 pos=find(md.mask.elementonfloatingice);
    12 md.friction.coefficient(md.mesh.elements(pos,:))=0.;
     10md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    1311md.friction.p=ones(md.mesh.numberofelements,1);
    1412md.friction.q=zeros(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/ISMIPD.py

    r15771 r15988  
    1111print "      creating drag"
    1212md.friction.coefficient=numpy.sqrt(md.constants.yts*(1000.+1000.*numpy.sin(md.mesh.x.reshape(-1,1)*2.*numpy.pi/numpy.max(md.mesh.x))))
    13 #Take care of iceshelves: no basal drag
    14 pos=numpy.nonzero(md.mask.elementonfloatingice)[0]
    15 md.friction.coefficient[md.mesh.elements[pos,:]-1]=0.
     13md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    1614md.friction.p=numpy.ones((md.mesh.numberofelements,1))
    1715md.friction.q=numpy.zeros((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/Pig.par

    r15771 r15988  
    2727
    2828%Friction
    29 pos=find(md.mask.elementonfloatingice);
    3029md.friction.coefficient=50*ones(md.mesh.numberofvertices,1);
    31 md.friction.coefficient(md.mesh.elements(pos,:))=0;
     30md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    3231md.friction.p=ones(md.mesh.numberofelements,1);
    3332md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/Pig.py

    r15771 r15988  
    3838
    3939#Friction
    40 pos=numpy.nonzero(md.mask.elementonfloatingice)
    4140md.friction.coefficient=50.*ones((md.mesh.numberofvertices,1))
    42 md.friction.coefficient[md.mesh.elements[pos].astype(int)-1]=0.
     41md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    4342md.friction.p=ones((md.mesh.numberofelements,1))
    4443md.friction.q=ones((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/RoundSheetEISMINT.par

    r15771 r15988  
    66
    77disp('      creating drag');
    8 md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1); %q=1. no drag is specified in the analytical solution
     8md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
    99md.friction.p=ones(md.mesh.numberofelements,1);
    1010md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/RoundSheetEISMINT.py

    r15771 r15988  
    99
    1010print "      creating drag"
    11 md.friction.coefficient=20.*numpy.ones((md.mesh.numberofvertices,1))    #q=1. no drag is specified in the analytical solution
     11md.friction.coefficient=20.*numpy.ones((md.mesh.numberofvertices,1))
    1212md.friction.p=numpy.ones((md.mesh.numberofelements,1))
    1313md.friction.q=numpy.ones((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/RoundSheetShelf.par

    r15950 r15988  
    6868
    6969%Friction
    70 pos=find(md.mask.elementonfloatingice);
    7170md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
    72 md.friction.coefficient(md.mesh.elements(pos,:))=0;
     71md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    7372md.friction.p=ones(md.mesh.numberofelements,1);
    7473md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/RoundSheetShelf.py

    r15948 r15988  
    7676
    7777#Friction
    78 pos=numpy.nonzero(md.mask.elementonfloatingice)
    7978md.friction.coefficient=20.*numpy.ones((md.mesh.numberofvertices,1))
    80 md.friction.coefficient[md.mesh.elements[pos,:].astype(int)-1]=0
     79md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    8180md.friction.p=numpy.ones((md.mesh.numberofelements,1))
    8281md.friction.q=numpy.ones((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/RoundSheetStaticEISMINT.par

    r15771 r15988  
    1010
    1111disp('      creating drag');
    12 md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1); %q=1. no drag is specified in the analytical solution
    13 %Take care of iceshelves: no basal drag
    14 pos=find(md.mask.elementonfloatingice);
    15 md.friction.coefficient(md.mesh.elements(pos,:))=0.;
     12md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
     13md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    1614md.friction.p=ones(md.mesh.numberofelements,1);
    1715md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/RoundSheetStaticEISMINT.py

    r15771 r15988  
    1313
    1414print "      creating drag"
    15 md.friction.coefficient=20.*numpy.ones((md.mesh.numberofvertices,1))    #q=1. no drag is specified in the analytical solution
    16 #Take care of iceshelves: no basal drag
    17 pos=numpy.nonzero(md.mask.elementonfloatingice)[0]
    18 md.friction.coefficient[md.mesh.elements[pos,:]-1]=0.
     15md.friction.coefficient=20.*numpy.ones((md.mesh.numberofvertices,1))
     16md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    1917md.friction.p=numpy.ones((md.mesh.numberofelements,1))
    2018md.friction.q=numpy.ones((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/SquareEISMINT.par

    r15767 r15988  
    99
    1010disp('      creating drag');
    11 md.friction.coefficient=200.*ones(md.mesh.numberofvertices,1); %q=1.
    12 %Take care of iceshelves: no basal drag
    13 pos=find(md.mask.elementonfloatingice);
    14 md.friction.coefficient(md.mesh.elements(pos,:))=0.;
     11md.friction.coefficient=200.*ones(md.mesh.numberofvertices,1);
     12md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    1513md.friction.p=ones(md.mesh.numberofelements,1);
    1614md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/SquareEISMINT.py

    r15767 r15988  
    1212
    1313print "      creating drag"
    14 md.friction.coefficient=200.*numpy.ones((md.mesh.numberofvertices,1))    #q=1.
    15 #Take care of iceshelves: no basal drag
    16 pos=numpy.nonzero(md.mask.elementonfloatingice)[0]
    17 md.friction.coefficient[md.mesh.elements[pos,:]-1]=0.
     14md.friction.coefficient=200.*numpy.ones((md.mesh.numberofvertices,1))
     15md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    1816md.friction.p=numpy.ones((md.mesh.numberofelements,1))
    1917md.friction.q=numpy.ones((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/SquareSheetConstrained.par

    r15771 r15988  
    3030
    3131%Friction
    32 pos=find(md.mask.elementonfloatingice);
    3332md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
    34 md.friction.coefficient(md.mesh.elements(pos,:))=0;
     33md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    3534md.friction.p=ones(md.mesh.numberofelements,1);
    3635md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/SquareSheetConstrained.py

    r15771 r15988  
    4040
    4141#Friction
    42 pos=numpy.nonzero(md.mask.elementonfloatingice)
    4342md.friction.coefficient=20.*numpy.ones((md.mesh.numberofvertices,1))
    44 md.friction.coefficient[md.mesh.elements[pos,:].astype(int)-1]=0.
     43md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    4544md.friction.p=numpy.ones((md.mesh.numberofelements,1))
    4645md.friction.q=numpy.ones((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/SquareSheetShelf.par

    r15771 r15988  
    3535
    3636%Friction
    37 pos=find(md.mask.elementonfloatingice);
    3837md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
    39 md.friction.coefficient(md.mesh.elements(pos,:))=0.;
     38md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    4039md.friction.p=ones(md.mesh.numberofelements,1);
    4140md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/SquareSheetShelf.py

    r15771 r15988  
    4646
    4747#Friction
    48 pos=numpy.nonzero(md.mask.elementonfloatingice)
    4948md.friction.coefficient=20.*ones((md.mesh.numberofvertices,1))
    50 md.friction.coefficient[md.mesh.elements[pos].astype(int)-1]=0.
     49md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    5150md.friction.p=ones((md.mesh.numberofelements,1))
    5251md.friction.q=ones((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/SquareShelf.par

    r15771 r15988  
    2828
    2929%Friction
    30 pos=find(md.mask.elementonfloatingice);
    3130md.friction.coefficient=20*ones(md.mesh.numberofvertices,1);
    32 md.friction.coefficient(md.mesh.elements(pos,:))=0;
     31md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    3332md.friction.p=ones(md.mesh.numberofelements,1);
    3433md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/SquareShelf.py

    r15771 r15988  
    5959
    6060#Friction
    61 pos=nonzero(md.mask.elementonfloatingice)
    6261md.friction.coefficient = 20.*ones((md.mesh.numberofvertices,1))
    63 md.friction.coefficient[md.mesh.elements[pos,:].astype(int)-1] = 0.
     62md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    6463md.friction.p = ones((md.mesh.numberofelements,1))
    6564md.friction.q = ones((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/SquareShelf2.par

    r15771 r15988  
    2828
    2929%Friction
    30 pos=find(md.mask.elementonfloatingice);
    3130md.friction.coefficient=20*ones(md.mesh.numberofvertices,1);
    32 md.friction.coefficient(md.mesh.elements(pos,:))=0;
     31md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    3332md.friction.p=ones(md.mesh.numberofelements,1);
    3433md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/test/Par/SquareShelf2.py

    r15771 r15988  
    5959
    6060#Friction
    61 pos=nonzero(md.mask.elementonfloatingice)
    6261md.friction.coefficient = 20.*ones((md.mesh.numberofvertices,1))
    63 md.friction.coefficient[md.mesh.elements[pos,:].astype(int)-1] = 0.
     62md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    6463md.friction.p = ones((md.mesh.numberofelements,1))
    6564md.friction.q = ones((md.mesh.numberofelements,1))

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

    r15771 r15988  
    3232
    3333%Friction
    34 pos=find(md.mask.elementonfloatingice);
    3534md.friction.coefficient=20*ones(md.mesh.numberofvertices,1);
    36 md.friction.coefficient(md.mesh.elements(pos,:))=0;
     35md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    3736md.friction.p=ones(md.mesh.numberofelements,1);
    3837md.friction.q=ones(md.mesh.numberofelements,1);

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

    r15771 r15988  
    4444
    4545#Friction
    46 pos=numpy.nonzero(md.mask.elementonfloatingice)
    4746md.friction.coefficient=20.*numpy.ones((md.mesh.numberofvertices,1))
    48 md.friction.coefficient[md.mesh.elements[pos,:].astype(int)-1]=0.
     47md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    4948md.friction.p=numpy.ones((md.mesh.numberofelements,1))
    5049md.friction.q=numpy.ones((md.mesh.numberofelements,1))

  • issm/trunk-jpl/test/Par/SquareThermal.par

    r15771 r15988  
    1515
    1616disp('      creating drag');
    17 md.friction.coefficient=200.*ones(md.mesh.numberofvertices,1); %q=1.
    18 %Take care of iceshelves: no basal drag
    19 pos=find(md.mask.elementonfloatingice);
    20 md.friction.coefficient(md.mesh.elements(pos,:))=0.;
     17md.friction.coefficient=200.*ones(md.mesh.numberofvertices,1);
     18md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
    2119md.friction.p=ones(md.mesh.numberofelements,1);
    2220md.friction.q=ones(md.mesh.numberofelements,1);
    … …  
    4139md.thermal.spctemperature(:)=md.initialization.temperature;
    4240md.basalforcings.geothermalflux=zeros(md.mesh.numberofvertices,1);
    43 pos=find(md.mask.elementongroundedice);md.basalforcings.geothermalflux(md.mesh.elements(pos,:))=1.*10^-3; %1 mW/m^2
     41md.basalforcings.geothermalflux(find(md.mask.groundedice_levelset>0.)=1.*10^-3; %1 mW/m^2

  • issm/trunk-jpl/test/Par/SquareThermal.py

    r15771 r15988  
    1919
    2020print "      creating drag"
    21 md.friction.coefficient=200.*numpy.ones((md.mesh.numberofvertices,1))    #q=1.
    22 #Take care of iceshelves: no basal drag
    23 pos=numpy.nonzero(md.mask.elementonfloatingice)[0]
    24 md.friction.coefficient[md.mesh.elements[pos,:]-1]=0.
     21md.friction.coefficient=200.*numpy.ones((md.mesh.numberofvertices,1))
     22md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
    2523md.friction.p=numpy.ones((md.mesh.numberofelements,1))
    2624md.friction.q=numpy.ones((md.mesh.numberofelements,1))
    … …  
    4543md.thermal.spctemperature[:]=md.initialization.temperature
    4644md.basalforcings.geothermalflux=numpy.zeros((md.mesh.numberofvertices,1))
    47 pos=numpy.nonzero(md.mask.elementongroundedice)[0]
    48 md.basalforcings.geothermalflux[md.mesh.elements[pos,:]-1]=1.*10**-3    #1 mW/m^2
     45md.basalforcings.geothermalflux[numpy.nonzero(md.mask.groundedice_levelset>0.)[0]]=1.*10**-3    #1 mW/m^2
Note: See TracChangeset for help on using the changeset viewer.