/[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.49 by heimbach, Fri Nov 20 22:33:21 2009 UTC revision 1.50 by jmc, Sun Jan 3 00:27:36 2010 UTC
# Line 18  C     === Global variables === Line 18  C     === Global variables ===
18  #include "PARAMS.h"  #include "PARAMS.h"
19  #include "GRID.h"  #include "GRID.h"
20  #include "DYNVARS.h"  #include "DYNVARS.h"
 #include "SEAICE.h"  
 #include "SEAICE_DIAGS.h"  
 #include "SEAICE_PARAMS.h"  
21  #include "FFIELDS.h"  #include "FFIELDS.h"
22    #include "SEAICE_PARAMS.h"
23    #include "SEAICE.h"
24    #include "SEAICE_TAVE.h"
25  #ifdef ALLOW_EXCH2  #ifdef ALLOW_EXCH2
26  # include "W2_EXCH2_SIZE.h"  # include "W2_EXCH2_SIZE.h"
27  # include "W2_EXCH2_TOPOLOGY.h"  # include "W2_EXCH2_TOPOLOGY.h"
# Line 39  CEndOfInterface Line 39  CEndOfInterface
39  C     === Local variables ===  C     === Local variables ===
40  C     i,j,k,bi,bj - Loop counters  C     i,j,k,bi,bj - Loop counters
41    
42        INTEGER i, j, k, bi, bj        INTEGER i, j, bi, bj
43        _RL PSTAR        _RL PSTAR
44        _RS  mask_uice        _RS  mask_uice
45        INTEGER myIter, myTile, kSurface        INTEGER myIter, kSurface
46    #ifdef SEAICE_MULTICATEGORY
47          INTEGER k
48    #endif
49  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
50        INTEGER I_obc, J_obc        INTEGER I_obc, J_obc
51  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
52    #ifdef ALLOW_EXCH2
53    #ifndef SEAICE_CGRID
54          INTEGER myTile
55    #endif
56    #endif
57    
58        IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN        IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN
59         kSurface = Nr         kSurface = Nr
# Line 140  C Line 148  C
148         ENDDO         ENDDO
149        ENDDO        ENDDO
150    
151    #ifdef ALLOW_TIMEAVE
152    C     Initialize averages to zero
153          DO bj = myByLo(myThid), myByHi(myThid)
154           DO bi = myBxLo(myThid), myBxHi(myThid)
155            CALL TIMEAVE_RESET( FUtave   , 1, bi, bj, myThid )
156            CALL TIMEAVE_RESET( FVtave   , 1, bi, bj, myThid )
157            CALL TIMEAVE_RESET( EmPmRtave, 1, bi, bj, myThid )
158            CALL TIMEAVE_RESET( QNETtave , 1, bi, bj, myThid )
159            CALL TIMEAVE_RESET( QSWtave  , 1, bi, bj, myThid )
160            CALL TIMEAVE_RESET( UICEtave , 1, bi, bj, myThid )
161            CALL TIMEAVE_RESET( VICEtave , 1, bi, bj, myThid )
162            CALL TIMEAVE_RESET( HEFFtave , 1, bi, bj, myThid )
163            CALL TIMEAVE_RESET( AREAtave , 1, bi, bj, myThid )
164            SEAICE_timeAve(bi,bj) = ZERO
165           ENDDO
166          ENDDO
167    #endif /* ALLOW_TIMEAVE */
168    
169  C--   Initialize (variable) grid info. As long as we allow masking of  C--   Initialize (variable) grid info. As long as we allow masking of
170  C--   velocities outside of ice covered areas (in seaice_dynsolver)  C--   velocities outside of ice covered areas (in seaice_dynsolver)
171  C--   we need to re-initialize seaiceMaskU/V here for TAF/TAMC  C--   we need to re-initialize seaiceMaskU/V here for TAF/TAMC
Legend:
Removed from v.1.49  
changed lines
  Added in v.1.50
  ViewVC Help
Powered by ViewVC 1.1.22