-
Notifications
You must be signed in to change notification settings - Fork 237
/
dic_read_co2_pickup.F
85 lines (73 loc) · 2.28 KB
/
dic_read_co2_pickup.F
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
#include "DIC_OPTIONS.h"
SUBROUTINE DIC_READ_CO2_PICKUP(
I myIter, myThid )
IMPLICIT NONE
C === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DIC_VARS.h"
#include "DIC_ATMOS.h"
C == Routine arguments ==
C myThid :: my Thread Id number
INTEGER myIter
INTEGER myThid
#ifdef ALLOW_DIC
C !FUNCTIONS
INTEGER MDS_RECLEN
EXTERNAL MDS_RECLEN
C !LOCAL VARIABLES:
C == Local variables ==
INTEGER fp, ioUnit
CHARACTER*(10) suff
CHARACTER*(MAX_LEN_FNAM) fn, filNam
CHARACTER*(MAX_LEN_MBUF) msgBuf
LOGICAL useCurrentDir, fileExist
REAL*8 tmpFld(2)
INTEGER length_of_rec
CEOP
IF ( dic_int1.EQ.3 ) THEN
IF (pickupSuff.EQ.' ') THEN
IF ( rwSuffixType.EQ.0 ) THEN
WRITE(fn,'(A,I10.10)') 'pickup_dic_co2atm.', myIter
ELSE
CALL RW_GET_SUFFIX( suff, startTime, myIter, myThid )
WRITE(fn,'(A,A)') 'pickup_dic_co2atm.', suff
ENDIF
ELSE
WRITE(fn,'(A,A10)') 'pickup_dic_co2atm.', pickupSuff
ENDIF
fp = precFloat64
C-- First check if pickup file exist
#ifdef ALLOW_MDSIO
useCurrentDir = .FALSE.
CALL MDS_CHECK4FILE(
I fn, '.data', 'DIC_READ_CO2_PICKUP',
O filNam, fileExist,
I useCurrentDir, myThid )
#else
STOP 'ABNORMAL END: S/R DIC_READ_CO2_PICKUP: Needs MDSIO pkg'
#endif
IF ( .NOT.fileExist ) THEN
WRITE(msgBuf,'(4A)') 'DIC_READ_CO2_PICKUP: ',
& 'with dic_int1=3, needs CO2-atmos pickup to restart'
CALL PRINT_ERROR( msgBuf, myThid )
STOP 'ABNORMAL END: S/R DIC_READ_CO2_PICKUP'
ENDIF
_BEGIN_MASTER(myThid)
CALL MDSFINDUNIT( ioUnit, myThid )
length_of_rec = MDS_RECLEN( fp, 2, myThid )
OPEN( ioUnit, file=filNam, status='old',
& access='direct', recl=length_of_rec )
READ(ioUnit,rec=1) tmpFld
#ifdef _BYTESWAPIO
CALL MDS_BYTESWAPR8( 2, tmpFld )
#endif
total_atmos_carbon = tmpFld(1)
atpco2 = tmpFld(2)
_END_MASTER(myThid)
c _BARRIER
ENDIF
#endif /* ALLOW_DIC */
RETURN
END