if (gencost_barfile(kgen)(1:15).EQ.'m_boxmean_theta') then

            elseif (gencost_barfile(kgen)(1:14).EQ.'m_boxmean_salt')

            elseif (gencost_barfile(kgen)(1:17).EQ.'m_boxmean_ptracer')

          if (gencost_barfile(kgen)(1:13).EQ.'m_boxmean_eta') then

            if (gencost_barfile(kgen)(1:13).EQ.'m_horflux_vol') then

      if (gencost_barfile(kgen)(1:9).EQ.'m_boxmean') then

     &         gencost_barfile,

         gencost_barfile(k)    = ' '

         elseif ( (gencost_barfile(k)(1:9).EQ.'m_boxmean').OR.

     &            (gencost_barfile(k)(1:9).EQ.'m_horflux') ) then

         if ( (gencost_barfile(k)(1:7).EQ.'m_theta').OR.

     &        (gencost_barfile(k)(1:6).EQ.'m_salt').OR.

     &        (gencost_barfile(k)(1:8).EQ.'m_diffkr').OR.

     &        (gencost_barfile(k)(1:7).EQ.'m_kapgm').OR.

     &        (gencost_barfile(k)(1:9).EQ.'m_kapredi').OR.

     &        (gencost_barfile(k)(1:7).EQ.'m_trVol').OR.

     &        (gencost_barfile(k)(1:9).EQ.'m_trHeat').OR.

     &        (gencost_barfile(k)(1:9).EQ.'m_trSalt')

        if (gencost_barfile(k).NE.' ') then

         if (gencost_barfile(k)(1:2).NE.'m_') then

          il=ilnblnk(gencost_barfile(k))

          WRITE(msgBuf,'(A,A)') 'm_',gencost_barfile(k)(1:il)

          gencost_barfile(k)=msgBuf(1:il+2)

         if (gencost_barfile(k)(1:8).EQ.'m_tauZon') then

          il=ilnblnk(gencost_barfile(k))

          WRITE(msgBuf,'(A,A)') 'm_ustress',gencost_barfile(k)(9:il)

          gencost_barfile(k)=msgBuf(1:il+1)

         if (gencost_barfile(k)(1:8).EQ.'m_tauMer') then

          il=ilnblnk(gencost_barfile(k))

          WRITE(msgBuf,'(A,A)') 'm_vstress',gencost_barfile(k)(9:il)

          gencost_barfile(k)=msgBuf(1:il+1)

         il = ilnblnk(gencost_barfile(k))

         write(msgbuf,'(a,a)') ' model file = ',gencost_barfile(k)(1:il)

         if ((gencost_barfile(k)(1:9).EQ.'m_ptracer').OR.

     &       (gencost_barfile(k)(1:17).EQ.'m_boxmean_ptracer')) then

     &       gencost_barfile(k),

     &     gencost_barfile(k),

     &     nnzbar, gencost_barfile(k), gencost_dummy(k),

      il = ilnblnk( gencost_barfile(kgen) )

     &    gencost_barfile(kgen)(1:il),'.',eccoiter

      il=ilnblnk( gencost_barfile(kgen) )

     &     gencost_barfile(kgen)(1:il),'.',eccoiter

         if (gencost_barfile(k)(1:5).EQ.'m_eta') then

         elseif (gencost_barfile(k)(1:9).EQ.'m_boxmean') then

         elseif (gencost_barfile(k)(1:9).EQ.'m_horflux') then

         elseif (gencost_barfile(k)(1:5).EQ.'m_sst') then

         elseif (gencost_barfile(k)(1:5).EQ.'m_sss') then

         elseif (gencost_barfile(k)(1:4).EQ.'m_bp') then

         elseif (gencost_barfile(k)(1:16).EQ.'m_geothermalflux') then

         elseif (gencost_barfile(k)(1:9).EQ.'m_ustress') then

         elseif (gencost_barfile(k)(1:9).EQ.'m_vstress') then

         elseif (gencost_barfile(k)(1:7).EQ.'m_uwind') then

         elseif (gencost_barfile(k)(1:7).EQ.'m_vwind') then

         elseif (gencost_barfile(k)(1:7).EQ.'m_atemp') then

         elseif (gencost_barfile(k)(1:5).EQ.'m_aqh') then

         elseif (gencost_barfile(k)(1:8).EQ.'m_precip') then

         elseif (gencost_barfile(k)(1:8).EQ.'m_swdown') then

         elseif (gencost_barfile(k)(1:8).EQ.'m_lwdown') then

         elseif (gencost_barfile(k)(1:8).EQ.'m_wspeed') then

         elseif (gencost_barfile(k)(1:12).EQ.'m_bottomdrag') then

     &            (gencost_barfile(k)(1:8).EQ.'m_siarea') ) then

     &            (gencost_barfile(k)(1:8).EQ.'m_siheff') ) then  

         elseif (gencost_barfile(k)(1:9).EQ.'m_sihsnow') then

         elseif (gencost_barfile(k)(1:7).EQ.'m_theta') then

         elseif (gencost_barfile(k)(1:6).EQ.'m_salt') then

         elseif (gencost_barfile(k)(1:9).EQ.'m_ptracer') then

         elseif (gencost_barfile(k)(1:4).EQ.'m_UE') then

         elseif (gencost_barfile(k)(1:4).EQ.'m_VN') then

         elseif (gencost_barfile(k)(1:7).EQ.'m_trVol') then

         elseif (gencost_barfile(k)(1:8).EQ.'m_trHeat') then

         elseif (gencost_barfile(k)(1:8).EQ.'m_trSalt') then

         elseif (gencost_barfile(k)(1:8).EQ.'m_diffkr') then

         elseif (gencost_barfile(k)(1:7).EQ.'m_kapgm') then

         elseif (gencost_barfile(k)(1:9).EQ.'m_kapredi') then

      il=ilnblnk( gencost_barfile(kgen) )

     &     gencost_barfile(kgen)(1:il),'.',eccoiter

      areabarfile=gencost_barfile(igen_conc)

      if(igen_deconc.ne.0) deconcbarfile=gencost_barfile(igen_deconc)

      if(igen_exconc.ne.0) exconcbarfile=gencost_barfile(igen_exconc)

       if (gencost_barfile(igen_tp)(1:5).EQ.'m_eta') 

       if (gencost_barfile(igen_ers)(1:5).EQ.'m_eta') 

       if (gencost_barfile(igen_gfo)(1:5).EQ.'m_eta') 

      ilps=ilnblnk( gencost_barfile(igen_etaday) )

     &     gencost_barfile(igen_etaday)(1:ilps),'.',eccoiter

      ilps=ilnblnk( gencost_barfile(kgen) )

     &     gencost_barfile(kgen)(1:ilps),'.',eccoiter

          il=ilnblnk(gencost_barfile(kgen(icount_transp)))

     &    ': ',gencost_barfile(kgen(icount_transp))(1:il)

          mybarfile=gencost_barfile(kgen(icount))

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

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

     &       gencost_barfile,

      character*(MAX_LEN_FNAM) gencost_barfile(NGENCOST)

     &          (gencost_barfile(k)(1:9).EQ.'m_boxmean')) then

     &          (gencost_barfile(k)(1:9).EQ.'m_horflux')) then

     &          (gencost_barfile(k)(1:15).EQ.'m_boxmean_theta').OR.

     &          (gencost_barfile(k)(1:13).EQ.'m_boxmean_eta').OR.

     &          (gencost_barfile(k)(1:14).EQ.'m_boxmean_salt').OR.

     &          (gencost_barfile(k)(1:17).EQ.'m_boxmean_ptracer').OR.

     &          (gencost_barfile(k)(1:14).EQ.'m_horflux_vol')

            il=ilnblnk(gencost_barfile(k))

     &      gencost_barfile(k)(1:il),': has no matched model var.'

          if(.not.( (gencost_barfile(k)(1:7).EQ.'m_trVol') .OR.

     &              (gencost_barfile(k)(1:8).EQ.'m_trHeat').OR.

     &              (gencost_barfile(k)(1:8).EQ.'m_trSalt') )) then

           il=ilnblnk(gencost_barfile(k))

     &     gencost_barfile(k)(1:il),': has no matched model var.'

        if ( (gencost_barfile(k)(1:5).EQ.'m_eta').AND.

      ilt=ilnblnk( gencost_barfile(k) )

     &     gencost_barfile(k)(1:ilt), '.', eccoiter

     &     yadprefix, gencost_barfile(k)(1:ilt), '.', eccoiter

       if ( gencost_barfile(k).EQ.' ' ) 

         if ( gencost_barfile(k2).EQ.gencost_barfile(k) )