--- MITgcm/lsopt/lswri.F	2002/02/05 20:34:34	1.1
+++ MITgcm/lsopt/lswri.F	2002/11/15 04:03:24	1.2
@@ -0,0 +1,54 @@
+
+      subroutine lswri( phniter, NN, xx, gg, lphprint )
+
+c     ==================================================================
+c     SUBROUTINE lswri
+c     ==================================================================
+c
+c     o writes formatted output of updated controls and 
+c       adjoint variables
+c
+c     o started: Patrick Heimbach, MIT/EAPS
+c                01-Apr-1999
+c
+c     ==================================================================
+c     SUBROUTINE lswri
+c     ==================================================================
+
+      implicit none
+
+C     ========= Routine arguments ======================================
+
+      INTEGER   phniter, NN
+      DOUBLE PRECISION      xx(NN), gg(NN)
+      LOGICAL   lphprint
+
+cph      CHARACTER yitnum*3
+c
+
+C     ========= Local variables ========================================
+
+
+C=======================================================================
+
+      if (lphprint) 
+     &   print *, 'pathei: writing in lswri: ', phniter
+c
+      open(95,
+     &     FILE='PH_X.dat',
+     &     FORM='UNFORMATTED',
+     &     ACCESS='DIRECT',
+     &     RECL=NN*8)
+      write(95,REC=phniter+1) xx
+      close(95)
+c
+      open(95,
+     &     FILE='PH_G.dat',
+     &     FORM='UNFORMATTED',
+     &     ACCESS='DIRECT',
+     &     RECL=NN*8)
+      write(95,REC=phniter+1) gg
+      close(95)
+c
+
+      end