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

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

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

revision 1.8 by heimbach, Tue Nov 16 05:42:12 2004 UTC revision 1.9 by heimbach, Mon May 14 23:06:13 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     globfld2d( snx,nsx,npx,sny,nsy,npy )        _RL     globfld2d( snx,sny,nsx,nsy )
55        _RL     globfld3d( snx,nsx,npx,sny,nsy,npy,nr )        _RL     globfld3d( snx,sny,nr,nsx,nsy )
56    
57  c     == external ==  c     == external ==
58    
# Line 68  c     == end of interface == Line 68  c     == end of interface ==
68        imax = snx        imax = snx
69    
70  c     Initialise temporary file  c     Initialise temporary file
       do jp = 1,nPy  
71           do bj = jtlo,jthi           do bj = jtlo,jthi
72              do j = jmin,jmax              do j = jmin,jmax
                do ip = 1,nPx  
73                    do bi = itlo,ithi                    do bi = itlo,ithi
74                       do i = imin,imax                       do i = imin,imax
75                          globfld2d(i,bi,ip,j,bj,jp) = 0. _d 0                          globfld2d(i,j,bi,bj) = 0. _d 0
76                       enddo                       enddo
                   enddo  
77                 enddo                 enddo
78              enddo              enddo
          enddo  
79        enddo        enddo
80  c     Initialise temporary file  c     Initialise temporary file
81        do k = 1,nr        do k = 1,nr
          do jp = 1,nPy  
82              do bj = jtlo,jthi              do bj = jtlo,jthi
83                 do j = jmin,jmax                 do j = jmin,jmax
                   do ip = 1,nPx  
84                       do bi = itlo,ithi                       do bi = itlo,ithi
85                          do i = imin,imax                          do i = imin,imax
86                             globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0                             globfld3d(i,j,k,bi,bj) = 0. _d 0
87                          enddo                          enddo
88                       enddo                       enddo
                   enddo  
89                 enddo                 enddo
90              enddo              enddo
          enddo  
91        enddo        enddo
92    
93  c--   Only the master thread will do I/O.  c--   Only the master thread will do I/O.
# Line 103  c--   Only the master thread will do I/O Line 95  c--   Only the master thread will do I/O
95    
96        nrec_nl=int(ncvarrecs(ivartype)/Nr)        nrec_nl=int(ncvarrecs(ivartype)/Nr)
97        do irec = 1,  nrec_nl        do irec = 1,  nrec_nl
98           call MDSWRITEFIELD_3D_GL( fname, ctrlprec, 'RL',           call MDSWRITEFIELD( fname, ctrlprec, .FALSE., 'RL',
99       &        Nr, globfld3d,       &        Nr, globfld3d,
100       &        irec,   optimcycle,  mythid)       &        irec,   optimcycle,  mythid)
101        enddo        enddo
102    
103        do irec = nrec_nl*Nr+1, ncvarrecs(ivartype)  CMM( add one for good measure
104           call MDSWRITEFIELD_2D_GL( fname, ctrlprec, 'RL',  CMM      do irec = nrec_nl*Nr+1, ncvarrecs(ivartype)
105          do irec = nrec_nl*Nr+1, ncvarrecs(ivartype)+1
106    CMM)
107             call MDSWRITEFIELD( fname, ctrlprec, .FALSE., 'RL',
108       &        1, globfld2d,       &        1, globfld2d,
109       &        irec,   optimcycle,  mythid)       &        irec,   optimcycle,  mythid)
110        enddo        enddo
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22