Index: /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 27260)
+++ /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 27261)
@@ -242,4 +242,5 @@
 			/*Retrieve autoregressive parameters*/
          parameters->AddObject(iomodel->CopyConstantObject("md.frontalforcings.ar_order",FrontalForcingsARMAarOrderEnum));
+         parameters->AddObject(iomodel->CopyConstantObject("md.frontalforcings.ma_order",FrontalForcingsARMAmaOrderEnum));
          parameters->AddObject(iomodel->CopyConstantObject("md.frontalforcings.arma_initialtime",FrontalForcingsARMAInitialTimeEnum));
          parameters->AddObject(iomodel->CopyConstantObject("md.frontalforcings.arma_timestep",FrontalForcingsARMATimestepEnum));
Index: /issm/trunk-jpl/src/c/modules/FrontalForcingsx/FrontalForcingsx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/FrontalForcingsx/FrontalForcingsx.cpp	(revision 27260)
+++ /issm/trunk-jpl/src/c/modules/FrontalForcingsx/FrontalForcingsx.cpp	(revision 27261)
@@ -53,5 +53,5 @@
    femmodel->parameters->FindParam(&numbasins,FrontalForcingsNumberofBasinsEnum);
    femmodel->parameters->FindParam(&arorder,FrontalForcingsARMAarOrderEnum);
-   femmodel->parameters->FindParam(&maorder,FrontalForcingsARMAarOrderEnum);
+   femmodel->parameters->FindParam(&maorder,FrontalForcingsARMAmaOrderEnum);
    IssmDouble tinit_arma;
    IssmDouble* termconstant  = NULL;
Index: /issm/trunk-jpl/src/m/classes/linearbasalforcingsarma.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/linearbasalforcingsarma.py	(revision 27260)
+++ /issm/trunk-jpl/src/m/classes/linearbasalforcingsarma.py	(revision 27261)
@@ -129,5 +129,4 @@
 
     def marshall(self, prefix, md, fid):  # {{{
-        print('ARMA is still under development')
         yts = md.constants.yts
 
Index: /issm/trunk-jpl/test/NightlyRun/test543.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test543.m	(revision 27260)
+++ /issm/trunk-jpl/test/NightlyRun/test543.m	(revision 27261)
@@ -49,5 +49,5 @@
 md.frontalforcings.basin_id             = idb_tf;
 md.frontalforcings.subglacial_discharge = 0.01*ones(md.mesh.numberofvertices,1);
-md.frontalforcings.const                = [0.005,0.001]; %intercept values of TF in basins [C]
+md.frontalforcings.const                = [0.05,0.01]; %intercept values of TF in basins [C]
 md.frontalforcings.trend                = [0.0001,0.00001]; %trend values of TF in basins [C/yr]
 md.frontalforcings.arma_initialtime     = md.timestepping.start_time; %initial time in the AR model parameterization [yr]
@@ -62,5 +62,5 @@
 
 % Covariance matrix
-covtf       = 1e-5*eye(nb_tf);
+covtf       = 1e-4*eye(nb_tf);
 covclv      = 1e-1*eye(nb_clv);
 covclv(1,1) = 1/10*covclv(1,1);
Index: /issm/trunk-jpl/test/NightlyRun/test543.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test543.py	(revision 27260)
+++ /issm/trunk-jpl/test/NightlyRun/test543.py	(revision 27261)
@@ -54,5 +54,5 @@
 md.frontalforcings.basin_id = idb_tf
 md.frontalforcings.subglacial_discharge = 0.01 * np.ones((md.mesh.numberofvertices,))
-md.frontalforcings.const = np.array([[0.005, 0.001]])  # intercept values of TF in basins [C]
+md.frontalforcings.const = np.array([[0.05, 0.01]])  # intercept values of TF in basins [C]
 md.frontalforcings.trend = np.array([[0.0001, 0.00001]])  # trend values of TF in basins [C/yr]
 md.frontalforcings.arma_initialtime = md.timestepping.start_time  # initial time in the AR model parameterization [yr]
@@ -76,6 +76,6 @@
 
 #Hard-coding covariance matrix because python is complaining
-covglob = np.array([[1e-5, 0., 0., 0., 0., 0.],
-                    [0., 1e-5, 0., 0., 0., 0.],
+covglob = np.array([[1e-4, 0., 0., 0., 0., 0.],
+                    [0., 1e-4, 0., 0., 0., 0.],
                     [0., 0., 1e-2, 0., 0., 0.],
                     [0., 0., 0., 1e-1, 0., 0.],
Index: /issm/trunk-jpl/test/NightlyRun/test544.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test544.py	(revision 27260)
+++ /issm/trunk-jpl/test/NightlyRun/test544.py	(revision 27261)
@@ -64,5 +64,5 @@
 md.basalforcings.arma_timestep  = 1.0  # timestep of the ARMA model [yr]
 md.basalforcings.arlag_coefs  = np.array([[0.0], [0.1]])  # autoregressive parameters
-md.basalforcings.malag_coefs  = np.array([[0.0], [0.0]])  # moving-average parameters
+md.basalforcings.malag_coefs  = np.array([[0.55], [0.34]])  # moving-average parameters
 md.basalforcings.deepwater_elevation = np.array([[-1000, -1520]])
 md.basalforcings.upperwater_elevation = np.array([[0, -50]])