/[MITgcm]/MITgcm/pkg/ctrl/ctrl_set_globfld_yz.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_set_globfld_yz.F

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

revision 1.3 by heimbach, Mon Sep 16 18:11:58 2002 UTC revision 1.4 by heimbach, Tue Jun 24 16:07:07 2003 UTC
# Line 12  c     o initialise field Line 12  c     o initialise field
12  c  c
13  c     started: heimbach@mit.edu, 16-Aug-2001  c     started: heimbach@mit.edu, 16-Aug-2001
14  c  c
15    c     changed: heimbach@mit.edu 17-Jun-2003
16    c              merged Armin's changes to replace write of
17    c              nr * globfld2d by 1 * globfld3d
18    c              (ad hoc fix to speed up global I/O)
19    c
20  c     ==================================================================  c     ==================================================================
21    
22        implicit none        implicit none
# Line 49  c     == local variables == Line 54  c     == local variables ==
54        integer jtlo,jthi        integer jtlo,jthi
55        integer jmin,jmax        integer jmin,jmax
56        integer imin,imax        integer imin,imax
57        integer irec        integer irec,nrec_nl
58    
59        _RL     globfldyz( nsx,npx,sny,nsy,npy,nr )        _RL     globfldyz( nsx,npx,sny,nsy,npy,nr )
60          _RL     globfld3d( snx,nsx,npx,sny,nsy,npy,nr )
61    
62  c     == external ==  c     == external ==
63    
# Line 73  c     Initialise temporary file Line 79  c     Initialise temporary file
79                 do j = jmin,jmax                 do j = jmin,jmax
80                    do ip = 1,nPx                    do ip = 1,nPx
81                       do bi = itlo,ithi                       do bi = itlo,ithi
 #ifdef ALLOW_TANGENTLINEAR_RUN  
                         globfldyz(bi,ip,j,bj,jp,k) = 0. _d 0  
 #else  
82                          globfldyz(bi,ip,j,bj,jp,k) = 0. _d 0                          globfldyz(bi,ip,j,bj,jp,k) = 0. _d 0
83  #endif                       enddo
84                      enddo
85                   enddo
86                enddo
87             enddo
88          enddo
89    c     Initialise temporary file
90          do k = 1,nr
91             do jp = 1,nPy
92                do bj = jtlo,jthi
93                   do j = jmin,jmax
94                      do ip = 1,nPx
95                         do bi = itlo,ithi
96                            do i = imin,imax
97                               globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0
98                            enddo
99                       enddo                       enddo
100                    enddo                    enddo
101                 enddo                 enddo
# Line 92  c     Initialise temporary file Line 110  c     Initialise temporary file
110  c--   Only the master thread will do I/O.  c--   Only the master thread will do I/O.
111        _BEGIN_MASTER( mythid )        _BEGIN_MASTER( mythid )
112    
113        do irec = 1, ncvarrecs(ivartype)        nrec_nl=int(ncvarrecs(ivartype)/snx)
114          do irec = 1, nrec_nl
115             call MDSWRITEFIELD_3D_GL( fname, ctrlprec, 'RL',
116         &        Nr, globfld3d,
117         &        irec,   optimcycle,  mythid)
118          enddo
119    
120          do irec = nrec_nl*snx+1, ncvarrecs(ivartype)
121           call MDSWRITEFIELD_YZ_GL( fname, ctrlprec, 'RL',           call MDSWRITEFIELD_YZ_GL( fname, ctrlprec, 'RL',
122       &        Nr, globfldyz,       &        Nr, globfldyz,
123       &        irec,   optimcycle,  mythid)       &        irec,   optimcycle,  mythid)
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22