Changeset 25172


Ignore:
Timestamp:
06/28/20 01:02:57 (5 years ago)
Author:
jdquinn
Message:

BUG: Corrected remaining solid earth failures (still need to identify source of errors)

Location:
issm/trunk-jpl
Files:
5 edited

Legend:

Unmodified
Added
Removed

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

    r25168 r25172  
    12641264                        md.frontalforcings  = frontalforcings();
    12651265                        md.gia              = giamme();
    1266                         md.esa              = esa();
    12671266                        md.love             = fourierlove();
    12681267                        md.esa              = esa();

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

    r25165 r25172  
    8787        self.mesh = mesh2d()
    8888        self.mask = mask()
     89        self.constants = constants()
    8990        self.geometry = geometry()
    90         self.constants = constants()
     91        self.initialization = initialization()
    9192        self.smb = SMBforcing()
    9293        self.basalforcings = basalforcings()
     94        self.friction = friction()
     95        self.rifts = rifts()
     96        self.solidearth = solidearth()
     97        self.dsl = dsl()
     98        self.timestepping = timestepping()
     99        self.groundingline = groundingline()
    93100        self.materials = matice()
    94101        self.damage = damage()
    95         self.friction = friction()
    96102        self.flowequation = flowequation()
    97         self.timestepping = timestepping()
    98         self.initialization = initialization()
    99         self.rifts = rifts()
    100         self.dsl = dsl()
    101         self.solidearth = solidearth()
    102 
    103103        self.debug = debug()
    104104        self.verbose = verbose()
    … …  
    106106        self.toolkits = toolkits()
    107107        self.cluster = generic()
    108 
    109108        self.balancethickness = balancethickness()
    110109        self.stressbalance = stressbalance()
    111         self.groundingline = groundingline()
    112110        self.hydrology = hydrologyshreve()
    113111        self.masstransport = masstransport()
    … …  
    125123        self.qmu = qmu()
    126124        self.amr = amr()
    127 
     125        self.radaroverlay = radaroverlay()
    128126        self.results = results()
    129127        self.outputdefinition = outputdefinition()
    130         self.radaroverlay = radaroverlay()
    131128        self.miscellaneous = miscellaneous()
    132129        self.private = private()
    … …  
    137134        return ['mesh',
    138135                'mask',
     136                'constants',
    139137                'geometry',
    140                 'constants',
     138                'initialization',
    141139                'smb',
    142140                'basalforcings',
     141                'friction',
     142                'rifts',
     143                'solidearth',
     144                'dsl',
     145                'timestepping',
     146                'groundingline',
    143147                'materials',
    144148                'damage',
    145                 'friction',
    146149                'flowequation',
    147                 'timestepping',
    148                 'initialization',
    149                 'rifts',
    150                 'dsl',
    151                 'solidearth',
    152150                'debug',
    153151                'verbose',
    … …  
    157155                'balancethickness',
    158156                'stressbalance',
    159                 'groundingline',
    160157                'hydrology',
    161158                'masstransport',
    … …  
    166163                'calving',
    167164                'frontalforcings',
     165                'gia',
    168166                'love',
    169                 'gia',
    170167                'esa',
    171168                'autodiff',
    … …  
    173170                'qmu',
    174171                'amr',
     172                'radaroverlay',
     173                'results',
    175174                'outputdefinition',
    176                 'results',
    177                 'radaroverlay',
    178175                'miscellaneous',
    179176                'private']

  • issm/trunk-jpl/test/NightlyRun/test2071.py

    r25131 r25172  
    1717
    1818#indicate what you want to compute
    19 md.gia.cross_section_shape = 1  # for square-edged x - section
     19md.gia.cross_section_shape = 1  # for square-edged x-section
    2020
    2121#define loading history
    2222md.timestepping.start_time = 0.3  # for t \approx 0 kyr : to get eleastic response!
    23 md.timestepping.final_time = 2500000  # 2, 500 kyr
     23md.timestepping.final_time = 2500000  # 2,500 kyr
    2424md.geometry.thickness = np.array([np.append(md.geometry.thickness * 0.0, 0.0),
    2525                                  np.append(md.geometry.thickness / 2.0, 0.1),

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

    r25129 r25172  
    1515nv = md.mesh.numberofvertices
    1616if (np.isnan(md.geometry.thickness)):
    17     md.geometry.thickness = np.zeros((md.mesh.numberofvertices, ))
     17    md.geometry.thickness = np.zeros((md.mesh.numberofvertices, 1))
    1818for i in range(nv):
    1919    dist = np.sqrt(md.mesh.x[i]**2 + md.mesh.y[i]**2)
    … …  
    2121        md.geometry.thickness[i] = 2000.0
    2222    else:
    23         md.geometry.thickness[i] = 1.0  # non - zero thickness
     23        md.geometry.thickness[i] = 1.0  # non-zero thickness
    2424
    2525md.geometry.thickness = md.geometry.thickness.reshape(-1, 1)
    26 md.geometry.base = np.zeros((md.mesh.numberofvertices, ))
     26md.geometry.base = np.zeros((md.mesh.numberofvertices, 1))
    2727md.geometry.surface = md.geometry.thickness + md.geometry.base.reshape(-1, 1)  #would otherwise create a 91x91 matrix
    2828
    29 #Ice density used for benchmarking, not 917 kg / m^3
    30 md.materials.rho_ice = 1000  #kg m^3
     29#Ice density used for benchmarking, not 917 kg/m^3
     30md.materials.rho_ice = 1000  #kg/m^3
    3131
    32 #GIA parameters specific to Experiments A  and B
     32#GIA parameters specific to Experiments A and B
    3333md.gia=giaivins();
    34 md.gia.mantle_viscosity = 1e21 * np.ones((md.mesh.numberofvertices, ))  #in Pa.s
    35 md.gia.lithosphere_thickness = 100 * np.ones((md.mesh.numberofvertices, ))  #in km
     34md.gia.mantle_viscosity = 1e21 * np.ones((md.mesh.numberofvertices, 1))  #in Pa.s
     35md.gia.lithosphere_thickness = 100 * np.ones((md.mesh.numberofvertices, 1))  #in km
    3636md.materials.lithosphere_shear_modulus = 6.7 * 1e10  #in Pa
    37 md.materials.lithosphere_density = 3.36  #in g / cm^3
     37md.materials.lithosphere_density = 3.36  #in g/cm^3
    3838md.materials.mantle_shear_modulus = 1.45 * 1e11  #in Pa
    39 md.materials.mantle_density = 3.38  #in g / cm^3
     39md.materials.mantle_density = 3.38  #in g/cm^3
    4040
    4141#Initial velocity
    … …  
    5353y = None
    5454index = None
    55 md.initialization.vz = np.zeros((md.mesh.numberofvertices, ))
    56 md.initialization.pressure = np.zeros((md.mesh.numberofvertices, ))
     55md.initialization.vz = np.zeros((md.mesh.numberofvertices, 1))
     56md.initialization.pressure = np.zeros((md.mesh.numberofvertices, 1))
    5757
    5858#Materials
    59 md.initialization.temperature = (273. - 20.) * np.ones((md.mesh.numberofvertices, ))
     59md.initialization.temperature = (273. - 20.) * np.ones((md.mesh.numberofvertices, 1))
    6060md.materials.rheology_B = paterson(md.initialization.temperature)
    61 md.materials.rheology_n = 3. * np.ones((md.mesh.numberofelements, ))
     61md.materials.rheology_n = 3. * np.ones((md.mesh.numberofelements, 1))
    6262
    6363#Friction
    64 md.friction.coefficient = 20. * np.ones((md.mesh.numberofvertices, ))
     64md.friction.coefficient = 20. * np.ones((md.mesh.numberofvertices, 1))
    6565md.friction.coefficient[np.where(md.mask.ocean_levelset < 0.)] = 0.
    66 md.friction.p = np.ones((md.mesh.numberofelements, ))
    67 md.friction.q = np.ones((md.mesh.numberofelements, ))
     66md.friction.p = np.ones((md.mesh.numberofelements, 1))
     67md.friction.q = np.ones((md.mesh.numberofelements, 1))
    6868
    6969#Numerical parameters
    … …  
    7676md.steadystate.reltol = 0.05
    7777md.stressbalance.reltol = 0.05
    78 md.stressbalance.abstol = float('NaN')
     78md.stressbalance.abstol = np.nan
    7979md.timestepping.time_step = 1.
    8080md.timestepping.final_time = 3.

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

    r25129 r25172  
    11#Geometry specific to Experiments C and D
    22
     3import inspect
    34import os.path
    4 import inspect
     5
    56import numpy as np
     7
    68from arch import *
     9from giaivins import giaivins
    710from InterpFromMeshToMesh2d import *
    811from paterson import *
    … …  
    1013from SetIceSheetBC import *
    1114
     15
    1216rad = 800000
    1317nv = md.mesh.numberofvertices
    1418if (np.isnan(md.geometry.thickness)):
    15     md.geometry.thickness = np.zeros((md.mesh.numberofvertices, ))
     19    md.geometry.thickness = np.zeros((md.mesh.numberofvertices, 1))
    1620for i in range(nv):
    1721    dist = np.sqrt(md.mesh.x[i]**2 + md.mesh.y[i]**2)
    … …  
    2226
    2327md.geometry.thickness = md.geometry.thickness.reshape(-1, 1)
    24 md.geometry.base = np.zeros((md.mesh.numberofvertices, ))
     28md.geometry.base = np.zeros((md.mesh.numberofvertices, 1))
    2529md.geometry.surface = md.geometry.thickness + md.geometry.base.reshape(-1, 1)  #would otherwise create a 91x91 matrix
    2630
    2731#Ice density used for benchmarking, not 917 kg / m^3
    28 md.materials.rho_ice = 1000  #kg m^3
     32md.materials.rho_ice = 1000 #kg m^3
    2933
    30 #GIA parameters specific to Experiments A  and B
    31 md.gia.mantle_viscosity = 1e21 * np.ones((md.mesh.numberofvertices, ))  #in Pa.s
    32 md.gia.lithosphere_thickness = 100 * np.ones((md.mesh.numberofvertices, ))  #in km
     34#GIA parameters specific to Experiments A and B
     35md.gia=giaivins();
     36md.gia.mantle_viscosity = 1e21 * np.ones((md.mesh.numberofvertices, 1))  #in Pa.s
     37md.gia.lithosphere_thickness = 100 * np.ones((md.mesh.numberofvertices, 1))  #in km
    3338md.materials.lithosphere_shear_modulus = 6.7e10  #in Pa
    3439md.materials.lithosphere_density = 3.32  #in g / cm^3
    … …  
    5055y = None
    5156index = None
    52 md.initialization.vz = np.zeros((md.mesh.numberofvertices, ))
    53 md.initialization.pressure = np.zeros((md.mesh.numberofvertices, ))
     57md.initialization.vz = np.zeros((md.mesh.numberofvertices, 1))
     58md.initialization.pressure = np.zeros((md.mesh.numberofvertices, 1))
    5459
    5560#Materials
    56 md.initialization.temperature = (273. - 20.) * np.ones((md.mesh.numberofvertices, ))
     61md.initialization.temperature = (273. - 20.) * np.ones((md.mesh.numberofvertices, 1))
    5762md.materials.rheology_B = paterson(md.initialization.temperature)
    58 md.materials.rheology_n = 3. * np.ones((md.mesh.numberofelements, ))
     63md.materials.rheology_n = 3. * np.ones((md.mesh.numberofelements, 1))
    5964
    6065#Friction
    61 md.friction.coefficient = 20. * np.ones((md.mesh.numberofvertices, ))
     66md.friction.coefficient = 20. * np.ones((md.mesh.numberofvertices, 1))
    6267md.friction.coefficient[np.where(md.mask.ocean_levelset < 0.)] = 0.
    63 md.friction.p = np.ones((md.mesh.numberofelements, ))
    64 md.friction.q = np.ones((md.mesh.numberofelements, ))
     68md.friction.p = np.ones((md.mesh.numberofelements, 1))
     69md.friction.q = np.ones((md.mesh.numberofelements, 1))
    6570
    6671#Numerical parameters
    … …  
    7378md.steadystate.reltol = 0.05
    7479md.stressbalance.reltol = 0.05
    75 md.stressbalance.abstol = float('NaN')
     80md.stressbalance.abstol = np.nan
    7681md.timestepping.time_step = 1.
    7782md.timestepping.final_time = 3.
Note: See TracChangeset for help on using the changeset viewer.