&           ProfDepthNo(num_file,bi,bj),prof_traj1D,num_var,

     &           profiles_dummy(num_file,num_var,bi,bj))

             call PROFILES_READVECTOR(num_file,num_var,

             call PROFILES_READVECTOR(num_file,-num_var,

               objf_profiles(num_file,num_var,bi,bj)=

     &             objf_profiles(num_file,num_var,bi,bj)

                 num_profiles(num_file,num_var,bi,bj)=

     &               num_profiles(num_file,num_var,bi,bj)

          do num_var=1,NVARMAX

     &      ProfNo(num_file,bi,bj),num_var

           write(msgbuf,'(a,D22.15,D22.15)') prof_names(num_var),

     &      objf_profiles(num_file,num_var,bi,bj),

     &      num_profiles(num_file,num_var,bi,bj)

      integer K,num_file,num_var,prof_num

           do num_var=1,NVARMAX

            if (vec_quantities(num_file,num_var,bi,bj).EQV..TRUE.) then

             if (num_var.EQ.6) then

     &         num_var,num_file,mytime,bi,bj,myThid)

     &         num_var,num_file,mytime,bi,bj,myThid)

      integer irec, num_var,num_file

      do num_var=1,NVARMAX

      err = NF_INQ_VARID(fid1,prof_names(num_var), varid0 )

      err = NF_DEF_VAR (fid2,prof_names(num_var), NF_DOUBLE,

     & 2,vecid,varid(2+(num_var-1)*2))

      err = NF_PUT_ATT_DOUBLE(fid2, varid(2+(num_var-1)*2),'_FillValue',

      err = NF_DEF_VAR (fid2,prof_namesmask(num_var),

     & NF_DOUBLE, 2,vecid, varid(3+(num_var-1)*2))

      err = NF_PUT_ATT_DOUBLE(fid2,varid(3+(num_var-1)*2),'_FillValue',

       do num_var=prof_num_var_tot(num_file,bi,bj),1,-1

     & +num_var-1)*2 +1) tmp_vec

     & +num_var-1)*2 +2) tmp_vec

      integer num_file,num_var

           do num_var=1,NVARMAX

            objf_profiles(num_file,num_var,bi,bj)= 0. _d 0

            num_profiles(num_file,num_var,bi,bj) = 0. _d 0

            profiles_dummy(num_file,num_var,bi,bj)=0. _d 0

      integer k,bi,bj,prof_num, num_file, num_var

      do num_var=1,NVARMAX

      if (vec_quantities(num_file,num_var,bi,bj).EQV..TRUE.) then

     & num_var,num_file,mytime,bi,bj,myThid)

     & num_var,num_file,mytime,bi,bj,myThid)

     &ProfDepthNo(num_file,bi,bj),prof_traj1D,num_var,

     &profiles_dummy(num_file,num_var,bi,bj))

      INTEGER iUnit, num_file, num_var

         do num_var=1,NVARMAX

          mult_profiles(num_file,num_var) = 0. _d 0

      INTEGER iUnit, num_file, num_var

      do num_var = 1, NGENCOST

      if ( gencost_barfile(num_var) .ne. ' ' ) then

     &     gencost_barfile(num_var), 

     &     gencost_barfld(1-Olx,1-Oly,1,1,num_var), 

     &     gencost_modfld(1-Olx,1-Oly,1,1,num_var),

     &     xx_genbar_dummy(num_var),

     &     startofgen(num_var), endofgen(num_var), ingen(num_var),

     &     sum1gen(num_var), genrec(num_var), 1, mythid )

      integer num_var

      integer num_var

      do num_var = 1, NGENCOST

         if ( gencost_avgperiod(num_var) .EQ. 'day' .OR.

     &        gencost_avgperiod(num_var) .EQ. 'DAY' ) then

            startofgen(num_var) = startofday

            endofgen(num_var)   = endofday

            ingen(num_var)      = inday

            sum1gen(num_var)    = sum1day

            genrec(num_var)     = dayrec

         else if ( gencost_avgperiod(num_var) .EQ. 'month' .OR.

     &        gencost_avgperiod(num_var) .EQ. 'MONTH' ) then

            startofgen(num_var) = startofmonth

            endofgen(num_var)   = endofmonth

            ingen(num_var)      = inmonth

            sum1gen(num_var)    = sum1mon

            genrec(num_var)     = monrec

         else if ( gencost_avgperiod(num_var) .EQ. 'year' .OR.

     &        gencost_avgperiod(num_var) .EQ. 'YEAR' ) then

            startofgen(num_var) = startofyear

            endofgen(num_var)   = endofyear

            ingen(num_var)      = inyear

            sum1gen(num_var)    = sum1year

            genrec(num_var)     = yearrec

      integer num_var

      integer num_var

      integer num_var

          do num_var=1,NSSHV4COST

     &                   + objf_sshv4cost(num_var,bi,bj)

     &                   + num_sshv4cost(num_var,bi,bj)

       do num_var=1,NVARMAX

        if (no_profiles(num_file,num_var).GT.0) then

     &    profilesfiles(num_file)(1:IL),' ',prof_names(num_var), ' = ',

     &    f_profiles(num_file,num_var),

     &    no_profiles(num_file,num_var)

       do num_var=1,NGENCOST

        if (no_gencost(num_var).GT.0) then

         IL  = ILNBLNK( gencost_name(num_var) )

     &    gencost_name(num_var)(1:IL),' (gencost ', num_var, ') = ',

     &    f_gencost(num_var),

     &    no_gencost(num_var)

       do num_var=1,NUSERCOST

        if (no_usercost(num_var).GT.0) then

     &    'usercost', num_var, ' = ',

     &    f_usercost(num_var),

     &    no_usercost(num_var)

       do num_var=1,NSSHV4COST

        if (no_sshv4cost(num_var).GT.0) then

     &    'sshv4cost', num_var, ' = ',

     &    f_sshv4cost(num_var),

     &    no_sshv4cost(num_var)

       do num_var=1,NVARMAX

        f_profiles(num_file,num_var)= 0. _d 0

       do num_var=1,NGENCOST

        f_gencost(num_var)= 0. _d 0

       do num_var=1,NUSERCOST

        f_usercost(num_var)= 0. _d 0

       do num_var=1,NSSHV4COST

        f_sshv4cost(num_var)= 0. _d 0

       do num_var=1,NVARMAX

        no_profiles(num_file,num_var)= 0. _d 0

       do num_var=1,NGENCOST

        no_gencost(num_var)= 0. _d 0

       do num_var=1,NUSERCOST

        no_usercost(num_var)= 0. _d 0

       do num_var=1,NSSHV4COST

        no_sshv4cost(num_var)= 0. _d 0

       do num_var=1,NVARMAX

     &            + mult_profiles(num_file,num_var)

     &            *objf_profiles(num_file,num_var,bi,bj)

       do num_var=1,NGENCOST

     &            + mult_gencost(num_var)

     &            *objf_gencost(bi,bj,num_var)

       do num_var=1,NUSERCOST

     &            + mult_usercost(num_var)

     &            *objf_usercost(num_var,bi,bj)

       do num_var=1,NSSHV4COST

     &            + mult_sshv4cost(num_var)

     &            *objf_sshv4cost(num_var,bi,bj)

       do num_var=1,NVARMAX

          f_profiles(num_file,num_var)=f_profiles(num_file,num_var)

     &            +objf_profiles(num_file,num_var,bi,bj)

       do num_var=1,NGENCOST

          f_gencost(num_var)=f_gencost(num_var)

     &            +objf_gencost(bi,bj,num_var)

       do num_var=1,NUSERCOST

          f_usercost(num_var)=f_usercost(num_var)

     &            +objf_usercost(num_var,bi,bj)

       do num_var=1,NSSHV4COST

          f_sshv4cost(num_var)=f_sshv4cost(num_var)

     &            +objf_sshv4cost(num_var,bi,bj)

      integer num_file,num_var

       do num_var=1,NVARMAX

          no_profiles(num_file,num_var)=no_profiles(num_file,num_var)

     &            +num_profiles(num_file,num_var,bi,bj)

       do num_var=1,NGENCOST

          no_gencost(num_var)=no_gencost(num_var)

     &            +num_gencost(bi,bj,num_var)

       do num_var=1,NUSERCOST

          no_usercost(num_var)=no_usercost(num_var)

     &            +num_usercost(num_var,bi,bj)

       do num_var=1,NSSHV4COST

          no_sshv4cost(num_var)=no_sshv4cost(num_var)

     &            +num_sshv4cost(num_var,bi,bj)

       do num_var=1,NVARMAX

        _GLOBAL_SUM_RL(f_profiles(num_file,num_var), myThid )

       do num_var=1,NGENCOST

        _GLOBAL_SUM_RL(f_gencost(num_var), myThid )

       do num_var=1,NUSERCOST

        _GLOBAL_SUM_RL(f_usercost(num_var), myThid )

       do num_var=1,NSSHV4COST

        _GLOBAL_SUM_RL(f_sshv4cost(num_var), myThid )

       do num_var=1,NVARMAX

        _GLOBAL_SUM_RL(no_profiles(num_file,num_var), myThid )

       do num_var=1,NGENCOST

        _GLOBAL_SUM_RL(no_gencost(num_var), myThid )

       do num_var=1,NUSERCOST

        _GLOBAL_SUM_RL(no_usercost(num_var), myThid )

       do num_var=1,NSSHV4COST

        _GLOBAL_SUM_RL(no_sshv4cost(num_var), myThid )

       do num_var=1,NVARMAX

        if (no_profiles(num_file,num_var).GT.0) then

     &     ' --> f_profiles =',f_profiles(num_file,num_var),

     &      num_file, num_var

       do num_var=1,NGENCOST

        if (no_gencost(num_var).GT.0) then

     &     ' --> f_gencost =',f_gencost(num_var),

     &      num_var

       do num_var=1,NUSERCOST

        if (no_usercost(num_var).GT.0) then

     &     ' --> f_usercost =',f_usercost(num_var),

     &      num_var

       do num_var=1,NSSHV4COST

        if (no_sshv4cost(num_var).GT.0) then

     &     ' --> f_sshv4cost =',f_sshv4cost(num_var),

     &      num_var

      integer num_file,num_var

              do num_var=1,NSSHV4COST

                wsshv4  (i,j,num_var,bi,bj) = 0. _d 0

      integer num_var

      do num_var=1,NSSHV4COST

      if ( sshv4cost_errfile(num_var) .NE. ' ' ) then

      call MDSREADFIELD( sshv4cost_errfile(num_var),

     &           *sshv4cost_errfactor(num_var)

        wsshv4(i,j,num_var,bi,bj)=wsshv4tmp(i,j,bi,bj)

     &   num_var, 0, mythid, dummy)

      INTEGER k, iUnit, num_file, num_var