diff --git a/var/da/da_minimisation/da_get_var_diagnostics.inc b/var/da/da_minimisation/da_get_var_diagnostics.inc
index c25236728d..aba293bcf1 100644
--- a/var/da/da_minimisation/da_get_var_diagnostics.inc
+++ b/var/da/da_minimisation/da_get_var_diagnostics.inc
@@ -13,7 +13,7 @@ subroutine da_get_var_diagnostics(it, iv, j)
    integer                      :: num_stats_tot
    integer                      :: i,k
    real                         :: jo_radiance
-   real                         :: temp(79)
+   real                         :: temp(82)
 
    if (trace_use) call da_trace_entry("da_get_var_diagnostics")
 
diff --git a/var/da/da_statistics/da_analysis_stats.inc b/var/da/da_statistics/da_analysis_stats.inc
index 652dbc9ba3..7ac2c831b2 100644
--- a/var/da/da_statistics/da_analysis_stats.inc
+++ b/var/da/da_statistics/da_analysis_stats.inc
@@ -85,7 +85,7 @@ use module_state_description, only : num_chem, PARAM_FIRST_SCALAR
                 'p     i    j          ', &
                 'q     i    j'
          case ( 1 )
-           write(unit=stats_unit, fmt='(8a/)') &
+           write(unit=stats_unit, fmt='(9a/)') &
                 ' Lvl         ', &
                 'u     i    j          ', &
                 'v     i    j          ', &
@@ -397,7 +397,7 @@ use module_state_description, only : num_chem, PARAM_FIRST_SCALAR
                 'p     i    j          ', &
                 'q     i    j'
          case ( 1 )
-           write(unit=stats_unit, fmt='(8a/)') &
+           write(unit=stats_unit, fmt='(9a/)') &
                 ' Lvl         ', &
                 'u     i    j          ', &
                 'v     i    j          ', &
@@ -1080,9 +1080,9 @@ use module_state_description, only : num_chem, PARAM_FIRST_SCALAR
             qcim=qcim+qciv(k)
             qsnm=qsnm+qsnv(k)
          end if
-		 if ( use_cv_w ) then
-		   wm=wm+wv(k)
-		 end if
+         if ( use_cv_w ) then
+            wm=wm+wv(k)
+         end if
       end do !k loop
    end if !rootproc