source: issm/oecreview/Archive/24307-24683/ISSM-24567-24568.diff@ 24684

Last change on this file since 24684 was 24684, checked in by Mathieu Morlighem, 5 years ago

CHG: added new review
File size: 7.2 KB
RevLine 
[24684]1Index: ../trunk-jpl/test/NightlyRun/test247.m
2===================================================================
3--- ../trunk-jpl/test/NightlyRun/test247.m (nonexistent)
4+++ ../trunk-jpl/test/NightlyRun/test247.m (revision 24568)
5@@ -0,0 +1,148 @@
6+%Test Name: SquareShelfTranIspddIsdeltaO18pdNoInterpSSA2d(
7+md=triangle(model(),'../Exp/Square.exp',60000.);
8+md=setmask(md,'all','');
9+md=parameterize(md,'../Par/SquareShelf.par');
10+
11+%md.verbose=verbose('all');
12+
13+% Use of ispdd and isdelta18o methods
14+md.smb = SMBd18opdd();
15+md.smb.isd18opd=1;
16+
17+% Add temperature, precipitation and delta18o needed to measure the surface mass balance
18+% creating delta18o
19+load '../Data/delta18o.data'
20+md.smb.delta18o=delta18o;
21+
22+% creating Present day temperatures
23+% Same temperature over the all region:
24+tmonth(1:12)=238.15+20.;
25+for imonth=0:11
26+ md.smb.temperatures_presentday(1:md.mesh.numberofvertices,imonth+1)=tmonth(imonth+1);
27+ % Time for the last line:
28+ md.smb.temperatures_presentday(md.mesh.numberofvertices+1,imonth+1)=((imonth+1)/12);
29+end
30+md.smb = initialize(md.smb,md);
31+
32+% creating initialization and spc temperatures initialization and
33+% spc
34+md.thermal.spctemperature=mean(md.smb.temperatures_presentday(1:md.mesh.numberofvertices,1:12),2)-10; %-10*ones(md.mesh.numberofvertices,1);
35+md.initialization.temperature=md.thermal.spctemperature; %md.smb.temperatures_presentday(1:md.mesh.numberofvertices,1);
36+
37+% creating precipitation
38+for imonth=0:11
39+ md.smb.precipitations_presentday(1:md.mesh.numberofvertices,imonth+1)=-0.4*10^(-6)*md.mesh.y+0.5;
40+ % Time for the last line:
41+ md.smb.precipitations_presentday(md.mesh.numberofvertices+1,imonth+1)=((imonth+1)/12);
42+end
43+
44+md.smb.s0p = (max(md.geometry.surface,0));
45+md.smb.s0t = (max(md.geometry.surface,0));
46+md.smb.issetpddfac = 1;
47+md.smb.pddfac_snow =4.3;
48+md.smb.pddfac_ice = 8.3;
49+
50+md=extrude(md,5,1.2);
51+md=setflowequation(md,'HO','all');
52+md.settings.results_on_nodes={'Temperature','Waterfraction','Enthalpy'};
53+
54+md.thermal.isenthalpy=1;
55+md.thermal.isdynamicbasalspc=1;
56+md.thermal.fe = 'P1xP2';
57+%md.thermal.spctemperature(find(md.mesh.vertexonbase)) = NaN;
58+
59+md.initialization.waterfraction= zeros(md.mesh.numberofvertices,1);
60+md.initialization.watercolumn = zeros(md.mesh.numberofvertices,1);
61+md.initialization.enthalpy = zeros(md.mesh.numberofvertices,1);
62+md.thermal.isdrainicecolumn=0;
63+
64+md=solve(md,'thermal');
65+
66+md.initialization.temperature=md.results.ThermalSolution.Temperature;
67+md.initialization.enthalpy=md.results.ThermalSolution.Enthalpy;
68+%x = find(md.initialization.temperature<210);
69+%md.initialization.temperature(x) = 268;
70+md.materials.rheology_B = cuffey(md.initialization.temperature);
71+
72+% Friction
73+TEMP = zeros(md.mesh.numberofvertices,1);
74+TEMP(md.mesh.elements) = md.initialization.temperature(:,1:6);
75+
76+temperature =TEMP;
77+pressure = md.initialization.pressure;%md.materials.rho_ice*md.constants.g*md.geometry.thickness;
78+Tm = md.materials.meltingpoint-md.materials.beta*pressure;
79+
80+md.friction.coefficient = md.friction.coefficient;%averaging(md,md.friction.coefficient,2);
81+md.friction=frictionjosh(md.friction);
82+md.friction.pressure_adjusted_temperature = temperature - Tm;
83+md.friction.gamma= 5;
84+
85+
86+% time steps and resolution
87+md.timestepping.time_step=0.5;
88+md.settings.output_frequency=1;
89+md.timestepping.final_time=2;
90+md.timestepping.interp_forcings=0;
91+
92+%md.transient.requested_outputs={'default','IceVolumeAboveFloatation','IceVolume','TemperaturePDD','SmbMonthlytemperatures','SmbPrecipitation'};
93+md.transient.requested_outputs={'default','IceVolumeAboveFloatation','IceVolume','TemperaturePDD'};
94+md=setflowequation(md,'SSA','all');
95+md.cluster=generic('name',oshostname(),'np',1); % 3 for the cluster
96+md=solve(md,'Transient');
97+
98+field_names = {...
99+ 'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','Enthalpy1','SmbMassBalance1', ...
100+ 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','Enthalpy2','SmbMassBalance2', ...
101+ 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','Enthalpy3','SmbMassBalance3', ...
102+ 'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','Temperature4','Enthalpy4','SmbMassBalance4'};
103+field_tolerances={...
104+ 1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,...
105+ 1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,...
106+ 1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,...
107+ 1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13};
108+field_values={...
109+ (md.results.TransientSolution(1).Vx),...
110+ (md.results.TransientSolution(1).Vy),...
111+ (md.results.TransientSolution(1).Vz),...
112+ (md.results.TransientSolution(1).Vel),...
113+ (md.results.TransientSolution(1).Pressure),...
114+ (md.results.TransientSolution(1).Base),...
115+ (md.results.TransientSolution(1).Surface),...
116+ (md.results.TransientSolution(1).Thickness),...
117+ (md.results.TransientSolution(1).Temperature),...
118+ (md.results.TransientSolution(1).Enthalpy),...
119+ (md.results.TransientSolution(1).SmbMassBalance),...
120+ (md.results.TransientSolution(2).Vx),...
121+ (md.results.TransientSolution(2).Vy),...
122+ (md.results.TransientSolution(2).Vz),...
123+ (md.results.TransientSolution(2).Vel),...
124+ (md.results.TransientSolution(2).Pressure),...
125+ (md.results.TransientSolution(2).Base),...
126+ (md.results.TransientSolution(2).Surface),...
127+ (md.results.TransientSolution(2).Thickness),...
128+ (md.results.TransientSolution(2).Temperature),...
129+ (md.results.TransientSolution(2).Enthalpy),...
130+ (md.results.TransientSolution(2).SmbMassBalance),...
131+ (md.results.TransientSolution(3).Vx),...
132+ (md.results.TransientSolution(3).Vy),...
133+ (md.results.TransientSolution(3).Vz),...
134+ (md.results.TransientSolution(3).Vel),...
135+ (md.results.TransientSolution(3).Pressure),...
136+ (md.results.TransientSolution(3).Base),...
137+ (md.results.TransientSolution(3).Surface),...
138+ (md.results.TransientSolution(3).Thickness),...
139+ (md.results.TransientSolution(3).Temperature),...
140+ (md.results.TransientSolution(3).Enthalpy),...
141+ (md.results.TransientSolution(3).SmbMassBalance),...
142+ (md.results.TransientSolution(4).Vx),...
143+ (md.results.TransientSolution(4).Vy),...
144+ (md.results.TransientSolution(4).Vz),...
145+ (md.results.TransientSolution(4).Vel),...
146+ (md.results.TransientSolution(4).Pressure),...
147+ (md.results.TransientSolution(4).Base),...
148+ (md.results.TransientSolution(4).Surface),...
149+ (md.results.TransientSolution(4).Thickness),...
150+ (md.results.TransientSolution(4).Temperature),...
151+ (md.results.TransientSolution(4).Enthalpy),...
152+ (md.results.TransientSolution(4).SmbMassBalance),...
153+ };
154Index: ../trunk-jpl/test/Archives/Archive247.arch
155===================================================================
156Cannot display: file marked as a binary type.
157svn:mime-type = application/octet-stream
158Index: ../trunk-jpl/test/Archives/Archive247.arch
159===================================================================
160--- ../trunk-jpl/test/Archives/Archive247.arch (nonexistent)
161+++ ../trunk-jpl/test/Archives/Archive247.arch (revision 24568)
162
163Property changes on: ../trunk-jpl/test/Archives/Archive247.arch
164___________________________________________________________________
165Added: svn:mime-type
166## -0,0 +1 ##
167+application/octet-stream
168\ No newline at end of property
Note: See TracBrowser for help on using the repository browser.