/[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.1 - (show annotations) (download)
Wed Apr 13 18:56:24 2011 UTC (14 years, 6 months ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt64k_20130723, ctrb_darwin2_ckpt63l_20120405, ctrb_darwin2_ckpt64h_20130528, ctrb_darwin2_ckpt62v_20110413, ctrb_darwin2_ckpt64m_20130820, ctrb_darwin2_ckpt64f_20130405, ctrb_darwin2_ckpt63f_20111201, ctrb_darwin2_ckpt64a_20121116, ctrb_darwin2_ckpt64n_20130826, ctrb_darwin2_ckpt62y_20110526, ctrb_darwin2_ckpt64i_20130622, ctrb_darwin2_ckpt62x_20110513, ctrb_darwin2_ckpt62w_20110426, ctrb_darwin2_ckpt63o_20120629, ctrb_darwin2_ckpt64e_20130305, ctrb_darwin2_ckpt63c_20111011, ctrb_darwin2_ckpt63i_20120124, ctrb_darwin2_ckpt63m_20120506, ctrb_darwin2_ckpt63s_20120908, ctrb_darwin2_ckpt63e_20111107, ctrb_darwin2_ckpt63b_20110830, ctrb_darwin2_ckpt63j_20120217, ctrb_darwin2_ckpt63r_20120817, ctrb_darwin2_ckpt64g_20130503, ctrb_darwin2_ckpt64l_20130806, ctrb_darwin2_ckpt63g_20111220, ctrb_darwin2_ckpt64c_20130120, ctrb_darwin2_ckpt63a_20110804, ctrb_darwin2_ckpt64j_20130704, ctrb_darwin2_ckpt63h_20111230, ctrb_darwin2_ckpt63p_20120707, ctrb_darwin2_ckpt63d_20111107, ctrb_darwin2_ckpt63q_20120731, ctrb_darwin2_ckpt63_20110728, ctrb_darwin2_ckpt64b_20121224, ctrb_darwin2_ckpt64d_20130219, ctrb_darwin2_ckpt64_20121012, ctrb_darwin2_baseline, ctrb_darwin2_ckpt63n_20120604, ctrb_darwin2_ckpt63k_20120317, ctrb_darwin2_ckpt62z_20110622 
darwin2 initial checkin
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 CALL READ_REC_3D_RL( fn, fp, 1, pH, 1, myIter, myThid )
64 pH_isLoaded = .TRUE.
65
66 _EXCH_XY_RL( pH, myThid )
67 ELSE
68 pH_isLoaded = .FALSE.
69 IF ( pickupStrictlyMatch ) THEN
70 WRITE(msgBuf,'(4A)') 'DIC_READ_PICKUP: ',
71 & 'try with " pickupStrictlyMatch=.FALSE.,"',
72 & ' in file: "data", NameList: "PARM03"'
73 CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
74 STOP 'ABNORMAL END: S/R DIC_READ_PICKUP'
75 ELSE
76 WRITE(msgBuf,'(2A)') 'WARNING >> DIC_READ_PICKUP: ',
77 & 'will restart from approximated pH'
78 CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
79 ENDIF
80 ENDIF
81
82
83 RETURN
84 END
85 #endif /*ALLOW_CARBON*/
86
87 #endif /*DARWIN*/
88 #endif /*ALLOW_PTRACERS*/
89 c ==================================================================
  ViewVC Help
Powered by ViewVC 1.1.22