/[MITgcm]/MITgcm/pkg/thsice/thsice_monitor.F
ViewVC logotype

Contents of /MITgcm/pkg/thsice/thsice_monitor.F

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

Revision 1.3 - (show annotations) (download)
Wed Apr 7 23:40:34 2004 UTC (20 years, 2 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint52n_post, checkpoint53d_post, checkpoint54a_pre, checkpoint54e_post, checkpoint54a_post, checkpoint53c_post, checkpoint53b_post, checkpoint54b_post, checkpoint53b_pre, checkpoint54d_post, checkpoint52m_post, checkpoint53a_post, checkpoint54, checkpoint53, checkpoint53g_post, checkpoint53f_post, checkpoint53d_pre, checkpoint54c_post
Changes since 1.2: +13 -5 lines 
major changes in pkg/thsice: allows atmospheric model (AIM) to use thsice.
- split thsice_therm.F in 2 S/R: thsice_solve4temp.F & thsice_calc_thickn.F
- move most of the ocean & bulk_force interface in thsice_main.F
- add a "slab ocean" component to be used with atmospheric model
1 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_monitor.F,v 1.2 2004/03/03 14:16:41 jmc Exp $
2 C $Name: $
3
4 #include "THSICE_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: THSICE_MONITOR
8 C !INTERFACE
9 SUBROUTINE THSICE_MONITOR( myTime, myIter, myThid )
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | S/R THSICE_MONITOR
14 C | o Do ICE global & Hemispheric diagnostic
15 C *==========================================================*
16 C *==========================================================*
17 C \ev
18
19 C !USES:
20 IMPLICIT NONE
21
22 C === Global variables ===
23 #include "SIZE.h"
24 #include "EEPARAMS.h"
25 #include "PARAMS.h"
26 #include "GRID.h"
27 #include "THSICE_PARAMS.h"
28 #include "THSICE_VARS.h"
29 #ifdef ALLOW_MONITOR
30 #include "MONITOR.h"
31 #endif
32
33 C !INPUT/OUTPUT PARAMETERS:
34 C == Routine arguments ==
35 C myTime - Current time of simulation ( s )
36 C myIter - Iteration number
37 C myThid - Number of this instance of INI_FORCING
38 _RL myTime
39 INTEGER myIter
40 INTEGER myThid
41 CEOP
42
43 #ifdef ALLOW_THSICE
44 #ifdef ALLOW_MONITOR
45
46 C == Local variables ==
47 LOGICAL DIFFERENT_MULTIPLE
48 EXTERNAL DIFFERENT_MULTIPLE
49 CHARACTER*(MAX_LEN_MBUF) msgBuf
50 CHARACTER*10 mon_var
51 CHARACTER*2 mon_sufx(0:2)
52 _RS yBand(2), locDr(1)
53 _RL theMin(2), theMax(2)
54 _RL theMean(2), theVar(2), theVol(2)
55 _RL theMeanG, theVarG, theVolG
56 _RL theMean0, theMean1, theMean2, theEnergy
57 _RL theMin0, theMax0, theSD, theDel2
58
59 DATA yBand / 0. , 0. /
60 DATA locDr / 1. /
61 DATA mon_sufx / '_G' , '_S' , '_N' /
62
63 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
64
65 IF ( DIFFERENT_MULTIPLE(thSIce_monFreq,myTime,myTime-deltaTclock)
66 & .OR. myIter.EQ.nIter0 ) THEN
67
68 _BEGIN_MASTER(myThid)
69 WRITE(msgBuf,'(A)')
70 &'// ======================================================='
71 CALL PRINT_MESSAGE( msgBuf, mon_ioUnit, SQUEEZE_RIGHT , 1)
72 WRITE(msgBuf,'(A)') '// Begin MONITOR Therm.SeaIce statistics'
73 CALL PRINT_MESSAGE( msgBuf, mon_ioUnit, SQUEEZE_RIGHT , 1)
74 WRITE(msgBuf,'(A)')
75 &'// ======================================================='
76 CALL PRINT_MESSAGE( msgBuf, mon_ioUnit, SQUEEZE_RIGHT , 1)
77 _END_MASTER(myThid)
78
79 CALL MON_SET_PREF('thSI_',myThid)
80 CALL MON_OUT_RL('time_sec', myTime,mon_string_none,myThid)
81
82 C-- Ice area and Ice thickness :
83 CALL MON_STATS_LATBND_RL(
84 I 1, 1, 1, 2, yBand,
85 I iceheight, iceMask, maskH, rA, yC, locDr,
86 O theMin, theMax, theMean, theVar, theVol,
87 I myThid )
88 theVolG= theVol(1)+theVol(2)
89 theMeanG= theMean(1)*theVol(1)+theMean(2)*theVol(2)
90 IF (theVolG.GT.0.) theMeanG = theMeanG / theVolG
91
92 mon_var='Ice_Area'
93 CALL MON_OUT_RL(mon_var, theVolG , mon_sufx(0), myThid)
94 CALL MON_OUT_RL(mon_var, theVol(1), mon_sufx(1), myThid)
95 CALL MON_OUT_RL(mon_var, theVol(2), mon_sufx(2), myThid)
96 mon_var='IceH_ave'
97 CALL MON_OUT_RL(mon_var,theMeanG , mon_sufx(0), myThid)
98 CALL MON_OUT_RL(mon_var,theMean(1), mon_sufx(1), myThid)
99 CALL MON_OUT_RL(mon_var,theMean(2), mon_sufx(2), myThid)
100 mon_var='IceH_max'
101 CALL MON_OUT_RL(mon_var, theMax(1), mon_sufx(1), myThid)
102 CALL MON_OUT_RL(mon_var, theMax(2), mon_sufx(2), myThid)
103
104 C-- Snow thickness :
105 CALL MON_STATS_LATBND_RL(
106 I 1, 1, 1, 2, yBand,
107 I snowheight, iceMask, maskH, rA, yC, locDr,
108 O theMin, theMax, theMean, theVar, theVol,
109 I myThid )
110 theVolG= theVol(1)+theVol(2)
111 theMeanG= theMean(1)*theVol(1)+theMean(2)*theVol(2)
112 theEnergy = -rhos*Lfresh*theMeanG
113 IF (theVolG.GT.0.) theMeanG = theMeanG / theVolG
114
115 mon_var='SnwH_ave'
116 CALL MON_OUT_RL(mon_var,theMeanG , mon_sufx(0), myThid)
117 CALL MON_OUT_RL(mon_var,theMean(1), mon_sufx(1), myThid)
118 CALL MON_OUT_RL(mon_var,theMean(2), mon_sufx(2), myThid)
119 mon_var='SnwH_max'
120 CALL MON_OUT_RL(mon_var, theMax(1), mon_sufx(1), myThid)
121 CALL MON_OUT_RL(mon_var, theMax(2), mon_sufx(2), myThid)
122
123 C-- Total Energy :
124 CALL MON_STATS_RL(
125 I 1, Qice1, iceMask, iceheight, rA, locDr,
126 O theMin0,theMax0,theMean1,theSD,theDel2,theVolG,
127 I myThid )
128 CALL MON_STATS_RL(
129 I 1, Qice2, iceMask, iceheight, rA, locDr,
130 O theMin0,theMax0,theMean2,theSD,theDel2,theVolG,
131 I myThid )
132 theEnergy = theEnergy -rhoi*(theMean1+theMean2)*theVolG/2
133 mon_var='TotEnerg'
134 CALL MON_OUT_RL(mon_var, theEnergy, mon_sufx(0), myThid)
135
136 C-- Surface Temp. :
137 CALL MON_STATS_LATBND_RL(
138 I 1, 1, 1, 2, yBand,
139 I Tsrf, iceMask, maskH, rA, yC, locDr,
140 O theMin, theMax, theMean, theVar, theVol,
141 I myThid )
142 theVolG= theVol(1)+theVol(2)
143 theMeanG= theMean(1)*theVol(1)+theMean(2)*theVol(2)
144 IF (theVolG.GT.0.) theMeanG = theMeanG / theVolG
145
146 mon_var='Tsrf_ave'
147 CALL MON_OUT_RL(mon_var,theMeanG , mon_sufx(0), myThid)
148 CALL MON_OUT_RL(mon_var,theMean(1), mon_sufx(1), myThid)
149 CALL MON_OUT_RL(mon_var,theMean(2), mon_sufx(2), myThid)
150 mon_var='Tsrf_min'
151 CALL MON_OUT_RL(mon_var, theMin(1), mon_sufx(1), myThid)
152 CALL MON_OUT_RL(mon_var, theMin(2), mon_sufx(2), myThid)
153 mon_var='Tsrf_max'
154 CALL MON_OUT_RL(mon_var, theMax(1), mon_sufx(1), myThid)
155 CALL MON_OUT_RL(mon_var, theMax(2), mon_sufx(2), myThid)
156
157 C-- 1rst level (volume-mean) Temp. :
158 CALL MON_STATS_LATBND_RL(
159 I 1, 1, 1, 2, yBand,
160 I Tice1, iceMask, iceheight, rA, yC, locDr,
161 O theMin, theMax, theMean, theVar, theVol,
162 I myThid )
163 theVolG = theVol(1)+theVol(2)
164 theMeanG= theMean(1)*theVol(1)+theMean(2)*theVol(2)
165 IF (theVolG.GT.0.) theMeanG = theMeanG / theVolG
166
167 c mon_var='IceVolum'
168 c CALL MON_OUT_RL(mon_var, theVolG , mon_sufx(0), myThid)
169 c CALL MON_OUT_RL(mon_var, theVol(1), mon_sufx(1), myThid)
170 c CALL MON_OUT_RL(mon_var, theVol(2), mon_sufx(2), myThid)
171 mon_var='Tic1_ave'
172 CALL MON_OUT_RL(mon_var,theMeanG , mon_sufx(0), myThid)
173 CALL MON_OUT_RL(mon_var,theMean(1), mon_sufx(1), myThid)
174 CALL MON_OUT_RL(mon_var,theMean(2), mon_sufx(2), myThid)
175 mon_var='Tic1_min'
176 CALL MON_OUT_RL(mon_var, theMin(1), mon_sufx(1), myThid)
177 CALL MON_OUT_RL(mon_var, theMin(2), mon_sufx(2), myThid)
178 mon_var='Tic1_max'
179 CALL MON_OUT_RL(mon_var, theMax(1), mon_sufx(1), myThid)
180 CALL MON_OUT_RL(mon_var, theMax(2), mon_sufx(2), myThid)
181
182 C-- 2nd level (volume-mean) Temp. :
183 CALL MON_STATS_LATBND_RL(
184 I 1, 1, 1, 2, yBand,
185 I Tice2, iceMask, iceheight, rA, yC, locDr,
186 O theMin, theMax, theMean, theVar, theVol,
187 I myThid )
188 theMeanG= theMean(1)*theVol(1)+theMean(2)*theVol(2)
189 IF (theVolG.GT.0.) theMeanG = theMeanG / theVolG
190
191 mon_var='Tic2_ave'
192 CALL MON_OUT_RL(mon_var,theMeanG , mon_sufx(0), myThid)
193 CALL MON_OUT_RL(mon_var,theMean(1), mon_sufx(1), myThid)
194 CALL MON_OUT_RL(mon_var,theMean(2), mon_sufx(2), myThid)
195 mon_var='Tic2_min'
196 CALL MON_OUT_RL(mon_var, theMin(1), mon_sufx(1), myThid)
197 CALL MON_OUT_RL(mon_var, theMin(2), mon_sufx(2), myThid)
198 mon_var='Tic2_max'
199 CALL MON_OUT_RL(mon_var, theMax(1), mon_sufx(1), myThid)
200 CALL MON_OUT_RL(mon_var, theMax(2), mon_sufx(2), myThid)
201
202 _BEGIN_MASTER(myThid)
203 WRITE(msgBuf,'(A)')
204 &'// ======================================================='
205 CALL PRINT_MESSAGE( msgBuf, mon_ioUnit, SQUEEZE_RIGHT , 1)
206 WRITE(msgBuf,'(A)') '// End MONITOR Therm.SeaIce statistics'
207 CALL PRINT_MESSAGE( msgBuf, mon_ioUnit, SQUEEZE_RIGHT , 1)
208 WRITE(msgBuf,'(A)')
209 &'// ======================================================='
210 CALL PRINT_MESSAGE( msgBuf, mon_ioUnit, SQUEEZE_RIGHT , 1)
211 _END_MASTER(myThid)
212
213 ENDIF
214
215 #endif /* ALLOW_MONITOR */
216 #endif /* ALLOW_THSICE */
217
218 RETURN
219 END
  ViewVC Help
Powered by ViewVC 1.1.22