printout now ok for kdm<1000 and idm,jdm<100,000

All large arrays are allocated at run time, and so HYCOM has no upper
limit on idm, jdm and kdm.  However, many print statements used i2, or
i3 with the expection that the 3rd digit is a space, for kdm and 2i4,
or 2i5 with the expection that the 3rd digit is a space, for idm,jdm.

All such print statements have been updated by adding a digit (i3,i4,
2i5,2i6) and so kdm of 1 to 999 and idm or jdm of 1 to 99999 are
fully supported.  There are extra spaces in .b archive files and in
the log output.

bigrid.F90

@@ -395,7 +395,7 @@ subroutine bigrid(depth, mapflg, util1,util2,util3)
           write (char3,'(i3)') ilast-ifrst
           fmt(8:10)=char3
           if     (mnproc.eq.1) then
-          write (lp,'(a,i5,a,i5)')  &
+          write (lp,'(a,i6,a,i6)')  &
             'ip array, cols',ifrst+1,' --',ilast
           endif
           do j= jtdm,1,-1
@@ -452,8 +452,8 @@ subroutine indxi(ipt,if,il,is)
             k = k+1
           elseif (last .eq. 0 .and. ipt(i,j) .eq. 1) then
             if     (k .gt. ms) then
-              write(lp,'(a,i5)')  'indxi problem on proc ',mnproc
-              write(lp,'(a,2i5)')  &
+              write(lp,'(a,i6)')  'indxi problem on proc ',mnproc
+              write(lp,'(a,2i6)')  &
                 ' error in indxi -- ms too small at i,j =',i0+i,j0+j
               call xchalt('(indxi)')
                      stop '(indxi)'
@@ -510,8 +510,8 @@ subroutine indxj(jpt,jf,jl,js)
             k = k+1
           elseif (last .eq. 0 .and. jpt(i,j) .eq. 1) then
             if     (k .gt. ms) then
-              write(lp,'(a,i5)')  'indxj problem on proc ',mnproc
-              write(lp,'(a,2i5)') &
+              write(lp,'(a,i6)')  'indxj problem on proc ',mnproc
+              write(lp,'(a,2i6)') &
                 ' error in indxj -- ms too small at i,j =',i0+i,j0+j
               call xchalt('(indxj)')
                      stop '(indxj)'
@@ -537,3 +537,4 @@ subroutine indxj(jpt,jf,jl,js)
 !> Oct  2008 - warning    on single-width inlets
 !> May  2014 - added ipim1,ipip1,ipjm1,ipjp1,ipim1x,ipip1x,ipjm1x,ipjp1x
 !> May  2014 - added allip,alliq,alliu,alliv
+!> Feb. 2025 - printout now ok for kdm<1000 and idm,jdm<100,000

blkdat.F90

@@ -69,7 +69,7 @@ subroutine blkdat(linit)
 !
       if     (itest.ne.itdm) then
         if (mnproc.eq.1) then
-        write(lp,'(/ a,i5 /)')  &
+        write(lp,'(/ a,i6 /)')  &
           'error - expected idm =',itdm
         call flush(lp)
         endif !1st tile
@@ -78,7 +78,7 @@ subroutine blkdat(linit)
       endif !error
       if     (jtest.ne.jtdm) then
         if (mnproc.eq.1) then
-        write(lp,'(/ a,i5 /)')  &
+        write(lp,'(/ a,i6 /)')  &
           'error - expected jdm =',jtdm
         call flush(lp)
         endif !1st tile
@@ -96,7 +96,7 @@ subroutine blkdat(linit)
 !
       if (ittest.gt.itdm) then
         if (mnproc.eq.1) then
-        write(lp,'(/ a,i5 /)')  &
+        write(lp,'(/ a,i6 /)')  &
           'error - maximum itest is',itdm
         call flush(lp)
         endif !1st tile
@@ -105,7 +105,7 @@ subroutine blkdat(linit)
       endif !error
       if (jttest.gt.jtdm) then
         if (mnproc.eq.1) then
-        write(lp,'(/ a,i5 /)')  &
+        write(lp,'(/ a,i6 /)')  &
           'error - maximum jtest is',jtdm
         call flush(lp)
         endif !1st tile
@@ -129,7 +129,7 @@ subroutine blkdat(linit)
 !     call xcsync(flush_lp)
 !     do k= 1,ijpr
 !       if     (mnproc.eq.k) then
-!         write(lp,'(a,3i5)') 'mnproc,[ij]test =',mnproc,itest,jtest
+!         write(lp,'(a,3i6)') 'mnproc,[ij]test =',mnproc,itest,jtest
 !       endif
 !       call xcsync(flush_lp)
 !     enddo !k
@@ -3087,3 +3087,4 @@ subroutine blkinl(lvar,cvar)
 !> Jan. 2025 - added tidnud to nudge towards the observed tides
 !> Feb. 2025 - Added cbtidc for adding tidal velocities to bottom speed
 !> Feb. 2025 - Negative cbar to input tidal amplitude flow speed
+!> Feb. 2025 - printout now ok for kdm<1000 and idm,jdm<100,000

cnuity.F90

@@ -389,7 +389,7 @@ subroutine cnuity(m,n)
 !diag   enddo
 !diag   enddo
 !diag endif
- 101  format (i9,2i5,i3,1p,e15.2,e30.2/a17,6e10.2/e37.2,e30.2)
+ 101  format (i9,2i6,i4,1p,e15.2,e30.2/a17,6e10.2/e37.2,e30.2)
 !
 ! --- at each grid point, determine the ratio of the largest permissible
 ! --- pos. (neg.) change in -dp- to the sum of all incoming (outgoing) fluxes
@@ -512,7 +512,7 @@ subroutine cnuity(m,n)
 !
 ! --- check for negative thicknesses.
 !
- 100  format (i9,' i,j,k=',2i5,i3,' neg. dp (m) in loop ',i2,g15.2,a)
+ 100  format (i9,' i,j,k=',2i6,i4,' neg. dp (m) in loop ',i2,g15.2,a)
 !
       if     (mod(nstep,3).eq.0) then  !skip some time steps for efficiency
         call xcminr(dpkmin(1:2*kk))
@@ -928,7 +928,7 @@ subroutine cnuity(m,n)
 !$OMP   END PARALLEL DO
 !
 !diag if (itest.gt.0.and.jtest.gt.0) &
-!diag write (lp,'(i9,2i5,i3,'' intfc.depth diffusion -- p_old,p_new ='', &
+!diag write (lp,'(i9,2i6,i4,'' intfc.depth diffusion -- p_old,p_new ='', &
 !diag 2f9.3)') nstep,itest,jtest,k,pold(itest,jtest)*qonem,p(itest, &
 !diag jtest,k)*qonem
 !
@@ -946,7 +946,7 @@ subroutine cnuity(m,n)
           do i=1-margin,ii+margin
             if (SEA_P) then
               if (p(i,j,k+1).lt.p(i,j,k)) then
-!diag           write (lp,'(i9,2i5,i3,a,g15.2,i4)') nstep,i+i0,j+j0,k, &
+!diag           write (lp,'(i9,2i6,i4,a,g15.2,i4)') nstep,i+i0,j+j0,k, &
 !diag           '  neg. dp after thknss smoothing', &
 !diag           qonem*(p(i,j,k+1)-p(i,j,k)),iflip
                 p(i,j,k+1)=p(i,j,k)
@@ -1080,7 +1080,7 @@ subroutine cnuity(m,n)
 !$OMP   END PARALLEL DO
 !
 !diag if (itest.gt.0.and.jtest.gt.0) &
-!diag write (lp,'(i9,2i5,i3," intfc.depth diffusion -- p_old,p_new =", &
+!diag write (lp,'(i9,2i6,i4," intfc.depth diffusion -- p_old,p_new =", &
 !diag 2f9.3)') nstep,itest+i0,jtest+j0,k,pold(itest,jtest)*qonem, &
 !diag p(itest,jtest,k)*qonem
 !
@@ -1098,7 +1098,7 @@ subroutine cnuity(m,n)
           do i=1-margin,ii+margin
             if (SEA_P) then
               if (p(i,j,k+1).lt.p(i,j,k)) then
-!diag           write (lp,'(i9,2i5,i3,a,g15.2,i4)') nstep,i+i0,j+j0,k, &
+!diag           write (lp,'(i9,2i6,i4,a,g15.2,i4)') nstep,i+i0,j+j0,k, &
 !diag           '  neg. dp after thknss smoothing', &
 !diag           qonem*(p(i,j,k+1)-p(i,j,k)),iflip
                 p(i,j,k+1)=p(i,j,k)
@@ -1449,3 +1449,4 @@ end subroutine cnuity
 !> Aug. 2018 - btrmas added, use onetamas to simplify logic
 !> Nov. 2018 - added oneta_u and oneta_v to correct and simplify logic
 !> Mar. 2023 - neg. dp in loop 19 is not fatal, might be corrected in loop 15
+!> Feb. 2025 - printout now ok for kdm<1000 and idm,jdm<100,000

convec.F90

@@ -34,7 +34,7 @@ subroutine convch(m,n)
 ! --- convective adjustment
 ! --- ---------------------
 !
- 103  format (i9,2i5,a/(33x,i3,2f8.3,f8.3,f8.2,f8.1))
+ 103  format (i9,2i6,a/(33x,i3,2f8.3,f8.3,f8.2,f8.1))
 !diag if (itest.gt.0 .and. jtest.gt.0) then
 !diag   write (lp,103) nstep,itest+i0,jtest+i0, &
 !diag   '  entering convec:  temp    saln    dens    thkns    dpth', &
@@ -242,7 +242,7 @@ subroutine convch(m,n)
               endif
               if     (llayer) then
 !$OMP           CRITICAL
-                write (lp,'(i9,2i5,i3,a,i3,a,i3,a,2f10.4)') &
+                write (lp,'(i9,2i6,i4,a,i4,a,i4,a,2f10.4)') &
                  nstep,i+i0,j+j0,k, &
                  ' colmn unstbl (wrt',ks,') after', &
                  iter-1,' its', &
@@ -276,7 +276,7 @@ subroutine convch(m,n)
 !diag     '  upr,lwr,final dens:',(sigup+thbase), &
 !diag     (siglo+thbase),(th3d(i,j,k,n)+thbase),q
 !diag endif
- 100    format (i9,2i5,i3,'  it',i2,a,3f8.3,f5.2)
+ 100    format (i9,2i6,i4,'  it',i2,a,3f8.3,f5.2)
 !
           end if
         end if
@@ -291,7 +291,7 @@ subroutine convch(m,n)
 !cc      colout(i)=colout(i)+temp(i,j,k,n)*dp(i,j,k,n)
 !cc      enddo !k
 !cc      if (abs((colout(i)-coluin(i))/coluin(i)).gt.1.e-6)
-!cc     .  write (lp,'(i9,2i5,a/1p,3e14.6)') nstep,i,j,
+!cc     .  write (lp,'(i9,2i6,a/1p,3e14.6)') nstep,i,j,
 !cc     .  '  column integral not conserved in convec:',
 !cc     .  coluin(i),colout(i),(colout(i)-coluin(i))/coluin(i)
       endif !ip
@@ -324,7 +324,7 @@ subroutine convcm(m,n)
 !
 # include "stmt_fns.h"
 !
- 103  format (i9,2i5,a/(33x,i3,2f8.3,3p,f8.3,0p,f8.2,f8.1))
+ 103  format (i9,2i6,a/(33x,i3,2f8.3,3p,f8.3,0p,f8.2,f8.1))
 !diag if (itest.gt.0 .and. jtest.gt.0) then
 !diag   write (lp,103) nstep,itest+i0,jtest+j0, &
 !diag   '  entering convec:  temp    saln    dens    thkns    dpth', &
@@ -381,8 +381,8 @@ subroutine convcm(m,n)
 !diag   +thbase,dp(i,j,1,n)*qonem,temp(i,j,1,n),saln(i,j,1,n),k, &
 !diag   th3d(i,j,k,n)+thbase,dp(i,j,k,n)*qonem,temp(i,j,k,n),saln(i,j,k,n)
 !diag endif
- 100    format (i9,2i5,a,i3,'  th3d,dp,t,s =',3pf7.3,0pf7.1,2f8.3 &
-          /26x,i3,15x,3pf7.3,0pf7.1,2f8.3)
+ 100    format (i9,2i6,a,i4,'  th3d,dp,t,s =',3pf7.3,0pf7.1,2f8.3 &
+          /26x,i4,15x,3pf7.3,0pf7.1,2f8.3)
 !
 ! --- layer -k- contains mass less dense than mixed layer. entrain it.
         delp=dp(i,j,1,n)+dp(i,j,k,n)
@@ -489,3 +489,4 @@ subroutine convcm(m,n)
 !> Oct  2010 - replaced two calls to dsiglocdX with one call at mid-point
 !> Aug  2011 - replaced dpold,dpoldm with dpo
 !> May  2014 - use land/sea masks (e.g. ip) to skip land
+!> Feb. 2025 - printout now ok for kdm<1000 and idm,jdm<100,000

diapfl.F90

@@ -482,11 +482,11 @@ subroutine diapf1aij(m,n, i,j)
 !
       return
 !
- 101  format(25x,'   thick      viscty    t diff    s diff  ' &
-           /(i9,2i5,i3,2x,4f10.2))
+ 101  format(28x,'   thick      viscty    t diff    s diff  ' &
+           /(i9,2i6,i4,2x,4f10.2))
  102  format(25x, &
       ' diff t  t old   t new   t chng  diff s  s old   s new   s chng' &
-           /(i9,2i5,i3,1x,8f8.3))
+           /(i9,2i6,i4,1x,8f8.3))
       end
       subroutine diapf1uij(m,n, i,j)
       use mod_xc         ! HYCOM communication interface
@@ -570,7 +570,7 @@ subroutine diapf1uij(m,n, i,j)
 !diag  (nstep,i+i0,j+j0,k,hm(k),u1do(k),u1dn(k),k=1,nlayer)
 !
       return
- 106  format(23x,'   thick   u old   u new'/(i9,2i5,i3,1x,f10.3,2f8.3))
+ 106  format(26x,'   thick   u old   u new'/(i9,2i6,i4,1x,f10.3,2f8.3))
       end
       subroutine diapf1vij(m,n, i,j)
       use mod_xc         ! HYCOM communication interface
@@ -655,7 +655,7 @@ subroutine diapf1vij(m,n, i,j)
 !diag  (nstep,i+i0,j+j0,k,hm(k),v1do(k),v1dn(k),k=1,nlayer)
 !
       return
- 107  format(23x,'   thick   v old   v new'/(i9,2i5,i3,1x,f10.3,2f8.3))
+ 107  format(26x,'   thick   v old   v new'/(i9,2i6,i4,1x,f10.3,2f8.3))
       end
 !
       subroutine diapf2(m,n)
@@ -747,7 +747,7 @@ subroutine diapf2j(m,n, j)
       enddo !ktr
 !
 !diag if (i.eq.itest.and.j.eq.jtest) &
-!diag   write (lp,'(i9,2i5,3x,a/(i36,4f10.3))') nstep,i+i0,j+j0, &
+!diag   write (lp,'(i9,2i6,3x,a/(i36,4f10.3))') nstep,i+i0,j+j0, &
 !diag   'before diapf2: thickness  salinity temperature density', &
 !diag   (k,dp(i,j,k,n)*qonem,saln(i,j,k,n), &
 !diag   temp(i,j,k,n),th3d(i,j,k,n)+thbase,k=1,kk)
@@ -768,7 +768,7 @@ subroutine diapf2j(m,n, j)
       enddo !k
 !
 !diag if (j.eq.jtest.and.itest.ge.ifp(j,l).and.itest.le.ilp(j,l)) &
-!diag   write (lp,'(i9,2i5,a,2i5)') &
+!diag   write (lp,'(i9,2i6,a,2i4)') &
 !diag     nstep,itest+i0,j+j0,' kmin,kmax =',kmin(itest),kmax(itest)
 !
 ! --- find buoyancy frequency for each layer
@@ -852,7 +852,7 @@ subroutine diapf2j(m,n, j)
         q=min(1.,.5*min(p(i,j,k)-p(i,j,k-1),p(i,j,k+2)-p(i,j,k+1))/ &
           max(flxu(i,k),flxl(i,k),epsil))
 !
-!diag   if (q.ne.1.) write (lp,'(i9,2i5,i3,a,1p,2e10.2,0p,2f7.2,f5.2)')  &
+!diag   if (q.ne.1.) write (lp,'(i9,2i6,i4,a,1p,2e10.2,0p,2f7.2,f5.2)')  &
 !diag     nstep,i+i0,j+j0,k,' flxu/l,dpu/l,q=',flxu(i,k),flxl(i,k), &
 !diag     (p(i,j,k)-p(i,j,k-1))*qonem,(p(i,j,k+2)-p(i,j,k+1))*qonem,q
 !
@@ -862,7 +862,7 @@ subroutine diapf2j(m,n, j)
       endif				!  kmin < k < kmax
 !
 !diag if (i.eq.itest.and.j.eq.jtest.and.k.ge.kmin(i).and.k.le.kmax(i)) &
-!diag    write (lp,'(i9,2i5,i3,3x,a/22x,f9.3,2f7.3,1p,3e10.3)') &
+!diag    write (lp,'(i9,2i6,i4,3x,a/22x,f9.3,2f7.3,1p,3e10.3)') &
 !diag    nstep,i+i0,j+j0,k, &
 !diag    'thknss   temp   saln    flngth      flxu      flxl', &
 !diag    dp(i,j,k,n)*qonem,temp(i,j,k,n),saln(i,j,k,n),flngth(i,k), &
@@ -999,10 +999,10 @@ subroutine diapf2j(m,n, j)
 !       if (abs(tndcys/tosal(i)).gt.1.e-11)
 !    .  write (lp,100)
 !    .  i+i0,j+i0,'  diapf2 saln.col.intgl.:',tosal(i),tndcys,clips(i)
-!100    format(2i5,a,1p,e16.8,2e13.5)
+!100    format(2i6,a,1p,e16.8,2e13.5)
 !
 !diag if (i.eq.itest.and.j.eq.jtest) &
-!diag   write (lp,'(i9,2i5,3x,a/(i36,0p,4f10.3))') &
+!diag   write (lp,'(i9,2i6,3x,a/(i36,0p,4f10.3))') &
 !diag   nstep,i+i0,j+j0, &
 !diag   'after  diapf2: thickness  salinity temperature density', &
 !diag   (k,dp(i,j,k,n)*qonem,saln(i,j,k,n), &
@@ -1090,7 +1090,7 @@ subroutine diapf3j(m,n, j)
       flxl(i,1)=0.
 !
 !diag if (i.eq.itest .and. j.eq.jtest) &
-!diag   write (lp,'(i9,2i5,3x,a/(i36,0p,3f10.3,3p,f10.3))') &
+!diag   write (lp,'(i9,2i6,3x,a/(i36,0p,3f10.3,3p,f10.3))') &
 !diag   nstep,i+i0,j+j0, &
 !diag   'before diapfl: thickness  salinity temperature density', &
 !diag   1,dp(i,j,1,n)*qonem,saln(i,j,1,n),temp(i,j,1,n), &
@@ -1119,7 +1119,7 @@ subroutine diapf3j(m,n, j)
       enddo !k
 !
 !diag if (j.eq.jtest.and.itest.ge.ifp(j,l).and.itest.le.ilp(j,l)) &
-!diag   write (lp,'(i9,2i5,a,2i5)') &
+!diag   write (lp,'(i9,2i6,a,2i4)') &
 !diag     nstep,itest+i0,j+j0,' kmin,kmax =',kmin(itest),kmax(itest)
 !
 ! --- temporarily swap layers  1  and  kmin-1
@@ -1229,7 +1229,7 @@ subroutine diapf3j(m,n, j)
           max(flxu(i,k),flxl(i,k),epsil))
 !
 !diag   if (i.eq.itest .and. j.eq.jtest .and. q.ne.1.) &
-!diag     write (lp,'(i9,2i5,i3,a,1p,2e10.2,0p,2f7.2,f4.2)')  &
+!diag     write (lp,'(i9,2i6,i4,a,1p,2e10.2,0p,2f7.2,f4.2)')  &
 !diag     nstep,i+i0,j+j0,k,' flxu/l,dpu/l,q=',flxu(i,k),flxl(i,k), &
 !diag     (p(i,j,k)-p(i,j,k-1))*qonem,(p(i,j,k+2)-p(i,j,k+1))*qonem,q
 !
@@ -1239,7 +1239,7 @@ subroutine diapf3j(m,n, j)
       end if				!  kmin < k < kmax
 !
 !diag if (i.eq.itest.and.j.eq.jtest.and.k.ge.kmin(i).and.k.le.kmax(i)) &
-!diag    write (lp,'(i9,2i5,i3,3x,a/22x,f9.3,2f7.3,1p,3e10.3)') &
+!diag    write (lp,'(i9,2i6,i4,3x,a/22x,f9.3,2f7.3,1p,3e10.3)') &
 !diag    nstep,i+i0,j+j0,k, &
 !diag    ' thknss   temp   saln   flngth      flxu      flxl', &
 !diag    dp(i,j,k,n)*qonem,temp(i,j,k,n),saln(i,j,k,n),flngth(i,k), &
@@ -1337,9 +1337,9 @@ subroutine diapf3j(m,n, j)
       if (k.eq.kmin(i) .and. p(i,j,k).lt..1*onemm) then
 !
 !diag   if (i.eq.itest .and. j.eq.jtest) then
-!diag   write (lp,'(i9,2i5,3x,a,i3,a)') nstep,i,j,'diapfl -- layer',k, &
+!diag   write (lp,'(i9,2i6,a,i4,a)') nstep,i,j,'diapfl -- layer',k, &
 !diag   ' erodes mixed layer'
-!diag   write (lp,'(i9,2i5,i3,3x,a/22x,f9.3,2f7.3,1p,3e10.3)') &
+!diag   write (lp,'(i9,2i6,i4,3x,a/22x,f9.3,2f7.3,1p,3e10.3)') &
 !diag    nstep,i+i0,j+j0,k, &
 !diag    ' thknss   temp   saln   flngth      flxu      flxl', &
 !diag    (p(i,j,k+1)-p(i,j,k))*qonem,temp(i,j,k,n), &
@@ -1364,7 +1364,7 @@ subroutine diapf3j(m,n, j)
 !
       dpmixl(i,j,n)=dp(i,j,1,n)
 !diag if (i.eq.itest.and.j.eq.jtest) &
-!diag   write (lp,'(i9,2i5,3x,a/(i36,0p,3f10.3,3p,f10.3))') &
+!diag   write (lp,'(i9,2i6,i4,a/(i36,0p,3f10.3,3p,f10.3))') &
 !diag   nstep,i+i0,j+j0, &
 !diag   'after  diapfl: thickness  salinity temperature density', &
 !diag   1,dp(i,j,1,n)*qonem,saln(i,j,1,n),temp(i,j,1,n), &
@@ -1403,3 +1403,4 @@ subroutine diapf3j(m,n, j)
 !> May  2014 - use land/sea masks (e.g. ip) to skip land
 !> July 2017 - added needed halo updates (xctilr)
 !> Dec. 2018 - add /* USE_NUOPC_CESMBETA */ macro and riv_input for coupled simulation
+!> Feb. 2025 - printout now ok for kdm<1000 and idm,jdm<100,000

forfun.F90

@@ -4175,21 +4175,21 @@ subroutine rdbaro_in(dtime,larchm,lslot)
         call xcsync(flush_lp)
         if     (i0.lt.ittest .and. i0+ii.ge.ittest .and. &
                 j0.lt.jttest .and. j0+jj.ge.jttest      ) then
-          write(lp,'(i5,i4,a,1p5e13.5)') &
+          write(lp,'(i6,i6,a,1p5e13.5)') &
              itest+i0,jtest+j0,' rdbaro: ub,vb,pb,ubp,vbp = ', &
              ubnest(itest,jtest,lslot), &
              vbnest(itest,jtest,lslot), &
              pbnest(itest,jtest,lslot), &
              ubpnst(itest,jtest,lslot), &
              vbpnst(itest,jtest,lslot)
-          write(lp,'(i5,i4,a,1p5e13.5)') &
+          write(lp,'(i6,i6,a,1p5e13.5)') &
              itest+1+i0,jtest+j0,' rdbaro: ub,vb,pb,ubp,vbp = ', &
              ubnest(itest+1,jtest,lslot), &
              vbnest(itest+1,jtest,lslot), &
              pbnest(itest+1,jtest,lslot), &
              ubpnst(itest+1,jtest,lslot), &
              vbpnst(itest+1,jtest,lslot)
-          write(lp,'(i5,i4,a,1p5e13.5)') &
+          write(lp,'(i6,i6,a,1p5e13.5)') &
              itest+i0,jtest+1+j0,' rdbaro: ub,vb,pb,ubp,vbp = ', &
              ubnest(itest,jtest+1,lslot), &
              vbnest(itest,jtest+1,lslot), &
@@ -4783,8 +4783,8 @@ subroutine rdnest_in(dtime,larchm,lslot)
         call xcsync(flush_lp)
         if     (i0.lt.ittest .and. i0+ii.ge.ittest .and. &
                 j0.lt.jttest .and. j0+jj.ge.jttest      ) then
- 103      format(i8,i5,i4,1x,a,a/ &
-                 (i8,5x,i4,1x,a,a,2f7.3,2f7.3,f8.4,f9.3,f9.2))
+ 103      format( i9,i6,i6,1x,a,a/ &
+                 (i9,6x,i6,1x,a,a,2f7.3,2f7.3,f8.4,f9.3,f9.2))
           write(lp,103) &
              nstep,itest+i0,jtest+j0,'rdnest', &
              ':   utot   vtot   temp   saln    dens    thkns     dpth', &
@@ -5047,3 +5047,4 @@ subroutine str2spd(wspd, tx,ty)
 !> Jan. 2025 - added forfunn for nudging towards the observed tides
 !> Jan. 2025 - salfac and hnudge in mod_tides
 !> Feb. 2025 - cbarmin (forfundf) ensures that BBL speed is not zero
+!> Feb. 2025 - printout now ok for kdm<1000 and idm,jdm<100,000