pkg/land/land_read_pickup.F

C $Header:  $
C $Name:  $

#include "LAND_OPTIONS.h"

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

C     !DESCRIPTION: \bv
C     *==========================================================*
C     | S/R LAND_READ_PICKUP
C     | o Reads current state of land model from a pickup file
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE

C     == Global variables ===
#include "LAND_SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "LAND_PARAMS.h"
#include "LAND_VARS.h"

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

#ifdef ALLOW_LAND

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
CEOP

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

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

C-    Read fields as consecutive records
      CALL MDSREADFIELD(fn,prec,'RL',land_nLev,land_groundT,1,myThid)
      CALL MDSREADFIELD(fn,prec,'RL',land_nLev,land_groundW,2,myThid)

C- jmc: exchange is not really necessary for land model
C-      and presently exchange S/R cannot work for 3d array where 3rd dim <> Nr

#endif /* ALLOW_LAND */

      RETURN
      END
1	jmc	1.1	C $Header: $
2			C $Name: $
3
4			#include "LAND_OPTIONS.h"
5
6			CBOP
7			C !ROUTINE: LAND_READ_PICKUP
8			C !INTERFACE:
9			SUBROUTINE LAND_READ_PICKUP( myIter, myThid )
10
11			C !DESCRIPTION: \bv
12			C ==========================================================
13			C \| S/R LAND_READ_PICKUP
14			C \| o Reads current state of land model from a pickup file
15			C ==========================================================
16			C \ev
17
18			C !USES:
19			IMPLICIT NONE
20
21			C == Global variables ===
22			#include "LAND_SIZE.h"
23			#include "EEPARAMS.h"
24			#include "PARAMS.h"
25			#include "LAND_PARAMS.h"
26			#include "LAND_VARS.h"
27
28			C == Routine Arguments ==
29			C myIter :: time-step number
30			C myThid :: Number of this instance
31			INTEGER myIter
32			INTEGER myThid
33
34			#ifdef ALLOW_LAND
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			CEOP
45
46			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
47
48			WRITE(fn,'(A,I10.10)') 'pickup_land.',myIter
49			prec = precFloat64
50
51			C- Read fields as consecutive records
52			CALL MDSREADFIELD(fn,prec,'RL',land_nLev,land_groundT,1,myThid)
53			CALL MDSREADFIELD(fn,prec,'RL',land_nLev,land_groundW,2,myThid)
54
55			C- jmc: exchange is not really necessary for land model
56			C- and presently exchange S/R cannot work for 3d array where 3rd dim <> Nr
57
58			#endif /* ALLOW_LAND */
59
60			RETURN
61			END