Changeset 20350

Timestamp:

03/24/16 14:52:03 (9 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixed test 243

Location:

issm/trunk-jpl/test

Files:

• Archives/Archive243.nc (modified) ( previous)
• Data/gemb_input.mat (modified) ( previous)
• NightlyRun/test243.m (modified) (4 diffs)

issm/trunk-jpl/test/NightlyRun/test243.m

@@ -1 +1 @@
 %Test Name: SquareShelfSMBGemb
-md=triangle(model(),'../Exp/Square.exp',500000.);
+md=triangle(model(),'../Exp/Square.exp',200000.);
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
 md=setflowequation(md,'SSA','all');
-md.cluster=generic('name',oshostname(),'np',4); % 3 for the cluster
-
-%md.verbose=verbose('all');
+md.materials.rho_ice=910;
+md.cluster=generic('name',oshostname(),'np',3);
 
 % Use of Gemb method for SMB computation
@@ -14 +13 @@
 inputs=load('../Data/gemb_input.mat');
 
-%setup the time line: 
-deltatindays=(inputs.dateN(2) - inputs.dateN(1));
-dateN=(1979:deltatindays/365.25:2011)';
-dateN(length(inputs.dateN)+1:end)=[];
-
 %setup the inputs: 
-md.smb.Ta=[repmat(inputs.Ta0',md.mesh.numberofelements,1);dateN'];
-md.smb.V=[repmat(inputs.V0',md.mesh.numberofelements,1);dateN'];
-md.smb.dswrf=[repmat(inputs.dsw0',md.mesh.numberofelements,1);dateN'];
-md.smb.dlwrf=[repmat(inputs.dlw0',md.mesh.numberofelements,1);dateN'];
-md.smb.P=[repmat(inputs.P0',md.mesh.numberofelements,1);dateN'];
-md.smb.eAir=[repmat(inputs.eAir0',md.mesh.numberofelements,1);dateN'];
-md.smb.pAir=[repmat(inputs.pAir0',md.mesh.numberofelements,1);dateN'];
+md.smb.Ta=[repmat(inputs.Ta0',md.mesh.numberofelements,1);inputs.dateN'];
+md.smb.V=[repmat(inputs.V0',md.mesh.numberofelements,1);inputs.dateN'];
+md.smb.dswrf=[repmat(inputs.dsw0',md.mesh.numberofelements,1);inputs.dateN'];
+md.smb.dlwrf=[repmat(inputs.dlw0',md.mesh.numberofelements,1);inputs.dateN'];
+md.smb.P=[repmat(inputs.P0',md.mesh.numberofelements,1);inputs.dateN'];
+md.smb.eAir=[repmat(inputs.eAir0',md.mesh.numberofelements,1);inputs.dateN'];
+md.smb.pAir=[repmat(inputs.pAir0',md.mesh.numberofelements,1);inputs.dateN'];
 md.smb.Vz=repmat(inputs.LP.Vz,md.mesh.numberofelements,1);
 md.smb.Tz=repmat(inputs.LP.Tz,md.mesh.numberofelements,1);
@@ -34 +28 @@
 %smb settings
 md.smb.requested_outputs={'SmbDz','SmbT','SmbD','SmbRe','SmbGdn','SmbGsp','SmbEC','SmbA','SmbMassBalance'};
-
-md.materials.rho_ice=910;
-md.timestepping.interp_forcings=0;
 
 %only run smb core: 
@@ -46 +37 @@
 md.timestepping.start_time=1979;
 md.timestepping.final_time=1980;
-md.timestepping.time_step=.1;
+md.timestepping.time_step=.5;
+md.timestepping.interp_forcings=0;
 
 %Run transient
-%md=solve(md,TransientSolutionEnum());
+md=solve(md,TransientSolutionEnum());
 
 %Fields and tolerances to track changes
-field_names          ={'SmbDz','SmbT' ,'SmbD' ,'SmbRe','SmbGdn','SmbGsp','SmbA' ,'SmbEC','SmbMassBalance'};
-field_tolerances     ={'1e-13','1e-13','1e-13','1e-13','1e-13' ,'1e-13' ,'1e-13','1e-13','1e-13'};
+field_names      ={'SmbDz','SmbT' ,'SmbD' ,'SmbRe','SmbGdn','SmbGsp','SmbA' ,'SmbEC','SmbMassBalance'};
+field_tolerances ={'1e-13','1e-13','1e-13','1e-13','1e-13' ,'1e-13' ,'1e-13','1e-13','1e-13'};
 
 field_values={...