/[MITgcm]/MITgcm/pkg/mnc/mnc_cw_readwrite.template
ViewVC logotype

Diff of /MITgcm/pkg/mnc/mnc_cw_readwrite.template

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

revision 1.3 by edhill, Thu Feb 5 00:13:47 2004 UTC revision 1.4 by edhill, Thu Feb 5 05:42:07 2004 UTC
# Line 79  C Line 79  C
79        DO lbj = bjs,bje        DO lbj = bjs,bje
80          DO lbi = bis,bie          DO lbi = bis,bie
81    
       write(*,*) 'lbi,lbj = ',lbi,lbj, '  ', vtype(nvf:nvl)  
   
82  C         Create the file name  C         Create the file name
83            CALL MNC_CW_GET_TILE_NUM(myThid, lbi,lbj, uniq_tnum)            CALL MNC_CW_GET_TILE_NUM(myThid, lbi,lbj, uniq_tnum)
84            fname(1:MNC_MAX_CHAR) = mnc_blank_name(1:MNC_MAX_CHAR)            fname(1:MNC_MAX_CHAR) = mnc_blank_name(1:MNC_MAX_CHAR)
# Line 141  C         Ensure that the variable is de Line 139  C         Ensure that the variable is de
139  C         Check for bi,bj indicies  C         Check for bi,bj indicies
140            bidim = mnc_cw_vbij(1,indv)            bidim = mnc_cw_vbij(1,indv)
141            bjdim = mnc_cw_vbij(2,indv)            bjdim = mnc_cw_vbij(2,indv)
142    CEH3      write(*,*) 'bidim,bjdim = ', bidim,bjdim
143    
144  C         Set the dimensions for the in-memory array  C         Set the dimensions for the in-memory array
145            ndim = mnc_cw_ndim(igrid)            ndim = mnc_cw_ndim(igrid)
# Line 199  C               Use the current unlimite Line 198  C               Use the current unlimite
198              s(bjdim) = lbj              s(bjdim) = lbj
199              e(bjdim) = lbj              e(bjdim) = lbj
200            ENDIF            ENDIF
201  CEH3          DO i = 1,7  CEH3      DO i = 1,9
202  CEH3            write(*,*) 'i,s(i),e(i) = ', i,s(i),e(i)  CEH3        write(*,*) 'i,p(i),s(i),e(i) = ', i,p(i),s(i),e(i)
203  CEH3          ENDDO  CEH3      ENDDO
204    
205  C         Add the global attributes  C         Add the global attributes
206            CALL MNC_CW_SET_GATTR(myThid, fname, lbi,lbj, uniq_tnum)            CALL MNC_CW_SET_GATTR(myThid, fname, lbi,lbj, uniq_tnum)
# Line 224  C         Add the per-variable attribute Line 223  C         Add the per-variable attribute
223    
224  C         Write the variable one vector at a time  C         Write the variable one vector at a time
225            DO j7 = s(7),e(7)            DO j7 = s(7),e(7)
226              k7 = (j7 - s(7))*p(6)              k7 = (j7 - 1)*p(6)
227              vstart(7) = udo(7) + j7 - s(7) + 1              vstart(7) = udo(7) + j7 - s(7) + 1
228              vcount(7) = 1              vcount(7) = 1
229              DO j6 = s(6),e(6)              DO j6 = s(6),e(6)
230                k6 = (j6 - s(6))*p(5) + k7                k6 = (j6 - 1)*p(5) + k7
231                vstart(6) = udo(6) + j6 - s(6) + 1                vstart(6) = udo(6) + j6 - s(6) + 1
232                vcount(6) = 1                vcount(6) = 1
233                DO j5 = s(5),e(5)                DO j5 = s(5),e(5)
234                  k5 = (j5 - s(5))*p(4) + k6                  k5 = (j5 - 1)*p(4) + k6
235                  vstart(5) = udo(5) + j5 - s(5) + 1                  vstart(5) = udo(5) + j5 - s(5) + 1
236                  vcount(5) = 1                  vcount(5) = 1
237                  DO j4 = s(4),e(4)                  DO j4 = s(4),e(4)
238                    k4 = (j4 - s(4))*p(3) + k5                    k4 = (j4 - 1)*p(3) + k5
239                    vstart(4) = udo(4) + j4 - s(4) + 1                    vstart(4) = udo(4) + j4 - s(4) + 1
240                    vcount(4) = 1                    vcount(4) = 1
241                    DO j3 = s(3),e(3)                    DO j3 = s(3),e(3)
242                      k3 = (j3 - s(3))*p(2) + k4                      k3 = (j3 - 1)*p(2) + k4
243                      vstart(3) = udo(3) + j3 - s(3) + 1                      vstart(3) = udo(3) + j3 - s(3) + 1
244                      vcount(3) = 1                      vcount(3) = 1
245                      DO j2 = s(2),e(2)                      DO j2 = s(2),e(2)
246                        k2 = (j2 - s(2))*p(1) + k3                        k2 = (j2 - 1)*p(1) + k3
247                        vstart(2) = udo(2) + j2 - s(2) + 1                        vstart(2) = udo(2) + j2 - s(2) + 1
248                        vcount(2) = 1                        vcount(2) = 1
249    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22