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

Annotation 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 - (hide 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 jahn 1.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 stephd 1.2 #ifdef pH_3D
64     CALL READ_REC_3D_RL( fn, fp, Nr, pH, 1, myIter, myThid )
65     #else
66 jahn 1.1 CALL READ_REC_3D_RL( fn, fp, 1, pH, 1, myIter, myThid )
67 stephd 1.2 #endif
68 jahn 1.1 pH_isLoaded = .TRUE.
69    
70 stephd 1.2 #ifdef pH_3D
71     _EXCH_XYZ_RL( pH, myThid )
72     #else
73 jahn 1.1 _EXCH_XY_RL( pH, myThid )
74 stephd 1.2 #endif
75 jahn 1.1 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