subroutine ADSTORE(chardum,int1,idow,int2,int3,icount)

      integer int1, int2, int3, idow, icount

      print *, 'adstore: ', chardum, int1, idow, int2, int3, icount

      if ( icount .gt. nidow ) then

      istoreidow(icount) = idow

      subroutine ADRESTO(chardum,int1,idow,int2,int3,icount)

      integer int1, int2, int3, idow, icount

      print *, 'adresto: ', chardum, int1, idow, int2, int3, icount

      if ( icount .gt. nidow ) then

      idow = istoreidow(icount)

      icount = 0

          icount=icount+1

      IF (icount .gt. 0) THEN

      INTEGER ii,jj,icount

      INTEGER ii, jj, icount, stdUnit

      icount = 0

           icount = icount + 1

      IF ( icount .gt. 0 ) THEN

     &      'THSICE_SOLVE4TEMP: there are ',icount,

       icount      = INT(tmp(5))

      print*, 'timesteps= ',icount

      allocate (time(icount))

      allocate (xpart(npart,icount))

      allocate (ypart(npart,icount))

      allocate (kpart(npart,icount))

      allocate (temp(npart,Nr,icount))

         allocate (uvel(npart,Nr,icount))

         allocate (vvel(npart,Nr,icount))

         allocate (salt(npart,Nr,icount))

      allocate (zeta(npart,icount))

         do n=1,icount

      do m=2,icount

      print*,icount,' x ',npart,' = ',icount*npart,' records expected,',

      print*,'==> ',icount*npart-itotalrecord,' float records missing'

      edges(1) = icount

      edges(2) = icount

      edges(3) = icount

      do icount=1,icount_transp

      if (kgen(icount).NE.0) then

            objf_gencost(bi,bj,kgen(icount))=0. _d 0

            num_gencost(bi,bj,kgen(icount))=0. _d 0

        nrecloc=gencost_nrec(kgen(icount))

          localperiod=gencost_period(kgen(icount))

          dummy=gencost_dummy(kgen(icount))

          spminloc=gencost_spmin(kgen(icount))

          spmaxloc=gencost_spmax(kgen(icount))

          spzeroloc=gencost_spzero(kgen(icount))

            preproc(k2)=gencost_preproc(k2,kgen(icount))

            preproc_i(k2)=gencost_preproc_i(k2,kgen(icount))

            preproc_c(k2)=gencost_preproc_c(k2,kgen(icount))

            preproc_r(k2)=gencost_preproc_r(k2,kgen(icount))

          il=ilnblnk(gencost_name(kgen(icount)))

     &    'Cost ',gencost_name(kgen(icount))(1:il),

          mybarfile=gencost_barfile(kgen(icount))

          if( .not. gencost_datafile(kgen(icount)).eq.' ') then

            call COST_GENCAL(mybarfile,gencost_errfile(kgen(icount)),

            il=ilnblnk(gencost_barfile(kgen(icount)))

     &        'globalsum transp ',gencost_barfile(kgen(icount))(1:il),

                objf_gencost(bi,bj,kgen(icount))=

     &             objf_gencost(bi,bj,kgen(icount))+tmpMeanTile(bi,bj)

              objf_gencost(bi,bj,kgen(icount))=

     &           objf_gencost(bi,bj,kgen(icount))*facCost

              num_gencost(bi,bj,kgen(icount))=nrecloc*facNum

      integer icount,icount_transp