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

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

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

revision 1.1 by torge, Wed Jun 13 19:33:51 2012 UTC revision 1.3 by heimbach, Mon Oct 22 20:06:27 2012 UTC
# Line 57  C     == Local variables == Line 57  C     == Local variables ==
57  #endif  #endif
58  #ifdef ALLOW_TIMEAVE  #ifdef ALLOW_TIMEAVE
59        INTEGER i, j, bi, bj        INTEGER i, j, bi, bj
60    #ifdef SEAICE_ITD
61          INTEGER k
62    #endif
63        LOGICAL dumpFiles        LOGICAL dumpFiles
64  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
65  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
# Line 98  CEOP Line 101  CEOP
101              CALL MNC_CW_RL_W(pf,'sice',0,0,'si_VICE',vIce,myThid)              CALL MNC_CW_RL_W(pf,'sice',0,0,'si_VICE',vIce,myThid)
102              CALL MNC_CW_RL_W(pf,'sice',0,0,'si_HEFF',hEff,myThid)              CALL MNC_CW_RL_W(pf,'sice',0,0,'si_HEFF',hEff,myThid)
103              CALL MNC_CW_RL_W(pf,'sice',0,0,'si_AREA',area,myThid)              CALL MNC_CW_RL_W(pf,'sice',0,0,'si_AREA',area,myThid)
104    #ifdef SEAICE_ITD
105                CALL MNC_CW_RL_W(pf,'sice',0,0,'si_HEFFITD',HEFFITD,myThid)
106                CALL MNC_CW_RL_W(pf,'sice',0,0,'si_AREAITD',AREAITD,myThid)
107    #endif
108              CALL MNC_CW_RL_W(pf,'sice',0,0,'si_UWIND',uwind,myThid)              CALL MNC_CW_RL_W(pf,'sice',0,0,'si_UWIND',uwind,myThid)
109              CALL MNC_CW_RL_W(pf,'sice',0,0,'si_VWIND',vwind,myThid)              CALL MNC_CW_RL_W(pf,'sice',0,0,'si_VWIND',vwind,myThid)
110              CALL MNC_CW_RS_W(pf,'sice',0,0,'fu',fu,myThid)              CALL MNC_CW_RS_W(pf,'sice',0,0,'fu',fu,myThid)
# Line 123  CEOP Line 130  CEOP
130              CALL WRITE_FLD_XY_RL( 'HEFF.',suff,hEff,myIter,myThid)              CALL WRITE_FLD_XY_RL( 'HEFF.',suff,hEff,myIter,myThid)
131              CALL WRITE_FLD_XY_RL( 'AREA.',suff,area,myIter,myThid)              CALL WRITE_FLD_XY_RL( 'AREA.',suff,area,myIter,myThid)
132              CALL WRITE_FLD_XY_RL( 'HSNOW.',suff,HSNOW,myIter,myThid)              CALL WRITE_FLD_XY_RL( 'HSNOW.',suff,HSNOW,myIter,myThid)
133    #ifdef SEAICE_ITD
134                CALL WRITE_FLD_XY_RL( 'HEFFITD.',suff,HEFFITD,myIter,myThid)
135                CALL WRITE_FLD_XY_RL( 'AREAITD.',suff,AREAITD,myIter,myThid)
136               CALL WRITE_FLD_XY_RL('HSNOWITD.',suff,HSNOWITD,myIter,myThid)
137    #endif
138  #ifdef SEAICE_VARIABLE_SALINITY  #ifdef SEAICE_VARIABLE_SALINITY
139              CALL WRITE_FLD_XY_RL( 'HSALT.',suff,HSALT,myIter,myThid)              CALL WRITE_FLD_XY_RL( 'HSALT.',suff,HSALT,myIter,myThid)
140  #endif  #endif
# Line 184  C- note(jmc): surf.Fluxes have not yet b Line 196  C- note(jmc): surf.Fluxes have not yet b
196       &         AREAtave(i,j,bi,bj) +AREA(i,j,bi,bj)*deltaTclock       &         AREAtave(i,j,bi,bj) +AREA(i,j,bi,bj)*deltaTclock
197            ENDDO            ENDDO
198           ENDDO           ENDDO
199    #ifdef SEAICE_ITD
200             DO k=1,nITD
201              DO j=1,sNy
202               DO i=1,sNx
203                HEFFITDtave(i,j,k,bi,bj) =
204         &        HEFFITDtave(i,j,k,bi,bj) +HEFFITD(i,j,k,bi,bj)*deltaTclock
205                AREAITDtave(i,j,k,bi,bj) =
206         &        AREAITDtave(i,j,k,bi,bj) +AREAITD(i,j,k,bi,bj)*deltaTclock
207               ENDDO
208              ENDDO
209             ENDDO
210    #endif
211           SEAICE_timeAve(bi,bj) = SEAICE_timeAve(bi,bj)+deltaTclock           SEAICE_timeAve(bi,bj) = SEAICE_timeAve(bi,bj)+deltaTclock
212          ENDDO          ENDDO
213         ENDDO         ENDDO
# Line 195  C     Dump files and restart average com Line 219  C     Dump files and restart average com
219       &     DIFFERENT_MULTIPLE(SEAICE_taveFreq,myTime,deltaTClock)       &     DIFFERENT_MULTIPLE(SEAICE_taveFreq,myTime,deltaTClock)
220  #ifdef ALLOW_CAL  #ifdef ALLOW_CAL
221          IF ( useCAL ) THEN          IF ( useCAL ) THEN
222            CALL CAL_TIME2DUMP( SEAICE_taveFreq, deltaTClock,            CALL CAL_TIME2DUMP( ZERO, SEAICE_taveFreq, deltaTClock,
223       U                        dumpFiles,       U                        dumpFiles,
224       I                        myTime, myIter, myThid )       I                        myTime, myIter, myThid )
225          ENDIF          ENDIF
# Line 224  C      Normalize by integrated time Line 248  C      Normalize by integrated time
248       &                            SEAICE_timeAve, 1, bi, bj, myThid )       &                            SEAICE_timeAve, 1, bi, bj, myThid )
249            CALL TIMEAVE_NORMALIZE( AREAtave,            CALL TIMEAVE_NORMALIZE( AREAtave,
250       &                            SEAICE_timeAve, 1, bi, bj, myThid )       &                            SEAICE_timeAve, 1, bi, bj, myThid )
251    #ifdef SEAICE_ITD
252              CALL TIMEAVE_NORMALIZE( HEFFITDtave,
253         &                            SEAICE_timeAve, nITD, bi, bj, myThid )
254              CALL TIMEAVE_NORMALIZE( AREAITDtave,
255         &                            SEAICE_timeAve, nITD, bi, bj, myThid )
256    #endif
257           ENDDO           ENDDO
258          ENDDO          ENDDO
259  c       IF (myIter.EQ.10) WRITE(0,*) myThid, dumpFiles  c       IF (myIter.EQ.10) WRITE(0,*) myThid, dumpFiles
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22