pkg/dic/dic_read_pickup.F

C $Header: /u/gcmpack/MITgcm/pkg/dic/dic_read_pickup.F,v 1.5 2008/04/10 00:49:25 jmc Exp $
C $Name:  $

#include "DIC_OPTIONS.h"

      SUBROUTINE DIC_READ_PICKUP(
     O                            pH_isLoaded,
     I                            myIter, myThid )

      IMPLICIT NONE
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DIC_VARS.h"

C     == Routine arguments ==
C     myThid  :: my Thread Id number
      LOGICAL pH_isLoaded
      INTEGER myIter
      INTEGER myThid

#ifdef ALLOW_DIC
#ifdef DIC_BIOTIC

C     !FUNCTIONS

C     !LOCAL VARIABLES:
C     == Local variables ==
      CHARACTER*(MAX_LEN_FNAM) fn, filNam
      CHARACTER*(MAX_LEN_MBUF) msgBuf
      LOGICAL fileExist
      INTEGER fp, ioUnit
CEOP

      pH_isLoaded =.FALSE.
      ioUnit = errorMessageUnit

C--   pickup file name :
      IF (pickupSuff.EQ.' ') THEN
        WRITE(fn,'(A,I10.10)') 'pickup_dic.', myIter
      ELSE
        WRITE(fn,'(A,A10)') 'pickup_dic.', pickupSuff
      ENDIF
      fp = precFloat64

C--   First check if pickup file exist
      CALL MDS_CHECK4FILE(
     I                     fn, '.data', 'DIC_READ_PICKUP',
     O                     filNam, fileExist,
     I                     myThid )

      IF ( fileExist ) THEN
C--   Read pickup file
        CALL READ_REC_3D_RL( fn, fp, 1, pH, 1, myIter, myThid )
        pH_isLoaded = .TRUE.

        _EXCH_XY_R8( pH, myThid )
      ELSE
        pH_isLoaded = .FALSE.
        IF ( pickupStrictlyMatch ) THEN
          WRITE(msgBuf,'(4A)') 'DIC_READ_PICKUP: ',
     &      'try with " pickupStrictlyMatch=.FALSE.,"',
     &      ' in file: "data", NameList: "PARM03"'
          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
          STOP 'ABNORMAL END: S/R DIC_READ_PICKUP'
        ELSE
          WRITE(msgBuf,'(2A)') 'WARNING >> DIC_READ_PICKUP: ',
     &      'will restart from approximated pH'
          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
        ENDIF
      ENDIF

#endif /*  DIC_BIOTIC  */
#endif /*  ALLOW_DIC  */

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/dic/dic_read_pickup.F,v 1.5 2008/04/10 00:49:25 jmc Exp $
2	C $Name: $
3
4	#include "DIC_OPTIONS.h"
5
6	SUBROUTINE DIC_READ_PICKUP(
7	O pH_isLoaded,
8	I myIter, myThid )
9
10	IMPLICIT NONE
11	C === Global variables ===
12	#include "SIZE.h"
13	#include "EEPARAMS.h"
14	#include "PARAMS.h"
15	#include "DIC_VARS.h"
16
17	C == Routine arguments ==
18	C myThid :: my Thread Id number
19	LOGICAL pH_isLoaded
20	INTEGER myIter
21	INTEGER myThid
22
23	#ifdef ALLOW_DIC
24	#ifdef DIC_BIOTIC
25
26	C !FUNCTIONS
27
28	C !LOCAL VARIABLES:
29	C == Local variables ==
30	CHARACTER*(MAX_LEN_FNAM) fn, filNam
31	CHARACTER*(MAX_LEN_MBUF) msgBuf
32	LOGICAL fileExist
33	INTEGER fp, ioUnit
34	CEOP
35
36	pH_isLoaded =.FALSE.
37	ioUnit = errorMessageUnit
38
39	C-- pickup file name :
40	IF (pickupSuff.EQ.' ') THEN
41	WRITE(fn,'(A,I10.10)') 'pickup_dic.', myIter
42	ELSE
43	WRITE(fn,'(A,A10)') 'pickup_dic.', pickupSuff
44	ENDIF
45	fp = precFloat64
46
47	C-- First check if pickup file exist
48	CALL MDS_CHECK4FILE(
49	I fn, '.data', 'DIC_READ_PICKUP',
50	O filNam, fileExist,
51	I myThid )
52
53	IF ( fileExist ) THEN
54	C-- Read pickup file
55	CALL READ_REC_3D_RL( fn, fp, 1, pH, 1, myIter, myThid )
56	pH_isLoaded = .TRUE.
57
58	_EXCH_XY_R8( pH, myThid )
59	ELSE
60	pH_isLoaded = .FALSE.
61	IF ( pickupStrictlyMatch ) THEN
62	WRITE(msgBuf,'(4A)') 'DIC_READ_PICKUP: ',
63	& 'try with " pickupStrictlyMatch=.FALSE.,"',
64	& ' in file: "data", NameList: "PARM03"'
65	CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
66	STOP 'ABNORMAL END: S/R DIC_READ_PICKUP'
67	ELSE
68	WRITE(msgBuf,'(2A)') 'WARNING >> DIC_READ_PICKUP: ',
69	& 'will restart from approximated pH'
70	CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
71	ENDIF
72	ENDIF
73
74	#endif /* DIC_BIOTIC */
75	#endif /* ALLOW_DIC */
76
77	RETURN
78	END