Ignore:
Timestamp:
10/11/19 00:27:00 (5 years ago)
Author:
bdef
Message:

CHG: syntax cahnge to meet most of Pep8 requirement

File:
1 edited

Legend:

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

    r23793 r24214  
    1515md.initialization.vy[:] = 1.
    1616md.geometry.thickness[:] = 500. - md.mesh.x / 10000.
    17 md.geometry.bed = -100. - md.mesh.x / 1000.
    18 md.geometry.base = -md.geometry.thickness * md.materials.rho_ice / md.materials.rho_water
     17md.geometry.bed = - 100. - md.mesh.x / 1000.
     18md.geometry.base = - md.geometry.thickness * md.materials.rho_ice / md.materials.rho_water
    1919md.mask.groundedice_levelset = md.geometry.thickness + md.materials.rho_water / md.materials.rho_ice * md.geometry.bed
    2020pos = np.where(md.mask.groundedice_levelset >= 0.)
     
    2525
    2626#Boundary conditions:
    27 md.mask.ice_levelset = -np.ones((md.mesh.numberofvertices,))
     27md.mask.ice_levelset = - np.ones((md.mesh.numberofvertices, ))
    2828md.mask.ice_levelset[np.where(md.mesh.x == max(md.mesh.x))] = 0.
    2929md.stressbalance.spcvx[:] = float('Nan')
     
    3131md.stressbalance.spcvz[:] = float('Nan')
    3232posA = 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)))
     33posB = np.intersect1d(np.array(np.where(md.mesh.y < 0.1)), np.array(np.where(md.mesh.y > - 0.1)))
    3434pos = np.unique(np.concatenate((posA, posB)))
    3535md.stressbalance.spcvy[pos] = 0.
    36 pos2 = np.intersect1d(np.array(np.where(md.mesh.x < 0.1)), np.array(np.where(md.mesh.x > -0.1)))
     36pos2 = np.intersect1d(np.array(np.where(md.mesh.x < 0.1)), np.array(np.where(md.mesh.x > - 0.1)))
    3737md.stressbalance.spcvx[pos2] = 0.
    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'
    42 md.friction.coefficient[:] = np.sqrt(1e7) * np.ones((md.mesh.numberofvertices,))
    43 md.friction.p = 3. * np.ones((md.mesh.numberofelements,))
     42md.friction.coefficient[:] = np.sqrt(1e7) * np.ones((md.mesh.numberofvertices, ))
     43md.friction.p = 3. * np.ones((md.mesh.numberofelements, ))
    4444md.smb.mass_balance[:] = 1.
    4545md.basalforcings.groundedice_melting_rate[:] = 0.
Note: See TracChangeset for help on using the changeset viewer.