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

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_init_fixed.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.2 by dgoldberg, Tue Sep 18 17:06:48 2012 UTC
# Line 29  CEOP Line 29  CEOP
29    
30  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
31  C     === Local variables ===  C     === Local variables ===
32        INTEGER bi, bj, i, j, Gi, Gj        INTEGER bi, bj, i, j, Gi, Gj, m
33        _RL x, y        _RL x, y
34        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
35    
# Line 62  C       The 0 is the "iteration" argumen Line 62  C       The 0 is the "iteration" argumen
62         ENDIF         ENDIF
63        ENDIF        ENDIF
64    
65    
66    !  INITIALIZE SIGMA COORD
67          IF (STREAMICEsigcoordInit.eq.'FILE') THEN
68           WRITE(msgBuf,'(A)') 'SIG FROM FILE - NOT IMPLENTED'
69           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
70         &                    SQUEEZE_RIGHT , 1)
71          ELSEIF (STREAMICEsigcoordInit.eq.'UNIFORM') THEN
72           do m=1,Nr
73            streamice_sigma_coord (m) = 1./Nr * (m-0.5)
74            streamice_delsigma (m) = 1./Nr
75           enddo
76          ELSE
77           WRITE(msgBuf,'(A)') 'SIG COORD INIT - NOT IMPLENTED'
78           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
79         &                    SQUEEZE_RIGHT , 1)
80          ENDIF
81    !!!!!!!!!!!!!!!!!!!!!!!!!
82    
83  C- fill in the overlap (+ BARRIER):  C- fill in the overlap (+ BARRIER):
84        _EXCH_XY_RL(STREAMICE_calve_mask, myThid )        _EXCH_XY_RL(STREAMICE_calve_mask, myThid )
85    
# Line 77  C         INIT VALUES FOR METRIC TERMS Line 95  C         INIT VALUES FOR METRIC TERMS
95            k2AtC_str(i,j,bi,bj) = recip_rA(i,j,bi,bj) *            k2AtC_str(i,j,bi,bj) = recip_rA(i,j,bi,bj) *
96       &     (dyG(i,j+1,bi,bj)-dyG(i,j,bi,bj))       &     (dyG(i,j+1,bi,bj)-dyG(i,j,bi,bj))
97    
98    #ifdef STREAMICE_HYBRID_STRESS
99              streamice_basal_geom (i,j,bi,bj) = 1.0
100    #endif          
101    
102  C         INIT BDRY CONDITIONS  C         INIT BDRY CONDITIONS
103    
104            Gi = (myXGlobalLo-1)+(bi-1)*sNx+i            Gi = (myXGlobalLo-1)+(bi-1)*sNx+i
# Line 87  C         INIT BDRY CONDITIONS Line 109  C         INIT BDRY CONDITIONS
109            STREAMICE_ufacemask_bdry (i,j,bi,bj) = -1.0            STREAMICE_ufacemask_bdry (i,j,bi,bj) = -1.0
110            STREAMICE_vfacemask_bdry (i,j,bi,bj) = -1.0            STREAMICE_vfacemask_bdry (i,j,bi,bj) = -1.0
111    
112            IF (Gi .EQ. 1) THEN            IF (Gi.EQ.1 .AND. .NOT. STREAMICE_EW_periodic) THEN
113  C         INITIALIZE BOUNDARY CONDS AT WEST  C         INITIALIZE BOUNDARY CONDS AT WEST
114             IF ((y .gt. min_y_noflow_WEST) .AND.             IF ((y .gt. min_y_noflow_WEST) .AND.
115       &         (y .le. max_y_noflow_WEST)) THEN       &         (y .le. max_y_noflow_WEST)) THEN
# Line 113  C         INITIALIZE BOUNDARY CONDS AT W Line 135  C         INITIALIZE BOUNDARY CONDS AT W
135             ENDIF             ENDIF
136            ENDIF            ENDIF
137    
138            IF (Gi .EQ. Nx) THEN            IF (Gi.EQ.Nx .AND. .NOT. STREAMICE_EW_periodic) THEN
139  C         INITIALIZE BOUNDARY CONDS AT EAST  C         INITIALIZE BOUNDARY CONDS AT EAST
140             IF ((y .gt. min_y_noflow_EAST) .AND.             IF ((y .gt. min_y_noflow_EAST) .AND.
141       &         (y .le. max_y_noflow_EAST)) THEN       &         (y .le. max_y_noflow_EAST)) THEN
# Line 139  C         INITIALIZE BOUNDARY CONDS AT E Line 161  C         INITIALIZE BOUNDARY CONDS AT E
161             ENDIF             ENDIF
162            ENDIF            ENDIF
163    
164            IF (Gj .EQ. 1) THEN            IF (Gj.EQ.1 .AND. .NOT. STREAMICE_NS_periodic) THEN
165  C         INITIALIZE BOUNDARY CONDS AT SOUTH  C         INITIALIZE BOUNDARY CONDS AT SOUTH
166             IF ((x .gt. min_x_noflow_SOUTH) .AND.             IF ((x .gt. min_x_noflow_SOUTH) .AND.
167       &         (x .le. max_x_noflow_SOUTH)) THEN       &         (x .le. max_x_noflow_SOUTH)) THEN
# Line 165  C         INITIALIZE BOUNDARY CONDS AT S Line 187  C         INITIALIZE BOUNDARY CONDS AT S
187             ENDIF             ENDIF
188            ENDIF            ENDIF
189    
190            IF (Gj .EQ. Ny) THEN            IF (Gj.EQ.Ny .AND. .NOT. STREAMICE_NS_periodic) THEN
191  C         INITIALIZE BOUNDARY CONDS AT NORTH  C         INITIALIZE BOUNDARY CONDS AT NORTH
192             IF ((x .gt. min_x_noflow_NORTH) .AND.             IF ((x .gt. min_x_noflow_NORTH) .AND.
193       &         (x .le. max_x_noflow_NORTH)) THEN       &         (x .le. max_x_noflow_NORTH)) THEN
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22