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

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

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

revision 1.12 by heimbach, Thu Jun 16 15:31:50 2005 UTC revision 1.13 by heimbach, Tue Nov 1 04:09:46 2005 UTC
# Line 170  c--   Only the master thread will do I/O Line 170  c--   Only the master thread will do I/O
170       &     Nr, globmsk, 1, mythid)       &     Nr, globmsk, 1, mythid)
171    
172        do irec = 1, ncvarrecs(ivartype)        do irec = 1, ncvarrecs(ivartype)
173    #ifndef ALLOW_ADMTLM
174           read(cunit) filencvarindex(ivartype)           read(cunit) filencvarindex(ivartype)
175           if (filencvarindex(ivartype) .NE. ncvarindex(ivartype))           if (filencvarindex(ivartype) .NE. ncvarindex(ivartype))
176       &        then       &        then
# Line 179  c--   Only the master thread will do I/O Line 180  c--   Only the master thread will do I/O
180           endif           endif
181           read(cunit) filej           read(cunit) filej
182           read(cunit) filei           read(cunit) filei
183    #endif /* ALLOW_ADMTLM */
184           do k = 1, Nr           do k = 1, Nr
185           irectrue = (irec-1)*nr + k           irectrue = (irec-1)*nr + k
186              if ( doZscaleUnpack ) then              if ( doZscaleUnpack ) then
# Line 189  cph               delZnorm = SQRT(delR(1 Line 191  cph               delZnorm = SQRT(delR(1
191              endif              endif
192              cbuffindex = nwetglobal(k)              cbuffindex = nwetglobal(k)
193              if ( cbuffindex .gt. 0 ) then              if ( cbuffindex .gt. 0 ) then
194    #ifndef ALLOW_ADMTLM
195                 read(cunit) filencbuffindex                 read(cunit) filencbuffindex
196                 if (filencbuffindex .NE. cbuffindex) then                 if (filencbuffindex .NE. cbuffindex) then
197                    print *, 'WARNING: wrong cbuffindex ',                    print *, 'WARNING: wrong cbuffindex ',
# Line 201  cph               delZnorm = SQRT(delR(1 Line 204  cph               delZnorm = SQRT(delR(1
204       &                 filek, k       &                 filek, k
205                    STOP 'in S/R ctrl_unpack'                    STOP 'in S/R ctrl_unpack'
206                 endif                 endif
207    #endif /* ALLOW_ADMTLM */
208                 read(cunit) (cbuff(ii), ii=1,cbuffindex)                 read(cunit) (cbuff(ii), ii=1,cbuffindex)
209              endif              endif
210    c
211              cbuffindex = 0              cbuffindex = 0
212              do jp = 1,nPy              do jp = 1,nPy
213               do bj = jtlo,jthi               do bj = jtlo,jthi
# Line 216  cph               delZnorm = SQRT(delR(1 Line 221  cph               delZnorm = SQRT(delR(1
221  cph(  cph(
222                       globfldtmp2(i,bi,ip,j,bj,jp) = cbuff(cbuffindex)                       globfldtmp2(i,bi,ip,j,bj,jp) = cbuff(cbuffindex)
223  cph)  cph)
224    #ifdef ALLOW_ADMTLM
225                         nveccount = nveccount + 1
226                         phtmpadmtlm(nveccount) = cbuff(cbuffindex)
227    #endif
228  #ifdef ALLOW_NONDIMENSIONAL_CONTROL_IO  #ifdef ALLOW_NONDIMENSIONAL_CONTROL_IO
229                       if ( lxxadxx ) then                       if ( lxxadxx ) then
230                          globfld3d(i,bi,ip,j,bj,jp,k) = delZnorm                          globfld3d(i,bi,ip,j,bj,jp,k) = delZnorm
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22