/[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.2 - (hide annotations) (download)
Wed Jun 20 19:58:09 2012 UTC (13 years, 1 month ago) by torge
Branch: MAIN
Changes since 1.1: +30 -0 lines
introducing 3-D ITD output fields SIheffN and SIareaN
containing ice thickness and fractinal area of each ice thickness category

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

  ViewVC Help
Powered by ViewVC 1.1.22