/[MITgcm]/MITgcm/model/src/ini_curvilinear_grid.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_curvilinear_grid.F

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

revision 1.5 by adcroft, Thu Dec 6 16:36:13 2001 UTC revision 1.6 by jmc, Sun Jul 13 18:44:51 2003 UTC
# Line 141  cs-   end block Line 141  cs-   end block
141  cs-   this block needed by cubed sphere until we write more useful I/O routines  cs-   this block needed by cubed sphere until we write more useful I/O routines
142  C !!! _EXCH_ZUV_XY_R4(DXV, DYU, unSigned, myThid)  C !!! _EXCH_ZUV_XY_R4(DXV, DYU, unSigned, myThid)
143  cs!   this is not correct <= need paired exchange for dxv,dyu  cs!   this is not correct <= need paired exchange for dxv,dyu
144        CALL EXCH_Z_XY_RS(DXV,myThid)  cs    CALL EXCH_Z_XY_RS(DXV,myThid)
145        CALL EXCH_Z_XY_RS(DYU,myThid)  cs    CALL EXCH_Z_XY_RS(DYU,myThid)
146        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
147         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
         DXV(sNx+1,1,bi,bj)=DXV(1,1,bi,bj)  
         DXV(1,sNy+1,bi,bj)=DXV(1,1,bi,bj)  
         DYU(sNx+1,1,bi,bj)=DYU(1,1,bi,bj)  
         DYU(1,sNy+1,bi,bj)=DYU(1,1,bi,bj)  
148  cs! fix overlaps:  cs! fix overlaps:
149          DO j=1,sNy          DO j=1,sNy
150           DO i=1,Olx           DO i=1,Olx
151            DXV(1-i,j,bi,bj)=DXV(1+i,j,bi,bj)            DXV(1-i,j,bi,bj)=DXV(1+i,j,bi,bj)
152            DXV(sNx+i,j,bi,bj)=DXV(sNx-i,j,bi,bj)            DXV(sNx+i,j,bi,bj)=DXV(i,j,bi,bj)
153            DYU(1-i,j,bi,bj)=DYU(1+i,j,bi,bj)            DYU(1-i,j,bi,bj)=DYU(1+i,j,bi,bj)
154            DYU(sNx+i,j,bi,bj)=DYU(sNx-i,j,bi,bj)            DYU(sNx+i,j,bi,bj)=DYU(i,j,bi,bj)
155           ENDDO           ENDDO
156          ENDDO          ENDDO
157          DO j=1,Oly          DO j=1,Oly
158           DO i=1,sNx           DO i=1-Olx,sNx+Olx
159            DXV(i,1-j,bi,bj)=DXV(i,1+j,bi,bj)            DXV(i,1-j,bi,bj)=DXV(i,1+j,bi,bj)
160            DXV(i,sNy+j,bi,bj)=DXV(i,sNy-j,bi,bj)            DXV(i,sNy+j,bi,bj)=DXV(i,j,bi,bj)
161            DYU(i,1-j,bi,bj)=DYU(i,1+j,bi,bj)            DYU(i,1-j,bi,bj)=DYU(i,1+j,bi,bj)
162            DYU(i,sNy+j,bi,bj)=DYU(i,sNy-j,bi,bj)            DYU(i,sNy+j,bi,bj)=DYU(i,j,bi,bj)
163           ENDDO           ENDDO
164          ENDDO          ENDDO
165         ENDDO         ENDDO
166        ENDDO        ENDDO
167  cs-   end block  cs-   end block
 C !!! _EXCH_ZUV_XY_R4(DXV, DYU, unSigned, myThid)  
 cs!   this is not correct <= need paired exchange for dxv,dyu  
 cs    CALL EXCH_Z_XY_RS(DXV,myThid)  
 cs    CALL EXCH_Z_XY_RS(DYU,myThid)  
168    
169        CALL MDSREADFIELD('RAZ.bin',readBinaryPrec,'RS',1,RAZ,  1,myThid)        CALL MDSREADFIELD('RAZ.bin',readBinaryPrec,'RS',1,RAZ,  1,myThid)
170  cs-   this block needed by cubed sphere until we write more useful I/O routines  cs-   this block needed by cubed sphere until we write more useful I/O routines
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22