source: issm/oecreview/Archive/13393-13976/ISSM-13517-13518.diff

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

preparing oecreview for 13393-13976'

File size: 5.4 KB
  • ../trunk-jpl/test/NightlyRun/test212.py

     
     1"""
     2== == == == == == == == == == == == == == == == == == ==
     3Auto generated python script for ISSM:   test212.m
     4Created on 2012-09-27 via translateToPy.py Ver 1.0 by mikep
     5== == == == == == == == == == == == == == == == == == ==
     6
     7Matlab script conversion into python
     8translateToPy.py Author: Michael Pellegrin
     9translateToPy.py Date: 09/24/12
     10== == == == == == == == == == == == == == == == == == ==
     11"""
     12
     13from MatlabFuncs import *
     14from model import *
     15from EnumDefinitions import *
     16from numpy import *
     17from triangle import *
     18from setmask import *
     19from parameterize import *
     20from setflowequation import *
     21from solve import *
     22
     23md=triangle(model(),'../Exp/Square.exp',200000)
     24md=setmask(md,'all','')
     25md=parameterize(md,'../Par/SquareShelf.py')
     26md=setflowequation(md,'macayeal','all')
     27
     28
     29# control parameters
     30
     31md.inversion.iscontrol=1
     32md.inversion.control_parameters=['MaterialsRheologyBbar']
     33md.inversion.min_parameters=10**6*ones((md.mesh.numberofvertices,1))
     34md.inversion.max_parameters=2*10**9*ones((md.mesh.numberofvertices,1))
     35md.inversion.nsteps=2
     36md.inversion.cost_functions=101*ones((md.inversion.nsteps,1))
     37md.inversion.cost_functions_coefficients=ones((md.mesh.numberofvertices,1))
     38md.inversion.gradient_scaling=10**7*ones((md.inversion.nsteps,1))
     39md.inversion.maxiter_per_step=2*ones((md.inversion.nsteps,1))
     40md.inversion.step_threshold=0.3*ones((md.inversion.nsteps,1))
     41md.inversion.vx_obs=md.initialization.vx
     42md.inversion.vy_obs=md.initialization.vy
     43
     44
     45md.cluster=generic('name',oshostname(),'np',3)
     46md=solve(md,DiagnosticSolutionEnum())
     47
     48
     49# Fields and tolerances to track changes
     50
     51field_names     =['Gradient','Misfits','MaterialsRheologyBbar','Pressure','Vel','Vx','Vy']
     52field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
     53field_values=[\
     54md.results['DiagnosticSolution'][1]['Gradient1'],\
     55md.results['DiagnosticSolution'][1]['J'],\
     56md.results['DiagnosticSolution'][1]['MaterialsRheologyBbar'],\
     57md.results['DiagnosticSolution'][1]['Pressure'],\
     58md.results['DiagnosticSolution'][1]['Vel'],\
     59md.results['DiagnosticSolution'][1]['Vx'],\
     60md.results['DiagnosticSolution'][1]['Vy'],\
     61]
  • ../trunk-jpl/test/NightlyRun/test212.m

     
    2424field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
    2525field_values={...
    2626        (md.results.DiagnosticSolution.Gradient1),...
    27         md.results.DiagnosticSolution.J,...
     27        (md.results.DiagnosticSolution.J),...
    2828        (md.results.DiagnosticSolution.MaterialsRheologyBbar),...
    2929        (md.results.DiagnosticSolution.Pressure),...
    3030        (md.results.DiagnosticSolution.Vel),...
  • ../trunk-jpl/test/NightlyRun/InNeedOfDebugging/test212.py

     
    1 """
    2 == == == == == == == == == == == == == == == == == == ==
    3 Auto generated python script for ISSM:   test212.m
    4 Created on 2012-09-27 via translateToPy.py Ver 1.0 by mikep
    5 == == == == == == == == == == == == == == == == == == ==
    6 
    7 Matlab script conversion into python
    8 translateToPy.py Author: Michael Pellegrin
    9 translateToPy.py Date: 09/24/12
    10 == == == == == == == == == == == == == == == == == == ==
    11 """
    12 
    13 from MatlabFuncs import *
    14 from model import *
    15 from EnumDefinitions import *
    16 from numpy import *
    17 from triangle import *
    18 from setmask import *
    19 from parameterize import *
    20 from setflowequation import *
    21 from solve import *
    22 
    23 md=triangle(model(),'../Exp/Square.exp',200000)
    24 md=setmask(md,'all','')
    25 md=parameterize(md,'../Par/SquareShelf.py')
    26 md=setflowequation(md,'macayeal','all')
    27 
    28 
    29 # control parameters
    30 
    31 md.inversion.iscontrol=1
    32 md.inversion.control_parameters=['MaterialsRheologyBbar']
    33 md.inversion.min_parameters=10^6*ones(md.mesh.numberofvertices,1)
    34 md.inversion.max_parameters=2*10^9*ones(md.mesh.numberofvertices,1)
    35 md.inversion.nsteps=2
    36 md.inversion.cost_functions=101*ones(md.inversion.nsteps,1)
    37 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1)
    38 md.inversion.gradient_scaling=10^7*ones(md.inversion.nsteps,1)
    39 md.inversion.maxiter_per_step=2*ones(md.inversion.nsteps,1)
    40 md.inversion.step_threshold=0.3*ones(md.inversion.nsteps,1)
    41 md.inversion.vx_obs=md.initialization.vx
    42 md.inversion.vy_obs=md.initialization.vy
    43 
    44 
    45 md.cluster=generic('name',oshostname(),'np',3)
    46 md=solve(md,DiagnosticSolutionEnum())
    47 
    48 
    49 # Fields and tolerances to track changes
    50 
    51 field_names     =['Gradient' 'Misfits' 'MaterialsRheologyBbar' 'Pressure' 'Vel' 'Vx' 'Vy']
    52 field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
    53 field_values=[\
    54 md.results['DiagnosticSolution'][1]['Gradient1'],\
    55 md.results['DiagnosticSolution.J'],\
    56 md.results['DiagnosticSolution'][1]['MaterialsRheologyBbar'],\
    57 md.results['DiagnosticSolution'][1]['Pressure'],\
    58 md.results['DiagnosticSolution'][1]['Vel'],\
    59 md.results['DiagnosticSolution'][1]['Vx'],\
    60 md.results['DiagnosticSolution'][1]['Vy'],\
    61 ]
Note: See TracBrowser for help on using the repository browser.