/[MITgcm]/MITgcm/pkg/seaice/seaice_init_fixed.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_init_fixed.F

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

revision 1.4 by mlosch, Wed Jun 17 07:33:30 2009 UTC revision 1.5 by heimbach, Thu Jun 25 14:36:15 2009 UTC
# Line 72  C--   Initialize grid info Line 72  C--   Initialize grid info
72          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
73           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
74            HEFFM(i,j,bi,bj)       = 0. _d 0            HEFFM(i,j,bi,bj)       = 0. _d 0
 #ifdef SEAICE_CGRID  
           seaiceMaskU(i,j,bi,bj) = 0. _d 0  
           seaiceMaskV(i,j,bi,bj) = 0. _d 0  
 #else  
           UVM(i,j,bi,bj)         = 0. _d 0  
 #endif /* not SEAICE_CGRID */  
75           ENDDO           ENDDO
76          ENDDO          ENDDO
77          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
# Line 89  C--   Initialize grid info Line 83  C--   Initialize grid info
83          ENDDO          ENDDO
84          DO j=1-OLy+1,sNy+OLy          DO j=1-OLy+1,sNy+OLy
85           DO i=1-OLx+1,sNx+OLx           DO i=1-OLx+1,sNx+OLx
86  #ifdef SEAICE_CGRID  #ifndef SEAICE_CGRID
           seaiceMaskU(i,j,bi,bj)=   0.0 _d 0  
           seaiceMaskV(i,j,bi,bj)=   0.0 _d 0  
           mask_uice=HEFFM(i,j,bi,bj)+HEFFM(i-1,j  ,bi,bj)  
           IF(mask_uice.GT.1.5 _d 0) seaiceMaskU(i,j,bi,bj)=1.0 _d 0  
           mask_uice=HEFFM(i,j,bi,bj)+HEFFM(i  ,j-1,bi,bj)  
           IF(mask_uice.GT.1.5 _d 0) seaiceMaskV(i,j,bi,bj)=1.0 _d 0  
 #else  
87            UVM(i,j,bi,bj)=0. _d 0            UVM(i,j,bi,bj)=0. _d 0
88            mask_uice=HEFFM(i,j,  bi,bj)+HEFFM(i-1,j-1,bi,bj)            mask_uice=HEFFM(i,j,  bi,bj)+HEFFM(i-1,j-1,bi,bj)
89       &             +HEFFM(i,j-1,bi,bj)+HEFFM(i-1,j,  bi,bj)       &             +HEFFM(i,j-1,bi,bj)+HEFFM(i-1,j,  bi,bj)
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22