![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 13582
- Timestamp:
- 10/10/12 10:28:32 (12 years ago)
- Location:
- issm/trunk-jpl/test/NightlyRun
- Files:
-
- 4 added
- 8 edited
-
test230.m (modified) (3 diffs)
-
test230.py (added)
-
test231.m (modified) (3 diffs)
-
test231.py (added)
-
test232.m (modified) (3 diffs)
-
test232.py (added)
-
test233.m (modified) (2 diffs)
-
test233.py (added)
-
test234.m (modified) (2 diffs)
-
test235.m (modified) (2 diffs)
-
test236.m (modified) (3 diffs)
-
test237.m (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test230.m
r12791 r13582 1 md=triangle(model ,'../Exp/Square.exp',350000);1 md=triangle(model(),'../Exp/Square.exp',350000); 2 2 md=setmask(md,'all',''); 3 3 md=parameterize(md,'../Par/SquareShelf.par'); … … 11 11 12 12 %Set up transient 13 smb =ones(md.mesh.numberofvertices,1)*3.6;14 smb=[ smb smb*-1 ];13 smb=ones(md.mesh.numberofvertices,1)*3.6; 14 smb=[ smb smb*-1. ]; 15 15 16 md.surfaceforcings.mass_balance= 17 md.surfaceforcings.mass_balance(end+1,:)=[1.5 3 ];16 md.surfaceforcings.mass_balance=smb; 17 md.surfaceforcings.mass_balance(end+1,:)=[1.5 3.]; 18 18 md.transient.isthermal=0; 19 19 20 md=solve(md,TransientSolutionEnum );20 md=solve(md,TransientSolutionEnum()); 21 21 22 22 %Fields and tolerances to track changes 23 23 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', ... 24 24 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', ... 25 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', ...25 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', ... 26 26 'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceMassbalance4'}; 27 27 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... … … 29 29 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... 30 30 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 31 32 31 field_values={... 33 32 (md.results.TransientSolution(1).Vx),... -
issm/trunk-jpl/test/NightlyRun/test231.m
r12791 r13582 1 md=triangle(model ,'../Exp/Square.exp',350000);1 md=triangle(model(),'../Exp/Square.exp',350000); 2 2 md=setmask(md,'all',''); 3 3 md=parameterize(md,'../Par/SquareShelf.par'); … … 11 11 12 12 %Set up transient 13 smb =ones(md.mesh.numberofvertices,1)*3.6;14 smb=[ smb smb*2 ];13 smb=ones(md.mesh.numberofvertices,1)*3.6; 14 smb=[ smb smb*2. ]; 15 15 16 md.surfaceforcings.mass_balance= 17 md.surfaceforcings.mass_balance(end+1,:)=[1.5 3 ];16 md.surfaceforcings.mass_balance=smb; 17 md.surfaceforcings.mass_balance(end+1,:)=[1.5 3.]; 18 18 md.transient.isthermal=0; 19 19 20 md=solve(md,TransientSolutionEnum );20 md=solve(md,TransientSolutionEnum()); 21 21 22 22 %Fields and tolerances to track changes 23 23 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', ... 24 24 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', ... 25 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', ...25 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', ... 26 26 'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceMassbalance4'}; 27 27 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... … … 29 29 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... 30 30 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 31 32 31 field_values={... 33 32 (md.results.TransientSolution(1).Vx),... -
issm/trunk-jpl/test/NightlyRun/test232.m
r12791 r13582 1 md=triangle(model ,'../Exp/Square.exp',180000);1 md=triangle(model(),'../Exp/Square.exp',180000); 2 2 md=setmask(md,'all',''); 3 3 md=parameterize(md,'../Par/SquareShelf.par'); … … 5 5 md=setflowequation(md,'macayeal','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md.thermal.spctemperature=[md.thermal.spctemperature, md.thermal.spctemperature+5 , md.thermal.spctemperature+10, md.thermal.spctemperature+15; 1.5 2.5 3.5 4];7 md.thermal.spctemperature=[md.thermal.spctemperature, md.thermal.spctemperature+5., md.thermal.spctemperature+10., md.thermal.spctemperature+15.; 1.5 2.5 3.5 4.]; 8 8 md.timestepping.time_step=1; 9 9 md.timestepping.final_time=4; … … 12 12 md.transient.isthermal=1; 13 13 md.transient.isgroundingline=0; 14 md=solve(md,TransientSolutionEnum );14 md=solve(md,TransientSolutionEnum()); 15 15 16 16 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test233.m
r12791 r13582 1 md=triangle(model ,'../Exp/Square.exp',200000);1 md=triangle(model(),'../Exp/Square.exp',200000); 2 2 md=setmask(md,'all',''); 3 3 md=parameterize(md,'../Par/SquareShelf.par'); … … 5 5 md=setflowequation(md,'pattyn','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md.thermal.spctemperature=[md.thermal.spctemperature, md.thermal.spctemperature+5 ; 1 2];7 md.thermal.spctemperature=[md.thermal.spctemperature, md.thermal.spctemperature+5.; 1. 2.]; 8 8 md.timestepping.time_step=0.5; 9 md.timestepping.final_time=2 ;10 md=solve(md,TransientSolutionEnum );9 md.timestepping.final_time=2.; 10 md=solve(md,TransientSolutionEnum()); 11 11 12 12 %Fields and tolerances to track changes 13 13 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsMeltingRate1', ... 14 15 16 14 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsMeltingRate2', ... 15 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsMeltingRate3', ... 16 'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','Temperature4','BasalforcingsMeltingRate4'}; 17 17 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09, ... 18 19 20 18 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-07, ... 19 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06, ... 20 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06}; 21 21 field_values={... 22 22 (md.results.TransientSolution(1).Vx),... -
issm/trunk-jpl/test/NightlyRun/test234.m
r13021 r13582 1 md=triangle(model ,'../Exp/Square.exp',180000);1 md=triangle(model(),'../Exp/Square.exp',180000); 2 2 md=setmask(md,'all',''); 3 3 md=parameterize(md,'../Par/SquareShelf.par'); … … 55 55 56 56 %solve 57 md=solve(md,TransientSolutionEnum ,'overwrite','y');57 md=solve(md,TransientSolutionEnum(),'overwrite','y'); 58 58 md.qmu.results=md.results.dakota; 59 59 -
issm/trunk-jpl/test/NightlyRun/test235.m
r13021 r13582 1 md=triangle(model ,'../Exp/Square.exp',180000);1 md=triangle(model(),'../Exp/Square.exp',180000); 2 2 md=setmask(md,'all',''); 3 3 md=parameterize(md,'../Par/SquareShelf.par'); … … 55 55 56 56 %solve 57 md=solve(md,TransientSolutionEnum ,'overwrite','y');57 md=solve(md,TransientSolutionEnum(),'overwrite','y'); 58 58 md.qmu.results=md.results.dakota; 59 59 -
issm/trunk-jpl/test/NightlyRun/test236.m
r12883 r13582 1 md=triangle(model ,'../Exp/Square.exp',150000);1 md=triangle(model(),'../Exp/Square.exp',150000); 2 2 md=setmask(md,'all',''); 3 3 … … 19 19 % Same temperature over the all region : 20 20 imonth=0:11; 21 tmonth(1:12)=238.15+20 ;21 tmonth(1:12)=238.15+20.; 22 22 for imonth=0:11 23 23 md.surfaceforcings.temperatures_presentday(1:md.mesh.numberofvertices,imonth+1)=tmonth(imonth+1); 24 md.surfaceforcings.temperatures_lgm(1:md.mesh.numberofvertices,imonth+1)=tmonth(imonth+1)-20 ;24 md.surfaceforcings.temperatures_lgm(1:md.mesh.numberofvertices,imonth+1)=tmonth(imonth+1)-20.; 25 25 % Time for the last line: 26 26 md.surfaceforcings.temperatures_presentday(md.mesh.numberofvertices+1,imonth+1)=((imonth+1)/12); … … 44 44 % time steps and resolution 45 45 md.timestepping.time_step=20; 46 md.timestepping.final_time= 46 md.timestepping.final_time=60; 47 47 48 48 % 49 49 md=setflowequation(md,'macayeal','all'); 50 50 md.cluster=generic('name',oshostname(),'np',3); 51 md=solve(md,TransientSolutionEnum); 52 51 md=solve(md,TransientSolutionEnum()); 53 52 54 53 %Fields and tolerances to track changes 55 54 field_names ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3'}; 56 55 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... 57 58 56 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... 57 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 59 58 field_values={... 60 59 (md.results.TransientSolution(1).Vx),... -
issm/trunk-jpl/test/NightlyRun/test237.m
r12891 r13582 1 md=triangle(model ,'../Exp/Square.exp',600000);%1800001 md=triangle(model(),'../Exp/Square.exp',600000);%180000 2 2 md=setmask(md,'all',''); 3 3 … … 19 19 % Same temperature over the all region : 20 20 imonth=0:11; 21 tmonth(1:12)=238.15+20 ;21 tmonth(1:12)=238.15+20.; 22 22 for imonth=0:11 23 23 md.surfaceforcings.temperatures_presentday(1:md.mesh.numberofvertices,imonth+1)=tmonth(imonth+1); 24 md.surfaceforcings.temperatures_lgm(1:md.mesh.numberofvertices,imonth+1)=tmonth(imonth+1)-20 ;24 md.surfaceforcings.temperatures_lgm(1:md.mesh.numberofvertices,imonth+1)=tmonth(imonth+1)-20.; 25 25 % Time for the last line: 26 26 md.surfaceforcings.temperatures_presentday(md.mesh.numberofvertices+1,imonth+1)=((imonth+1)/12); … … 45 45 md.timestepping.time_step=20; 46 46 md.settings.output_frequency=1; 47 md.timestepping.final_time= 47 md.timestepping.final_time=60; 48 48 49 49 % … … 51 51 md=setflowequation(md,'macayeal','all'); 52 52 md.cluster=generic('name',oshostname(),'np',3); 53 md=solve(md,TransientSolutionEnum );53 md=solve(md,TransientSolutionEnum()); 54 54 55 55 %Fields and tolerances to track changes 56 56 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsMeltingRate1', ... 57 58 57 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsMeltingRate2', ... 58 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsMeltingRate3'}; 59 59 field_tolerances={... 60 60 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-8,... 61 61 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-8,... 62 62 1e-13,1e-13,1e-08,1e-13,1e-13,1e-10,1e-10,1e-10,1e-13,1e-8}; 63 64 63 field_values={... 65 64 (md.results.TransientSolution(1).Vx),...
Note:
See TracChangeset
for help on using the changeset viewer.
