source: trunk/SOURCES/vtoutput-0.2.f @ 23

!> \file vtoutput-0.2.f
!! TOOOOOOOO DOOOOOOOO
!<

!> SUBROUTINE: VTOUTPUT()
!! \author ...
!! \date ...
!! @note TOOOOOOO DOOOOOOO
!! @note Used modules:
!! @note    - use module3D_phy
!! \param FILENAME
!<

c  *****************************************************************************
      subroutine VTOUTPUT(FILENAME)

       USE module3D_phy

      implicit none
!      real,dimension(nz) :: E   ! vertical coordinate in ice, scaled to H zeta
      character FILENAME*13
c     char SccsId[] = "@(#)vtoutput.f  1.3   date delta= 01/06/95    heure delta= 14:56:26    date lecture= 01/09/95"
      
      E(1)=0.
      E(NZ)=1.
      do K=1,NZ
         if ((K.ne.1).and.(K.ne.NZ)) E(K)=(K-1.)/(NZ-1.)
      end do
      
      open(num_file2,file=FILENAME)
      
      do K=NZ,1,-1

c     diverses erreurs corrigees (semble t'il)

        W16(K)=UZR(16,16,K)+(UX(16,16,K)+UX(17,16,K))/2
     &    *((E(K)*(HMX(17,16)-HMX(16,16))/DX)
     &    -(SDX(17,16)+SDX(16,16))/2)
     &    - CDE(K)*HDOT(16,16)

        W24(K)=UZR(24,16,K)+(UX(24,16,K)+UX(25,16,K))/2
     &    *((E(K)*(HMX(25,16)-HMX(24,16))/DX)
     &    -(SDX(25,16)+SDX(24,16))/2)
     &    - CDE(K)*HDOT(24,16)

        write(num_file2,999) CDE(K),T(16,16,K)-TPMP(16,16,K)+
     &     273.15,-W16(K),
     &     CDE(K),(UX(24,16,K)+UX(25,16,K))/2.,
     &     -W24(K),T(24,16,K)-TPMP(24,16,K)+273.15

      end do
      close(num_file2)

999   format(7(1X,f9.3))
      end