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

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

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


Revision 1.4 - (hide annotations) (download)
Mon Dec 10 22:19:49 2012 UTC (12 years, 7 months ago) by torge
Branch: MAIN
Changes since 1.3: +13 -8 lines
include updates from main branch

1 torge 1.4 C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_output.F,v 1.15 2012/11/15 00:07:59 jmc Exp $
2 torge 1.1 C $Name: $
3    
4     #include "SEAICE_OPTIONS.h"
5 torge 1.4 #ifdef ALLOW_EXF
6     # include "EXF_OPTIONS.h"
7     #endif
8 torge 1.1
9     CBOP 0
10     C !ROUTINE: SEAICE_OUTPUT
11    
12     C !INTERFACE:
13     SUBROUTINE SEAICE_OUTPUT( myTime, myIter, myThid )
14    
15     C !DESCRIPTION:
16     C *==========================================================*
17     C | SUBROUTINE SEAICE_OUTPUT
18     C | o Do SEAICE output and TimeAve averaging and output.
19     C *==========================================================*
20    
21     C !USES:
22     IMPLICIT NONE
23    
24     C === Global variables ===
25     #include "SIZE.h"
26     #include "EEPARAMS.h"
27     #include "PARAMS.h"
28     #include "FFIELDS.h"
29     #include "SEAICE_SIZE.h"
30     #include "SEAICE_PARAMS.h"
31     #include "SEAICE.h"
32     #include "SEAICE_TRACER.h"
33     #include "SEAICE_TAVE.h"
34     #ifdef ALLOW_EXF
35     # include "EXF_FIELDS.h"
36     #endif
37    
38     C !INPUT PARAMETERS:
39     C == Routine arguments ==
40     C myTime :: my time in simulation ( s )
41     C myIter :: my Iteration number
42     C myThid :: my Thread Id number
43     _RL myTime
44     INTEGER myIter
45     INTEGER myThid
46    
47     C !FUNCTIONS:
48     LOGICAL DIFFERENT_MULTIPLE
49     EXTERNAL DIFFERENT_MULTIPLE
50    
51     C !LOCAL VARIABLES:
52     C == Local variables in common block :
53    
54     C == Local variables ==
55     CHARACTER*(10) suff
56     #ifdef ALLOW_SITRACER
57     CHARACTER*(13) pref
58     INTEGER iTracer
59     #endif
60     #ifdef ALLOW_TIMEAVE
61     INTEGER i, j, bi, bj
62 torge 1.2 #ifdef SEAICE_ITD
63     INTEGER k
64     #endif
65 torge 1.1 LOGICAL dumpFiles
66     #endif /* ALLOW_TIMEAVE */
67     #ifdef ALLOW_MNC
68     CHARACTER*(1) pf
69     #endif /* ALLOW_MNC */
70     CEOP
71    
72     IF (SEAICEwriteState) THEN
73    
74     IF ( DIFFERENT_MULTIPLE(SEAICE_dumpFreq,myTime,deltaTClock)
75     & .OR. dumpInitAndLast.AND.( myTime.EQ.endTime .OR.
76     & myTime.EQ.startTime )
77     & ) THEN
78    
79     #ifdef ALLOW_MNC
80     IF ( useMNC .AND. SEAICE_dump_mnc ) THEN
81     IF ( writeBinaryPrec .EQ. precFloat64 ) THEN
82     pf(1:1) = 'D'
83     ELSE
84     pf(1:1) = 'R'
85     ENDIF
86     CALL MNC_CW_SET_UDIM('sice', -1, myThid)
87     CALL MNC_CW_RL_W_S('D','sice',0,0,'T', myTime, myThid)
88     CALL MNC_CW_SET_UDIM('sice', 0, myThid)
89     CALL MNC_CW_I_W_S('I','sice',0,0,'iter', myIter, myThid)
90     CALL MNC_CW_RL_W_S('D','sice',0,0,'model_time',
91     & myTime,myThid)
92     #if ( defined(SEAICE_CGRID) && defined(SEAICE_ALLOW_EVP) )
93     IF ( SEAICEuseEVP ) THEN
94     CALL MNC_CW_RL_W(pf,'sice',0,0,'si_sigma1',
95     & seaice_sigma1,myThid)
96     CALL MNC_CW_RL_W(pf,'sice',0,0,'si_sigma2',
97     & seaice_sigma2,myThid)
98     CALL MNC_CW_RL_W(pf,'sice',0,0,'si_sigma12',
99     & seaice_sigma12,myThid)
100     ENDIF
101     #endif /* SEAICE_CGRID and SEAICE_ALLOW_EVP */
102     CALL MNC_CW_RL_W(pf,'sice',0,0,'si_UICE',uIce,myThid)
103     CALL MNC_CW_RL_W(pf,'sice',0,0,'si_VICE',vIce,myThid)
104     CALL MNC_CW_RL_W(pf,'sice',0,0,'si_HEFF',hEff,myThid)
105     CALL MNC_CW_RL_W(pf,'sice',0,0,'si_AREA',area,myThid)
106 torge 1.2 #ifdef SEAICE_ITD
107 torge 1.4 c CALL MNC_CW_RL_W(pf,'sice',0,0,'si_HEFFITD',HEFFITD,myThid)
108     c CALL MNC_CW_RL_W(pf,'sice',0,0,'si_AREAITD',AREAITD,myThid)
109 torge 1.2 #endif
110 torge 1.1 CALL MNC_CW_RL_W(pf,'sice',0,0,'si_UWIND',uwind,myThid)
111     CALL MNC_CW_RL_W(pf,'sice',0,0,'si_VWIND',vwind,myThid)
112     CALL MNC_CW_RS_W(pf,'sice',0,0,'fu',fu,myThid)
113     CALL MNC_CW_RS_W(pf,'sice',0,0,'fv',fv,myThid)
114     CALL MNC_CW_RS_W(pf,'sice',0,0,'EmPmR',EmPmR,myThid)
115     CALL MNC_CW_RS_W(pf,'sice',0,0,'Qnet',Qnet,myThid)
116     CALL MNC_CW_RS_W(pf,'sice',0,0,'Qsw',Qsw,myThid)
117     ENDIF
118     #endif /* ALLOW_MNC */
119     IF (SEAICE_dump_mdsio) THEN
120     WRITE(suff,'(I10.10)') myIter
121     IF ( myIter.NE.nIter0 ) THEN
122     CALL WRITE_FLD_XY_RL('UWIND.',suff,uwind,myIter,myThid)
123     CALL WRITE_FLD_XY_RL('VWIND.',suff,vwind,myIter,myThid)
124     CALL WRITE_FLD_XY_RS('FU.',suff,fu,myIter,myThid)
125     CALL WRITE_FLD_XY_RS('FV.',suff,fv,myIter,myThid)
126     CALL WRITE_FLD_XY_RS('EmPmR.',suff,EmPmR,myIter,myThid)
127     CALL WRITE_FLD_XY_RS('Qnet.',suff,Qnet,myIter,myThid)
128     CALL WRITE_FLD_XY_RS('Qsw.',suff,Qsw,myIter,myThid)
129     ENDIF
130     CALL WRITE_FLD_XY_RL( 'UICE.',suff,uIce,myIter,myThid)
131     CALL WRITE_FLD_XY_RL( 'VICE.',suff,vIce,myIter,myThid)
132     CALL WRITE_FLD_XY_RL( 'HEFF.',suff,hEff,myIter,myThid)
133     CALL WRITE_FLD_XY_RL( 'AREA.',suff,area,myIter,myThid)
134     CALL WRITE_FLD_XY_RL( 'HSNOW.',suff,HSNOW,myIter,myThid)
135 torge 1.2 #ifdef SEAICE_ITD
136 torge 1.4 CALL WRITE_FLD_3D_RL( 'HEFFITD.', suff, nITD, HEFFITD,
137     & myIter, myThid )
138     CALL WRITE_FLD_3D_RL( 'AREAITD.', suff, nITD, AREAITD,
139     & myIter, myThid )
140     CALL WRITE_FLD_3D_RL('HSNOWITD.', suff, nITD, HSNOWITD,
141     & myIter, myThid )
142 torge 1.2 #endif
143 torge 1.1 #ifdef SEAICE_VARIABLE_SALINITY
144     CALL WRITE_FLD_XY_RL( 'HSALT.',suff,HSALT,myIter,myThid)
145     #endif
146     #ifdef ALLOW_SITRACER
147     DO iTracer = 1, SItrNumInUse
148     WRITE(pref,'(A8,I2.2,A1)') 'SITRACER', iTracer, '.'
149     CALL WRITE_FLD_XY_RL(
150     & pref,suff,
151     & SItracer(1-OLx,1-OLy,1,1,iTracer),
152     & myIter,myThid)
153     ENDDO
154     #endif
155     #if ( defined(SEAICE_CGRID) && defined(SEAICE_ALLOW_EVP) )
156     IF ( SEAICEuseEVP ) THEN
157     CALL WRITE_FLD_XY_RL('SIGMA1.',suff,seaice_sigma1,
158     & myIter,myThid)
159     CALL WRITE_FLD_XY_RL('SIGMA2.',suff,seaice_sigma2,
160     & myIter,myThid)
161     CALL WRITE_FLD_XY_RL('SIGMA12.',suff,seaice_sigma12,
162     & myIter,myThid)
163     ENDIF
164     #endif /* SEAICE_CGRID and SEAICE_ALLOW_EVP */
165     C-- end SEAICE_dump_mdsio block
166     ENDIF
167    
168     ENDIF
169     ENDIF
170    
171     C----------------------------------------------------------------
172     C Do SEAICE time averaging.
173     C----------------------------------------------------------------
174    
175     #ifdef ALLOW_TIMEAVE
176     IF ( SEAICE_taveFreq.GT.0. _d 0 ) THEN
177    
178     C-- Time-cumulations
179     DO bj = myByLo(myThid), myByHi(myThid)
180     DO bi = myBxLo(myThid), myBxHi(myThid)
181     DO j=1,sNy
182     DO i=1,sNx
183     C- note(jmc): surf.Fluxes have not yet been computed when called @ nIter0
184     FUtave(i,j,bi,bj) =
185     & FUtave(i,j,bi,bj) +FU(i,j,bi,bj) *deltaTclock
186     FVtave(i,j,bi,bj) =
187     & FVtave(i,j,bi,bj) +FV(i,j,bi,bj) *deltaTclock
188     EmPmRtave(i,j,bi,bj)=
189     & EmPmRtave(i,j,bi,bj)+EmPmR(i,j,bi,bj) *deltaTclock
190     QNETtave(i,j,bi,bj) =
191     & QNETtave(i,j,bi,bj) +QNET(i,j,bi,bj) *deltaTclock
192     QSWtave(i,j,bi,bj) =
193     & QSWtave(i,j,bi,bj) +QSW(i,j,bi,bj) *deltaTclock
194     UICEtave(i,j,bi,bj) =
195     & UICEtave(i,j,bi,bj) +UICE(i,j,bi,bj)*deltaTclock
196     VICEtave(i,j,bi,bj) =
197     & VICEtave(i,j,bi,bj) +VICE(i,j,bi,bj)*deltaTclock
198     HEFFtave(i,j,bi,bj) =
199     & HEFFtave(i,j,bi,bj) +HEFF(i,j,bi,bj)*deltaTclock
200     AREAtave(i,j,bi,bj) =
201     & AREAtave(i,j,bi,bj) +AREA(i,j,bi,bj)*deltaTclock
202     ENDDO
203     ENDDO
204 torge 1.2 #ifdef SEAICE_ITD
205 torge 1.4 DO k=1,nITD
206 torge 1.2 DO j=1,sNy
207     DO i=1,sNx
208     HEFFITDtave(i,j,k,bi,bj) =
209     & HEFFITDtave(i,j,k,bi,bj) +HEFFITD(i,j,k,bi,bj)*deltaTclock
210     AREAITDtave(i,j,k,bi,bj) =
211     & AREAITDtave(i,j,k,bi,bj) +AREAITD(i,j,k,bi,bj)*deltaTclock
212     ENDDO
213     ENDDO
214     ENDDO
215     #endif
216 torge 1.1 SEAICE_timeAve(bi,bj) = SEAICE_timeAve(bi,bj)+deltaTclock
217     ENDDO
218     ENDDO
219    
220     C Dump files and restart average computation if needed
221     dumpFiles = .FALSE.
222     IF ( myIter .NE. nIter0 ) THEN
223     dumpFiles =
224     & DIFFERENT_MULTIPLE(SEAICE_taveFreq,myTime,deltaTClock)
225     #ifdef ALLOW_CAL
226     IF ( useCAL ) THEN
227 heimbach 1.3 CALL CAL_TIME2DUMP( ZERO, SEAICE_taveFreq, deltaTClock,
228 torge 1.1 U dumpFiles,
229     I myTime, myIter, myThid )
230     ENDIF
231     #endif
232     ENDIF
233    
234     IF (dumpFiles) THEN
235     C Normalize by integrated time
236     DO bj = myByLo(myThid), myByHi(myThid)
237     DO bi = myBxLo(myThid), myBxHi(myThid)
238     CALL TIMEAVE_NORMALIZE( FUtave,
239     & SEAICE_timeAve, 1, bi, bj, myThid )
240     CALL TIMEAVE_NORMALIZE( FVtave,
241     & SEAICE_timeAve, 1, bi, bj, myThid )
242     CALL TIMEAVE_NORMALIZE( EmPmRtave,
243     & SEAICE_timeAve, 1, bi, bj, myThid )
244     CALL TIMEAVE_NORMALIZE( QNETtave,
245     & SEAICE_timeAve, 1, bi, bj, myThid )
246     CALL TIMEAVE_NORMALIZE( QSWtave,
247     & SEAICE_timeAve, 1, bi, bj, myThid )
248     CALL TIMEAVE_NORMALIZE( UICEtave,
249     & SEAICE_timeAve, 1, bi, bj, myThid )
250     CALL TIMEAVE_NORMALIZE( VICEtave,
251     & SEAICE_timeAve, 1, bi, bj, myThid )
252     CALL TIMEAVE_NORMALIZE( HEFFtave,
253     & SEAICE_timeAve, 1, bi, bj, myThid )
254     CALL TIMEAVE_NORMALIZE( AREAtave,
255     & SEAICE_timeAve, 1, bi, bj, myThid )
256 torge 1.2 #ifdef SEAICE_ITD
257     CALL TIMEAVE_NORMALIZE( HEFFITDtave,
258     & SEAICE_timeAve, nITD, bi, bj, myThid )
259     CALL TIMEAVE_NORMALIZE( AREAITDtave,
260     & SEAICE_timeAve, nITD, bi, bj, myThid )
261     #endif
262 torge 1.1 ENDDO
263     ENDDO
264     c IF (myIter.EQ.10) WRITE(0,*) myThid, dumpFiles
265    
266     #ifdef ALLOW_MNC
267     IF (useMNC .AND. SEAICE_tave_mnc) THEN
268     IF ( writeBinaryPrec .EQ. precFloat64 ) THEN
269     pf(1:1) = 'D'
270     ELSE
271     pf(1:1) = 'R'
272     ENDIF
273     CALL MNC_CW_SET_UDIM('sice_tave', -1, myThid)
274     CALL MNC_CW_RL_W_S('D','sice_tave',0,0,'T', myTime, myThid)
275     CALL MNC_CW_SET_UDIM('sice_tave', 0, myThid)
276     CALL MNC_CW_I_W_S('I','sice_tave',0,0,'iter', myIter, myThid)
277     C CALL MNC_CW_RL_W_S('D','sice_tave',0,0,'model_time',
278     C & myTime,myThid)
279     CALL MNC_CW_RL_W(pf,'sice_tave',0,0,
280     & 'si_UICEtave',UICEtave,myThid)
281     CALL MNC_CW_RL_W(pf,'sice_tave',0,0,
282     & 'si_VICEtave',VICEtave,myThid)
283     CALL MNC_CW_RL_W(pf,'sice_tave',0,0,
284     & 'si_FUtave',FUtave,myThid)
285     CALL MNC_CW_RL_W(pf,'sice_tave',0,0,
286     & 'si_FVtave',FVtave,myThid)
287     CALL MNC_CW_RL_W(pf,'sice_tave',0,0,
288     & 'si_EmPmRtave',EmPmRtave,myThid)
289     CALL MNC_CW_RL_W(pf,'sice_tave',0,0,
290     & 'si_QNETtave',QNETtave,myThid)
291     CALL MNC_CW_RL_W(pf,'sice_tave',0,0,
292     & 'si_QSWtave',QSWtave,myThid)
293     CALL MNC_CW_RL_W(pf,'sice_tave',0,0,
294     & 'si_HEFFtave',HEFFtave,myThid)
295     CALL MNC_CW_RL_W(pf,'sice_tave',0,0,
296     & 'si_AREAtave',AREAtave,myThid)
297     ENDIF
298     #endif
299     IF (SEAICE_tave_mdsio) THEN
300     WRITE(suff,'(I10.10)') myIter
301     CALL WRITE_FLD_XY_RL('FUtave.' ,suff,FUtave ,myIter,myThid)
302     CALL WRITE_FLD_XY_RL('FVtave.' ,suff,FVtave ,myIter,myThid)
303     CALL WRITE_FLD_XY_RL('EmPmRtave.',suff,EmPmRtave,myIter,myThid)
304     CALL WRITE_FLD_XY_RL('QNETtave.' ,suff,QNETtave ,myIter,myThid)
305     CALL WRITE_FLD_XY_RL('QSWtave.' ,suff,QSWtave ,myIter,myThid)
306     CALL WRITE_FLD_XY_RL('UICEtave.' ,suff,UICEtave ,myIter,myThid)
307     CALL WRITE_FLD_XY_RL('VICEtave.' ,suff,VICEtave ,myIter,myThid)
308     CALL WRITE_FLD_XY_RL('HEFFtave.' ,suff,HEFFtave ,myIter,myThid)
309     CALL WRITE_FLD_XY_RL('AREAtave.' ,suff,AREAtave ,myIter,myThid)
310     ENDIF
311    
312     C Reset averages to zero
313     DO bj = myByLo(myThid), myByHi(myThid)
314     DO bi = myBxLo(myThid), myBxHi(myThid)
315     CALL TIMEAVE_RESET( FUtave , 1, bi, bj, myThid )
316     CALL TIMEAVE_RESET( FVtave , 1, bi, bj, myThid )
317     CALL TIMEAVE_RESET( EmPmRtave, 1, bi, bj, myThid )
318     CALL TIMEAVE_RESET( QNETtave , 1, bi, bj, myThid )
319     CALL TIMEAVE_RESET( QSWtave , 1, bi, bj, myThid )
320     CALL TIMEAVE_RESET( UICEtave , 1, bi, bj, myThid )
321     CALL TIMEAVE_RESET( VICEtave , 1, bi, bj, myThid )
322     CALL TIMEAVE_RESET( HEFFtave , 1, bi, bj, myThid )
323     CALL TIMEAVE_RESET( AREAtave , 1, bi, bj, myThid )
324     SEAICE_timeAve(bi,bj) = ZERO
325     ENDDO
326     ENDDO
327    
328     C-- end dumpFiles block
329     ENDIF
330    
331     C-- end if SEAICE_taveFreq > 0
332     ENDIF
333     #endif /* ALLOW_TIMEAVE */
334    
335     C-- do SEAICE monitor output : print some statistics about seaice fields
336     CALL SEAICE_MONITOR( myTime, myIter, myThid )
337    
338     C-- do SEAICE Open-Boundary output
339     IF ( useOBCS ) CALL SEAICE_OBCS_OUTPUT( myTime, myIter, myThid )
340    
341     RETURN
342     END

  ViewVC Help
Powered by ViewVC 1.1.22