/[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.2 by dimitri, Fri Apr 27 22:25:23 2012 UTC revision 1.3 by torge, Mon Jun 4 16:35:57 2012 UTC
# Line 64  C#if defined(SEAICE_GROWTH_LEGACY) || de Line 64  C#if defined(SEAICE_GROWTH_LEGACY) || de
64  C>>>ToM  C>>>ToM
65        INTEGER i, j, bi, bj        INTEGER i, j, bi, bj
66  #endif  #endif
 CToM<<<  
 #ifdef SEAICE_ITD  
       INTEGER k  
 #endif  
 C>>>ToM  
67  #ifdef ALLOW_SITRACER  #ifdef ALLOW_SITRACER
68        INTEGER iTr        INTEGER iTr
69  #endif  #endif
# Line 212  C NOW DO ADVECTION and DIFFUSION Line 207  C NOW DO ADVECTION and DIFFUSION
207       &        .OR. SEAICEadvSalt ) THEN       &        .OR. SEAICEadvSalt ) THEN
208  CToM<<<  CToM<<<
209  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
 c just for testing until seaice_growth includes itd  
 c       DO bj=myByLo(myThid),myByHi(myThid)  
 c        DO bi=myBxLo(myThid),myBxHi(myThid)  
 c         DO j=1-OLy,sNy+OLy  
 c          DO i=1-OLx,sNx+OLx  
 c           AREAITD (i,j,1,bi,bj)=AREA (i,j,bi,bj)  
 c           HEFFITD (i,j,1,bi,bj)=HEFF (i,j,bi,bj)  
 c           HSNOWITD(i,j,1,bi,bj)=HSNOW(i,j,bi,bj)  
 c          ENDDO  
 c         ENDDO  
 c         DO k=2,nITD  
 c          DO j=1-OLy,sNy+OLy  
 c           DO i=1-OLx,sNx+OLx  
 c            AREAITD (i,j,k,bi,bj) = ZERO  
 c            HEFFITD (i,j,k,bi,bj) = ZERO  
 c            HSNOWITD(i,j,k,bi,bj) = ZERO  
 c           ENDDO  
 c          ENDDO  
 c         ENDDO  
 c        ENDDO  
 c       ENDDO  
 c end this test  
210  C     test if REDIST works (1/2):  C     test if REDIST works (1/2):
211         WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F6.2)'         WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F6.2)'
212         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
# Line 242  c         DO j=1-OLy,sNy+OLy Line 215  c         DO j=1-OLy,sNy+OLy
215  c          DO i=1-OLx,sNx+OLx  c          DO i=1-OLx,sNx+OLx
216  ccc           WRITE(msgBuf,HlimitMsgFormat)  ccc           WRITE(msgBuf,HlimitMsgFormat)
217             WRITE(msgBuf,'(A,F6.2,x,F6.2)')             WRITE(msgBuf,'(A,F6.2,x,F6.2)')
218       &       ' SEAICE_MODEL: AREA and HEFF before advection: ',       &       ' SEAICE_MODEL: AREA and HEFF/AREA before advection: ',
219       &       AREA(20,20,bi,bj), HEFF(20,20,bi,bj)       &       AREA(20,20,bi,bj), HEFF(20,20,bi,bj)/AREA(20,20,bi,bj)
220             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
221       &       SQUEEZE_RIGHT , myThid)       &       SQUEEZE_RIGHT , myThid)
222             WRITE(msgBuf,HlimitMsgFormat)             WRITE(msgBuf,HlimitMsgFormat)
223       &       ' SEAICE_MODEL: HEFFITD before advection: ',       &       ' SEAICE_MODEL: HEFFITD/AREAITD before advection: ',
224       &       HEFFITD(20,20,:,bi,bj)       &       HEFFITD(20,20,:,bi,bj) / AREAITD(20,20,:,bi,bj)
225             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
226       &       SQUEEZE_RIGHT , myThid)       &       SQUEEZE_RIGHT , myThid)
227             WRITE(msgBuf,HlimitMsgFormat)             WRITE(msgBuf,HlimitMsgFormat)
228       &       ' SEAICE_MODEL: AREAITD before advection: ',       &       ' SEAICE_MODEL: AREAITD         before advection: ',
229       &       AREAITD(20,20,:,bi,bj)       &       AREAITD(20,20,:,bi,bj)
230             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
231       &       SQUEEZE_RIGHT , myThid)       &       SQUEEZE_RIGHT , myThid)
# Line 273  C     test if REDIST works (1/2): Line 246  C     test if REDIST works (1/2):
246         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
247          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
248             WRITE(msgBuf,HlimitMsgFormat)             WRITE(msgBuf,HlimitMsgFormat)
249       &       ' SEAICE_MODEL: HEFFITD before 1. sorting: ',       &       ' SEAICE_MODEL: HEFFITD/AREAITD before 1. sorting: ',
250       &       HEFFITD(20,20,:,bi,bj)       &       HEFFITD(20,20,:,bi,bj) / AREAITD(20,20,:,bi,bj)
251             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
252       &       SQUEEZE_RIGHT , myThid)       &       SQUEEZE_RIGHT , myThid)
253             WRITE(msgBuf,HlimitMsgFormat)             WRITE(msgBuf,HlimitMsgFormat)
254       &       ' SEAICE_MODEL: AREAITD before 1. sorting: ',       &       ' SEAICE_MODEL: AREAITD before         1. sorting: ',
255       &       AREAITD(20,20,:,bi,bj)       &       AREAITD(20,20,:,bi,bj)
256             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
257       &       SQUEEZE_RIGHT , myThid)       &       SQUEEZE_RIGHT , myThid)
# Line 329  C     test if GROWTH works: Line 302  C     test if GROWTH works:
302          DO bj=myByLo(myThid),myByHi(myThid)          DO bj=myByLo(myThid),myByHi(myThid)
303           DO bi=myBxLo(myThid),myBxHi(myThid)           DO bi=myBxLo(myThid),myBxHi(myThid)
304              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
305       &        ' SEAICE_MODEL: HEFFITD after growth: ',       &       ' SEAICE_MODEL: HEFFITD/AREAITD after growth: ',
306       &        HEFFITD(20,20,:,bi,bj)       &       HEFFITD(20,20,:,bi,bj) / AREAITD(20,20,:,bi,bj)
307              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
308       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
309              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
310       &        ' SEAICE_MODEL: AREAITD after growth: ',       &        ' SEAICE_MODEL: AREAITD        after growth: ',
311       &        AREAITD(20,20,:,bi,bj)       &        AREAITD(20,20,:,bi,bj)
312              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
313       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
314              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
315       &        ' SEAICE_MODEL: HSNOWITD after growth: ',       &        ' SEAICE_MODEL: HSNOWITD       after growth: ',
316       &        HSNOWITD(20,20,:,bi,bj)       &        HSNOWITD(20,20,:,bi,bj)
317              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
318       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
# Line 361  C     test if REDIST works (2/2): Line 334  C     test if REDIST works (2/2):
334          DO bj=myByLo(myThid),myByHi(myThid)          DO bj=myByLo(myThid),myByHi(myThid)
335           DO bi=myBxLo(myThid),myBxHi(myThid)           DO bi=myBxLo(myThid),myBxHi(myThid)
336              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
337       &        ' SEAICE_MODEL: HEFFITD after 2. sorting: ',       &        ' SEAICE_MODEL: HEFFITD/AREAITD after 2. sorting: ',
338       &        HEFFITD(20,20,:,bi,bj)       &        HEFFITD(20,20,:,bi,bj) / AREAITD(20,20,:,bi,bj)
339              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
340       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
341              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
342       &        ' SEAICE_MODEL: AREAITD after 2. sorting: ',       &        ' SEAICE_MODEL: AREAITD         after 2. sorting: ',
343       &        AREAITD(20,20,:,bi,bj)       &        AREAITD(20,20,:,bi,bj)
344              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
345       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
346              WRITE(msgBuf,HlimitMsgFormat)              WRITE(msgBuf,HlimitMsgFormat)
347       &        ' SEAICE_MODEL: HSNOWITD after 2. sorting: ',       &        ' SEAICE_MODEL: HSNOWITD        after 2. sorting: ',
348       &        HSNOWITD(20,20,:,bi,bj)       &        HSNOWITD(20,20,:,bi,bj)
349              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
350       &        SQUEEZE_RIGHT , myThid)       &        SQUEEZE_RIGHT , myThid)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22