C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/plumes/plumes_read_pickup.F,v 1.2 2004/05/25 18:08:51 molod Exp $
C $Name:  $

CBOP
C     !ROUTINE: PLUMES_READ_PICKUP
C     !INTERFACE:
      SUBROUTINE PLUMES_READ_PICKUP( myIter, myThid )

C     !DESCRIPTION: \bv
C     *==========================================================*
C     | S/R PLUMES_READ_PICKUP
C     | o Reads current plumes state from a pickup file
C     *==========================================================*
C     \ev

C     !USES:
CEOP
      IMPLICIT NONE

C     == Global variables ===
#include "CPP_OPTIONS.h"
#include "SIZE.h"
#include "plumes_SIZE.h"
#include "land_SIZE.h"
#include "plumes_coms.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

C     == Routine Arguments ==
C     myIter  :: time-step number
C     myThid  :: Number of this instance
      INTEGER myIter
      INTEGER myThid

C     !LOCAL VARIABLES:
C     fn      :: character buffer for creating filename
C     prec    :: precision of pickup files
C     lgf     :: flag to write "global" files
C     k       :: loop index
      INTEGER prec, iChar, lChar, k
      CHARACTER*(MAX_LEN_FNAM) fn
      LOGICAL lgf

      WRITE(fn,'(A,I10.10)') 'pickup_plumes.',myIter
      prec = precFloat64

C-    Read fields as consecutive records
      CALL MDSREADFIELD(fn,prec,'RL',Nrplume,uplume,1,myThid)
      CALL MDSREADFIELD(fn,prec,'RL',Nrplume,thplume,2,myThid)

      RETURN
      END