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

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

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

revision 1.24 by gforget, Tue Sep 4 14:58:53 2012 UTC revision 1.25 by heimbach, Fri Feb 1 19:25:32 2013 UTC
# Line 101  c     == end of interface == Line 101  c     == end of interface ==
101    
102  c     Initialise temporary file  c     Initialise temporary file
103        do k = 1,nr        do k = 1,nr
104           do jp = 1,nPy         do jp = 1,nPy
105              do bj = jtlo,jthi          do bj = jtlo,jthi
106                 do j = jmin,jmax           do j = jmin,jmax
107                    do ip = 1,nPx            do ip = 1,nPx
108                       do bi = itlo,ithi             do bi = itlo,ithi
109                          do i = imin,imax              do i = imin,imax
110                             globfld3d  (i,bi,ip,j,bj,jp,k) = 0. _d 0               globfld3d  (i,bi,ip,j,bj,jp,k) = 0. _d 0
111                             globmsk    (i,bi,ip,j,bj,jp,k) = 0. _d 0               globmsk    (i,bi,ip,j,bj,jp,k) = 0. _d 0
112                             globfldtmp2(i,bi,ip,j,bj,jp)   = 0.               globfldtmp2(i,bi,ip,j,bj,jp)   = 0. _d 0
113                             globfldtmp3(i,bi,ip,j,bj,jp)   = 0.               globfldtmp3(i,bi,ip,j,bj,jp)   = 0. _d 0
                         enddo  
                      enddo  
                   enddo  
                enddo  
114              enddo              enddo
115               enddo
116              enddo
117           enddo           enddo
118            enddo
119           enddo
120        enddo        enddo
121    
122  c--   Only the master thread will do I/O.  c--   Only the master thread will do I/O.
# Line 126  c--   Only the master thread will do I/O Line 126  c--   Only the master thread will do I/O
126           write(cfile2(1:80),'(80a)') ' '           write(cfile2(1:80),'(80a)') ' '
127           write(cfile3(1:80),'(80a)') ' '           write(cfile3(1:80),'(80a)') ' '
128           if ( lxxadxx ) then           if ( lxxadxx ) then
129              write(cfile2(1:80),'(a,I2.2,a,I4.4,a)')              write(cfile2(1:80),'(a,I3.3,a,I4.4,a)')
130       &           'diag_pack_nonout_ctrl_',       &           'diag_pack_nonout_ctrl_',
131       &           ivartype, '_', optimcycle, '.bin'       &           ivartype, '_', optimcycle, '.bin'
132              write(cfile3(1:80),'(a,I2.2,a,I4.4,a)')              write(cfile3(1:80),'(a,I3.3,a,I4.4,a)')
133       &           'diag_pack_dimout_ctrl_',       &           'diag_pack_dimout_ctrl_',
134       &           ivartype, '_', optimcycle, '.bin'       &           ivartype, '_', optimcycle, '.bin'
135           else           else
136              write(cfile2(1:80),'(a,I2.2,a,I4.4,a)')              write(cfile2(1:80),'(a,I3.3,a,I4.4,a)')
137       &           'diag_pack_nonout_grad_',       &           'diag_pack_nonout_grad_',
138       &           ivartype, '_', optimcycle, '.bin'       &           ivartype, '_', optimcycle, '.bin'
139              write(cfile3(1:80),'(a,I2.2,a,I4.4,a)')              write(cfile3(1:80),'(a,I3.3,a,I4.4,a)')
140       &           'diag_pack_dimout_grad_',       &           'diag_pack_dimout_grad_',
141       &           ivartype, '_', optimcycle, '.bin'       &           ivartype, '_', optimcycle, '.bin'
142           endif           endif
# Line 151  c--   Only the master thread will do I/O Line 151  c--   Only the master thread will do I/O
151        endif        endif
152    
153  #ifdef CTRL_PACK_PRECISE  #ifdef CTRL_PACK_PRECISE
154        il=ilnblnk( weighttype)        if (weighttype.NE.' ') then
155        write(weightname(1:80),'(80a)') ' '         il=ilnblnk( weighttype)
156        write(weightname(1:80),'(a)') weighttype(1:il)         write(weightname(1:80),'(80a)') ' '
157           write(weightname(1:80),'(a)') weighttype(1:il)
158        call MDSREADFIELD_3D_GL(         call MDSREADFIELD_3D_GL(
159       &     weightname, ctrlprec, 'RL',       &     weightname, ctrlprec, 'RL',
160       &     Nr, weightfld3d, 1, mythid)       &     Nr, weightfld3d, 1, mythid)
161          else
162           do k = 1,nr
163            do jp = 1,nPy
164             do bj = jtlo,jthi
165              do j = jmin,jmax
166               do ip = 1,nPx
167                do bi = itlo,ithi
168                 do i = imin,imax
169                  weightfld3d(i,bi,ip,j,bj,jp,k) = 1. _d 0
170                 enddo
171                enddo
172               enddo
173              enddo
174             enddo
175            enddo
176           enddo
177          endif
178  #endif  #endif
179    
180        call MDSREADFIELD_3D_GL(        call MDSREADFIELD_3D_GL(
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22