Ignore:
Timestamp:
03/13/19 03:17:46 (6 years ago)
Author:
bdef
Message:

pep8 compliance of NTs

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/test/NightlyRun/test442.py

    r23004 r23793  
    99from solve import *
    1010
    11 md = triangle(model(),'../Exp/Square.exp',100000.)
    12 md = setmask(md,'../Exp/SquareShelf.exp','')
    13 md = parameterize(md,'../Par/SquareSheetShelf.py')
     11md = triangle(model(), '../Exp/Square.exp', 100000.)
     12md = setmask(md, '../Exp/SquareShelf.exp', '')
     13md = parameterize(md, '../Par/SquareSheetShelf.py')
    1414md.initialization.vx[:] = 1.
    1515md.initialization.vy[:] = 1.
     
    2121md.geometry.base[pos] = md.geometry.bed[pos]
    2222md.geometry.surface = md.geometry.base + md.geometry.thickness
    23 md = md.extrude(4,1.)
    24 md = setflowequation(md,'HO','all')
     23md = md.extrude(4, 1.)
     24md = setflowequation(md, 'HO', 'all')
    2525
    2626#Boundary conditions:
     
    3030md.stressbalance.spcvy[:] = float('Nan')
    3131md.stressbalance.spcvz[:] = float('Nan')
    32 posA = np.intersect1d(np.array(np.where(md.mesh.y < 1000000.1)),np.array(np.where(md.mesh.y > 999999.9)))
    33 posB = np.intersect1d(np.array(np.where(md.mesh.y < 0.1)),np.array(np.where(md.mesh.y > -0.1)))
    34 pos = np.unique(np.concatenate((posA,posB)))
     32posA = np.intersect1d(np.array(np.where(md.mesh.y < 1000000.1)), np.array(np.where(md.mesh.y > 999999.9)))
     33posB = np.intersect1d(np.array(np.where(md.mesh.y < 0.1)), np.array(np.where(md.mesh.y > -0.1)))
     34pos = np.unique(np.concatenate((posA, posB)))
    3535md.stressbalance.spcvy[pos] = 0.
    3636pos2 = np.intersect1d(np.array(np.where(md.mesh.x < 0.1)), np.array(np.where(md.mesh.x > -0.1)))
     
    3838md.stressbalance.spcvy[pos2] = 0.
    3939
    40 md.materials.rheology_B = 1. / ((10**-25)**(1./3.)) * np.ones((md.mesh.numberofvertices,))
     40md.materials.rheology_B = 1. / ((10**-25)**(1. / 3.)) * np.ones((md.mesh.numberofvertices,))
    4141md.materials.rheology_law = 'None'
    4242md.friction.coefficient[:] = np.sqrt(1e7) * np.ones((md.mesh.numberofvertices,))
     
    5050md.transient.ismasstransport = 1
    5151md.transient.issmb = 1
    52 md.transient.requested_outputs = ['default','BasalforcingsFloatingiceMeltingRate']
     52md.transient.requested_outputs = ['default', 'BasalforcingsFloatingiceMeltingRate']
    5353md.groundingline.migration = 'SubelementMigration'
    54 md.groundingline.friction_interpolation= 'SubelementFriction2'
    55 md.groundingline.melt_interpolation= 'SubelementMelt2'
     54md.groundingline.friction_interpolation = 'SubelementFriction2'
     55md.groundingline.melt_interpolation = 'SubelementMelt2'
    5656md.timestepping.final_time = 30
    5757md.timestepping.time_step = 10
    5858
    59 md.cluster = generic('name',gethostname(),'np',3)
    60 md = solve(md,'Transient')
     59md.cluster = generic('name', gethostname(), 'np', 3)
     60md = solve(md, 'Transient')
    6161
    6262#Fields and tolerances to track changes
    63 field_names     = [
    64         'Bed1','Surface1','Thickness1','Floatingice1','Vx1','Vy1','Vz1','Pressure1','FloatingiceMeltingrate1',
    65         'Bed2','Surface2','Thickness2','Floatingice2','Vx2','Vy2','Vz2','Pressure2','FloatingiceMeltingrate2',
    66         'Bed3','Surface3','Thickness3','Floatingice3','Vx3','Vy3','Vz3','Pressure3','FloatingiceMeltingrate3']
    67 field_tolerances = [
    68         2e-11,5e-12,2e-11,1e-11,5e-10,3e-08,6e-10,1e-13,1e-13,
    69         3e-11,3e-11,9e-10,7e-11,7e-09,1e-07,1e-09,1e-10,1e-13,
    70         1e-8,2e-08,7e-09,2e-7 ,1e-03,8e-04,2e-09,1e-10,1e-13]
    71 field_values = [
    72         md.results.TransientSolution[0].Base,
    73         md.results.TransientSolution[0].Surface,
    74         md.results.TransientSolution[0].Thickness,
    75         md.results.TransientSolution[0].MaskGroundediceLevelset,
    76         md.results.TransientSolution[0].Vx,
    77         md.results.TransientSolution[0].Vy,
    78         md.results.TransientSolution[0].Vz,
    79         md.results.TransientSolution[0].Pressure,
    80         md.results.TransientSolution[0].BasalforcingsFloatingiceMeltingRate,
    81         md.results.TransientSolution[1].Base,
    82         md.results.TransientSolution[1].Surface,
    83         md.results.TransientSolution[1].Thickness,
    84         md.results.TransientSolution[1].MaskGroundediceLevelset,
    85         md.results.TransientSolution[1].Vx,
    86         md.results.TransientSolution[1].Vy,
    87         md.results.TransientSolution[1].Vz,
    88         md.results.TransientSolution[1].Pressure,
    89         md.results.TransientSolution[1].BasalforcingsFloatingiceMeltingRate,
    90         md.results.TransientSolution[2].Base,
    91         md.results.TransientSolution[2].Surface,
    92         md.results.TransientSolution[2].Thickness,
    93         md.results.TransientSolution[2].MaskGroundediceLevelset,
    94         md.results.TransientSolution[2].Vx,
    95         md.results.TransientSolution[2].Vy,
    96         md.results.TransientSolution[2].Vz,
    97         md.results.TransientSolution[2].Pressure,
    98         md.results.TransientSolution[2].BasalforcingsFloatingiceMeltingRate,
    99         ]
     63field_names = ['Bed1', 'Surface1', 'Thickness1', 'Floatingice1', 'Vx1', 'Vy1', 'Vz1', 'Pressure1', 'FloatingiceMeltingrate1',
     64               'Bed2', 'Surface2', 'Thickness2', 'Floatingice2', 'Vx2', 'Vy2', 'Vz2', 'Pressure2', 'FloatingiceMeltingrate2',
     65               'Bed3', 'Surface3', 'Thickness3', 'Floatingice3', 'Vx3', 'Vy3', 'Vz3', 'Pressure3', 'FloatingiceMeltingrate3']
     66field_tolerances = [2e-11, 5e-12, 2e-11, 1e-11, 5e-10, 3e-08, 6e-10, 1e-13, 1e-13,
     67                    3e-11, 3e-11, 9e-10, 7e-11, 7e-09, 1e-07, 1e-09, 1e-10, 1e-13,
     68                    1e-8, 2e-08, 7e-09, 2e-7, 1e-03, 8e-04, 2e-09, 1e-10, 1e-13]
     69field_values = [md.results.TransientSolution[0].Base,
     70                md.results.TransientSolution[0].Surface,
     71                md.results.TransientSolution[0].Thickness,
     72                md.results.TransientSolution[0].MaskGroundediceLevelset,
     73                md.results.TransientSolution[0].Vx,
     74                md.results.TransientSolution[0].Vy,
     75                md.results.TransientSolution[0].Vz,
     76                md.results.TransientSolution[0].Pressure,
     77                md.results.TransientSolution[0].BasalforcingsFloatingiceMeltingRate,
     78                md.results.TransientSolution[1].Base,
     79                md.results.TransientSolution[1].Surface,
     80                md.results.TransientSolution[1].Thickness,
     81                md.results.TransientSolution[1].MaskGroundediceLevelset,
     82                md.results.TransientSolution[1].Vx,
     83                md.results.TransientSolution[1].Vy,
     84                md.results.TransientSolution[1].Vz,
     85                md.results.TransientSolution[1].Pressure,
     86                md.results.TransientSolution[1].BasalforcingsFloatingiceMeltingRate,
     87                md.results.TransientSolution[2].Base,
     88                md.results.TransientSolution[2].Surface,
     89                md.results.TransientSolution[2].Thickness,
     90                md.results.TransientSolution[2].MaskGroundediceLevelset,
     91                md.results.TransientSolution[2].Vx,
     92                md.results.TransientSolution[2].Vy,
     93                md.results.TransientSolution[2].Vz,
     94                md.results.TransientSolution[2].Pressure,
     95                md.results.TransientSolution[2].BasalforcingsFloatingiceMeltingRate]
Note: See TracChangeset for help on using the changeset viewer.