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/test344.py

    r22267 r23793  
    1010from SMBgradientsela import *
    1111
    12 md = triangle(model(),'../Exp/Square.exp',150000.)
    13 md = setmask(md,'','')
    14 md = parameterize(md,'../Par/SquareSheetConstrained.py')
     12md = triangle(model(), '../Exp/Square.exp', 150000.)
     13md = setmask(md, '', '')
     14md = parameterize(md, '../Par/SquareSheetConstrained.py')
    1515
    1616#Change geometry
     
    1818md.geometry.surface = md.geometry.base + md.geometry.thickness
    1919
    20 md = md.extrude(3,1.)
    21 md = setflowequation(md,'HO','all')
     20md = md.extrude(3, 1.)
     21md = setflowequation(md, 'HO', 'all')
    2222md.smb = SMBgradientsela()
    23 md.smb.ela = 1500. * np.ones((md.mesh.numberofvertices+1,))
    24 md.smb.b_pos = 0.002 * np.ones((md.mesh.numberofvertices+1,))
    25 md.smb.b_neg = 0.005 * np.ones((md.mesh.numberofvertices+1,))
    26 md.smb.b_max = 4. * (md.materials.rho_freshwater / md.materials.rho_ice) * np.ones((md.mesh.numberofvertices+1,))
    27 md.smb.b_min = -4. * (md.materials.rho_freshwater / md.materials.rho_ice) * np.ones((md.mesh.numberofvertices+1,))
     23md.smb.ela = 1500. * np.ones((md.mesh.numberofvertices + 1,))
     24md.smb.b_pos = 0.002 * np.ones((md.mesh.numberofvertices + 1,))
     25md.smb.b_neg = 0.005 * np.ones((md.mesh.numberofvertices + 1,))
     26md.smb.b_max = 4. * (md.materials.rho_freshwater / md.materials.rho_ice) * np.ones((md.mesh.numberofvertices + 1,))
     27md.smb.b_min = -4. * (md.materials.rho_freshwater / md.materials.rho_ice) * np.ones((md.mesh.numberofvertices + 1,))
    2828
    2929
    3030#Transient options
    31 md.transient.requested_outputs = ['default','TotalSmb']
    32 md.cluster = generic('name',gethostname(),'np',3)
    33 md = solve(md,'Transient')
     31md.transient.requested_outputs = ['default', 'TotalSmb']
     32md.cluster = generic('name', gethostname(), 'np', 3)
     33md = solve(md, 'Transient')
    3434
    3535#Fields and tolerances to track changes
    36 field_names     = ['Vx1','Vy1','Vz1','Vel1','Bed1','Surface1','Thickness1','Temperature1','SMB1','TotalSmb1',
    37  'Vx2','Vy2','Vz2','Vel2','Bed2','Surface2','Thickness2','Temperature2','SMB2','TotalSmb2',
    38  'Vx3','Vy3','Vz3','Vel3','Bed3','Surface3','Thickness3','Temperature3','SMB3','TotalSmb3']
    39 field_tolerances = [1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,
    40         1e-09,1e-09,1e-10,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,
    41         1e-09,5e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10]
    42 field_values = [
    43         md.results.TransientSolution[0].Vx,
    44         md.results.TransientSolution[0].Vy,
    45         md.results.TransientSolution[0].Vz,
    46         md.results.TransientSolution[0].Vel,
    47         md.results.TransientSolution[0].Base,
    48         md.results.TransientSolution[0].Surface,
    49         md.results.TransientSolution[0].Thickness,
    50         md.results.TransientSolution[0].Temperature,
    51         md.results.TransientSolution[0].SmbMassBalance,
    52         md.results.TransientSolution[0].TotalSmb,
    53         md.results.TransientSolution[1].Vx,
    54         md.results.TransientSolution[1].Vy,
    55         md.results.TransientSolution[1].Vz,
    56         md.results.TransientSolution[1].Vel,
    57         md.results.TransientSolution[1].Base,
    58         md.results.TransientSolution[1].Surface,
    59         md.results.TransientSolution[1].Thickness,
    60         md.results.TransientSolution[1].Temperature,
    61         md.results.TransientSolution[1].SmbMassBalance,
    62         md.results.TransientSolution[1].TotalSmb,
    63         md.results.TransientSolution[2].Vx,
    64         md.results.TransientSolution[2].Vy,
    65         md.results.TransientSolution[2].Vz,
    66         md.results.TransientSolution[2].Vel,
    67         md.results.TransientSolution[2].Base,
    68         md.results.TransientSolution[2].Surface,
    69         md.results.TransientSolution[2].Thickness,
    70         md.results.TransientSolution[2].Temperature,
    71         md.results.TransientSolution[2].SmbMassBalance,
    72         md.results.TransientSolution[2].TotalSmb,
    73         ]
     36field_names = ['Vx1', 'Vy1', 'Vz1', 'Vel1', 'Bed1', 'Surface1', 'Thickness1', 'Temperature1', 'SMB1', 'TotalSmb1',
     37               'Vx2', 'Vy2', 'Vz2', 'Vel2', 'Bed2', 'Surface2', 'Thickness2', 'Temperature2', 'SMB2', 'TotalSmb2',
     38               'Vx3', 'Vy3', 'Vz3', 'Vel3', 'Bed3', 'Surface3', 'Thickness3', 'Temperature3', 'SMB3', 'TotalSmb3']
     39field_tolerances = [1e-09, 1e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10,
     40                    1e-09, 1e-09, 1e-10, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10,
     41                    1e-09, 5e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10]
     42field_values = [md.results.TransientSolution[0].Vx,
     43                md.results.TransientSolution[0].Vy,
     44                md.results.TransientSolution[0].Vz,
     45                md.results.TransientSolution[0].Vel,
     46                md.results.TransientSolution[0].Base,
     47                md.results.TransientSolution[0].Surface,
     48                md.results.TransientSolution[0].Thickness,
     49                md.results.TransientSolution[0].Temperature,
     50                md.results.TransientSolution[0].SmbMassBalance,
     51                md.results.TransientSolution[0].TotalSmb,
     52                md.results.TransientSolution[1].Vx,
     53                md.results.TransientSolution[1].Vy,
     54                md.results.TransientSolution[1].Vz,
     55                md.results.TransientSolution[1].Vel,
     56                md.results.TransientSolution[1].Base,
     57                md.results.TransientSolution[1].Surface,
     58                md.results.TransientSolution[1].Thickness,
     59                md.results.TransientSolution[1].Temperature,
     60                md.results.TransientSolution[1].SmbMassBalance,
     61                md.results.TransientSolution[1].TotalSmb,
     62                md.results.TransientSolution[2].Vx,
     63                md.results.TransientSolution[2].Vy,
     64                md.results.TransientSolution[2].Vz,
     65                md.results.TransientSolution[2].Vel,
     66                md.results.TransientSolution[2].Base,
     67                md.results.TransientSolution[2].Surface,
     68                md.results.TransientSolution[2].Thickness,
     69                md.results.TransientSolution[2].Temperature,
     70                md.results.TransientSolution[2].SmbMassBalance,
     71                md.results.TransientSolution[2].TotalSmb]
Note: See TracChangeset for help on using the changeset viewer.