Changeset 22857


Ignore:
Timestamp:
06/21/18 06:00:43 (7 years ago)
Author:
bdef
Message:

CHG: modifying naming related to hydrology sub-step

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

Legend:

Unmodified
Added
Removed

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

    r21965 r22857  
    3838md=solve(md,'Transient');
    3939
    40 %fields to track, results can also be found in 
     40%fields to track, results can also be found in
    4141%Wang 2009 Fig 6b (jouranl of Hydrology)
    4242field_names={'SedimentWaterHead1',...
    … …  
    4444field_tolerances={1e-13,...
    4545                                                                        1e-13};
    46 field_values={md.results.TransientSolution(11).SedimentHead,...
    47                                                         md.results.TransientSolution(31).SedimentHead};
     46field_values={md.results.TransientSolution(11).SedimentHeadHydrostep,...
     47                                                        md.results.TransientSolution(31).SedimentHeadHydrostep};

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

    r21965 r22857  
    4949md=solve(md,'Transient')
    5050
    51 #fields to track, results can also be found in 
     51#fields to track, results can also be found in
    5252#Wang 2009 Fig 6b (jouranl of Hydrology)
    5353field_names=['SedimentWaterHead1',
    5454                                                 'SedimentWaterHead2']
    55 field_tolerances=[1e-13, 
     55field_tolerances=[1e-13,
    5656                                                                        1e-13]
    57 field_values=[md.results.TransientSolution[10].SedimentHead,
    58                                                         md.results.TransientSolution[30].SedimentHead]
     57field_values=[md.results.TransientSolution[10].SedimentHeadHydrostep,
     58                                                        md.results.TransientSolution[30].SedimentHeadHydrostep]

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

    r21411 r22857  
    6666for i in range(0,4000):
    6767        sedvol[i]=np.mean(md.results.TransientSolution[i].SedimentHead)*sedstore
    68         eplvol[i]=np.mean(md.results.TransientSolution[i].EplHead)*store*np.mean(md.results.TransientSolution[i].HydrologydcEplThickness)
     68        eplvol[i]=np.mean(md.results.TransientSolution[i].EplHead)*store*np.mean(md.results.TransientSolution[i].HydrologydcEplThicknessHydrostep)
    6969        totvol[i+1]=totvol[i]+md.basalforcings.groundedice_melting_rate[0,i]*0.002
    7070
    7171field_names     =['SedimentWaterHead5','EplWaterHead5','SedimentWaterHead40','EplWaterHead40']
    7272field_tolerances=[1e-13, 1e-13, 1e-13, 1e-13, 1e-13]
    73 field_values=[md.results.TransientSolution[5].SedimentHead,
    74                                                         md.results.TransientSolution[5].EplHead,
    75                                                         md.results.TransientSolution[40].SedimentHead,
    76                                                         md.results.TransientSolution[40].EplHead]
     73field_values=[md.results.TransientSolution[5].SedimentHeadHydrostep,
     74                                                        md.results.TransientSolution[5].EplHeadHydrostep,
     75                                                        md.results.TransientSolution[40].SedimentHeadHydrostep,
     76                                                        md.results.TransientSolution[40].EplHeadHydrostep]

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

    r21480 r22857  
    2929field_names     ={'SedimentWaterHead','SedimentHeadResidual'};
    3030field_tolerances={1e-13, 3e-10};
    31 field_values={md.results.HydrologySolution.SedimentHead,md.results.HydrologySolution.SedimentHeadResidual};
     31field_values={md.results.HydrologySolution.SedimentHeadHydrostep,md.results.HydrologySolution.SedimentHeadResidual};

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

    r21480 r22857  
    4545field_names     =['SedimentWaterHead','SedimentHeadResidual']
    4646field_tolerances=[1e-13, 3e-10]
    47 field_values=[md.results.HydrologySolution.SedimentHead,md.results.HydrologySolution.SedimentHeadResidual]
     47field_values=[md.results.HydrologySolution.SedimentHeadHydrostep,md.results.HydrologySolution.SedimentHeadResidual]

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

    r22289 r22857  
    2828md.hydrology.epl_colapse_thickness=1.0e-3;
    2929md.hydrology.epl_max_thickness=1;
    30 md.timestepping.time_step=0.2;
     30md.hydrology.steps_per_step=10;
     31md.timestepping.time_step=2.0;
    3132md.timestepping.final_time=2.0;
    3233
    … …  
    3940mdfine.hydrology.steps_per_step=1;
    4041mdfine.timestepping.time_step=0.2;
    41 mdfine=solve(md,'Transient');
     42mdfine=solve(mdfine,'Transient');
    4243
    4344%store=md.constants.g*md.hydrology.sediment_porosity*md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility);
    … …  
    5253                  'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5',...
    5354                  'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9',...
    54                   'EplWaterHead9'};
     55                  'EplWaterHead10'};
    5556field_tolerances={...
    5657        1e-13, 1e-13, 1e-13,...
    … …  
    5960        1e-13, 5e-12, 1e-11,...
    6061  1e-13};
    61 field_values={mdfine.results.TransientSolution(1).SedimentHead, ...
    62                                                         mdfine.results.TransientSolution(1).EplHead,...
     62field_values={mdfine.results.TransientSolution(1).SedimentHeadHydrostep, ...
     63                                                        mdfine.results.TransientSolution(1).EplHeadHydrostep,...
    6364                                                        mdfine.results.TransientSolution(1).SedimentHeadResidual,...
    64                                                         mdfine.results.TransientSolution(4).SedimentHead,...
    65                                                         mdfine.results.TransientSolution(4).EplHead,...
     65                                                        mdfine.results.TransientSolution(4).SedimentHeadHydrostep,...
     66                                                        mdfine.results.TransientSolution(4).EplHeadHydrostep,...
    6667                                                        mdfine.results.TransientSolution(4).SedimentHeadResidual, ...
    67                                                         mdfine.results.TransientSolution(5).SedimentHead,...
    68                                                         mdfine.results.TransientSolution(5).EplHead,...
     68                                                        mdfine.results.TransientSolution(5).SedimentHeadHydrostep,...
     69                                                        mdfine.results.TransientSolution(5).EplHeadHydrostep,...
    6970                                                        mdfine.results.TransientSolution(5).SedimentHeadResidual, ...
    70                                                         mdfine.results.TransientSolution(9).SedimentHead,...
    71                                                         mdfine.results.TransientSolution(9).EplHead,...
     71                                                        mdfine.results.TransientSolution(9).SedimentHeadHydrostep,...
     72                                                        mdfine.results.TransientSolution(9).EplHeadHydrostep,...
    7273                                                        mdfine.results.TransientSolution(9).SedimentHeadResidual,...
    7374              md.results.TransientSolution(1).EplHead};

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

    r22287 r22857  
    4949
    5050#re-run with no substeps
    51 mdfine=md
     51mdfine=copy.deepcopy(md)
    5252mdfine.result=[]
    5353mdfine.hydrology.steps_per_step=1
    5454mdfine.timestepping.time_step=0.2
    55 mdfine=solve(md,'Transient')
     55mdfine=solve(mdfine,'Transient')
    5656
    5757
    58 # sedvol=np.zeros(10)
    59 # eplvol=np.zeros(10)
    60 # time=np.arange(0.2,2.1,0.2)
    61 # store=md.constants.g*md.hydrology.sediment_porosity*md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility)
    62 # sedstore=20.0*store
    63 # for i in range(0,10):
    64 #       sedvol[i]=np.mean(md.results.TransientSolution[i].SedimentHead)*sedstore
    65 #       eplvol[i]=np.mean(md.results.TransientSolution[i].EplHead)*store*np.mean(md.results.TransientSolution[i].HydrologydcEplThickness)
     58sedmean=mdfine.results.TransientSolution[0].SedimentHeadHydrostep
     59eplmean=mdfine.results.TransientSolution[0].EplHeadHydrostep
     60for i in range(1,10):
     61        sedmean=sedmean+(mdfine.results.TransientSolution[i].SedimentHeadHydrostep)
     62        eplmean=eplmean+(mdfine.results.TransientSolution[i].EplHeadHydrostep)
    6663
    6764field_names=['SedimentWaterHead1','EplWaterHead1','SedimentHeadResidual1',
    6865                                                 'SedimentWaterHead4','EplWaterHead4','SedimentHeadResidual4',
    6966                                                 'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5',
    70                                                  'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9', 
    71                                                  'EplWaterHead9']
     67                                                 'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9',
     68                                                 'EplWaterHead10']
    7269field_tolerances=[1e-13, 1e-13, 1e-13,
    7370                                                                        1e-13, 1e-13, 1e-13,
    … …  
    7572                                                                        1e-13, 5e-12, 1e-11,
    7673                                                                        1e-13]
    77 field_values=[mdfine.results.TransientSolution[0].SedimentHead,
    78                                                         mdfine.results.TransientSolution[0].EplHead,
     74field_values=[mdfine.results.TransientSolution[0].SedimentHeadHydrostep,
     75                                                        mdfine.results.TransientSolution[0].EplHeadHydrostep,
    7976                                                        mdfine.results.TransientSolution[0].SedimentHeadResidual,
    80                                                         mdfine.results.TransientSolution[3].SedimentHead,
    81                                                         mdfine.results.TransientSolution[3].EplHead,
     77                                                        mdfine.results.TransientSolution[3].SedimentHeadHydrostep,
     78                                                        mdfine.results.TransientSolution[3].EplHeadHydrostep,
    8279                                                        mdfine.results.TransientSolution[3].SedimentHeadResidual,
    83                                                         mdfine.results.TransientSolution[4].SedimentHead,
    84                                                         mdfine.results.TransientSolution[4].EplHead,
     80                                                        mdfine.results.TransientSolution[4].SedimentHeadHydrostep,
     81                                                        mdfine.results.TransientSolution[4].EplHeadHydrostep,
    8582                                                        mdfine.results.TransientSolution[4].SedimentHeadResidual,
    86                                                         mdfine.results.TransientSolution[8].SedimentHead,
    87                                                         mdfine.results.TransientSolution[8].EplHead,
     83                                                        mdfine.results.TransientSolution[8].SedimentHeadHydrostep,
     84                                                        mdfine.results.TransientSolution[8].EplHeadHydrostep,
    8885                                                        mdfine.results.TransientSolution[8].SedimentHeadResidual,
    89                                                         md.results.TransientSolution[0].EplHead]
    90 
     86                                                        md.results.TransientSolution[-1].EplHead]

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

    r21480 r22857  
    3030field_names     ={'SedimentWaterHead','SedimentHeadResidual'};
    3131field_tolerances={1e-13, 3e-10};
    32 field_values={md.results.HydrologySolution.SedimentHead,md.results.HydrologySolution.SedimentHeadResidual};
     32field_values={md.results.HydrologySolution.SedimentHeadHydrostep,md.results.HydrologySolution.SedimentHeadResidual};

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

    r21480 r22857  
    4545field_names     =['SedimentWaterHead','SedimentHeadResidual']
    4646field_tolerances=[1e-13, 3e-10]
    47 field_values=[md.results.HydrologySolution.SedimentHead,
     47field_values=[md.results.HydrologySolution.SedimentHeadHydrostep,
    4848                                                        md.results.HydrologySolution.SedimentHeadResidual]

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

    r21480 r22857  
    5858                                                1e-13, 5e-12, 2e-11,...
    5959                                                1e-13, 5e-12, 2e-11};
    60 field_values={md.results.TransientSolution(1).SedimentHead, ...
    61                                                         md.results.TransientSolution(1).EplHead,...
     60field_values={md.results.TransientSolution(1).SedimentHeadHydrostep, ...
     61                                                        md.results.TransientSolution(1).EplHeadHydrostep,...
    6262                                                        md.results.TransientSolution(1).SedimentHeadResidual,...
    63                                                         md.results.TransientSolution(4).SedimentHead,...
    64                                                         md.results.TransientSolution(4).EplHead,...
     63                                                        md.results.TransientSolution(4).SedimentHeadHydrostep,...
     64                                                        md.results.TransientSolution(4).EplHeadHydrostep,...
    6565                                                        md.results.TransientSolution(4).SedimentHeadResidual, ...
    66                                                         md.results.TransientSolution(5).SedimentHead,...
    67                                                         md.results.TransientSolution(5).EplHead,...
     66                                                        md.results.TransientSolution(5).SedimentHeadHydrostep,...
     67                                                        md.results.TransientSolution(5).EplHeadHydrostep,...
    6868                                                        md.results.TransientSolution(5).SedimentHeadResidual, ...
    69                                                         md.results.TransientSolution(9).SedimentHead,...
    70                                                         md.results.TransientSolution(9).EplHead,...
     69                                                        md.results.TransientSolution(9).SedimentHeadHydrostep,...
     70                                                        md.results.TransientSolution(9).EplHeadHydrostep,...
    7171                                                        md.results.TransientSolution(9).SedimentHeadResidual};

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

    r21480 r22857  
    11#Test Name: SquareSheetExtrudedHydrologyDCTwoLayers
    2 import numpy as np 
     2import numpy as np
    33from model import *
    44from socket import gethostname
    … …  
    5757                                                                        1e-13, 5e-12, 2e-11,
    5858                                                                        1e-13, 5e-12, 2e-11]
    59 field_values=[md.results.TransientSolution[0].SedimentHead,
    60                                                         md.results.TransientSolution[0].EplHead,
     59field_values=[md.results.TransientSolution[0].SedimentHeadHydrostep,
     60                                                        md.results.TransientSolution[0].EplHeadHydrostep,
    6161                                                        md.results.TransientSolution[0].SedimentHeadResidual,
    62                                                         md.results.TransientSolution[3].SedimentHead,
    63                                                         md.results.TransientSolution[3].EplHead,
     62                                                        md.results.TransientSolution[3].SedimentHeadHydrostep,
     63                                                        md.results.TransientSolution[3].EplHeadHydrostep,
    6464                                                        md.results.TransientSolution[3].SedimentHeadResidual,
    65                                                         md.results.TransientSolution[4].SedimentHead,
    66                                                         md.results.TransientSolution[4].EplHead,
     65                                                        md.results.TransientSolution[4].SedimentHeadHydrostep,
     66                                                        md.results.TransientSolution[4].EplHeadHydrostep,
    6767                                                        md.results.TransientSolution[4].SedimentHeadResidual,
    68                                                         md.results.TransientSolution[8].SedimentHead,
    69                                                         md.results.TransientSolution[8].EplHead,
     68                                                        md.results.TransientSolution[8].SedimentHeadHydrostep,
     69                                                        md.results.TransientSolution[8].EplHeadHydrostep,
    7070                                                        md.results.TransientSolution[8].SedimentHeadResidual]

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

    r22274 r22857  
    22import numpy as np
    33from scipy.interpolate import interp1d
    4 from model import  * 
    5 from setflowequation import  * 
     4from model import  *
     5from setflowequation import  *
    66from solve import  *
    77from NowickiProfile import *
    … …  
    5656
    5757#boundary conditions
    58 md.stressbalance.spcvx = float('NaN') * np.ones((md.mesh.numberofvertices,))
    59 md.stressbalance.spcvy = float('NaN') * np.ones((md.mesh.numberofvertices,))
    60 md.stressbalance.spcvz = float('NaN') * np.ones((md.mesh.numberofvertices,))
    61 md.stressbalance.referential = float('NaN') * np.ones((md.mesh.numberofvertices,6))
     58md.stressbalance.spcvx = np.nan * np.ones((md.mesh.numberofvertices,))
     59md.stressbalance.spcvy = np.nan * np.ones((md.mesh.numberofvertices,))
     60md.stressbalance.spcvz = np.nan * np.ones((md.mesh.numberofvertices,))
     61md.stressbalance.referential = np.nan * np.ones((md.mesh.numberofvertices,6))
    6262md.stressbalance.loadingforce = 0 * np.ones((md.mesh.numberofvertices,3))
    6363md.stressbalance.spcvx[np.where(md.mesh.vertexflags(4))] = 800.
    … …  
    7878md = setflowequation(md,'FS','all')
    7979md.flowequation.fe_FS = 'TaylorHood'
    80 md.stressbalance.abstol = float('NaN')
     80md.stressbalance.abstol = np.nan
    8181md.miscellaneous.name = 'test703'
    8282
    … …  
    8585md.timestepping.final_time = 0.000005
    8686md.stressbalance.shelf_dampening = 1.
    87 md.smb.mass_balance = np.zeros((md.mesh.numberofvertices,))
    88 md.basalforcings.groundedice_melting_rate = np.zeros((md.mesh.numberofvertices,))
    89 md.basalforcings.floatingice_melting_rate = np.zeros((md.mesh.numberofvertices,))
    90 md.basalforcings.geothermalflux = np.zeros((md.mesh.numberofvertices,))
     87md.smb.mass_balance = np.zeros((md.mesh.numberofvertices))
     88md.basalforcings.groundedice_melting_rate = np.zeros((md.mesh.numberofvertices))
     89md.basalforcings.floatingice_melting_rate = np.zeros((md.mesh.numberofvertices))
     90md.basalforcings.geothermalflux = np.zeros((md.mesh.numberofvertices))
    9191posb = np.intersect1d(np.where(md.mesh.x > 0.), np.where(md.mesh.vertexonbase))
    9292md.basalforcings.groundedice_melting_rate[posb] = 18.
    9393md.basalforcings.floatingice_melting_rate[posb] = 18.
    94 md.initialization.vx = np.zeros((md.mesh.numberofvertices,))
    95 md.initialization.vy = np.zeros((md.mesh.numberofvertices,))
    96 md.initialization.pressure = np.zeros((md.mesh.numberofvertices,))
    97 md.masstransport.spcthickness = float('NaN') * np.ones((md.mesh.numberofvertices,))
    98 md.thermal.spctemperature = float('NaN') * np.ones((md.mesh.numberofvertices,))
     94md.initialization.vx = np.zeros((md.mesh.numberofvertices))
     95md.initialization.vy = np.zeros((md.mesh.numberofvertices))
     96md.initialization.pressure = np.zeros((md.mesh.numberofvertices))
     97md.masstransport.spcthickness = np.nan * np.ones((md.mesh.numberofvertices))
     98md.thermal.spctemperature = np.nan * np.ones((md.mesh.numberofvertices))
    9999md.transient.isthermal = 0
    100100md.masstransport.isfreesurface = 1
Note: See TracChangeset for help on using the changeset viewer.