Changeset 21827


Ignore:
Timestamp:
07/19/17 17:40:37 (8 years ago)
Author:
schlegel
Message:

CHG: add a python test for 443, regional output

Location:
issm/trunk-jpl
Files:
1 added
6 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/classes/model.m

    r21819 r21827  
    177177                        %Start with changing all the fields from the 3d mesh
    178178
    179                         %dealing with the friciton law
     179                        %dealing with the friction law
    180180                        %drag is limited to nodes that are on the bedrock.
    181181                        if isa(md.friction,'friction'),
    … …  
    206206                        else
    207207                                disp('friction type not supported');
    208             end
     208                        end
    209209
    210210                        %observations

  • issm/trunk-jpl/src/m/classes/model.py

    r21808 r21827  
    689689                        raise StandardError("only a 3D model can be collapsed")
    690690               
     691                #dealing with the friction law
    691692                #drag is limited to nodes that are on the bedrock.
    692                 md.friction.coefficient=project2d(md,md.friction.coefficient,1)
     693                if hasattr(md.friction,'coefficient'): md.friction.coefficient=project2d(md,md.friction.coefficient,1)
    693694
    694695                #p and q (same deal, except for element that are on the bedrock: )
    695                 md.friction.p=project2d(md,md.friction.p,1)
    696                 md.friction.q=project2d(md,md.friction.q,1)
     696                if hasattr(md.friction,'p'): md.friction.p=project2d(md,md.friction.p,1)
     697                if hasattr(md.friction,'q'): md.friction.q=project2d(md,md.friction.q,1)
     698               
     699                if hasattr(md.friction,'coefficientcoulomb'): md.friction.q=project2d(md,md.friction.coefficientcoulomb,1)
     700                if hasattr(md.friction,'C'): md.friction.q=project2d(md,md.friction.C,1)
     701                if hasattr(md.friction,'As'): md.friction.q=project2d(md,md.friction.As,1)
     702                if hasattr(md.friction,'effective_pressure'): md.friction.q=project2d(md,md.friction.effective_pressure,1)
     703                if hasattr(md.friction,'water_layer'): md.friction.q=project2d(md,md.friction.water_layer,1)
     704                if hasattr(md.friction,'m'): md.friction.q=project2d(md,md.friction.m,1)
    697705
    698706                #observations
    … …  
    731739                        md.flowequation.borderFS=project2d(md,md.flowequation.borderFS,1)
    732740
    733 
    734       # Hydrologydc variables
     741                # Hydrologydc variables
    735742                if hasattr(md.hydrology,'hydrologydc'):
    736743                        md.hydrology.spcsediment_head=project2d(md,md.hydrology.spcsediment_head,1)
    … …  
    794801                                                                elif np.size(subfield)==md.mesh.numberofelements:
    795802                                                                        setattr(fieldr,solutionsubfield,project2d(md,subfield,1))
    796                                                                 else:
    797                                                                         setattr(fieldr,solutionsubfield,subfield)
     803
    798804
    799805                #Initialize with the 2d mesh
    … …  
    812818                mesh.epsg=md.mesh.epsg
    813819                md.mesh=mesh
     820                md.mesh.vertexconnectivity=NodeConnectivity(md.mesh.elements,md.mesh.numberofvertices)[0]
     821                md.mesh.elementconnectivity=ElementConnectivity(md.mesh.elements,md.mesh.vertexconnectivity)[0]
     822                md.mesh.segments=contourenvelope(md)
    814823
    815824                return md

  • issm/trunk-jpl/src/m/classes/regionaloutput.py

    r21808 r21827  
     1from project3d import project3d
    12from fielddisplay import fielddisplay
    23from pairoptions import pairoptions
    … …  
    45from WriteData import WriteData
    56from MeshProfileIntersection import MeshProfileIntersection
     7from ContourToMesh import ContourToMesh
     8import numpy as np
    69import os
    710
    … …  
    3740
    3841                #get name
    39                 modelname=getfieldvalue(options,'model');
    40                 self.maskexpstring=getfieldvalue(options,'maskexpstring');
    41                 self.setmaskfromexp(self,md)
    42                 error('regionaloutput error message: ''mask'' field or ''maskexpstring'' and ''model'' fields should be defined!');
     42                if options.getfieldvalue('model',0):
     43                        if options.getfieldvalue('maskexpstring',0):
     44                                modelname=options.getfieldvalue('model')
     45                                self.maskexpstring=options.getfieldvalue('maskexpstring')
     46                                self.setmaskfromexp(modelname)
     47                       
     48                if (len(self.mask)<=1 & np.any(np.isnan(self.mask))):
     49                        error('regionaloutput error message: ''mask'' field or ''maskexpstring'' and ''model'' fields should be defined!');
    4350
    4451                #}}}
    … …  
    6168        #}}}
    6269        def setmaskfromexp(self,md):    # {{{
    63                 if len(self.maskexpstring,0) > 0:
     70                if len(self.maskexpstring) > 0:
    6471                        self.mask=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,self.maskexpstring,'node',1)
    6572                       
    … …  
    7481                        raise RuntimeError("regionaloutput error message: 'outputnamestring' field should be a string!")
    7582               
    76                 if len(self.maskexpstring,0) > 0:
    77                         if not os.path.isfile(self.profilename):
     83                if len(self.maskexpstring) > 0:
     84                        if not os.path.isfile(self.maskexpstring):
    7885                                raise RuntimeError("regionaloutput error message: file name for mask exp does not point to a legitimate file on disk!")
    7986                        else:
    80                                 self.setmaskfromexp(self,md)
     87                                self.setmaskfromexp(md)
    8188
    8289                OutputdefinitionStringArray=[]
    … …  
    8592                        OutputdefinitionStringArray.append(x)
    8693
    87                 md = checkfield(md,'field',self.definitionstring,'values',OutputdefinitionStringArray])
    88       md = checkfield(md,'field',self.mask,'size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1)
    89                
     94                md = checkfield(md,'field',self.definitionstring,'values',OutputdefinitionStringArray)
     95                md = checkfield(md,'field',self.mask,'size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    9096                return md
    9197        # }}}
    … …  
    9399
    94100                #before marshalling, make sure mask is set:
    95                 self.setmaskfromexp(self,md)
     101                self.setmaskfromexp(md)
    96102
    97103                #ok, marshall strings and mask:

  • issm/trunk-jpl/test/NightlyRun/test408.py

    r21613 r21827  
    2424md.transient.isthermal=False
    2525md.transient.isgroundingline=True
    26 md.extrude(3,1.)
    2726md=setflowequation(md,'SSA','all')
    2827md.cluster=generic('name',gethostname(),'np',3)
    2928md.transient.requested_outputs=['default','GroundedArea','FloatingArea','TotalFloatingBmb','TotalGroundedBmb','TotalSmb']
     29md.extrude(3,1.)
    3030md=solve(md,'Transient')
    3131

  • issm/trunk-jpl/test/NightlyRun/test443.m

    r21820 r21827  
    4444
    4545md=extrude(md,3,1);
     46md2=collapse(md);
    4647md=solve(md,'Transient');
    4748
    … …  
    5051        'Thickness1','GroundedArea1','FloatingArea1','TotalFloatingBmb1','TotalGroundedBmb1','TotalSmb1',...
    5152        'IceMass3','IceVolume3','IceVolumeAboveFloatation3','IceVolumeAboveFloatation23',...
    52         'Thickness3','GroundedArea3','FloatingArea3','TotalFloatingBmb3','TotalGroundedBmb3','TotalSmb3'};
     53        'Thickness3','GroundedArea3','FloatingArea3','TotalFloatingBmb3','TotalGroundedBmb3','TotalSmb3',...
     54        'ExtrudedMask','CollapsedMask'};
    5355field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
    54         1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     56        1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
    5557field_values={...
    5658        (md.results.TransientSolution(1).IceMass1),...
    … …  
    7476        (md.results.TransientSolution(3).TotalGroundedBmb1),...
    7577        (md.results.TransientSolution(3).TotalSmb1),...
     78        length(md.outputdefinition.definitions{1}.mask),...
     79        length(md2.outputdefinition.definitions{1}.mask),...
    7680        };
Note: See TracChangeset for help on using the changeset viewer.