/[MITgcm]/MITgcm/pkg/seaice/seaice_do_diags.F
ViewVC logotype

Contents of /MITgcm/pkg/seaice/seaice_do_diags.F

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

Revision 1.7 - (show annotations) (download)
Thu Oct 9 04:19:20 2003 UTC (20 years, 7 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint53f_post, checkpoint54a_pre, checkpoint55c_post, checkpoint53b_pre, checkpoint51l_post, checkpoint52e_pre, hrcube4, checkpoint52n_post, checkpoint52j_post, checkpoint51o_pre, checkpoint52e_post, checkpoint51n_pre, checkpoint57b_post, checkpoint52d_pre, checkpoint53c_post, checkpoint53d_post, checkpoint55d_pre, checkpoint55j_post, checkpoint56b_post, checkpoint51q_post, checkpoint52j_pre, checkpoint54a_post, branch-netcdf, checkpoint52l_post, checkpoint55h_post, checkpoint51r_post, checkpoint52k_post, checkpoint52b_pre, checkpoint54b_post, checkpoint51i_post, checkpoint54d_post, checkpoint52l_pre, checkpoint56c_post, checkpoint54e_post, checkpoint55b_post, checkpoint51l_pre, checkpoint52m_post, checkpoint55, checkpoint53a_post, checkpoint55a_post, checkpoint57a_post, hrcube5, checkpoint55g_post, checkpoint51o_post, checkpoint55f_post, checkpoint57a_pre, checkpoint54, checkpoint57, checkpoint56, checkpoint53, checkpoint52, checkpoint52d_post, checkpoint52a_post, checkpoint52b_post, checkpoint53g_post, checkpoint52f_post, checkpoint52c_post, ecco_c52_e35, checkpoint54f_post, checkpoint55e_post, checkpoint53b_post, checkpoint52a_pre, checkpoint52i_post, checkpoint51m_post, checkpoint51t_post, checkpoint53d_pre, checkpoint51p_post, checkpoint51n_post, checkpoint55i_post, checkpoint54c_post, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint56a_post, checkpoint51s_post, checkpoint55d_post
Branch point for: branch-nonh, netcdf-sm0, tg2-branch, checkpoint51n_branch
Changes since 1.6: +2 -1 lines 
 o first check-in for the "branch-genmake2" merge
 o verification suite as run on shelley (gcc 3.2.2):

Wed Oct  8 23:42:29 EDT 2003
                T           S           U           V
G D M    c        m  s        m  s        m  s        m  s
E p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
N n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
2 d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .

OPTFILE=NONE

Y Y Y Y 13 16 16 16  0 16 16 16 16 16 16 16 16 13 12  0  0 pass  adjustment.128x64x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16  0  0 16 16  0  0 pass  adjustment.cs-32x32x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16 22  0 16 16 22  0 pass  adjust_nlfs.cs-32x32x1
Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O   advect_cs
Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O   advect_xy
Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O   advect_xz
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  aim.5l_cs
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass  aim.5l_Equatorial_Channel
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass  aim.5l_LatLon
Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass  exp0
Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass  exp1
Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass  exp2
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  exp4
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass  exp5
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  front_relax
Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass  global_ocean.90x40x15
Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL  global_ocean.cs32x15
Y Y Y Y  6 11 12 13 13 12 13 16 13  9  9  9  9 10  9  9 11 FAIL  global_ocean_pressure
Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass  global_with_exf
Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16  9 16 pass  hs94.128x64x5
Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass  hs94.1x64x5
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass  hs94.cs-32x32x5
Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL  ideal_2D_oce
Y Y Y Y  8 16 16 16 16 16 16 16 16 13 13  8 16 16 16 16 16 FAIL  internal_wave
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass  inverted_barometer
Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL  lab_sea
Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL  natl_box
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  plume_on_slope
Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass  solid-body.cs-32x32x1
1 C $Header: /u/u3/gcmpack/MITgcm/pkg/seaice/seaice_do_diags.F,v 1.6.2.1 2003/10/02 18:18:34 adcroft Exp $
2 C $Name: $
3
4 #include "SEAICE_OPTIONS.h"
5
6 SUBROUTINE SEAICE_DO_DIAGS( myTime, myIter, myThid )
7 C /==========================================================\
8 C | SUBROUTINE SEAICE_DO_DIAGS |
9 C | o Do SEAICE diagnostic output. |
10 C \==========================================================/
11 IMPLICIT NONE
12
13 C === Global variables ===
14 #include "SIZE.h"
15 #include "EEPARAMS.h"
16 #include "PARAMS.h"
17 #include "FFIELDS.h"
18 #include "SEAICE_DIAGS.h"
19 #include "SEAICE_PARAMS.h"
20 #include "SEAICE_FFIELDS.h"
21 #include "SEAICE.h"
22
23 C == Routine arguments ==
24 C myTime - Current time of simulation ( s )
25 C myIter - Iteration number
26 C myThid - Number of this instance of SEAICE_DO_DIAGS
27 _RL myTime
28 INTEGER myIter
29 INTEGER myThid
30
31 #ifdef ALLOW_SEAICE
32
33 C == Local variables ==
34 CHARACTER*(MAX_LEN_MBUF) suff
35 LOGICAL DIFFERENT_MULTIPLE
36 EXTERNAL DIFFERENT_MULTIPLE
37 INTEGER i, j, k, bi, bj
38 _RS arr(1-oLx:sNx+oLx,1-oLy:sNy+oLy,nSx,nSy)
39
40 IF (SEAICEwriteState) THEN
41
42 IF ( DIFFERENT_MULTIPLE(SEAICE_dumpFreq,myTime,
43 & myTime-deltaTClock) ) THEN
44
45 WRITE(suff,'(I10.10)') myIter
46 _BARRIER
47 _BEGIN_MASTER( myThid )
48 CALL WRITE_FLD_XY_RS( 'UWIND.',suff,uwind,myIter,myThid)
49 CALL WRITE_FLD_XY_RS( 'VWIND.',suff,vwind,myIter,myThid)
50 CALL WRITE_FLD_XY_RS( 'FU.',suff,fu,myIter,myThid)
51 CALL WRITE_FLD_XY_RS( 'FV.',suff,fv,myIter,myThid)
52 CALL WRITE_FLD_XY_RS( 'EmPmR.',suff,EmPmR,myIter,myThid)
53 CALL WRITE_FLD_XY_RS( 'Qnet.',suff,Qnet,myIter,myThid)
54 CALL WRITE_FLD_XY_RS( 'Qsw.',suff,Qsw,myIter,myThid)
55 _END_MASTER( myThid )
56 _BARRIER
57
58 #ifdef SEAICE_DEBUG
59 CALL PLOT_FIELD_XYRS( uwind , 'Current uwind ', myIter, myThid )
60 CALL PLOT_FIELD_XYRS( vwind , 'Current vwind ', myIter, myThid )
61 CALL PLOT_FIELD_XYRS( atemp , 'Current atemp ', myIter, myThid )
62 CALL PLOT_FIELD_XYRS( aqh , 'Current aqh ', myIter, myThid )
63 CALL PLOT_FIELD_XYRS( lwdown, 'Current lwdown', myIter, myThid )
64 CALL PLOT_FIELD_XYRS( swdown, 'Current swdown', myIter, myThid )
65 CALL PLOT_FIELD_XYRS( precip, 'Current precip', myIter, myThid )
66 CALL PLOT_FIELD_XYRL( evap , 'Current evap ', myIter, myThid )
67 CALL PLOT_FIELD_XYRS( runoff, 'Current runoff', myIter, myThid )
68 CALL PLOT_FIELD_XYRS( SSS , 'Current SSS ', myIter, myThid )
69 CALL PLOT_FIELD_XYRS( SST , 'Current SST ', myIter, myThid )
70 CALL PLOT_FIELD_XYRL( fu , 'Current fu ', myIter, myThid )
71 CALL PLOT_FIELD_XYRL( fv , 'Current fv ', myIter, myThid )
72 CALL PLOT_FIELD_XYRL( EmPmR , 'Current EmPmR ', myIter, myThid )
73 CALL PLOT_FIELD_XYRL( Qnet , 'Current Qnet ', myIter, myThid )
74 CALL PLOT_FIELD_XYRL( Qsw , 'Current Qsw ', myIter, myThid )
75 #endif
76
77 DO bj=myByLo(myThid),myByHi(myThid)
78 DO bi=myBxLo(myThid),myBxHi(myThid)
79 DO j=1,sNy
80 DO i=1,sNx
81 arr(i,j,bi,bj)=UICE(i,j,1,bi,bj)
82 ENDDO
83 ENDDO
84 ENDDO
85 ENDDO
86 _BARRIER
87 _BEGIN_MASTER( myThid )
88 CALL WRITE_FLD_XY_RS( 'UICE.',suff,arr,myIter,myThid)
89 _END_MASTER( myThid )
90 _BARRIER
91 #ifdef SEAICE_DEBUG
92 _EXCH_XY_R4( arr, myThid )
93 CALL PLOT_FIELD_XYRS( arr , 'Current uice ', myIter, myThid )
94 #endif
95
96 DO bj=myByLo(myThid),myByHi(myThid)
97 DO bi=myBxLo(myThid),myBxHi(myThid)
98 DO j=1,sNy
99 DO i=1,sNx
100 arr(i,j,bi,bj)=VICE(i,j,1,bi,bj)
101 ENDDO
102 ENDDO
103 ENDDO
104 ENDDO
105 _BARRIER
106 _BEGIN_MASTER( myThid )
107 CALL WRITE_FLD_XY_RS( 'VICE.',suff,arr,myIter,myThid)
108 _END_MASTER( myThid )
109 _BARRIER
110 #ifdef SEAICE_DEBUG
111 _EXCH_XY_R4( arr, myThid )
112 CALL PLOT_FIELD_XYRS( arr , 'Current vice ', myIter, myThid )
113 #endif
114
115 DO bj=myByLo(myThid),myByHi(myThid)
116 DO bi=myBxLo(myThid),myBxHi(myThid)
117 DO j=1,sNy
118 DO i=1,sNx
119 arr(i,j,bi,bj)=HEFF(i,j,1,bi,bj)
120 ENDDO
121 ENDDO
122 ENDDO
123 ENDDO
124 _BARRIER
125 _BEGIN_MASTER( myThid )
126 CALL WRITE_FLD_XY_RS( 'HEFF.',suff,arr,myIter,myThid)
127 _END_MASTER( myThid )
128 _BARRIER
129 #ifdef SEAICE_DEBUG
130 _EXCH_XY_R4( arr, myThid )
131 CALL PLOT_FIELD_XYRS( arr , 'Current heff ', myIter, myThid )
132 #endif
133
134 DO bj=myByLo(myThid),myByHi(myThid)
135 DO bi=myBxLo(myThid),myBxHi(myThid)
136 DO j=1,sNy
137 DO i=1,sNx
138 arr(i,j,bi,bj)=AREA(i,j,1,bi,bj)
139 ENDDO
140 ENDDO
141 ENDDO
142 ENDDO
143 _BARRIER
144 _BEGIN_MASTER( myThid )
145 CALL WRITE_FLD_XY_RS( 'AREA.',suff,arr,myIter,myThid)
146 _END_MASTER( myThid )
147 _BARRIER
148 #ifdef SEAICE_DEBUG
149 _EXCH_XY_R4( arr, myThid )
150 CALL PLOT_FIELD_XYRS( arr , 'Current area ', myIter, myThid )
151 #endif
152
153 ENDIF
154 ENDIF
155
156 C----------------------------------------------------------------
157 C Do SEAICE time averaging.
158 C----------------------------------------------------------------
159
160 #ifdef ALLOW_TIMEAVE
161
162 C-- Time-cumulations
163 DO bj = myByLo(myThid), myByHi(myThid)
164 DO bi = myBxLo(myThid), myBxHi(myThid)
165 DO j=1,sNy
166 DO i=1,sNx
167 FUtave(i,j,1,bi,bj) =
168 & FUtave(i,j,1,bi,bj) +FU(i,j,bi,bj) *deltaTclock
169 FVtave(i,j,1,bi,bj) =
170 & FVtave(i,j,1,bi,bj) +FV(i,j,bi,bj) *deltaTclock
171 EmPmRtave(i,j,1,bi,bj)=
172 & EmPmRtave(i,j,1,bi,bj)+EmPmR(i,j,bi,bj) *deltaTclock
173 QNETtave(i,j,1,bi,bj) =
174 & QNETtave(i,j,1,bi,bj) +QNET(i,j,bi,bj) *deltaTclock
175 QSWtave(i,j,1,bi,bj) =
176 & QSWtave(i,j,1,bi,bj) +QSW(i,j,bi,bj) *deltaTclock
177 UICEtave(i,j,1,bi,bj) =
178 & UICEtave(i,j,1,bi,bj) +UICE(i,j,1,bi,bj)*deltaTclock
179 VICEtave(i,j,1,bi,bj) =
180 & VICEtave(i,j,1,bi,bj) +VICE(i,j,1,bi,bj)*deltaTclock
181 HEFFtave(i,j,1,bi,bj) =
182 & HEFFtave(i,j,1,bi,bj) +HEFF(i,j,1,bi,bj)*deltaTclock
183 AREAtave(i,j,1,bi,bj) =
184 & AREAtave(i,j,1,bi,bj) +AREA(i,j,1,bi,bj)*deltaTclock
185 ENDDO
186 ENDDO
187 DO k=1,Nr
188 SEAICE_TimeAve(k,bi,bj)=SEAICE_TimeAve(k,bi,bj)+deltaTclock
189 ENDDO
190 ENDDO
191 ENDDO
192
193 C Dump files and restart average computation if needed
194 IF ( myIter.NE.nIter0 .AND.
195 & DIFFERENT_MULTIPLE(SEAICE_taveFreq,myTime,myTime-deltaTClock)
196 & ) THEN
197
198 C Normalize by integrated time
199 DO bj = myByLo(myThid), myByHi(myThid)
200 DO bi = myBxLo(myThid), myBxHi(myThid)
201 CALL TIMEAVE_NORMALIZ(FUtave ,SEAICE_timeave, 1,
202 & bi,bj,myThid)
203 CALL TIMEAVE_NORMALIZ(FVtave ,SEAICE_timeave, 1,
204 & bi,bj,myThid)
205 CALL TIMEAVE_NORMALIZ(EmPmRtave,SEAICE_timeave, 1,
206 & bi,bj,myThid)
207 CALL TIMEAVE_NORMALIZ(QNETtave ,SEAICE_timeave, 1,
208 & bi,bj,myThid)
209 CALL TIMEAVE_NORMALIZ(QSWtave ,SEAICE_timeave, 1,
210 & bi,bj,myThid)
211 CALL TIMEAVE_NORMALIZ(UICEtave ,SEAICE_timeave, 1,
212 & bi,bj,myThid)
213 CALL TIMEAVE_NORMALIZ(VICEtave ,SEAICE_timeave, 1,
214 & bi,bj,myThid)
215 CALL TIMEAVE_NORMALIZ(HEFFtave ,SEAICE_timeave, 1,
216 & bi,bj,myThid)
217 CALL TIMEAVE_NORMALIZ(AREAtave ,SEAICE_timeave, 1,
218 & bi,bj,myThid)
219 ENDDO
220 ENDDO
221
222 WRITE(suff,'(I10.10)') myIter
223 _BARRIER
224 _BEGIN_MASTER( myThid )
225 CALL WRITE_FLD_XY_RL('FUtave.' ,suff,FUtave ,myIter,myThid)
226 CALL WRITE_FLD_XY_RL('FVtave.' ,suff,FVtave ,myIter,myThid)
227 CALL WRITE_FLD_XY_RL('EmPmRtave.',suff,EmPmRtave,myIter,myThid)
228 CALL WRITE_FLD_XY_RL('QNETtave.' ,suff,QNETtave ,myIter,myThid)
229 CALL WRITE_FLD_XY_RL('QSWtave.' ,suff,QSWtave ,myIter,myThid)
230 CALL WRITE_FLD_XY_RL('UICEtave.' ,suff,UICEtave ,myIter,myThid)
231 CALL WRITE_FLD_XY_RL('VICEtave.' ,suff,VICEtave ,myIter,myThid)
232 CALL WRITE_FLD_XY_RL('HEFFtave.' ,suff,HEFFtave ,myIter,myThid)
233 CALL WRITE_FLD_XY_RL('AREAtave.' ,suff,AREAtave ,myIter,myThid)
234 _END_MASTER( myThid )
235 _BARRIER
236
237 C Reset averages to zero
238 DO bj = myByLo(myThid), myByHi(myThid)
239 DO bi = myBxLo(myThid), myBxHi(myThid)
240 CALL TIMEAVE_RESET(FUtave ,1,bi,bj,myThid)
241 CALL TIMEAVE_RESET(FVtave ,1,bi,bj,myThid)
242 CALL TIMEAVE_RESET(EmPmRtave,1,bi,bj,myThid)
243 CALL TIMEAVE_RESET(QNETtave ,1,bi,bj,myThid)
244 CALL TIMEAVE_RESET(QSWtave ,1,bi,bj,myThid)
245 CALL TIMEAVE_RESET(UICEtave ,1,bi,bj,myThid)
246 CALL TIMEAVE_RESET(VICEtave ,1,bi,bj,myThid)
247 CALL TIMEAVE_RESET(HEFFtave ,1,bi,bj,myThid)
248 CALL TIMEAVE_RESET(AREAtave ,1,bi,bj,myThid)
249 DO k=1,Nr
250 SEAICE_TimeAve(k,bi,bj)=ZERO
251 ENDDO
252 ENDDO
253 ENDDO
254
255 ENDIF
256
257 #endif /* ALLOW_TIMEAVE */
258
259 #endif /* ALLOW_SEAICE */
260
261 RETURN
262 END
  ViewVC Help
Powered by ViewVC 1.1.22