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

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

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

revision 1.2 by dgoldberg, Thu Jul 19 18:46:56 2012 UTC revision 1.3 by dgoldberg, Tue Sep 4 21:11:44 2012 UTC
# Line 54  C     ZERO OUT FLOATING POINT ARRAYS Line 54  C     ZERO OUT FLOATING POINT ARRAYS
54            surf_el_streamice(i,j,bi,bj) = 0. _d 0            surf_el_streamice(i,j,bi,bj) = 0. _d 0
55            area_shelf_streamice(i,j,bi,bj) = 0. _d 0            area_shelf_streamice(i,j,bi,bj) = 0. _d 0
56            mass_ice_streamice(i,j,bi,bj) = 0. _d 0            mass_ice_streamice(i,j,bi,bj) = 0. _d 0
57              BDOT_streamice(i,j,bi,bj) = 0. _d 0
58            C_basal_friction(i,j,bi,bj) = C_basal_fric_const            C_basal_friction(i,j,bi,bj) = C_basal_fric_const
59              A_glen(i,j,bi,bj) = A_glen_isothermal
60  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
61            ru_old_si(i,j,bi,bj) = 0. _d 0            ru_old_si(i,j,bi,bj) = 0. _d 0
62            rv_old_si(i,j,bi,bj) = 0. _d 0            rv_old_si(i,j,bi,bj) = 0. _d 0
# Line 123  C     INIT. INTEGER ARRAYS Line 125  C     INIT. INTEGER ARRAYS
125            DO i=1,sNx            DO i=1,sNx
126             Gi = (myXGlobalLo-1)+(bi-1)*sNx+i             Gi = (myXGlobalLo-1)+(bi-1)*sNx+i
127             Gj = (myYGlobalLo-1)+(bj-1)*sNy+j             Gj = (myYGlobalLo-1)+(bj-1)*sNy+j
   
             
128                        
129             IF ((Gi.lt.Nx).and.(Gj.lt.Ny)) THEN             IF ((Gi.lt.Nx).and.(Gj.lt.Ny)) THEN
130                            
# Line 158  C             IF (flow_dir .EQ. 2.0) THE Line 158  C             IF (flow_dir .EQ. 2.0) THE
158                 STREAMICE_hmask(i,j,bi,bj) = 1.0                 STREAMICE_hmask(i,j,bi,bj) = 1.0
159                ENDIF                ENDIF
160    
 c              WRITE(msgBuf,'(A,I3,I3,E9.2)') 'thickness',  
 c     &          i,j,shelf_edge_pos  
 c              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
 c     &                    SQUEEZE_RIGHT , 1)  
161    
162               ENDIF               ENDIF
163              ENDIF              ENDIF
# Line 171  c     &                    SQUEEZE_RIGHT Line 167  c     &                    SQUEEZE_RIGHT
167          ENDDO          ENDDO
168         ENDDO         ENDDO
169    
170  !        DO bj = myByLo(myThid), myByHi(myThid)  
 !         DO bi = myBxLo(myThid), myBxHi(myThid)  
 !          DO j=1,sNy  
 !           DO i=1,sNx  
 !            STREAMICE_dummy_array(i,j,bi,bj) =  
 !      &      REAL(STREAMICE_hmask(i,j,bi,bj))  
 !            WRITE(msgBuf,'(F3.0)') STREAMICE_dummy_array(i,j,bi,bj)  
 !            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
 !      &                    SQUEEZE_RIGHT , 1)  
 !           ENDDO  
 !          ENDDO  
 !         ENDDO  
 !        ENDDO  
 !          
 !        CALL WRITE_FLD_XY_RL("STREAMICE_hmask"," ",  
 !      &    STREAMICE_dummy_array,0,myThid)      
 !        CALL WRITE_FLD_XY_RL("h_streamice"," ",  
 !      &    H_streamice,0,myThid)  
171    
172        ELSE IF ( STREAMICEthickInit.EQ.'FILE' ) THEN        ELSE IF ( STREAMICEthickInit.EQ.'FILE' ) THEN
173    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22