/[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.10 by heimbach, Mon May 14 23:06:13 2007 UTC revision 1.11 by heimbach, Tue May 15 01:28:57 2007 UTC
# Line 51  c     == local variables == Line 51  c     == local variables ==
51        integer imin,imax        integer imin,imax
52        integer irec,nrec_nl        integer irec,nrec_nl
53    
54        _RL     globfldyz( snx,nr,snx,nsx )        _RL     globfldyz( nsx,npx,sny,nsy,npy,nr )
55        _RL     globfld3d( snx,sny,nr,nsx,nsy )        _RL     globfld3d( snx,nsx,npx,sny,nsy,npy,nr )
   
56    
57  c     == external ==  c     == external ==
58    
# Line 70  c     == end of interface == Line 69  c     == end of interface ==
69    
70  c     Initialise temporary file  c     Initialise temporary file
71        do k = 1,nr        do k = 1,nr
72             do jp = 1,nPy
73              do bj = jtlo,jthi              do bj = jtlo,jthi
74                 do j = jmin,jmax                 do j = jmin,jmax
75                      do ip = 1,nPx
76                       do bi = itlo,ithi                       do bi = itlo,ithi
77                          globfldyz(j,k,bi,bj) = 0. _d 0                          globfldyz(bi,ip,j,bj,jp,k) = 0. _d 0
78                       enddo                       enddo
79                    enddo                    enddo
80                   enddo
81              enddo              enddo
82             enddo
83        enddo        enddo
84  c     Initialise temporary file  c     Initialise temporary file
85        do k = 1,nr        do k = 1,nr
86             do jp = 1,nPy
87              do bj = jtlo,jthi              do bj = jtlo,jthi
88                 do j = jmin,jmax                 do j = jmin,jmax
89                      do ip = 1,nPx
90                       do bi = itlo,ithi                       do bi = itlo,ithi
91                          do i = imin,imax                          do i = imin,imax
92                             globfld3d(i,j,k,bi,bj) = 0. _d 0                             globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0
93                          enddo                          enddo
94                       enddo                       enddo
95                      enddo
96                 enddo                 enddo
97              enddo              enddo
98             enddo
99        enddo        enddo
100    
101  c--   Only the master thread will do I/O.  c--   Only the master thread will do I/O.
# Line 96  c--   Only the master thread will do I/O Line 103  c--   Only the master thread will do I/O
103    
104        nrec_nl=int(ncvarrecs(ivartype)/Nx)        nrec_nl=int(ncvarrecs(ivartype)/Nx)
105        do irec = 1, nrec_nl        do irec = 1, nrec_nl
106                    call MDSWRITEFIELD( fname, ctrlprec,.FALSE., 'RL',           call MDSWRITEFIELD_3D_GL( fname, ctrlprec, 'RL',
107       &        Nr, globfld3d,       &        Nr, globfld3d,
108       &        irec,   optimcycle,  mythid)       &        irec,   optimcycle,  mythid)
   
109        enddo        enddo
110    
111  CMM( addone for good measure        do irec = nrec_nl*Nx+1, ncvarrecs(ivartype)
112  CMM      do irec = nrec_nl*Nx+1, ncvarrecs(ivartype)           call MDSWRITEFIELD_YZ_GL( fname, ctrlprec, 'RL',
       do irec = nrec_nl*Nx+1, ncvarrecs(ivartype)+1  
 CMM)  
          call MDSWRITEFIELDYZ( fname, ctrlprec,.FALSE., 'RL',  
113       &        Nr, globfldyz,       &        Nr, globfldyz,
114       &        irec,   optimcycle,  mythid)       &        irec,   optimcycle,  mythid)
115        enddo        enddo
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22