/[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.4 by adcroft, Wed Nov 28 16:45:16 2001 UTC revision 1.5 by adcroft, Thu Dec 6 16:36:13 2001 UTC
# Line 87  C !!! _EXCH_OUV_XY_R4(DXF, DYF, unSigned Line 87  C !!! _EXCH_OUV_XY_R4(DXF, DYF, unSigned
87  cs!   this is not correct! <= need paired exchange for DXF,DYF  cs!   this is not correct! <= need paired exchange for DXF,DYF
88        _EXCH_XY_R4(DXF,myThid)        _EXCH_XY_R4(DXF,myThid)
89        _EXCH_XY_R4(DYF,myThid)        _EXCH_XY_R4(DYF,myThid)
90    cs! fix overlaps:
91          DO bj = myByLo(myThid), myByHi(myThid)
92           DO bi = myBxLo(myThid), myBxHi(myThid)
93            DO j=1,sNy
94             DO i=1,Olx
95              DXF(1-i,j,bi,bj)=DXF(i,j,bi,bj)
96              DXF(sNx+i,j,bi,bj)=DXF(sNx+1-i,j,bi,bj)
97              DYF(1-i,j,bi,bj)=DYF(i,j,bi,bj)
98              DYF(sNx+i,j,bi,bj)=DYF(sNx+1-i,j,bi,bj)
99             ENDDO
100            ENDDO
101            DO j=1,Oly
102             DO i=1,sNx
103              DXF(i,1-j,bi,bj)=DXF(i,j,bi,bj)
104              DXF(i,sNy+j,bi,bj)=DXF(i,sNy+1-j,bi,bj)
105              DYF(i,1-j,bi,bj)=DYF(i,j,bi,bj)
106              DYF(i,sNy+j,bi,bj)=DYF(i,sNy+1-j,bi,bj)
107             ENDDO
108            ENDDO
109           ENDDO
110          ENDDO
111    cs
112    
113        CALL MDSREADFIELD('RA.bin',readBinaryPrec,'RS',1,RA,  1,myThid)        CALL MDSREADFIELD('RA.bin',readBinaryPrec,'RS',1,RA,  1,myThid)
114        _EXCH_XY_R4(RA,myThid )        _EXCH_XY_R4(RA,myThid )
# Line 127  cs!   this is not correct <= need paired Line 149  cs!   this is not correct <= need paired
149          DXV(1,sNy+1,bi,bj)=DXV(1,1,bi,bj)          DXV(1,sNy+1,bi,bj)=DXV(1,1,bi,bj)
150          DYU(sNx+1,1,bi,bj)=DYU(1,1,bi,bj)          DYU(sNx+1,1,bi,bj)=DYU(1,1,bi,bj)
151          DYU(1,sNy+1,bi,bj)=DYU(1,1,bi,bj)          DYU(1,sNy+1,bi,bj)=DYU(1,1,bi,bj)
152    cs! fix overlaps:
153            DO j=1,sNy
154             DO i=1,Olx
155              DXV(1-i,j,bi,bj)=DXV(1+i,j,bi,bj)
156              DXV(sNx+i,j,bi,bj)=DXV(sNx-i,j,bi,bj)
157              DYU(1-i,j,bi,bj)=DYU(1+i,j,bi,bj)
158              DYU(sNx+i,j,bi,bj)=DYU(sNx-i,j,bi,bj)
159             ENDDO
160            ENDDO
161            DO j=1,Oly
162             DO i=1,sNx
163              DXV(i,1-j,bi,bj)=DXV(i,1+j,bi,bj)
164              DXV(i,sNy+j,bi,bj)=DXV(i,sNy-j,bi,bj)
165              DYU(i,1-j,bi,bj)=DYU(i,1+j,bi,bj)
166              DYU(i,sNy+j,bi,bj)=DYU(i,sNy-j,bi,bj)
167             ENDDO
168            ENDDO
169         ENDDO         ENDDO
170        ENDDO        ENDDO
171  cs-   end block  cs-   end block
172  C !!! _EXCH_ZUV_XY_R4(DXV, DYU, unSigned, myThid)  C !!! _EXCH_ZUV_XY_R4(DXV, DYU, unSigned, myThid)
173  cs!   this is not correct <= need paired exchange for dxv,dyu  cs!   this is not correct <= need paired exchange for dxv,dyu
174        CALL EXCH_Z_XY_RS(DXV,myThid)  cs    CALL EXCH_Z_XY_RS(DXV,myThid)
175        CALL EXCH_Z_XY_RS(DYU,myThid)  cs    CALL EXCH_Z_XY_RS(DYU,myThid)
176    
177        CALL MDSREADFIELD('RAZ.bin',readBinaryPrec,'RS',1,RAZ,  1,myThid)        CALL MDSREADFIELD('RAZ.bin',readBinaryPrec,'RS',1,RAZ,  1,myThid)
178  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.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22