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

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

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

revision 1.32 by heimbach, Sat Oct 11 20:37:05 2008 UTC revision 1.33 by dimitri, Wed Dec 17 03:33:30 2008 UTC
# Line 118  C--   Initialise all variables in common Line 118  C--   Initialise all variables in common
118  #ifdef SEAICE_SALINITY  #ifdef SEAICE_SALINITY
119            HSALT(i,j,bi,bj)  = 0. _d 0            HSALT(i,j,bi,bj)  = 0. _d 0
120  #endif  #endif
121    #ifdef SEAICE_AGE
122              ICEAGE(i,j,bi,bj) = 0. _d 0
123    #endif
124            HEFFM(i,j,bi,bj)  = 0. _d 0            HEFFM(i,j,bi,bj)  = 0. _d 0
125            YNEG (i,j,bi,bj)  = 0. _d 0            YNEG (i,j,bi,bj)  = 0. _d 0
126            RIVER(i,j,bi,bj)  = 0. _d 0            RIVER(i,j,bi,bj)  = 0. _d 0
# Line 421  C--   Read initial sea ice salinity from Line 424  C--   Read initial sea ice salinity from
424         ENDIF         ENDIF
425  #endif /* SEAICE_SALINITY */  #endif /* SEAICE_SALINITY */
426    
427    #ifdef SEAICE_AGE
428    C--   Read initial sea ice age from file if available.
429           IF ( IceAgeFile .NE. ' ' ) THEN
430            CALL READ_FLD_XY_RL( IceAgeFile, ' ', ZETA, 0, myThid )
431            _EXCH_XY_R8(ZETA,myThid)
432            DO bj=myByLo(myThid),myByHi(myThid)
433             DO bi=myBxLo(myThid),myBxHi(myThid)
434              DO j=1-OLy,sNy+OLy
435               DO i=1-OLx,sNx+OLx
436                ICEAGE(i,j,bi,bj) = ZETA(i,j,bi,bj)
437               ENDDO
438              ENDDO
439             ENDDO
440            ENDDO
441           ENDIF
442    #endif /* SEAICE_AGE */
443    
444        ENDIF        ENDIF
445    
446  C---  Complete initialization  C---  Complete initialization
Legend:
Removed from v.1.32  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22