/[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.11 by heimbach, Tue May 15 01:28:57 2007 UTC revision 1.12 by heimbach, Wed May 16 00:03:37 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( nsx,npx,sny,nsy,npy,nr )        _RL globfld3d(1-oLx:sNx+oLx,1-oLy:sNy+oLy,nr,nSx,nSy)
55        _RL     globfld3d( snx,nsx,npx,sny,nsy,npy,nr )        _RL globfldyz(1-oLy:sNy+oLy,nr,nSx,nSy)
56    
57    
58  c     == external ==  c     == external ==
59    
# Line 62  c     == end of interface == Line 63  c     == end of interface ==
63        jthi = nsy        jthi = nsy
64        itlo = 1        itlo = 1
65        ithi = nsx        ithi = nsx
66        jmin = 1        jmin = 1-oLy
67        jmax = sny        jmax = sny+oLy
68        imin = 1        imin = 1-oLx
69        imax = snx        imax = snx+oLx
70    
71  c     Initialise temporary file  c     Initialise temporary file
72        do k = 1,nr        do bj = jtlo,jthi
73           do jp = 1,nPy           do bi = itlo,ithi
74              do bj = jtlo,jthi              do k = 1,nr
75                 do j = jmin,jmax                 do j = jmin,jmax
76                    do ip = 1,nPx                    globfldyz(j,k,bi,bj) = 0. _d 0
                      do bi = itlo,ithi  
                         globfldyz(bi,ip,j,bj,jp,k) = 0. _d 0  
                      enddo  
                   enddo  
77                 enddo                 enddo
78              enddo              enddo
79           enddo           enddo
80        enddo        enddo
81  c     Initialise temporary file  c     Initialise temporary file
82        do k = 1,nr        do bj = jtlo,jthi
83           do jp = 1,nPy           do bi = itlo,ithi
84              do bj = jtlo,jthi              do k = 1,nr
85                 do j = jmin,jmax                 do j = jmin,jmax
86                    do ip = 1,nPx                    do i = imin,imax
87                       do bi = itlo,ithi                       globfld3d(i,j,k,bi,bj) = 0. _d 0
                         do i = imin,imax  
                            globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0  
                         enddo  
                      enddo  
88                    enddo                    enddo
89                 enddo                 enddo
90              enddo              enddo
# Line 103  c--   Only the master thread will do I/O Line 96  c--   Only the master thread will do I/O
96    
97        nrec_nl=int(ncvarrecs(ivartype)/Nx)        nrec_nl=int(ncvarrecs(ivartype)/Nx)
98        do irec = 1, nrec_nl        do irec = 1, nrec_nl
99           call MDSWRITEFIELD_3D_GL( fname, ctrlprec, 'RL',          call MDSWRITEFIELD( fname, ctrlprec,.FALSE., 'RL',
100       &        Nr, globfld3d,       &        Nr, globfld3d,
101       &        irec,   optimcycle,  mythid)       &        irec,   optimcycle,  mythid)
102    
103        enddo        enddo
104    
105        do irec = nrec_nl*Nx+1, ncvarrecs(ivartype)  CMM( addone for good measure
106           call MDSWRITEFIELD_YZ_GL( fname, ctrlprec, 'RL',  CMM      do irec = nrec_nl*Nx+1, ncvarrecs(ivartype)
107          do irec = nrec_nl*Nx+1, ncvarrecs(ivartype)+1
108    CMM)
109            call MDSWRITEFIELDYZ( fname, ctrlprec,.FALSE., 'RL',
110       &        Nr, globfldyz,       &        Nr, globfldyz,
111       &        irec,   optimcycle,  mythid)       &        irec,   optimcycle,  mythid)
112        enddo        enddo
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22