[Dart-dev] [4023] DART/trunk/DART_LAB/matlab/advance_oned.m: Added a model bias capability.

[nancy at ucar.edu]

An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20090828/5e0dbd22/attachment.html 
-------------- next part --------------
Modified: DART/trunk/DART_LAB/matlab/advance_oned.m
===================================================================
--- DART/trunk/DART_LAB/matlab/advance_oned.m	2009-08-28 15:49:03 UTC (rev 4022)
+++ DART/trunk/DART_LAB/matlab/advance_oned.m	2009-08-28 22:44:56 UTC (rev 4023)
@@ -1,4 +1,4 @@
-function x_new = advance_oned(x, alpha)
+function x_new = advance_oned(x, alpha, model_bias)
 
 % Data Assimilation Research Testbed -- DART
 % Copyright 2004-2009, Data Assimilation Research Section
@@ -11,15 +11,16 @@
 % $Revision$
 % $Date$
 
-x_new = x + comp_dt(x, alpha);
+x_new = x + comp_dt(x, alpha, model_bias);
 end
 
 %---------------------------------------------------
 
 % Internal function comp_dt
-function dx = comp_dt(x, alpha)
+function dx = comp_dt(x, alpha, model_bias)
 
 % Compute the time tendency; alpha controls nonlinearity
-dx = x + alpha .* x .* abs(x);
+% model_bias controls a shift in the model dynamics
+dx = (x + model_bias) + alpha .* x .* abs(x);
 
 end