/[MITgcm]/MITgcm/pkg/mnc/mnc_cwrapper.F
ViewVC logotype

Diff of /MITgcm/pkg/mnc/mnc_cwrapper.F

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

revision 1.16 by edhill, Tue Jul 6 22:38:31 2004 UTC revision 1.17 by edhill, Fri Jul 9 02:41:52 2004 UTC
# Line 237  C     !ROUTINE: MNC_CW_INIT Line 237  C     !ROUTINE: MNC_CW_INIT
237    
238  C     !INTERFACE:  C     !INTERFACE:
239        SUBROUTINE MNC_CW_INIT(        SUBROUTINE MNC_CW_INIT(
240       I     sNx,sNy, OLx,OLy, nSx,nSy, nPx,nPy, Nr, NrPhys,       I     sNx,sNy, OLx,OLy, nSx,nSy, nPx,nPy, Nr,
241       I     myThid )       I     myThid )
242    
243  C     !DESCRIPTION:  C     !DESCRIPTION:
# Line 273  C     !USES: Line 273  C     !USES:
273  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
274        integer myThid        integer myThid
275        integer sNx,sNy, OLx,OLy, nSx,nSy, nPx,nPy, Nr        integer sNx,sNy, OLx,OLy, nSx,nSy, nPx,nPy, Nr
       integer NrPhys  
276  CEOP  CEOP
277    
278  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
279        integer CW_MAX_LOC        integer CW_MAX_LOC
280        parameter ( CW_MAX_LOC = 5 )        parameter ( CW_MAX_LOC = 5 )
281        integer i, ihorz,ihsub,ivert,itime,ihalo, is,ih, n,ntot        integer i, ihorz,ihsub,ivert,itime,ihalo, is,ih, n,ntot
282        integer ndim, ncomb, nvch, NrPh        integer ndim, ncomb, nvch
283        character*(MNC_MAX_CHAR) name        character*(MNC_MAX_CHAR) name
284        character*(MNC_MAX_CHAR) dn(CW_MAX_LOC)        character*(MNC_MAX_CHAR) dn(CW_MAX_LOC)
285        character*(5) horz_dat(CW_MAX_LOC), hsub_dat(CW_MAX_LOC),        character*(5) horz_dat(CW_MAX_LOC), hsub_dat(CW_MAX_LOC),
# Line 300  C     ......12345....12345....12345....1 Line 299  C     ......12345....12345....12345....1
299        data halo_dat /        data halo_dat /
300       &     'Hn   ', 'Hy   ', '--   ', '     ', '     '  /       &     'Hn   ', 'Hy   ', '--   ', '     ', '     '  /
301        data vert_dat /        data vert_dat /
302       &     '-    ', 'C    ', 'I    ', 'Phys ', 'PhysI'  /       &     '-    ', 'C    ', 'I    ', '     ', '     '  /
303        data time_dat /        data time_dat /
304       &     '-    ', 't    ', '     ', '     ', '     '  /       &     '-    ', 't    ', '     ', '     ', '     '  /
305    
       if (NrPhys .lt. 1) then  
         NrPh = Nr  
       else  
         NrPh = NrPhys  
       endif  
   
306        ncomb = 0        ncomb = 0
307        DO ihorz = 1,5        DO ihorz = 1,5
308          DO is = 1,3          DO is = 1,3
# Line 327  C           Loop just ONCE if the Horiz Line 320  C           Loop just ONCE if the Horiz
320                ENDIF                ENDIF
321              ENDIF              ENDIF
322                            
323              DO ivert = 1,5              DO ivert = 1,3
324                DO itime = 1,2                DO itime = 1,2
325                                    
326  C               horiz and hsub  C               horiz and hsub
# Line 431  C               Vertical dimension Line 424  C               Vertical dimension
424                    ib(ndim)  = 1                    ib(ndim)  = 1
425                    ie(ndim)  = Nr + 1                    ie(ndim)  = Nr + 1
426                  ENDIF                  ENDIF
                 IF (vert_dat(ivert)(1:5) .EQ. 'Phys ') THEN  
                   ndim = ndim + 1  
                   dn(ndim)(1:5) = 'Zphys'  
                   dim(ndim) = NrPh  
                   ib(ndim)  = 1  
                   ie(ndim)  = NrPh  
                 ENDIF  
                 IF (vert_dat(ivert)(1:5) .EQ. 'PhysI') THEN  
                   ndim = ndim + 1  
                   dn(ndim)(1:7) = 'Zphysm1'  
                   dim(ndim) = NrPh - 1  
                   ib(ndim)  = 1  
                   ie(ndim)  = NrPh - 1  
                 ENDIF  
427    
428  C               Time dimension  C               Time dimension
429                  IF (time_dat(itime)(1:1) .EQ. 't') THEN                  IF (time_dat(itime)(1:1) .EQ. 't') THEN
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22