Changeset 22331


Ignore:
Timestamp:
01/04/18 16:06:16 (7 years ago)
Author:
seroussi
Message:

NEW: added test for new scaled routines

Location:
issm/trunk-jpl/test
Files:
6 edited

Legend:

Unmodified
Added
Removed

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

    r21056 r22331  
    55md=setflowequation(md,'SSA','all');
    66md.cluster=generic('name',oshostname(),'np',3);
     7md.geometry.bed=md.geometry.base;
     8[md.mesh.lat,md.mesh.long] = xy2ll(md.mesh.x,md.mesh.y,-1);
     9md.mesh.scale_factor=0.9*ones(md.mesh.numberofvertices,1);
     10md.transient.requested_outputs={'default','IceVolume','IceVolumeScaled','GroundedArea','GroundedAreaScaled','FloatingArea','FloatingAreaScaled','TotalSmb','TotalSmbScaled','TotalFloatingBmb','TotalFloatingBmbScaled'};
    711md=solve(md,'Transient');
    812
    913%Fields and tolerances to track changes
    10 field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2'};
    11 field_tolerances={1e-12,2e-12,2e-12,1e-13,1e-13,1e-13,1e-13,1e-12,1e-12,1e-12,1e-13,1e-13,1e-13,1e-13};
     14field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','IceVolume1','IceVolumeScaled1','GroundedArea1','GroundedAreaScaled1','FloatingArea1','FloatingAreaScaled1','TotalSmb1','TotalSmbScaled1','TotalFloatingBmb1','TotalFloatingBmbScaled1',...
     15        'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','IceVolume2','IceVolumeScaled2','GroundedArea2','GroundedAreaScaled2','FloatingArea2','FloatingAreaScaled2','TotalSmb2','TotalSmbScaled2','TotalFloatingBmb2','TotalFloatingBmbScaled2'};
     16field_tolerances={1e-12,2e-12,2e-12,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,...
     17        1e-12,1e-12,1e-12,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};
    1218field_values={...
    1319        (md.results.TransientSolution(1).Vx),...
    … …  
    1824        (md.results.TransientSolution(1).Surface),...
    1925        (md.results.TransientSolution(1).Thickness),...
     26        (md.results.TransientSolution(1).IceVolume),...
     27        (md.results.TransientSolution(1).IceVolumeScaled),...
     28        (md.results.TransientSolution(1).GroundedArea),...
     29        (md.results.TransientSolution(1).GroundedAreaScaled),...
     30        (md.results.TransientSolution(1).FloatingArea),...
     31        (md.results.TransientSolution(1).FloatingAreaScaled),...
     32        (md.results.TransientSolution(1).TotalSmb),...
     33        (md.results.TransientSolution(1).TotalSmbScaled),...
     34        (md.results.TransientSolution(1).TotalFloatingBmb),...
     35        (md.results.TransientSolution(1).TotalFloatingBmbScaled),...
    2036        (md.results.TransientSolution(2).Vx),...
    2137        (md.results.TransientSolution(2).Vy),...
    … …  
    2541        (md.results.TransientSolution(2).Surface),...
    2642        (md.results.TransientSolution(2).Thickness),...
     43        (md.results.TransientSolution(2).IceVolume),...
     44        (md.results.TransientSolution(2).IceVolumeScaled),...
     45        (md.results.TransientSolution(2).GroundedArea),...
     46        (md.results.TransientSolution(2).GroundedAreaScaled),...
     47        (md.results.TransientSolution(2).FloatingArea),...
     48        (md.results.TransientSolution(2).FloatingAreaScaled),...
     49        (md.results.TransientSolution(2).TotalSmb),...
     50        (md.results.TransientSolution(2).TotalSmbScaled),...
     51        (md.results.TransientSolution(2).TotalFloatingBmb),...
     52        (md.results.TransientSolution(2).TotalFloatingBmbScaled),...
    2753        };

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

    r21408 r22331  
    1313md=parameterize(md,'../Par/Pig.py')
    1414md=setflowequation(md,'SSA','all')
     15md.mesh.scale_factor=0.9*np.ones((md.mesh.numberofvertices))
     16md.transient.requested_outputs=['default','IceVolume','IceVolumeScaled','GroundedArea','GroundedAreaScaled','FloatingArea','FloatingAreaScaled','TotalSmb','TotalSmbScaled','TotalFloaringBmb','TotalFloatingBmbScaled']
    1517md.cluster=generic('name',gethostname(),'np',3)
    1618md=solve(md,'Transient')
    1719
    1820# Fields and tolerances to track changes
    19 field_names     =['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2']
    20 field_tolerances=[1e-12,2e-12,2e-12,1e-13,1e-13,1e-13,1e-13,1e-12,1e-12,1e-12,1e-13,1e-13,1e-13,1e-13]
     21field_names     =['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','IceVolume1','IceVolumeScaled1','GroundedArea1','GroundedAreaScaled1','FloatingArea1','FloatingAreaScaled1','TotalSmb1','TotalSmbScaled1','TotalFloatingBmb1','TotalFloatingBmbScaled1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','IceVolume2','IceVolumeScaled2','GroundedArea2','GroundedAreaScaled2','FloatingArea2','FloatingAreaScaled2','TotalSmb2','TotalSmbScaled2','TotalFloatingBmb2','TotalFloatingBmbScaled2']
     22field_tolerances=[1e-12,2e-12,2e-12,1e-13,1e-13,1e-13,1e-13,1e-12,1e-12,1e-12,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,1e-13,1e-13],
    2123field_values=[\
    2224        md.results.TransientSolution[0].Vx,\
    … …  
    2729        md.results.TransientSolution[0].Surface,\
    2830        md.results.TransientSolution[0].Thickness,\
     31        md.results.TransientSolution[0].IceVolume,\
     32        md.results.TransientSolution[0].IceVolumeScaled,\
     33        md.results.TransientSolution[0].GroundedArea,\
     34        md.results.TransientSolution[0].GroundedAreaScaled,\
     35        md.results.TransientSolution[0].FloatingArea,\
     36        md.results.TransientSolution[0].FloatingAreaScaled,\
     37        md.results.TransientSolution[0].TotalSmb,\
     38        md.results.TransientSolution[0].TotalSmbScaled,\
     39        md.results.TransientSolution[0].TotalFloatingBmb,\
     40        md.results.TransientSolution[0].TotalFloatingBmbScaled,\
    2941        md.results.TransientSolution[1].Vx,\
    3042        md.results.TransientSolution[1].Vy,\
    … …  
    3446        md.results.TransientSolution[1].Surface,\
    3547        md.results.TransientSolution[1].Thickness,\
     48        md.results.TransientSolution[1].IceVolume,\
     49        md.results.TransientSolution[1].IceVolumeScaled,\
     50        md.results.TransientSolution[1].GroundedArea,\
     51        md.results.TransientSolution[1].GroundedAreaScaled,\
     52        md.results.TransientSolution[1].FloatingArea,\
     53        md.results.TransientSolution[1].FloatingAreaScaled,\
     54        md.results.TransientSolution[1].TotalSmb,\
     55        md.results.TransientSolution[1].TotalSmbScaled,\
     56        md.results.TransientSolution[1].TotalFloatingBmb,\
     57        md.results.TransientSolution[1].TotalFloatingBmbScaled,\
    3658        ]

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

    r21056 r22331  
    33md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
    44md=parameterize(md,'../Par/Pig.par');
     5md.mesh.scale_factor=0.9*ones(md.mesh.numberofvertices,1);
    56md=extrude(md,2,1.);
    67md=setflowequation(md,'HO','all');
    78md.cluster=generic('name',oshostname(),'np',3);
     9md.transient.requested_outputs={'default','IceVolume','IceVolumeScaled','TotalSmb','TotalSmbScaled','TotalFloatingBmb','TotalFloatingBmbScaled','FloatingArea','FloatingAreaScaled','GroundedArea','GroundedAreaScaled'};
    810md=solve(md,'Transient');
    911
    1012%Fields and tolerances to track changes
    11 field_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1', ...
    12                                       'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2'};
    13 field_tolerances={1e-10,1e-10,1e-10,1e-10,1e-12,1e-11,2e-12,1e-11,1e-12,1e-09,...
    14                                                 1e-11,1e-11,1e-09,1e-11,1e-11,1e-10,1e-11,1e-10,1e-11,2e-08};
     13field_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1','IceVolume1','IceVolumeScaled1','GroundedArea1','GroundedAreaScaled1','FloatingArea1','FloatingAreaScaled1','TotalSmb1','TotalSmbScaled1','TotalFloatingBmb1','TotalFloatingBmbScaled1', ...
     14                                      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2','IceVolume2','IceVolumeScaled2','GroundedArea2','GroundedAreaScaled2','FloatingArea2','FloatingAreaScaled2','TotalSmb2','TotalSmbScaled2','TotalFloatingBmb2','TotalFloatingBmbScaled2'};
     15field_tolerances={1e-10,1e-10,1e-10,1e-10,1e-12,1e-11,2e-12,1e-11,1e-12,1e-09,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
     16                                                1e-11,1e-11,1e-09,1e-11,1e-11,1e-10,1e-11,1e-10,1e-11,2e-08,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
    1517field_values={...
    1618        (md.results.TransientSolution(1).Vx),...
    … …  
    2426        (md.results.TransientSolution(1).Temperature),...
    2527        (md.results.TransientSolution(1).BasalforcingsGroundediceMeltingRate),...
     28        (md.results.TransientSolution(1).IceVolume),...
     29        (md.results.TransientSolution(1).IceVolumeScaled),...
     30        (md.results.TransientSolution(1).GroundedArea),...
     31        (md.results.TransientSolution(1).GroundedAreaScaled),...
     32        (md.results.TransientSolution(1).FloatingArea),...
     33        (md.results.TransientSolution(1).FloatingAreaScaled),...
     34        (md.results.TransientSolution(1).TotalSmb),...
     35        (md.results.TransientSolution(1).TotalSmbScaled),...
     36        (md.results.TransientSolution(1).TotalFloatingBmb),...
     37        (md.results.TransientSolution(1).TotalFloatingBmbScaled),...
    2638        (md.results.TransientSolution(2).Vx),...
    2739        (md.results.TransientSolution(2).Vy),...
    … …  
    3446        (md.results.TransientSolution(2).Temperature),...
    3547        (md.results.TransientSolution(2).BasalforcingsGroundediceMeltingRate),...
     48        (md.results.TransientSolution(2).IceVolume),...
     49        (md.results.TransientSolution(2).IceVolumeScaled),...
     50        (md.results.TransientSolution(2).GroundedArea),...
     51        (md.results.TransientSolution(2).GroundedAreaScaled),...
     52        (md.results.TransientSolution(2).FloatingArea),...
     53        (md.results.TransientSolution(2).FloatingAreaScaled),...
     54        (md.results.TransientSolution(2).TotalSmb),...
     55        (md.results.TransientSolution(2).TotalSmbScaled),...
     56        (md.results.TransientSolution(2).TotalFloatingBmb),...
     57        (md.results.TransientSolution(2).TotalFloatingBmbScaled),...
    3658        };

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

    r21408 r22331  
    1212md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp')
    1313md=parameterize(md,'../Par/Pig.py')
     14md.mesh.scale_factor=0.9*np.ones((md.mesh.numberofvertices))
    1415md.extrude(2,1.)
    1516md=setflowequation(md,'HO','all')
     17md.transient.requested_outputs=['default','IceVolume','IceVolumeScaled','GroundedArea','GroundedAreaScaled','FloatingArea','FloatingAreaScaled','TotalSmb','TotalSmbScaled','TotalFloaringBmb','TotalFloatingBmbScaled']
    1618md.cluster=generic('name',gethostname(),'np',3)
    1719md=solve(md,'Transient')
    1820
    1921# Fields and tolerances to track changes
    20 field_names     =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1', \
    21                                       'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2']
    22 field_tolerances=[1e-10,1e-10,1e-10,1e-10,1e-12,1e-11,2e-12,1e-11,1e-12,1e-09,\
    23                                                 1e-11,1e-11,1e-09,1e-11,1e-11,1e-10,1e-11,1e-10,1e-11,2e-08]
     22field_names     =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1','IceVolume1','IceVolumeScaled1','GroundedArea1','GroundedAreaScaled1','FloatingArea1','FloatingAreaScaled1','TotalSmb1','TotalSmbScaled1','TotalFloatingBmb1','TotalFloatingBmbScaled1', \
     23                                      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2','IceVolume2','IceVolumeScaled2','GroundedArea2','GroundedAreaScaled2','FloatingArea2','FloatingAreaScaled2','TotalSmb2','TotalSmbScaled2','TotalFloatingBmb2','TotalFloatingBmbScaled2']
     24field_tolerances=[1e-10,1e-10,1e-10,1e-10,1e-12,1e-11,2e-12,1e-11,1e-12,1e-09,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
     25                                                1e-11,1e-11,1e-09,1e-11,1e-11,1e-10,1e-11,1e-10,1e-11,2e-08,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
    2426field_values=[\
    2527        md.results.TransientSolution[0].Vx,\
    … …  
    3335        md.results.TransientSolution[0].Temperature,\
    3436        md.results.TransientSolution[0].BasalforcingsGroundediceMeltingRate,\
     37        md.results.TransientSolution[0].IceVolume,\
     38        md.results.TransientSolution[0].IceVolumeScaled,\
     39        md.results.TransientSolution[0].GroundedArea,\
     40        md.results.TransientSolution[0].GroundedAreaScaled,\
     41        md.results.TransientSolution[0].FloatingArea,\
     42        md.results.TransientSolution[0].FloatingAreaScaled,\
     43        md.results.TransientSolution[0].TotalSmb,\
     44        md.results.TransientSolution[0].TotalSmbScaled,\
     45        md.results.TransientSolution[0].TotalFloatingBmb,\
     46        md.results.TransientSolution[0].TotalFloatingBmbScaled,\
    3547        md.results.TransientSolution[1].Vx,\
    3648        md.results.TransientSolution[1].Vy,\
    … …  
    4355        md.results.TransientSolution[1].Temperature,\
    4456        md.results.TransientSolution[1].BasalforcingsGroundediceMeltingRate,\
     57        md.results.TransientSolution[1].IceVolume,\
     58        md.results.TransientSolution[1].IceVolumeScaled,\
     59        md.results.TransientSolution[1].GroundedArea,\
     60        md.results.TransientSolution[1].GroundedAreaScaled,\
     61        md.results.TransientSolution[1].FloatingArea,\
     62        md.results.TransientSolution[1].FloatingAreaScaled,\
     63        md.results.TransientSolution[1].TotalSmb,\
     64        md.results.TransientSolution[1].TotalSmbScaled,\
     65        md.results.TransientSolution[1].TotalFloatingBmb,\
     66        md.results.TransientSolution[1].TotalFloatingBmbScaled,\
    4567        ]
Note: See TracChangeset for help on using the changeset viewer.