anchor.diff:          swradf(1:im,1:jm) = swr
         end if
         ! read heat forcing corresponding to next theat
  -      if (iint.eq.1 .or. mod(iint,iheat).eq.0.) then
  +      if (iint.eq.1 .or. mod(iint+cont_bry,iheat).eq.0.) then
           do i=1,im
             do j=1,jm
               wtsurfb(i,j)=wtsurff(i,j)
  
positive.code:         wvsurff(1:im,1:jm) = wv
        end if
        ! read wind file corresponding to next time
        if (iint.eq.1 .or. mod(iint,iwind).eq.0.) then
          do i=1,im
            do j=1,jm
              wusurfb(i,j)=wusurff(i,j)
  
negative.code: 
        end if
        ! read bc file corresponding to next time
        if (iint.eq.1 .or. mod(iint+cont_bry,ibc).eq.0.) then
          tbwb = tbwf
          sbwb = sbwf
          tbeb = tbef
  