Changeset 9607 for issm/trunk/test/NightlyRun/test1502.m

Timestamp:

09/06/11 10:08:46 (14 years ago)

Author:

Mathieu Morlighem

Message:

moved fields to surfaceforcings object

File:

• issm/trunk/test/NightlyRun/test1502.m (modified) (2 diffs)

issm/trunk/test/NightlyRun/test1502.m

@@ -15 +15 @@
 %Solve for thinning rate -> -1 * surface mass balance
 smb= 2*ones(md.numberofnodes,1);   
-md.surface_mass_balance= smb;
+md.surfaceforcings.mass_balance= smb;
 md.basal_melting_rate= smb;
 
@@ -22 +22 @@
 for i=1:10
          md=solve(md,PrognosticSolutionEnum);
-         md.surface_mass_balance= md.surface_mass_balance - (PatchToVec(md.results.PrognosticSolution.Thickness)-md.thickness);
+         md.surfaceforcings.mass_balance= md.surfaceforcings.mass_balance - (PatchToVec(md.results.PrognosticSolution.Thickness)-md.thickness);
 end
 
 %Set up transient
-smb = md.surface_mass_balance;
+smb = md.surfaceforcings.mass_balance;
 
 tooth= [ [ones(400,1)*(smb') - 10]' [ones(400,1)*(smb')]' ];
 smb=[ [ones(399,1)*(smb')]' smb  tooth tooth];
 
-md.surface_mass_balance= smb;
-md.surface_mass_balance(end+1,:)=[1:2000];
+md.surfaceforcings.mass_balance= smb;
+md.surfaceforcings.mass_balance(end+1,:)=[1:2000];
 
 md=solve(md,TransientSolutionEnum);