/[MITgcm]/MITgcm/pkg/ptracers/ptracers_read_checkpoint.F
ViewVC logotype

Contents of /MITgcm/pkg/ptracers/ptracers_read_checkpoint.F

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

Revision 1.4 - (show annotations) (download)
Wed Jun 25 23:47:05 2003 UTC (21 years ago) by stephd
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint53f_post, checkpoint54a_pre, checkpoint53b_pre, checkpoint51l_post, checkpoint51j_post, checkpoint52l_pre, checkpoint52e_pre, hrcube4, checkpoint52n_post, checkpoint52j_post, checkpoint51o_pre, checkpoint52e_post, checkpoint51n_pre, checkpoint52d_pre, checkpoint53c_post, checkpoint53d_post, checkpoint51f_pre, branchpoint-genmake2, checkpoint52j_pre, checkpoint54a_post, branch-netcdf, checkpoint51r_post, checkpoint52b_pre, checkpoint51i_post, checkpoint51e_post, checkpoint51b_post, checkpoint51l_pre, checkpoint52m_post, checkpoint51c_post, checkpoint53a_post, hrcube5, checkpoint53b_post, checkpoint51o_post, checkpoint51q_post, checkpoint52l_post, checkpoint52k_post, checkpoint54, checkpoint53, checkpoint52, checkpoint52d_post, checkpoint51b_pre, checkpoint52a_post, checkpoint52b_post, checkpoint53g_post, checkpoint52f_post, checkpoint52c_post, checkpoint51h_pre, checkpoint51g_post, ecco_c52_e35, checkpoint51f_post, checkpoint52a_pre, checkpoint51d_post, checkpoint51m_post, checkpoint51t_post, checkpoint53d_pre, checkpoint52i_post, checkpoint51p_post, checkpoint51n_post, checkpoint51i_pre, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post
Branch point for: netcdf-sm0, branch-genmake2, branch-nonh, tg2-branch, checkpoint51n_branch
Changes since 1.3: +3 -3 lines 
fix bug in pickup file name
1 C $Header: /u/gcmpack/MITgcm/pkg/ptracers/ptracers_read_checkpoint.F,v 1.3 2003/06/12 18:41:46 jmc Exp $
2 C $Name: $
3
4 #include "PTRACERS_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: PTRACERS_READ_CHECKPOINT
8
9 C !INTERFACE: ==========================================================
10 SUBROUTINE PTRACERS_READ_CHECKPOINT( myIter,myThid )
11
12 C !DESCRIPTION:
13 C Reads current state of passive tracers from a pickup file
14
15 C !USES: ===============================================================
16 IMPLICIT NONE
17 #include "SIZE.h"
18 #include "EEPARAMS.h"
19 #include "PARAMS.h"
20 #include "PTRACERS.h"
21
22 C !INPUT PARAMETERS: ===================================================
23 C myIter :: time-step number
24 C suff :: filename suffix
25 C myThid :: thread number
26 INTEGER myIter
27 CHARACTER*(10) suff
28 INTEGER myThid
29
30 C !OUTPUT PARAMETERS: ==================================================
31 C none
32
33 #ifdef ALLOW_PTRACERS
34
35 C !LOCAL VARIABLES: ====================================================
36 C iTracer :: loop indices
37 C iRec :: record number
38 C fn :: character buffer for creating filename
39 C prec :: precision of pickup files
40 C fn :: suffix for pickup file (eg. 0000000010)
41 INTEGER iTracer,prec,iRec
42 CHARACTER*(MAX_LEN_MBUF) fn
43 CEOP
44
45 WRITE(fn,'(A,I10.10)') 'pickup_ptracers.',myIter
46 prec = precFloat64
47
48 C Write fields as consecutive records
49 DO iTracer=1,PTRACERS_num
50 iRec = iTracer
51 CALL MDSREADFIELD(fn,prec,'RL',Nr,
52 & pTracer(1-Olx,1-Oly,1,1,1,iTracer),iRec,myThid)
53 _EXCH_XYZ_R8(pTracer(1-Olx,1-Oly,1,1,1,iTracer),myThid)
54 ENDDO
55
56 C Write historical tendancies as consecutive records
57 DO iTracer=1,PTRACERS_num
58 iRec = iTracer + PTRACERS_num
59 CALL MDSREADFIELD(fn,prec,'RL',Nr,
60 & gPtr(1-Olx,1-Oly,1,1,1,iTracer),iRec,myThid)
61 _EXCH_XYZ_R8(gPtr(1-Olx,1-Oly,1,1,1,iTracer),myThid)
62 ENDDO
63 DO iTracer=1,PTRACERS_num
64 iRec = iTracer + PTRACERS_num*2
65 CALL MDSREADFIELD(fn,prec,'RL',Nr,
66 & gPtrNm1(1-Olx,1-Oly,1,1,1,iTracer),iRec,myThid)
67 _EXCH_XYZ_R8(gPtrNm1(1-Olx,1-Oly,1,1,1,iTracer),myThid)
68 ENDDO
69
70 #endif /* ALLOW_PTRACERS */
71
72 RETURN
73 END
  ViewVC Help
Powered by ViewVC 1.1.22