c
      subroutine wrxdr(flnm,f,nx,ny,xx,yy,xlab,ylab,lab1,
     +  lab2,lab3,lab4,mtime,iclose)
c
c Write xdr data file containing field f(nx,ny)
c
c On input:
c   flnm = xdr file name
c   f(nx,ny) = data field to write (nx or ny can be 1 for 1d fields)
c   xx(nx) = x-axis values
c   yy(ny) = y-axis values
c   xlab,ylab = x,y axis labels
c   lab1,lab2,lab3,lab4 = information labels
c   mtime(3) = model time (day,hr,min) (set 0,0,0 for multi-ut plots)
c On output:
c   4-line header, xx, yy, and f, are written to lu (6e13.5 format for data)
c btf 4/8/97: added mtime (see also wrxdrc.c)
c
      character*(*) flnm,lab1,lab2,lab3,lab4,xlab,ylab
      real f(nx,ny),xx(nx),yy(ny)
      character*80 hdr1,hdr2,hdr3,hdr4,xlabel,ylabel
      integer mtime(3)
c
      if (iclose.gt.0) goto 100
c
c lab1:
c
      call clearstr(hdr1)
      lenlab = lenstr(lab1)
      if (lenlab.le.80) then
        write(hdr1,"(a)") lab1(1:lenlab)
      else
        write(hdr1,"(a)") lab1(1:80)
      endif
      lhdr1 = lenstr(hdr1)
c
c lab2:
c
      call clearstr(hdr2)
      lenlab = lenstr(lab2)
      if (lenlab.le.80) then
        write(hdr2,"(a)") lab2(1:lenlab)
      else
        write(hdr2,"(a)") lab2(1:80)
      endif
c
c lab3:
c
      call clearstr(hdr3)
      lenlab = lenstr(lab3)
      if (lenlab.le.80) then
        write(hdr3,"(a)") lab3(1:lenlab)
      else
        write(hdr3,"(a)") lab3(1:80)
      endif
c
c lab4:
c
      call clearstr(hdr4)
      lenlab = lenstr(lab4)
      if (lenlab.le.80) then
        write(hdr4,"(a)") lab4(1:lenlab)
      else
        write(hdr4,"(a)") lab4(1:80)
      endif
c
c x and y labels:
c
      call clearstr(xlabel)
      lenlab = lenstr(xlab)
      if (lenlab.le.80) then
        write(xlabel,"(a)") xlab(1:lenlab)
      else
        write(xlabel,"(a)") xlab(1:80)
      endif
c
      call clearstr(ylabel)
      lenlab = lenstr(ylab)
      if (lenlab.le.80) then
        write(ylabel,"(a)") ylab(1:lenlab)
      else
        write(ylabel,"(a)") ylab(1:80)
      endif
c
c int WRXDR(_fcd flnm, _fcd hdr1, _fcd hdr2, _fcd hdr3, _fcd hdr4,
c	  float *f2d, short *nx, short *ny, float *xx, float *yy,
c	  _fcd xlab, _fcd ylab, short *iclose)
c
 100  continue
      istat = wrxdrc(flnm(1:lenstr(flnm)),hdr1,hdr2,hdr3,hdr4,
     +  f,nx,ny,xx,yy,xlabel,ylabel,mtime,iclose)
      return
      end