MITgcm_contrib/plumes/plumes_read_pickup.F

C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_read_pickup.F,v 1.2 2004/05/05 00:39:21 edhill Exp $
C $Name:  $

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

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

C     !USES:
CEOP
      IMPLICIT NONE

C     == Global variables ===
#include "CPP_OPTIONS.h"
#include "SIZE.h"
#include "fizhi_SIZE.h"
#include "land_SIZE.h"
#include "fizhi_coms.h"
#include "land_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_fizhi.',myIter
      prec = precFloat64

C-    Read fields as consecutive records
      CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)
      CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)
      CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)
      CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_read_pickup.F,v 1.2 2004/05/05 00:39:21 edhill Exp $
2	C $Name: $
3
4	CBOP
5	C !ROUTINE: FIZHI_READ_PICKUP
6	C !INTERFACE:
7	SUBROUTINE FIZHI_READ_PICKUP( myIter, myThid )
8
9	C !DESCRIPTION: \bv
10	C ==========================================================
11	C \| S/R FIZHI_READ_PICKUP
12	C \| o Reads current fizhi state from a pickup file
13	C ==========================================================
14	C \ev
15
16	C !USES:
17	CEOP
18	IMPLICIT NONE
19
20	C == Global variables ===
21	#include "CPP_OPTIONS.h"
22	#include "SIZE.h"
23	#include "fizhi_SIZE.h"
24	#include "land_SIZE.h"
25	#include "fizhi_coms.h"
26	#include "land_coms.h"
27	#include "EEPARAMS.h"
28	#include "PARAMS.h"
29
30	C == Routine Arguments ==
31	C myIter :: time-step number
32	C myThid :: Number of this instance
33	INTEGER myIter
34	INTEGER myThid
35
36	C !LOCAL VARIABLES:
37	C fn :: character buffer for creating filename
38	C prec :: precision of pickup files
39	C lgf :: flag to write "global" files
40	C k :: loop index
41	INTEGER prec, iChar, lChar, k
42	CHARACTER*(MAX_LEN_FNAM) fn
43	LOGICAL lgf
44
45	WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
46	prec = precFloat64
47
48	C- Read fields as consecutive records
49	CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)
50	CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)
51	CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)
52	CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)
53
54	RETURN
55	END