/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/streamice_velmask_upd.F
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_velmask_upd.F

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

revision 1.1 by heimbach, Thu Mar 29 15:59:21 2012 UTC revision 1.3 by dgoldberg, Wed Jan 9 21:56:18 2013 UTC
# Line 31  C     === Global variables === Line 31  C     === Global variables ===
31        INTEGER i, j, bi, bj, ki, kj        INTEGER i, j, bi, bj, ki, kj
32        INTEGER maskFlag        INTEGER maskFlag
33    
34          _EXCH_XY_RL( H_streamice, myThid )
35          _EXCH_XY_RL( area_shelf_streamice, myThid )
36          _EXCH_XY_RL( STREAMICE_hmask, myThid )
37    
38        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
39         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
40          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
# Line 62  C     === Global variables === Line 66  C     === Global variables ===
66               IF (maskFlag.EQ.3) THEN               IF (maskFlag.EQ.3) THEN
67                DO kj=0,1                DO kj=0,1
68                 STREAMICE_umask(i+ki,j+kj,bi,bj) = 3.0                 STREAMICE_umask(i+ki,j+kj,bi,bj) = 3.0
69                 STREAMICE_vmask(i+ki,j+kj,bi,bj) = 0.0                 STREAMICE_vmask(i+ki,j+kj,bi,bj) = 3.0
70                ENDDO                ENDDO
71                STREAMICE_ufacemask(i+ki,j,bi,bj) = 3.0                STREAMICE_ufacemask(i+ki,j,bi,bj) = 3.0
72               ELSE IF (maskFlag.EQ.2) THEN               ELSE IF (maskFlag.EQ.2) THEN
73                DO kj=0,1                !DO kj=0,1
74                 STREAMICE_ufacemask(i+ki,j+kj,bi,bj) = 2.0                 STREAMICE_ufacemask(i+ki,j,bi,bj) = 2.0
75                ENDDO                !ENDDO
76               ELSE IF (maskFlag.EQ.4) THEN               ELSE IF (maskFlag.EQ.4) THEN
77                DO kj=0,1                DO kj=0,1
78                 STREAMICE_umask(i+ki,j+kj,bi,bj) = 0.0                 STREAMICE_umask(i+ki,j+kj,bi,bj) = 0.0
# Line 93  C     === Global variables === Line 97  C     === Global variables ===
97              IF (maskFlag.EQ.3) THEN              IF (maskFlag.EQ.3) THEN
98                DO ki=0,1                DO ki=0,1
99                 STREAMICE_vmask(i+ki,j+kj,bi,bj) = 3.0                 STREAMICE_vmask(i+ki,j+kj,bi,bj) = 3.0
100                 STREAMICE_umask(i+ki,j+kj,bi,bj) = 0.0                 STREAMICE_umask(i+ki,j+kj,bi,bj) = 3.0
101                ENDDO                ENDDO
102                STREAMICE_vfacemask(i+ki,j,bi,bj) = 3.0                STREAMICE_vfacemask(i,j+kj,bi,bj) = 3.0
103             ELSE IF (maskFlag.EQ.2) THEN             ELSE IF (maskFlag.EQ.2) THEN
104                DO ki=0,1                !DO ki=0,1
105                 STREAMICE_vfacemask(i+ki,j+kj,bi,bj) = 2.0                 STREAMICE_vfacemask(i,j+kj,bi,bj) = 2.0
106                ENDDO                !ENDDO
107             ELSE IF (maskFlag.EQ.4) THEN             ELSE IF (maskFlag.EQ.4) THEN
108                DO ki=0,1                DO ki=0,1
109                 STREAMICE_umask(i+ki,j+kj,bi,bj) = 0.0                 STREAMICE_umask(i+ki,j+kj,bi,bj) = 0.0
110                 STREAMICE_vmask(i+ki,j+kj,bi,bj) = 0.0                             STREAMICE_vmask(i+ki,j+kj,bi,bj) = 0.0            
111                ENDDO                ENDDO
112                STREAMICE_vfacemask(i+ki,j,bi,bj) = 4.0                STREAMICE_vfacemask(i,j+kj,bi,bj) = 4.0
113             ELSE IF (maskFlag.EQ.0) THEN             ELSE IF (maskFlag.EQ.0) THEN
114                DO ki=0,1                DO ki=0,1
115                 STREAMICE_umask(i+ki,j+kj,bi,bj) = 0.0                 STREAMICE_umask(i+ki,j+kj,bi,bj) = 0.0
# Line 164  C     === Global variables === Line 168  C     === Global variables ===
168    
169  !       CALL WRITE_FULLARRAY_RL ("umask",STREAMICE_umask,  !       CALL WRITE_FULLARRAY_RL ("umask",STREAMICE_umask,
170  !      c   1,0,0,1,0,myThid)  !      c   1,0,0,1,0,myThid)
171  !       CALL WRITE_FLD_XY_RL ("umask2","",STREAMICE_umask,0,myThid)        CALL WRITE_FLD_XY_RL ("umask","",STREAMICE_umask,0,myThid)
172          CALL WRITE_FLD_XY_RL ("vmask","",STREAMICE_vmask,0,myThid)
173          CALL WRITE_FLD_XY_RL ("ufacemask","",STREAMICE_ufacemask,0,myThid)
174          CALL WRITE_FLD_XY_RL ("vfacemask","",STREAMICE_vfacemask,0,myThid)
175    
176    
177    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22