/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/dic_read_pickup.F
ViewVC logotype

Contents of /MITgcm_contrib/darwin2/pkg/darwin/dic_read_pickup.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.2 - (show annotations) (download)
Wed Oct 9 17:14:37 2013 UTC (11 years, 9 months ago) by stephd
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt65w_20160512, ctrb_darwin2_ckpt65j_20150225, ctrb_darwin2_ckpt66g_20170424, ctrb_darwin2_ckpt66k_20171025, ctrb_darwin2_ckpt66n_20180118, ctrb_darwin2_ckpt65v_20160409, ctrb_darwin2_ckpt65s_20160114, ctrb_darwin2_ckpt65_20140718, ctrb_darwin2_ckpt66d_20170214, ctrb_darwin2_ckpt64r_20131210, ctrb_darwin2_ckpt65m_20150615, ctrb_darwin2_ckpt65q_20151118, ctrb_darwin2_ckpt65o_20150914, ctrb_darwin2_ckpt65p_20151023, ctrb_darwin2_ckpt65e_20140929, ctrb_darwin2_ckpt64o_20131024, ctrb_darwin2_ckpt64v_20140411, ctrb_darwin2_ckpt64z_20140711, ctrb_darwin2_ckpt65l_20150504, ctrb_darwin2_ckpt65z_20160929, ctrb_darwin2_ckpt65n_20150729, ctrb_darwin2_ckpt64y_20140622, ctrb_darwin2_ckpt65d_20140915, ctrb_darwin2_ckpt64t_20140202, ctrb_darwin2_ckpt66h_20170602, ctrb_darwin2_ckpt64s_20140105, ctrb_darwin2_ckpt64x_20140524, ctrb_darwin2_ckpt65x_20160612, ctrb_darwin2_ckpt66f_20170407, ctrb_darwin2_ckpt65g_20141120, ctrb_darwin2_ckpt65k_20150402, ctrb_darwin2_ckpt64w_20140502, ctrb_darwin2_ckpt66a_20161020, ctrb_darwin2_ckpt65f_20141014, ctrb_darwin2_ckpt66b_20161219, ctrb_darwin2_ckpt64u_20140308, ctrb_darwin2_ckpt65i_20150123, ctrb_darwin2_ckpt66j_20170815, ctrb_darwin2_ckpt65y_20160801, ctrb_darwin2_ckpt66c_20170121, ctrb_darwin2_ckpt65a_20140728, ctrb_darwin2_ckpt65b_20140812, ctrb_darwin2_ckpt65t_20160221, ctrb_darwin2_ckpt64p_20131118, ctrb_darwin2_ckpt66o_20180209, ctrb_darwin2_ckpt66e_20170314, ctrb_darwin2_ckpt64q_20131118, ctrb_darwin2_ckpt64p_20131024, ctrb_darwin2_ckpt65u_20160315, ctrb_darwin2_ckpt65r_20151221, ctrb_darwin2_ckpt66i_20170718, ctrb_darwin2_ckpt65c_20140830, ctrb_darwin2_ckpt66l_20171025, ctrb_darwin2_ckpt65h_20141217, ctrb_darwin2_ckpt66m_20171213, HEAD
Changes since 1.1: +8 -0 lines 
o changes so that pH and pCO2 can be calculated for full water column
  by defining pH_3D in DARWIN_OPTIONS.h. Includes 3D diags and pickup
1 #include "CPP_OPTIONS.h"
2 #include "PTRACERS_OPTIONS.h"
3 #include "DARWIN_OPTIONS.h"
4
5 #ifdef ALLOW_PTRACERS
6 #ifdef ALLOW_DARWIN
7
8 #ifdef ALLOW_CARBON
9
10
11 SUBROUTINE DIC_READ_PICKUP(
12 O pH_isLoaded,
13 I myIter, myThid )
14
15 IMPLICIT NONE
16 C === Global variables ===
17 #include "SIZE.h"
18 #include "EEPARAMS.h"
19 #include "PARAMS.h"
20 #include "DARWIN_FLUX.h"
21
22 C == Routine arguments ==
23 C myThid :: my Thread Id number
24 LOGICAL pH_isLoaded
25 INTEGER myIter
26 INTEGER myThid
27
28
29 C !FUNCTIONS
30
31 C !LOCAL VARIABLES:
32 C == Local variables ==
33 CHARACTER*(MAX_LEN_FNAM) fn, filNam
34 CHARACTER*(MAX_LEN_MBUF) msgBuf
35 LOGICAL useCurrentDir, fileExist
36 INTEGER fp, ioUnit
37 CEOP
38
39 pH_isLoaded =.FALSE.
40 ioUnit = errorMessageUnit
41
42 C-- pickup file name :
43 IF (pickupSuff.EQ.' ') THEN
44 WRITE(fn,'(A,I10.10)') 'pickup_dic.', myIter
45 ELSE
46 WRITE(fn,'(A,A10)') 'pickup_dic.', pickupSuff
47 ENDIF
48 fp = precFloat64
49
50 C-- First check if pickup file exist
51 #ifdef ALLOW_MDSIO
52 useCurrentDir = .FALSE.
53 CALL MDS_CHECK4FILE(
54 I fn, '.data', 'DIC_READ_PICKUP',
55 O filNam, fileExist,
56 I useCurrentDir, myThid )
57 #else
58 STOP 'ABNORMAL END: S/R DIC_READ_PICKUP: Needs MDSIO pkg'
59 #endif
60
61 IF ( fileExist ) THEN
62 C-- Read pickup file
63 #ifdef pH_3D
64 CALL READ_REC_3D_RL( fn, fp, Nr, pH, 1, myIter, myThid )
65 #else
66 CALL READ_REC_3D_RL( fn, fp, 1, pH, 1, myIter, myThid )
67 #endif
68 pH_isLoaded = .TRUE.
69
70 #ifdef pH_3D
71 _EXCH_XYZ_RL( pH, myThid )
72 #else
73 _EXCH_XY_RL( pH, myThid )
74 #endif
75 ELSE
76 pH_isLoaded = .FALSE.
77 IF ( pickupStrictlyMatch ) THEN
78 WRITE(msgBuf,'(4A)') 'DIC_READ_PICKUP: ',
79 & 'try with " pickupStrictlyMatch=.FALSE.,"',
80 & ' in file: "data", NameList: "PARM03"'
81 CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
82 STOP 'ABNORMAL END: S/R DIC_READ_PICKUP'
83 ELSE
84 WRITE(msgBuf,'(2A)') 'WARNING >> DIC_READ_PICKUP: ',
85 & 'will restart from approximated pH'
86 CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
87 ENDIF
88 ENDIF
89
90
91 RETURN
92 END
93 #endif /*ALLOW_CARBON*/
94
95 #endif /*DARWIN*/
96 #endif /*ALLOW_PTRACERS*/
97 c ==================================================================
  ViewVC Help
Powered by ViewVC 1.1.22