source: issm/oecreview/Archive/13393-13976/ISSM-13661-13662.diff

Last change on this file was 13980, checked in by Mathieu Morlighem, 12 years ago

preparing oecreview for 13393-13976'
File size: 9.9 KB

  • ../trunk-jpl/test/NightlyRun/test512.py

     
     1import numpy
     2from model import *
     3from triangle import *
     4from setmask import *
     5from parameterize import *
     6from setflowequation import *
     7from EnumDefinitions import *
     8from solve import *
     9from MatlabFuncs import *
     10
     11md=triangle(model(),'../Exp/Pig.exp',20000.)
     12md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp')
     13md=parameterize(md,'../Par/Pig.py')
     14md.extrude(3,1.)
     15md=setflowequation(md,'pattyn','all')
     16
     17# control parameters
     18md.inversion.iscontrol=1
     19md.inversion.control_parameters=['FrictionCoefficient']
     20md.inversion.min_parameters=1.*numpy.ones((md.mesh.numberofvertices,1))
     21md.inversion.max_parameters=200.*numpy.ones((md.mesh.numberofvertices,1))
     22md.inversion.nsteps=2
     23md.inversion.cost_functions=numpy.hstack((103*numpy.ones((md.inversion.nsteps,1)),501*numpy.ones((md.inversion.nsteps,1))))
     24md.inversion.cost_functions_coefficients=numpy.ones((md.mesh.numberofvertices,2))
     25md.inversion.cost_functions_coefficients[:,1]=2.*10**-7
     26md.inversion.gradient_scaling=3.*numpy.ones((md.inversion.nsteps,1))
     27md.inversion.maxiter_per_step=2.*numpy.ones((md.inversion.nsteps,1))
     28md.inversion.step_threshold=2.99*numpy.ones((md.inversion.nsteps,1))
     29md.inversion.vx_obs=md.initialization.vx
     30md.inversion.vy_obs=md.initialization.vy
     31
     32md.cluster=generic('name',oshostname(),'np',3)
     33md=solve(md,DiagnosticSolutionEnum())
     34
     35# Fields and tolerances to track changes
     36field_names     =['Gradient','Misfits','FrictionCoefficient','Pressure','Vel','Vx','Vy']
     37field_tolerances=[1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11]
     38field_values=[\
     39        md.results['DiagnosticSolution'][1]['Gradient1'],\
     40        md.results['DiagnosticSolution'][1]['J'],\
     41        md.results['DiagnosticSolution'][1]['FrictionCoefficient'],\
     42        md.results['DiagnosticSolution'][1]['Pressure'],\
     43        md.results['DiagnosticSolution'][1]['Vel'],\
     44        md.results['DiagnosticSolution'][1]['Vx'],\
     45        md.results['DiagnosticSolution'][1]['Vy']
     46]

  • ../trunk-jpl/test/NightlyRun/test510.m

     
    99
    1010%Fields and tolerances to track changes
    1111field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalforcingsMeltingRate'};
    12 field_tolerances={1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-06
    13 };
     12field_tolerances={1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-06};
    1413field_values={...
    1514        (md.results.SteadystateSolution.Vx),...
    1615        (md.results.SteadystateSolution.Vy),...

  • ../trunk-jpl/test/NightlyRun/test513.py

     
     1from model import *
     2from triangle import *
     3from setmask import *
     4from parameterize import *
     5from setflowequation import *
     6from EnumDefinitions import *
     7from solve import *
     8from MatlabFuncs import *
     9
     10md=triangle(model(),'../Exp/Pig.exp',30000.)
     11md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp')
     12md=parameterize(md,'../Par/Pig.py')
     13md.extrude(3,1)
     14md=setflowequation(md,'macayeal','all')
     15
     16# control parameters
     17md.inversion.iscontrol=1
     18md.inversion.control_parameters=['FrictionCoefficient']
     19md.inversion.min_parameters=1.*numpy.ones((md.mesh.numberofvertices,1))
     20md.inversion.max_parameters=200.*numpy.ones((md.mesh.numberofvertices,1))
     21md.inversion.nsteps=2
     22md.inversion.cost_functions=numpy.hstack((103*numpy.ones((md.inversion.nsteps,1)),501*numpy.ones((md.inversion.nsteps,1))))
     23md.inversion.cost_functions_coefficients=numpy.ones((md.mesh.numberofvertices,2))
     24md.inversion.cost_functions_coefficients[:,1]=2.*10**-7
     25md.inversion.gradient_scaling=3.*numpy.ones((md.inversion.nsteps,1))
     26md.inversion.maxiter_per_step=2.*numpy.ones((md.inversion.nsteps,1))
     27md.inversion.step_threshold=0.99*numpy.ones((md.inversion.nsteps,1))
     28md.inversion.vx_obs=md.initialization.vx
     29md.inversion.vy_obs=md.initialization.vy
     30md.timestepping.time_step=0
     31
     32md.thermal.penalty_lock=5
     33md.cluster=generic('name',oshostname(),'np',3)
     34md=solve(md,SteadystateSolutionEnum())
     35
     36# Fields and tolerances to track changes
     37field_names     =['Gradient','Misfits','FrictionCoefficient','Pressure','Vel','Vx','Vy','Vz','Temperature','BasalforcingsMeltingRate']
     38field_tolerances=[1e-10,1e-10,1e-10,1e-10,1e-6,1e-6,1e-6,1e-6,1e-6,1e-06]
     39field_values=[\
     40        md.results['SteadystateSolution'][1]['Gradient1'],\
     41        md.results['SteadystateSolution'][1]['J'],\
     42        md.results['SteadystateSolution'][1]['FrictionCoefficient'],\
     43        md.results['SteadystateSolution'][1]['Pressure'],\
     44        md.results['SteadystateSolution'][1]['Vel'],\
     45        md.results['SteadystateSolution'][1]['Vx'],\
     46        md.results['SteadystateSolution'][1]['Vy'],\
     47        md.results['SteadystateSolution'][1]['Vz'],\
     48        md.results['SteadystateSolution'][1]['Temperature'],\
     49        md.results['SteadystateSolution'][1]['BasalforcingsMeltingRate']
     50]

  • ../trunk-jpl/test/NightlyRun/test512.m

     
    2121md=solve(md,DiagnosticSolutionEnum);
    2222
    2323%Fields and tolerances to track changes
    24 field_names     ={'Gradient' 'Misfits' 'FrictionCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy'};
     24field_names     ={'Gradient','Misfits','FrictionCoefficient','Pressure','Vel','Vx','Vy'};
    2525field_tolerances={1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11};
    2626field_values={...
    2727        (md.results.DiagnosticSolution.Gradient1),...

  • ../trunk-jpl/test/NightlyRun/test513.m

     
    2323md=solve(md,SteadystateSolutionEnum);
    2424
    2525%Fields and tolerances to track changes
    26 field_names     ={'Gradient' 'Misfits' 'FrictionCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'BasalforcingsMeltingRate'};
     26field_names     ={'Gradient','Misfits','FrictionCoefficient','Pressure','Vel','Vx','Vy','Vz','Temperature','BasalforcingsMeltingRate'};
    2727field_tolerances={1e-10,1e-10,1e-10,1e-10,1e-6,1e-6,1e-6,1e-6,1e-6,1e-06};
    2828field_values={...
    2929        (md.results.SteadystateSolution.Gradient1),...

  • ../trunk-jpl/test/NightlyRun/test515.py

     
     1from model import *
     2from triangle import *
     3from setmask import *
     4from parameterize import *
     5from setflowequation import *
     6from EnumDefinitions import *
     7from solve import *
     8from MatlabFuncs import *
     9
     10md=triangle(model(),'../Exp/Pig.exp',30000)
     11md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp')
     12md=parameterize(md,'../Par/Pig.py')
     13md.extrude(3,1)
     14md=setflowequation(md,'pattyn','all')
     15md.thermal.stabilization=2
     16md.cluster=generic('name',oshostname(),'np',3)
     17md.transient.isdiagnostic=0
     18md.transient.isprognostic=0
     19md.transient.isthermal=1
     20md.transient.isgroundingline=0
     21md=solve(md,TransientSolutionEnum())
     22
     23# Fields and tolerances to track changes
     24field_names     =['Temperature1','BasalforcingsMeltingRate1', \
     25                                      'Temperature2','BasalforcingsMeltingRate2']
     26field_tolerances=[1e-13,1e-8,1e-13,1e-8]
     27field_values=[\
     28        md.results['TransientSolution'][1]['Temperature'],\
     29        md.results['TransientSolution'][1]['BasalforcingsMeltingRate'],\
     30        md.results['TransientSolution'][2]['Temperature'],\
     31        md.results['TransientSolution'][2]['BasalforcingsMeltingRate'],\
     32        ]

  • ../trunk-jpl/test/NightlyRun/test510.py

     
     1from model import *
     2from triangle import *
     3from setmask import *
     4from parameterize import *
     5from setflowequation import *
     6from EnumDefinitions import *
     7from solve import *
     8from MatlabFuncs import *
     9
     10md=triangle(model(),'../Exp/Pig.exp',20000)
     11md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp')
     12md=parameterize(md,'../Par/Pig.py')
     13md.extrude(2,1)
     14md=setflowequation(md,'stokes','all')
     15md.cluster=generic('name',oshostname(),'np',3)
     16md.timestepping.time_step=0
     17md=solve(md,SteadystateSolutionEnum())
     18
     19# Fields and tolerances to track changes
     20field_names     =['Vx','Vy','Vz','Vel','Pressure','Temperature','BasalforcingsMeltingRate']
     21field_tolerances=[1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-06]
     22field_values=[\
     23        md.results['SteadystateSolution'][1]['Vx'],\
     24        md.results['SteadystateSolution'][1]['Vy'],\
     25        md.results['SteadystateSolution'][1]['Vz'],\
     26        md.results['SteadystateSolution'][1]['Vel'],\
     27        md.results['SteadystateSolution'][1]['Pressure'],\
     28        md.results['SteadystateSolution'][1]['Temperature'],\
     29        md.results['SteadystateSolution'][1]['BasalforcingsMeltingRate'],\
     30        ]

  • ../trunk-jpl/test/NightlyRun/test514.m

     
    4343y7=md2.mesh.y;
    4444
    4545%Fields and tolerances to track changes
    46 field_names     ={'x1' 'y1' 'x2' 'y2' 'x3' 'y3' 'x4' 'y4' 'x5' 'y5' 'x6' 'y6'  'x7' 'y7' };
    47 field_tolerances={1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13};
     46field_names     ={'x1y1x2y2x3y3x4y4x5y5x6y6'  'x7y7' };
     47field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
    4848field_values={...
    4949        x1, y1,...
    5050        y2, y2,...
Note: See TracBrowser for help on using the repository browser.