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

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

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

Revision 1.7 - (hide 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 edhill 1.7 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 heimbach 1.2
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 dimitri 1.6 CALL WRITE_FLD_XY_RS( 'UWIND.',suff,uwind,myIter,myThid)
49     CALL WRITE_FLD_XY_RS( 'VWIND.',suff,vwind,myIter,myThid)
50 heimbach 1.2 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 dimitri 1.6 #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 heimbach 1.2 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 dimitri 1.6 #ifdef SEAICE_DEBUG
92     _EXCH_XY_R4( arr, myThid )
93     CALL PLOT_FIELD_XYRS( arr , 'Current uice ', myIter, myThid )
94     #endif
95    
96 heimbach 1.2 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 dimitri 1.6 #ifdef SEAICE_DEBUG
111     _EXCH_XY_R4( arr, myThid )
112     CALL PLOT_FIELD_XYRS( arr , 'Current vice ', myIter, myThid )
113     #endif
114 heimbach 1.2
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 dimitri 1.6 #ifdef SEAICE_DEBUG
130     _EXCH_XY_R4( arr, myThid )
131     CALL PLOT_FIELD_XYRS( arr , 'Current heff ', myIter, myThid )
132     #endif
133 heimbach 1.2
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 dimitri 1.6 _EXCH_XY_R4( arr, myThid )
150     CALL PLOT_FIELD_XYRS( arr , 'Current area ', myIter, myThid )
151 dimitri 1.5 #endif
152 heimbach 1.2
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 dimitri 1.5 _BARRIER
224     _BEGIN_MASTER( myThid )
225 heimbach 1.2 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 dimitri 1.5 _END_MASTER( myThid )
235     _BARRIER
236 heimbach 1.2
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 dimitri 1.3 SEAICE_TimeAve(k,bi,bj)=ZERO
251 heimbach 1.2 ENDDO
252     ENDDO
253     ENDDO
254    
255     ENDIF
256    
257     #endif /* ALLOW_TIMEAVE */
258    
259 dimitri 1.5 #endif /* ALLOW_SEAICE */
260 heimbach 1.2
261     RETURN
262     END
  ViewVC Help
Powered by ViewVC 1.1.22