/[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.3 by molod, Fri Feb 18 16:20:26 2005 UTC revision 1.5 by molod, Wed Feb 23 16:01:54 2005 UTC
# Line 149  C-      Check for consistency with Nb of Line 149  C-      Check for consistency with Nb of
149           undef = getcon('UNDEF')           undef = getcon('UNDEF')
150           do i= 1,Lena           do i= 1,Lena
151            jindx = 1 + int((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
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
155             qdiag(newindx,1,kd,bi,bj) = undef             qdiag(newindx,1,kd,bi,bj) = undef
156            else            else
157             qdiag(newindx,1,kd,bi,bj)=qdiag(newindx,1,kd,bi,bj)+             qdiag(newindx,1,kd,bi,bj)=qdiag(newindx,1,kd,bi,bj)+
158       .                                             field(i,k)*chfr(i)       .                                                field(i,k)*chfr(i)
159            endif            endif
160           enddo           enddo
161          else          else
162           do i= 1,Lena           do i= 1,Lena
163            jindx = 1 + int((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
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)
167           enddo           enddo
168          endif          endif
169         ENDDO         ENDDO
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22