/[MITgcm]/MITgcm_contrib/torge/itd/code/seaice_model.F
ViewVC logotype

Diff of /MITgcm_contrib/torge/itd/code/seaice_model.F

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

revision 1.3 by torge, Mon Jun 4 16:35:57 2012 UTC revision 1.12 by torge, Fri May 3 18:59:40 2013 UTC
# Line 55  CEndOfInterface Line 55  CEndOfInterface
55    
56  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
57  C     i,j,bi,bj :: Loop counters  C     i,j,bi,bj :: Loop counters
58  CToM<<<  #ifdef SEAICE_DEBUG
59    C     coordinates for debugging output
60          INTEGER i_dbOut, j_dbOut
61  C     msgBuf      :: Informational/error message buffer  C     msgBuf      :: Informational/error message buffer
62        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
63        CHARACTER*10 HlimitMsgFormat        CHARACTER*10 HlimitMsgFormat
64  C#if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC)  #endif
65  #if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC) || defined(SEAICE_ITD)  #ifdef SEAICE_ITD
66  C>>>ToM        INTEGER IT
67    #endif
68    #if defined(ALLOW_AUTODIFF_TAMC) || defined(SEAICE_ITD)
69        INTEGER i, j, bi, bj        INTEGER i, j, bi, bj
70  #endif  #endif
71  #ifdef ALLOW_SITRACER  #ifdef ALLOW_SITRACER
72        INTEGER iTr        INTEGER iTr
73  #endif  #endif
74  CEOP  CEOP
75    #ifdef SEAICE_DEBUG
76          i_dbOut = 1
77          j_dbOut = 1
78    #endif
79    
80  #ifdef ALLOW_DEBUG  #ifdef ALLOW_DEBUG
81        IF (debugMode) CALL DEBUG_ENTER( 'SEAICE_MODEL', myThid )        IF (debugMode) CALL DEBUG_ENTER( 'SEAICE_MODEL', myThid )
# Line 83  C--   Map thSice-variables to HEFF and A Line 91  C--   Map thSice-variables to HEFF and A
91        ENDIF        ENDIF
92  #endif /* ALLOW_THSICE */  #endif /* ALLOW_THSICE */
93    
 #ifdef SEAICE_GROWTH_LEGACY  
       IF ( .NOT.useThSice ) THEN  
 #ifdef ALLOW_AUTODIFF_TAMC  
 CADJ STORE heff  = comlev1, key=ikey_dynamics, kind=isbyte  
 CADJ STORE heffm = comlev1, key=ikey_dynamics, kind=isbyte  
 CADJ STORE area  = comlev1, key=ikey_dynamics, kind=isbyte  
 CADJ STORE hsnow = comlev1, key=ikey_dynamics, kind=isbyte  
 CADJ STORE tice  = comlev1, key=ikey_dynamics, kind=isbyte  
 #ifdef SEAICE_VARIABLE_SALINITY  
 CADJ STORE hsalt = comlev1, key=ikey_dynamics, kind=isbyte  
 #endif  
 #endif  
       DO bj=myByLo(myThid),myByHi(myThid)  
        DO bi=myBxLo(myThid),myBxHi(myThid)  
         DO j=1-OLy,sNy+OLy  
          DO i=1-OLx,sNx+OLx  
           IF ( (heff(i,j,bi,bj).EQ.0.)  
      &     .OR.(area(i,j,bi,bj).EQ.0.)  
      &     ) THEN  
            HEFF(i,j,bi,bj) = 0. _d 0  
            AREA(i,j,bi,bj) = 0. _d 0  
            HSNOW(i,j,bi,bj) = 0. _d 0  
            TICE(i,j,bi,bj) = celsius2K  
 #ifdef SEAICE_VARIABLE_SALINITY  
            HSALT(i,j,bi,bj) = 0. _d 0  
 #endif  
           ENDIF  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
       ENDIF  
 #endif  
   
94  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
95        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
96         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
97          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
98           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
 # ifdef SEAICE_GROWTH_LEGACY  
           areaNm1(i,j,bi,bj) = 0. _d 0  
           hEffNm1(i,j,bi,bj) = 0. _d 0  
 # endif  
99            uIceNm1(i,j,bi,bj) = 0. _d 0            uIceNm1(i,j,bi,bj) = 0. _d 0
100            vIceNm1(i,j,bi,bj) = 0. _d 0            vIceNm1(i,j,bi,bj) = 0. _d 0
101  # ifdef ALLOW_SITRACER  # ifdef ALLOW_SITRACER
# Line 195  C--   Apply ice velocity open boundary c Line 165  C--   Apply ice velocity open boundary c
165  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
166    
167  #ifdef ALLOW_THSICE  #ifdef ALLOW_THSICE
168        IF ( .NOT.useThSice ) THEN        IF ( useThSice ) THEN
169    #ifndef OLD_THSICE_CALL_SEQUENCE
170    #ifdef ALLOW_DEBUG
171            IF (debugMode) CALL DEBUG_CALL( 'THSICE_DO_ADVECT', myThid )
172    #endif
173            CALL THSICE_DO_ADVECT( 0, 0, myTime, myIter, myThid )
174    #endif /* OLD_THSICE_CALL_SEQUENCE */
175          ELSE
176  #endif  #endif
177  C--   Only call advection of heff, area, snow, and salt and  C--   Only call advection of heff, area, snow, and salt and
178  C--   growth for the generic 0-layer thermodynamics of seaice  C--   growth for the generic 0-layer thermodynamics of seaice
# Line 205  C--   (called from DO_OCEANIC_PHYSICS) t Line 182  C--   (called from DO_OCEANIC_PHYSICS) t
182  C NOW DO ADVECTION and DIFFUSION  C NOW DO ADVECTION and DIFFUSION
183        IF ( SEAICEadvHeff .OR. SEAICEadvArea .OR. SEAICEadvSnow        IF ( SEAICEadvHeff .OR. SEAICEadvArea .OR. SEAICEadvSnow
184       &        .OR. SEAICEadvSalt ) THEN       &        .OR. SEAICEadvSalt ) THEN
 CToM<<<  
185  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
186  C     test if REDIST works (1/2):  #ifdef SEAICE_DEBUG
187         WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F6.2)'         WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F8.4)'
188         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
189          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
190  c         DO j=1-OLy,sNy+OLy             WRITE(msgBuf,'(A,F8.4,x,F8.4)')
191  c          DO i=1-OLx,sNx+OLx       &       ' SEAICE_MODEL: AREA and HEFF before advection: ',
192  ccc           WRITE(msgBuf,HlimitMsgFormat)       &       AREA(i_dbOut,j_dbOut,bi,bj), HEFF(i_dbOut,j_dbOut,bi,bj)
193             WRITE(msgBuf,'(A,F6.2,x,F6.2)')  c     &       ' SEAICE_MODEL: AREA and HEFF/AREA before advection: ',
194       &       ' SEAICE_MODEL: AREA and HEFF/AREA before advection: ',  c     &       AREA(i_dbOut,j_dbOut,bi,bj), HEFF(i_dbOut,j_dbOut,bi,bj)
195       &       AREA(20,20,bi,bj), HEFF(20,20,bi,bj)/AREA(20,20,bi,bj)  c     &                                   /AREA(i_dbOut,j_dbOut,bi,bj)
196             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
197       &       SQUEEZE_RIGHT , myThid)       &       SQUEEZE_RIGHT , myThid)
198             WRITE(msgBuf,HlimitMsgFormat)             WRITE(msgBuf,HlimitMsgFormat)
199       &       ' SEAICE_MODEL: HEFFITD/AREAITD before advection: ',       &       ' SEAICE_MODEL: HEFFITD       before advection: ',
200       &       HEFFITD(20,20,:,bi,bj) / AREAITD(20,20,:,bi,bj)       &       HEFFITD(i_dbOut,j_dbOut,:,bi,bj)
201    c     &       ' SEAICE_MODEL: HEFFITD/AREAITD before advection: ',
202    c     &       HEFFITD(i_dbOut,j_dbOut,:,bi,bj)
203    c     &      /AREAITD(i_dbOut,j_dbOut,:,bi,bj)
204             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
205       &       SQUEEZE_RIGHT , myThid)       &       SQUEEZE_RIGHT , myThid)
206             WRITE(msgBuf,HlimitMsgFormat)             WRITE(msgBuf,HlimitMsgFormat)
207       &       ' SEAICE_MODEL: AREAITD         before advection: ',       &       ' SEAICE_MODEL: AREAITD       before advection: ',
208       &       AREAITD(20,20,:,bi,bj)       &       AREAITD(i_dbOut,j_dbOut,:,bi,bj)
209             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
210       &       SQUEEZE_RIGHT , myThid)       &       SQUEEZE_RIGHT , myThid)
 c          ENDDO  
 c         ENDDO  
211          ENDDO          ENDDO
212         ENDDO         ENDDO
213  #endif  #endif
214  C>>>ToM  #endif
215  #ifdef ALLOW_DEBUG  #ifdef ALLOW_DEBUG
216         IF (debugMode) CALL DEBUG_CALL( 'SEAICE_ADVDIFF', myThid )         IF (debugMode) CALL DEBUG_CALL( 'SEAICE_ADVDIFF', myThid )
217  #endif  #endif
218         CALL SEAICE_ADVDIFF( myTime, myIter, myThid )         CALL SEAICE_ADVDIFF( myTime, myIter, myThid )
 CToM<<<  
219  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
220  C     test if REDIST works (1/2):  #ifdef SEAICE_DEBUG
221         WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F6.2)'         WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F8.4)'
222         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
223          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
224             WRITE(msgBuf,HlimitMsgFormat)             WRITE(msgBuf,HlimitMsgFormat)
225       &       ' SEAICE_MODEL: HEFFITD/AREAITD before 1. sorting: ',       &       ' SEAICE_MODEL: HEFFITD        after advection: ',
226       &       HEFFITD(20,20,:,bi,bj) / AREAITD(20,20,:,bi,bj)       &       HEFFITD(i_dbOut,j_dbOut,:,bi,bj)
227             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
228       &       SQUEEZE_RIGHT , myThid)       &       SQUEEZE_RIGHT , myThid)
229             WRITE(msgBuf,HlimitMsgFormat)             WRITE(msgBuf,HlimitMsgFormat)
230       &       ' SEAICE_MODEL: AREAITD before         1. sorting: ',       &       ' SEAICE_MODEL: AREAITD        after advection: ',
231       &       AREAITD(20,20,:,bi,bj)       &       AREAITD(i_dbOut,j_dbOut,:,bi,bj)
232             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
233       &       SQUEEZE_RIGHT , myThid)       &       SQUEEZE_RIGHT , myThid)
234                WRITE(msgBuf,'(A)')
235         &        ' --------------------------------------------- '
236                CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
237         &        SQUEEZE_RIGHT , myThid)
238          ENDDO          ENDDO
239         ENDDO         ENDDO
240    #endif
241  C  C
242  C     check that all ice thickness categories meet their limits  C     check that all ice thickness categories meet their limits
243  C     (includes Hibler-type ridging)  C     (includes Hibler-type ridging)
244  #ifdef ALLOW_DEBUG  #ifdef ALLOW_DEBUG
245          IF (debugMode) CALL DEBUG_CALL( 'SEAICE_ITD_REDIST', myThid )          IF (debugMode) CALL DEBUG_CALL( 'SEAICE_ITD_REDIST', myThid )
246  #endif  #endif
247         CALL SEAICE_ITD_REDIST(myTime, myIter, myThid)         DO bj=myByLo(myThid),myByHi(myThid)
248  C     store the mean ice thickness in HEFF (for dynamic solver and diagnostics)          DO bi=myBxLo(myThid),myBxHi(myThid)
249         CALL SEAICE_ITD_SUM(myTime, myIter, myThid)           CALL SEAICE_ITD_REDIST(bi, bj, myTime, myIter, myThid)
250            ENDDO
251           ENDDO
252    C     update mean ice thickness HEFF and total ice concentration AREA
253    C     to match single category values
254    #ifdef ALLOW_DEBUG
255            IF (debugMode) CALL DEBUG_CALL( 'SEAICE_ITD_SUM', myThid )
256  #endif  #endif
 C>>>ToM  
 #ifdef SEAICE_GROWTH_LEGACY  
       ELSE  
257         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
258          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
259           DO j=1-OLy,sNy+OLy           CALL SEAICE_ITD_SUM(bi, bj, myTime, myIter, myThid)
           DO i=1-OLx,sNx+OLx  
            areaNm1(i,j,bi,bj) = AREA(i,j,bi,bj)  
            hEffNm1(i,j,bi,bj) = HEFF(i,j,bi,bj)  
           ENDDO  
          ENDDO  
260          ENDDO          ENDDO
261         ENDDO         ENDDO
262  #endif /* SEAICE_GROWTH_LEGACY */  #endif
263        ENDIF        ENDIF
264  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
265  CADJ STORE heffm  = comlev1, key=ikey_dynamics, kind=isbyte  CADJ STORE heffm  = comlev1, key=ikey_dynamics, kind=isbyte
266  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
267    
268    #ifdef SEAICE_ITD
269    #ifdef SEAICE_DEBUG
270           WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F8.4)'
271           DO bj=myByLo(myThid),myByHi(myThid)
272            DO bi=myBxLo(myThid),myBxHi(myThid)
273               WRITE(msgBuf,HlimitMsgFormat)
274         &       ' SEAICE_MODEL: HEFFITD          before growth: ',
275         &       HEFFITD(i_dbOut,j_dbOut,:,bi,bj)
276               CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
277         &       SQUEEZE_RIGHT , myThid)
278               WRITE(msgBuf,HlimitMsgFormat)
279         &       ' SEAICE_MODEL: AREAITD          before growth: ',
280         &       AREAITD(i_dbOut,j_dbOut,:,bi,bj)
281               CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
282         &       SQUEEZE_RIGHT , myThid)
283               WRITE(msgBuf,HlimitMsgFormat)
284         &       ' SEAICE_MODEL: HSNOWITD         before growth: ',
285         &       HSNOWITD(i_dbOut,j_dbOut,:,bi,bj)
286               CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
287         &       SQUEEZE_RIGHT , myThid)
288            ENDDO
289           ENDDO
290    #endif
291    #endif
292    
293  #ifndef DISABLE_SEAICE_GROWTH  #ifndef DISABLE_SEAICE_GROWTH
294  C     thermodynamics growth  C     thermodynamics growth
295  C     must call growth after calling advection  C     must call growth after calling advection
# Line 295  C     because of ugly time level busines Line 299  C     because of ugly time level busines
299          IF (debugMode) CALL DEBUG_CALL( 'SEAICE_GROWTH', myThid )          IF (debugMode) CALL DEBUG_CALL( 'SEAICE_GROWTH', myThid )
300  #endif  #endif
301          CALL SEAICE_GROWTH( myTime, myIter, myThid )          CALL SEAICE_GROWTH( myTime, myIter, myThid )
 CToM<<<  
302  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
303  C     test if GROWTH works:  #ifdef SEAICE_DEBUG
304          WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F6.2)'          WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F8.4)'
305          DO bj=myByLo(myThid),myByHi(myThid)          DO bj=myByLo(myThid),myByHi(myThid)
306           DO bi=myBxLo(myThid),myBxHi(myThid)           DO bi=myBxLo(myThid),myBxHi(myThid)
307              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
308       &       ' SEAICE_MODEL: HEFFITD/AREAITD after growth: ',       &        ' SEAICE_MODEL: HEFFITD           after growth: ',
309       &       HEFFITD(20,20,:,bi,bj) / AREAITD(20,20,:,bi,bj)       &        HEFFITD(i_dbOut,j_dbOut,:,bi,bj)
310              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
311       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
312              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
313       &        ' SEAICE_MODEL: AREAITD        after growth: ',       &        ' SEAICE_MODEL: AREAITD           after growth: ',
314       &        AREAITD(20,20,:,bi,bj)       &        AREAITD(i_dbOut,j_dbOut,:,bi,bj)
315              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
316       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
317              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
318       &        ' SEAICE_MODEL: HSNOWITD       after growth: ',       &        ' SEAICE_MODEL: HSNOWITD          after growth: ',
319       &        HSNOWITD(20,20,:,bi,bj)       &        HSNOWITD(i_dbOut,j_dbOut,:,bi,bj)
320                CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
321         &        SQUEEZE_RIGHT , myThid)
322                WRITE(msgBuf,'(A)')
323         &        ' --------------------------------------------- '
324              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
325       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
326           ENDDO           ENDDO
327          ENDDO          ENDDO
328    #endif
329  C  C
330  C     redistribute sea ice into proper sea ice category after growth/melt  C     redistribute sea ice into proper sea ice category after growth/melt
331  C     in case model runs with ice thickness distribution  C     in case model runs with ice thickness distribution
# Line 325  C---+-|--1----+----2----+----3----+----4 Line 333  C---+-|--1----+----2----+----3----+----4
333  #ifdef ALLOW_DEBUG  #ifdef ALLOW_DEBUG
334          IF (debugMode) CALL DEBUG_CALL( 'SEAICE_ITD_REDIST', myThid )          IF (debugMode) CALL DEBUG_CALL( 'SEAICE_ITD_REDIST', myThid )
335  #endif  #endif
336         CALL SEAICE_ITD_REDIST(myTime, myIter, myThid)         DO bj=myByLo(myThid),myByHi(myThid)
337            DO bi=myBxLo(myThid),myBxHi(myThid)
338             CALL SEAICE_ITD_REDIST(bi, bj, myTime, myIter, myThid)
339            ENDDO
340           ENDDO
341  C     store the mean ice thickness in HEFF (for dynamic solver and diagnostics)  C     store the mean ice thickness in HEFF (for dynamic solver and diagnostics)
342         CALL SEAICE_ITD_SUM(myTime, myIter, myThid)         DO bj=myByLo(myThid),myByHi(myThid)
343            DO bi=myBxLo(myThid),myBxHi(myThid)
344             CALL SEAICE_ITD_SUM(bi, bj, myTime, myIter, myThid)
345            ENDDO
346           ENDDO
347    
348  C     test if REDIST works (2/2):  #ifdef SEAICE_DEBUG
349          WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F6.2)'          WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F8.4)'
350          DO bj=myByLo(myThid),myByHi(myThid)          DO bj=myByLo(myThid),myByHi(myThid)
351           DO bi=myBxLo(myThid),myBxHi(myThid)           DO bi=myBxLo(myThid),myBxHi(myThid)
352              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
353       &        ' SEAICE_MODEL: HEFFITD/AREAITD after 2. sorting: ',       &        ' SEAICE_MODEL: HEFFITD    after final sorting: ',
354       &        HEFFITD(20,20,:,bi,bj) / AREAITD(20,20,:,bi,bj)       &        HEFFITD(i_dbOut,j_dbOut,:,bi,bj)
355              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
356       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
357              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
358       &        ' SEAICE_MODEL: AREAITD         after 2. sorting: ',       &        ' SEAICE_MODEL: AREAITD    after final sorting: ',
359       &        AREAITD(20,20,:,bi,bj)       &        AREAITD(i_dbOut,j_dbOut,:,bi,bj)
360              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
361       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
362              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
363       &        ' SEAICE_MODEL: HSNOWITD        after 2. sorting: ',       &        ' SEAICE_MODEL: HSNOWITD   after final sorting: ',
364       &        HSNOWITD(20,20,:,bi,bj)       &        HSNOWITD(i_dbOut,j_dbOut,:,bi,bj)
365                CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
366         &        SQUEEZE_RIGHT , myThid)
367                WRITE(msgBuf,'(A)')
368         &        ' ============================================= '
369              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
370       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
371           ENDDO           ENDDO
372          ENDDO          ENDDO
373  #endif  #endif
374    #endif
375  C  C
 C>>>ToM  
376         ENDIF         ENDIF
377  #endif /* DISABLE_SEAICE_GROWTH */  #endif /* DISABLE_SEAICE_GROWTH */
378    
# Line 427  C     the seaice model Line 447  C     the seaice model
447          CALL DIAGNOSTICS_FILL_RS(EmPmR,'SIempmr ',0,1 ,0,1,1,myThid)          CALL DIAGNOSTICS_FILL_RS(EmPmR,'SIempmr ',0,1 ,0,1,1,myThid)
448          CALL DIAGNOSTICS_FILL_RS(Qnet ,'SIqnet  ',0,1 ,0,1,1,myThid)          CALL DIAGNOSTICS_FILL_RS(Qnet ,'SIqnet  ',0,1 ,0,1,1,myThid)
449          CALL DIAGNOSTICS_FILL_RS(Qsw  ,'SIqsw   ',0,1 ,0,1,1,myThid)          CALL DIAGNOSTICS_FILL_RS(Qsw  ,'SIqsw   ',0,1 ,0,1,1,myThid)
450    #ifdef SEAICE_ITD
451            CALL DIAGNOSTICS_FILL(HEFFITD ,'SIheffN ',0,nITD,0,1,1,myThid)
452            CALL DIAGNOSTICS_FILL(AREAITD ,'SIareaN ',0,nITD,0,1,1,myThid)
453    #endif
454         ENDIF         ENDIF
455  #endif /* ALLOW_DIAGNOSTICS */  #endif /* ALLOW_DIAGNOSTICS */
456    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22