Changeset 26906


Ignore:
Timestamp:
02/25/22 00:30:49 (3 years ago)
Author:
bdef
Message:

CHG: syntax fix

Location:
issm/trunk-jpl/test/NightlyRun
Files:
7 edited

Legend:

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

    r26402 r26906  
    1 #Test Name: SquareShelfTranIspddIsdeltaO18pdNoInterpSSA2d 
     1#Test Name: SquareShelfTranIspddIsdeltaO18pdNoInterpSSA2d
    22import numpy as np
    33from model import *
     
    5757tim1 = np.linspace(1,12,12)/12
    5858
    59 md.smb.temperatures_reconstructed[md.mesh.numberofvertices,0:12] = tim1 
     59md.smb.temperatures_reconstructed[md.mesh.numberofvertices,0:12] = tim1
    6060md.smb.temperatures_reconstructed[md.mesh.numberofvertices,12:24] = tim1+3
    6161md.smb.temperatures_reconstructed[md.mesh.numberofvertices,24:36] = tim1+5
    6262
    6363md.smb.precipitations_reconstructed[md.mesh.numberofvertices,0:12] = tim1
    64 md.smb.precipitations_reconstructed[md.mesh.numberofvertices,12:24] = tim1+3 
     64md.smb.precipitations_reconstructed[md.mesh.numberofvertices,12:24] = tim1+3
    6565md.smb.precipitations_reconstructed[md.mesh.numberofvertices,24:36] = tim1+5
    6666
     
    6868# spc
    6969md.thermal.spctemperature=np.mean(md.smb.temperatures_presentday[0:md.mesh.numberofvertices,0:12],axis=1)-10.
    70 md.initialization.temperature=md.thermal.spctemperature 
     70md.initialization.temperature=md.thermal.spctemperature
    7171
    7272md.smb.s0p = np.maximum(md.geometry.surface.reshape(-1, 1),np.zeros((md.mesh.numberofvertices,1)))
     
    8989md.thermal.isdrainicecolumn = 0
    9090
    91 md=solve(md,'thermal')
     91md = solve(md, 'thermal')
    9292
    9393md.initialization.temperature = md.results.ThermalSolution.Temperature
     
    121121
    122122#Fields and tolerances to track changes
    123 field_names = ['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','Enthalpy1','SmbMassBalance1',
    124    'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','Enthalpy2','SmbMassBalance2',
    125    'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','Enthalpy3','SmbMassBalance3',
    126    'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','Temperature4','Enthalpy4','SmbMassBalance4']
    127 field_tolerances=[1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,
    128    1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,
    129    1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,
    130    1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13]
     123field_names = ['Vx1', 'Vy1', 'Vz1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'Temperature1', 'Enthalpy1', 'SmbMassBalance1',
     124               'Vx2', 'Vy2', 'Vz2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'Temperature2', 'Enthalpy2', 'SmbMassBalance2',
     125               'Vx3', 'Vy3', 'Vz3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'Temperature3', 'Enthalpy3', 'SmbMassBalance3',
     126               'Vx4', 'Vy4', 'Vz4', 'Vel4', 'Pressure4', 'Bed4', 'Surface4', 'Thickness4', 'Temperature4', 'Enthalpy4', 'SmbMassBalance4']
     127field_tolerances=[1e-09, 1e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-13,
     128   1e-09, 1e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-13,
     129   1e-09, 1e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-13,
     130   1e-09, 1e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-13]
    131131field_values = [md.results.TransientSolution[0].Vx,
    132    md.results.TransientSolution[0].Vy,
    133    md.results.TransientSolution[0].Vz,
    134    md.results.TransientSolution[0].Vel,
    135    md.results.TransientSolution[0].Pressure,
    136    md.results.TransientSolution[0].Base,
    137    md.results.TransientSolution[0].Surface,
    138    md.results.TransientSolution[0].Thickness,
    139    md.results.TransientSolution[0].Temperature,
    140    md.results.TransientSolution[0].Enthalpy,
    141    md.results.TransientSolution[0].SmbMassBalance,
    142    md.results.TransientSolution[1].Vx,
    143    md.results.TransientSolution[1].Vy,
    144    md.results.TransientSolution[1].Vz,
    145    md.results.TransientSolution[1].Vel,
    146    md.results.TransientSolution[1].Pressure,
    147    md.results.TransientSolution[1].Base,
    148    md.results.TransientSolution[1].Surface,
    149    md.results.TransientSolution[1].Thickness,
    150    md.results.TransientSolution[1].Temperature,
    151    md.results.TransientSolution[1].Enthalpy,
    152    md.results.TransientSolution[1].SmbMassBalance,
    153    md.results.TransientSolution[2].Vx,
    154    md.results.TransientSolution[2].Vy,
    155    md.results.TransientSolution[2].Vz,
    156    md.results.TransientSolution[2].Vel,
    157    md.results.TransientSolution[2].Pressure,
    158    md.results.TransientSolution[2].Base,
    159    md.results.TransientSolution[2].Surface,
    160    md.results.TransientSolution[2].Thickness,
    161    md.results.TransientSolution[2].Temperature,
    162    md.results.TransientSolution[2].Enthalpy,
    163    md.results.TransientSolution[2].SmbMassBalance,
    164    md.results.TransientSolution[3].Vx,
    165    md.results.TransientSolution[3].Vy,
    166    md.results.TransientSolution[3].Vz,
    167    md.results.TransientSolution[3].Vel,
    168    md.results.TransientSolution[3].Pressure,
    169    md.results.TransientSolution[3].Base,
    170    md.results.TransientSolution[3].Surface,
    171    md.results.TransientSolution[3].Thickness,
    172    md.results.TransientSolution[3].Temperature,
    173    md.results.TransientSolution[3].Enthalpy,
    174    md.results.TransientSolution[3].SmbMassBalance]
     132                md.results.TransientSolution[0].Vy,
     133                md.results.TransientSolution[0].Vz,
     134                md.results.TransientSolution[0].Vel,
     135                md.results.TransientSolution[0].Pressure,
     136                md.results.TransientSolution[0].Base,
     137                md.results.TransientSolution[0].Surface,
     138                md.results.TransientSolution[0].Thickness,
     139                md.results.TransientSolution[0].Temperature,
     140                md.results.TransientSolution[0].Enthalpy,
     141                md.results.TransientSolution[0].SmbMassBalance,
     142                md.results.TransientSolution[1].Vx,
     143                md.results.TransientSolution[1].Vy,
     144                md.results.TransientSolution[1].Vz,
     145                md.results.TransientSolution[1].Vel,
     146                md.results.TransientSolution[1].Pressure,
     147                md.results.TransientSolution[1].Base,
     148                md.results.TransientSolution[1].Surface,
     149                md.results.TransientSolution[1].Thickness,
     150                md.results.TransientSolution[1].Temperature,
     151                md.results.TransientSolution[1].Enthalpy,
     152                md.results.TransientSolution[1].SmbMassBalance,
     153                md.results.TransientSolution[2].Vx,
     154                md.results.TransientSolution[2].Vy,
     155                md.results.TransientSolution[2].Vz,
     156                md.results.TransientSolution[2].Vel,
     157                md.results.TransientSolution[2].Pressure,
     158                md.results.TransientSolution[2].Base,
     159                md.results.TransientSolution[2].Surface,
     160                md.results.TransientSolution[2].Thickness,
     161                md.results.TransientSolution[2].Temperature,
     162                md.results.TransientSolution[2].Enthalpy,
     163                md.results.TransientSolution[2].SmbMassBalance,
     164                md.results.TransientSolution[3].Vx,
     165                md.results.TransientSolution[3].Vy,
     166                md.results.TransientSolution[3].Vz,
     167                md.results.TransientSolution[3].Vel,
     168                md.results.TransientSolution[3].Pressure,
     169                md.results.TransientSolution[3].Base,
     170                md.results.TransientSolution[3].Surface,
     171                md.results.TransientSolution[3].Thickness,
     172                md.results.TransientSolution[3].Temperature,
     173                md.results.TransientSolution[3].Enthalpy,
     174                md.results.TransientSolution[3].SmbMassBalance]
  • issm/trunk-jpl/test/NightlyRun/test317.py

    r23793 r26906  
    2020
    2121#Fields and tolerances to track changes
    22 field_names = ['Vx1', 'Vy1', 'Vz1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'Temperature1', 'BasalforcingsGroundediceMeltingRate1', 'GroundedArea1', 'FloatingArea1', 'TotalFloatingBmb1', 'TotalGroundedBmb1'
     22field_names = ['Vx1', 'Vy1', 'Vz1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'Temperature1', 'BasalforcingsGroundediceMeltingRate1', 'GroundedArea1', 'FloatingArea1', 'TotalFloatingBmb1', 'TotalGroundedBmb1',
    2323               'Vx2', 'Vy2', 'Vz2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'Temperature2', 'BasalforcingsGroundediceMeltingRate2', 'GroundedArea2', 'FloatingArea2', 'TotalFloatingBmb2', 'TotalGroundedBmb2',
    2424               'Vx3', 'Vy3', 'Vz3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'Temperature3', 'BasalforcingsGroundediceMeltingRate3', 'GroundedArea3', 'FloatingArea3', 'TotalFloatingBmb2', 'TotalGroundedBmb2']
  • issm/trunk-jpl/test/NightlyRun/test427.py

    r24862 r26906  
    22from model import *
    33from socket import gethostname
    4 
    5 from setmask import *
    6 from parameterize import *
    7 from setflowequation import *
    8 from solve import *
     4from triangle import triangle
     5from setmask import setmask
     6from parameterize import parameterize
     7from setflowequation import setflowequation
     8from solve import solve
    99
    1010md = triangle(model(), '../Exp/Square.exp', 150000.)
  • issm/trunk-jpl/test/NightlyRun/test436.py

    r24214 r26906  
    33from model import *
    44from socket import gethostname
    5 from triangle import *
    6 from setmask import *
    7 from parameterize import *
    8 from setflowequation import *
    9 from solve import *
     5from triangle import triangle
     6from setmask import setmask
     7from parameterize import parameterize
     8from setflowequation import setflowequation
     9from solve import solve
    1010
    1111md = triangle(model(), '../Exp/Square.exp', 150000.)
     
    2626for i in ['LliboutryDuval', 'CuffeyTemperate']:
    2727    print(' ')
    28     print(' == == == Testing rheology law: ' + i + ' = ')
     28    print(' ====== Testing rheology law: ' + i + ' = ')
    2929
    3030    md.materials.rheology_law = i
  • issm/trunk-jpl/test/NightlyRun/test443.py

    r24862 r26906  
    33from model import *
    44from socket import gethostname
    5 
    6 from triangle import *
    7 from setmask import *
    8 from parameterize import *
    9 from setflowequation import *
    10 from solve import *
     5from triangle import triangle
     6from setmask import setmask
     7from parameterize import parameterize
     8from setflowequation import setflowequation
     9from solve import solve
    1110from frictioncoulomb import frictioncoulomb
    1211from ContourToMesh import ContourToMesh
  • issm/trunk-jpl/test/NightlyRun/test612.py

    r25023 r26906  
    22import numpy as np
    33from model import *
    4 from setmask import *
     4from setmask import setmask
     5from setflowequation import setflowequation
    56from triangle import triangle
    67from parameterize import parameterize
  • issm/trunk-jpl/test/NightlyRun/test701.py

    r24862 r26906  
    6666    md.flowequation.fe_FS = i
    6767    md = solve(md, 'Stressbalance')
    68     field_names = field_names + [['Vx' + i], ['Vy' + i], ['Vel' + i], ['Pressure' + i]]
     68    field_names = field_names + ['Vx' + i, 'Vy' + i, 'Vel' + i, 'Pressure' + i]
    6969    field_tolerances = field_tolerances + [9e-5, 9e-5, 9e-5, 1e-10]
    7070    field_values = field_values + [md.results.StressbalanceSolution.Vx,
Note: See TracChangeset for help on using the changeset viewer.