/[MITgcm]/MITgcm/pkg/diagnostics/diag_vegtile_fill.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diag_vegtile_fill.F

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

revision 1.2 by molod, Thu Feb 17 22:40:16 2005 UTC revision 1.3 by molod, Fri Feb 18 16:20:26 2005 UTC
# Line 148  C-      Check for consistency with Nb of Line 148  C-      Check for consistency with Nb of
148          if( check ) then          if( check ) then
149           undef = getcon('UNDEF')           undef = getcon('UNDEF')
150           do i= 1,Lena           do i= 1,Lena
151            jindx = 1.0001 + (indx(i+offset-1)-1/sNx)            jindx = 1 + int((indx(i+offset-1)-1)/sNx)
152            newindx = indx(i+offset-1)+(jindx-1)*2*Olx+(sNx+2*Olx)*Oly+Olx            newindx = indx(i+offset-1)+(jindx-1)*2*Olx+(sNx+2*Olx)*Oly+Olx
153            if(qdiag(newindx,1,kd,bi,bj).eq.undef            if(qdiag(newindx,1,kd,bi,bj).eq.undef
154       .                                  .or.field(i,k).eq.undef)then       .                                  .or.field(i,k).eq.undef)then
# Line 160  C-      Check for consistency with Nb of Line 160  C-      Check for consistency with Nb of
160           enddo           enddo
161          else          else
162           do i= 1,Lena           do i= 1,Lena
163            jindx = 1.0001 + (indx(i+offset-1)-1/sNx)            jindx = 1 + int((indx(i+offset-1)-1)/sNx)
164            newindx = indx(i+offset-1)+(jindx-1)*2*Olx+(sNx+2*Olx)*Oly+Olx            newindx = indx(i+offset-1)+(jindx-1)*2*Olx+(sNx+2*Olx)*Oly+Olx
165            qdiag(newindx,1,kd,bi,bj)=qdiag(newindx,1,kd,bi,bj)+            qdiag(newindx,1,kd,bi,bj)=qdiag(newindx,1,kd,bi,bj)+
166       .                                             field(i,k)*chfr(i)       .                                             field(i,k)*chfr(i)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22