copyright (C) 2001  MSC-RPN COMM  %%%MC2%%%
*

      subroutine wrvref ( unf ) 1
      implicit none
*
      integer unf
*
#include "lcldim.cdk"
#include "levels.cdk"
#include "yomdyn1.cdk"
*
      integer  fstecr
      external fstecr
*
      integer k,err
      real wke (gnk+3)
*
      do k=1,gnk
         wke(k)= zt(k)
      end do
      wke(gnk+1)= htop
      wke(gnk+2)= iscal(1)
      wke(gnk+3)= iscal(2)
      err = fstecr (wke,wke,-32,unf,0,0,0,gnk+3,1,1,0,0,0,
     $          'C','HT','THERMO_H','X',0,0,0,0,5,.false.)
*
      return
      end