pkg/fizhi/fizhi_read_pickup.F

C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_read_pickup.F,v 1.6 2004/08/06 21:12:14 edhill 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 "fizhi_io_comms.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#ifdef ALLOW_MNC
#include "MNC_PARAMS.h"
#endif

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 i, prec
      CHARACTER*(MAX_LEN_FNAM) fn

      DO i = 1,MAX_LEN_FNAM
        fn(i:i) = ' '
      ENDDO
      WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
      prec = precFloat64

C#ifdef ALLOW_MDSIO
      IF ( fizhi_mdsio_read_pickup ) THEN

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)

      ENDIF

#ifdef ALLOW_MNC
      IF ( useMNC .and. fizhi_mnc_read_pickup ) THEN

        _BEGIN_MASTER(myThid)

        WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter

        CALL MNC_CW_RL_R('D',fn,0,0,'uphy',uphy,myThid)
        CALL MNC_CW_RL_R('D',fn,0,0,'vphy',vphy,myThid)
        CALL MNC_CW_RL_R('D',fn,0,0,'thphy',thphy,myThid)
        CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)


        _END_MASTER()

      ENDIF
#endif

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_read_pickup.F,v 1.6 2004/08/06 21:12:14 edhill Exp $
2	C $Name: $
3
4	#include "FIZHI_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: FIZHI_READ_PICKUP
8	C !INTERFACE:
9	SUBROUTINE FIZHI_READ_PICKUP( myIter, myThid )
10
11	C !DESCRIPTION: \bv
12	C ==========================================================
13	C \| S/R FIZHI_READ_PICKUP
14	C \| o Reads current fizhi state from a pickup file
15	C ==========================================================
16	C \ev
17
18	C !USES:
19	CEOP
20	IMPLICIT NONE
21
22	C == Global variables ===
23	#include "SIZE.h"
24	#include "fizhi_SIZE.h"
25	#include "fizhi_land_SIZE.h"
26	#include "fizhi_coms.h"
27	#include "fizhi_land_coms.h"
28	#include "fizhi_io_comms.h"
29	#include "EEPARAMS.h"
30	#include "PARAMS.h"
31	#ifdef ALLOW_MNC
32	#include "MNC_PARAMS.h"
33	#endif
34
35	C == Routine Arguments ==
36	C myIter :: time-step number
37	C myThid :: Number of this instance
38	INTEGER myIter
39	INTEGER myThid
40
41	C !LOCAL VARIABLES:
42	C fn :: character buffer for creating filename
43	C prec :: precision of pickup files
44	C k :: loop index
45	INTEGER i, prec
46	CHARACTER*(MAX_LEN_FNAM) fn
47
48	DO i = 1,MAX_LEN_FNAM
49	fn(i:i) = ' '
50	ENDDO
51	WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
52	prec = precFloat64
53
54	C#ifdef ALLOW_MDSIO
55	IF ( fizhi_mdsio_read_pickup ) THEN
56
57	C Read fields as consecutive records
58	CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)
59	CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)
60	CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)
61	CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)
62
63	ENDIF
64
65	#ifdef ALLOW_MNC
66	IF ( useMNC .and. fizhi_mnc_read_pickup ) THEN
67
68	_BEGIN_MASTER(myThid)
69
70	WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
71
72	CALL MNC_CW_RL_R('D',fn,0,0,'uphy',uphy,myThid)
73	CALL MNC_CW_RL_R('D',fn,0,0,'vphy',vphy,myThid)
74	CALL MNC_CW_RL_R('D',fn,0,0,'thphy',thphy,myThid)
75	CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)
76
77
78	_END_MASTER()
79
80	ENDIF
81	#endif
82
83	RETURN
84	END