pkg/fizhi/fizhi_read_pickup.F

C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_read_pickup.F,v 1.4 2004/07/14 17:31:57 molod Exp $
C $Name:  $

#include "FIZHI_OPTIONS.h"
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 "SIZE.h"
#include "fizhi_SIZE.h"
#include "fizhi_land_SIZE.h"
#include "fizhi_coms.h"
#include "fizhi_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     k       :: loop index
      INTEGER prec
      CHARACTER*(MAX_LEN_FNAM) fn

      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	molod	1.5	C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_read_pickup.F,v 1.4 2004/07/14 17:31:57 molod Exp $
2	edhill	1.2	C $Name: $
3
4	molod	1.5	#include "FIZHI_OPTIONS.h"
5	molod	1.1	CBOP
6			C !ROUTINE: FIZHI_READ_PICKUP
7			C !INTERFACE:
8			SUBROUTINE FIZHI_READ_PICKUP( myIter, myThid )
9
10			C !DESCRIPTION: \bv
11			C ==========================================================
12			C \| S/R FIZHI_READ_PICKUP
13			C \| o Reads current fizhi state from a pickup file
14			C ==========================================================
15			C \ev
16
17			C !USES:
18			CEOP
19			IMPLICIT NONE
20
21			C == Global variables ===
22			#include "SIZE.h"
23			#include "fizhi_SIZE.h"
24	molod	1.3	#include "fizhi_land_SIZE.h"
25	molod	1.1	#include "fizhi_coms.h"
26	molod	1.3	#include "fizhi_land_coms.h"
27	molod	1.1	#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 k :: loop index
40	molod	1.4	INTEGER prec
41	molod	1.1	CHARACTER*(MAX_LEN_FNAM) fn
42
43			WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
44			prec = precFloat64
45
46			C- Read fields as consecutive records
47			CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)
48			CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)
49			CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)
50			CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)
51
52			RETURN
53			END