diff --git a/parallel/src/mcipparm_mod.f90 b/parallel/src/mcipparm_mod.f90
index 151ae54..ab0e703 100644
--- a/parallel/src/mcipparm_mod.f90
+++ b/parallel/src/mcipparm_mod.f90
@@ -367,8 +367,8 @@ MODULE mcipparm
 !-------------------------------------------------------------------------------
  
   CHARACTER(LEN=16),  PARAMETER     :: progname_nacc   = 'NACC'
-  CHARACTER(LEN=10),  PARAMETER     :: vdate_nacc      = '06/01/2023'
-  CHARACTER(LEN=8),   PARAMETER     :: ver_nacc        = 'V2.1.2'
+  CHARACTER(LEN=10),  PARAMETER     :: vdate_nacc      = '12/17/2023'
+  CHARACTER(LEN=8),   PARAMETER     :: ver_nacc        = 'V2.1.3'
 
 
   real projparm(6)  ! map projection type, lat1,lat2,projlon, center_lon, center_lat
diff --git a/parallel/src/rdfv3.f90 b/parallel/src/rdfv3.f90
index 86a3b03..02924f3 100644
--- a/parallel/src/rdfv3.f90
+++ b/parallel/src/rdfv3.f90
@@ -829,7 +829,7 @@ END SUBROUTINE windrotation
           yoff = -2.0  ! V-face: no offset in Y from dot-point center value
 
           DO j = 1, nrows_x+1
-            DO i = 1, nrows_x+1  ! use all X to fill array; last col outside domain
+            DO i = 1, ncols_x+1  ! use all X to fill array; last col outside domain
 
               xxin = met_xxctr + (FLOAT(i) + xoff) * met_resoln
               yyin = met_yyctr + (FLOAT(j) + yoff) * met_resoln
@@ -962,7 +962,7 @@ END SUBROUTINE windrotation
           yoff = -2.0  ! V-face: no offset in Y from dot-point center value
 
           DO j = 1, nrows_x+1
-            DO i = 1, nrows_x+1  ! use all X to fill array; last col outside domain
+            DO i = 1, ncols_x+1  ! use all X to fill array; last col outside domain
 
               xxin = met_xxctr + (FLOAT(i) + xoff) * met_resoln
               yyin = met_yyctr + (FLOAT(j) + yoff) * met_resoln
diff --git a/parallel/src/rdfv3_lam.f90 b/parallel/src/rdfv3_lam.f90
index c2dda7a..1cf9bcb 100644
--- a/parallel/src/rdfv3_lam.f90
+++ b/parallel/src/rdfv3_lam.f90
@@ -861,7 +861,7 @@ END SUBROUTINE windrotation
           yoff = -2.0  ! V-face: no offset in Y from dot-point center value
 
           DO j = 1, nrows_x+1
-            DO i = 1, nrows_x+1  ! use all X to fill array; last col outside domain
+            DO i = 1, ncols_x+1  ! use all X to fill array; last col outside domain
 
               xxin = met_xxctr + (FLOAT(i) + xoff) * met_resoln
               yyin = met_yyctr + (FLOAT(j) + yoff) * met_resoln
diff --git a/serial/src/mcipparm_mod.f90 b/serial/src/mcipparm_mod.f90
index 65e0586..b593e85 100644
--- a/serial/src/mcipparm_mod.f90
+++ b/serial/src/mcipparm_mod.f90
@@ -367,8 +367,8 @@ MODULE mcipparm
 !-------------------------------------------------------------------------------
  
   CHARACTER(LEN=16),  PARAMETER     :: progname_nacc   = 'NACC'
-  CHARACTER(LEN=10),  PARAMETER     :: vdate_nacc      = '06/01/2023'
-  CHARACTER(LEN=8),   PARAMETER     :: ver_nacc        = 'V2.1.2'
+  CHARACTER(LEN=10),  PARAMETER     :: vdate_nacc      = '12/17/2023'
+  CHARACTER(LEN=8),   PARAMETER     :: ver_nacc        = 'V2.1.3'
 
 
   real projparm(6)  ! map projection type, lat1,lat2,projlon, center_lon, center_lat
diff --git a/serial/src/rdfv3.f90 b/serial/src/rdfv3.f90
index fb9d9c6..b2590d3 100644
--- a/serial/src/rdfv3.f90
+++ b/serial/src/rdfv3.f90
@@ -826,7 +826,7 @@ END SUBROUTINE windrotation
           yoff = -2.0  ! V-face: no offset in Y from dot-point center value
 
           DO j = 1, nrows_x+1
-            DO i = 1, nrows_x+1  ! use all X to fill array; last col outside domain
+            DO i = 1, ncols_x+1  ! use all X to fill array; last col outside domain
 
               xxin = met_xxctr + (FLOAT(i) + xoff) * met_resoln
               yyin = met_yyctr + (FLOAT(j) + yoff) * met_resoln
@@ -957,7 +957,7 @@ END SUBROUTINE windrotation
           yoff = -2.0  ! V-face: no offset in Y from dot-point center value
 
           DO j = 1, nrows_x+1
-            DO i = 1, nrows_x+1  ! use all X to fill array; last col outside domain
+            DO i = 1, ncols_x+1  ! use all X to fill array; last col outside domain
 
               xxin = met_xxctr + (FLOAT(i) + xoff) * met_resoln
               yyin = met_yyctr + (FLOAT(j) + yoff) * met_resoln
diff --git a/serial/src/rdfv3_lam.f90 b/serial/src/rdfv3_lam.f90
index 9ef2ada..0277fad 100644
--- a/serial/src/rdfv3_lam.f90
+++ b/serial/src/rdfv3_lam.f90
@@ -860,7 +860,7 @@ END SUBROUTINE windrotation
           yoff = -2.0  ! V-face: no offset in Y from dot-point center value
 
           DO j = 1, nrows_x+1
-            DO i = 1, nrows_x+1  ! use all X to fill array; last col outside domain
+            DO i = 1, ncols_x+1  ! use all X to fill array; last col outside domain
 
               xxin = met_xxctr + (FLOAT(i) + xoff) * met_resoln
               yyin = met_yyctr + (FLOAT(j) + yoff) * met_resoln