Changeset 19554 for issm/trunk-jpl/test/NightlyRun/gemb.m

Timestamp:

09/16/15 17:15:08 (10 years ago)

Author:

Eric.Larour

Message:

CHG: finished first prototype implementation of the GEMB model inside the Gembx module of the SurfaceMassBalancex module,
and inside the SMBgemb.m class, as well as the SMBGemb method of the Element.cpp class.
Had to introduce the concept also of a DoubleArrayInput to keep track of the vertical gridded datasets within each element.

File:

• issm/trunk-jpl/test/NightlyRun/gemb.m (modified) (1 diff)

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

@@ -23 +23 @@
 md.smb.pAir=[repmat(inputs.pAir0',md.mesh.numberofelements,1);dateN'];
 
-md.smb.Vz=inputs.LP.Vz;
-md.smb.Tz=inputs.LP.Tz;
-md.smb.Tmean=inputs.LP.Tmean;
-md.smb.C=inputs.LP.C;
+md.smb.Vz=repmat(inputs.LP.Vz,md.mesh.numberofelements,1);
+md.smb.Tz=repmat(inputs.LP.Tz,md.mesh.numberofelements,1);
+md.smb.Tmean=repmat(inputs.LP.Tmean,md.mesh.numberofelements,1);
+md.smb.C=repmat(inputs.LP.C,md.mesh.numberofelements,1);
 
 %settings
-md.smb.spinUp=2;
+md.smb.ismelt=0;
+md.smb.isthermal=0;
+md.smb.isalbedo=0;
+md.transient.isstressbalance=0;
+md.transient.ismasstransport=0;
+md.transient.isthermal=0;
+md.timestepping.start_time=2000;
+md.timestepping.final_time=2001;
+md.timestepping.time_step=1/12; %monthly
+
+md.verbose=verbose('11111111');
 
 %Run transient