/[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.8 by jmc, Tue Feb 5 15:31:19 2008 UTC revision 1.9 by jmc, Tue Nov 18 21:41:07 2008 UTC
# Line 65  C =============== Line 65  C ===============
65        INTEGER kd, kd0, ksgn, kStore        INTEGER kd, kd0, ksgn, kStore
66        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
67        integer i,offset,Lena,newindx,jindx        integer i,offset,Lena,newindx,jindx
68        _RL undef,getcon        _RL undef
69        INTEGER iSp, ndId, j,l        INTEGER iSp, ndId, j,l
70        INTEGER region2fill(0:nRegions)        INTEGER region2fill(0:nRegions)
71        _RL     scaleFact        _RL     scaleFact
72        _RL     gridField(sNx*sNy,nlevs), gridFrac(sNx*sNy)        _RL     gridField(sNx*sNy,nlevs), gridFrac(sNx*sNy)
73    
74    #ifdef ALLOW_FIZHI
75          _RL   getcon
76          EXTERNAL getcon
77    #endif
78    
79  C Run through list of active diagnostics to make sure  C Run through list of active diagnostics to make sure
80  C we are trying to fill a valid diagnostic  C we are trying to fill a valid diagnostic
81    
82          undef = UNSET_RL
83        ndiagnum = 0        ndiagnum = 0
84        ipointer = 0        ipointer = 0
85        DO n=1,nlists        DO n=1,nlists
# Line 143  C-      Check for consistency with Nb of Line 149  C-      Check for consistency with Nb of
149         DO k = kFirst,kLast         DO k = kFirst,kLast
150          kd = kd0 + ksgn*k          kd = kd0 + ksgn*k
151          if( check ) then          if( check ) then
152    #ifdef ALLOW_FIZHI
153           undef = getcon('UNDEF')           undef = getcon('UNDEF')
154    #endif
155           do i= 1,Lena           do i= 1,Lena
156            jindx = 1 + int((indx(i+offset-1)-1)/sNx)            jindx = 1 + int((indx(i+offset-1)-1)/sNx)
157            newindx = indx(i+offset-1)+(jindx-1)*2*Olx            newindx = indx(i+offset-1)+(jindx-1)*2*Olx
# Line 228  C-    Fill local array with grid-space f Line 236  C-    Fill local array with grid-space f
236           gridField(i,k)= 0.           gridField(i,k)= 0.
237          ENDDO          ENDDO
238          if( check ) then          if( check ) then
239    #ifdef ALLOW_FIZHI
240           undef = getcon('UNDEF')           undef = getcon('UNDEF')
241    #endif
242           do i= 1,Lena           do i= 1,Lena
243            newindx = indx(i+offset-1)            newindx = indx(i+offset-1)
244            if(gridField(newindx,k).eq.undef            if(gridField(newindx,k).eq.undef
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22