1 | MODULE limdia |
---|
2 | !!====================================================================== |
---|
3 | !! *** MODULE limdia *** |
---|
4 | !! diagnostics of ice model |
---|
5 | !!====================================================================== |
---|
6 | !! To add a new field |
---|
7 | !! 1) in lim_dia : add its definition for both hemispheres if wished |
---|
8 | !! 2) add the new titles in lim_dia_init |
---|
9 | !!---------------------------------------------------------------------- |
---|
10 | #if defined key_lim3 |
---|
11 | !!---------------------------------------------------------------------- |
---|
12 | !! 'key_lim3' : LIM3 sea-ice model |
---|
13 | !!---------------------------------------------------------------------- |
---|
14 | !! lim_dia : computation of the time evolution of keys var. |
---|
15 | !! lim_dia_init : initialization and namelist read |
---|
16 | !!---------------------------------------------------------------------- |
---|
17 | !! * Modules used |
---|
18 | USE phycst |
---|
19 | USE in_out_manager |
---|
20 | USE par_ice ! ice parameters |
---|
21 | USE daymod |
---|
22 | USE dom_ice |
---|
23 | USE ice |
---|
24 | USE iceini |
---|
25 | USE limistate |
---|
26 | USE dom_oce |
---|
27 | USE sbc_oce ! Surface boundary condition: ocean fields |
---|
28 | |
---|
29 | IMPLICIT NONE |
---|
30 | PRIVATE |
---|
31 | |
---|
32 | !! * Routine accessibility |
---|
33 | PUBLIC lim_dia ! called by ice_step |
---|
34 | |
---|
35 | !! * Shared module variables |
---|
36 | INTEGER, PUBLIC :: & |
---|
37 | ntmoy = 1 , & !: instantaneous values of ice evolution or averaging ntmoy |
---|
38 | ninfo = 1 !: frequency of ouputs on file ice_evolu in case of averaging |
---|
39 | |
---|
40 | !! * Module variables |
---|
41 | INTEGER, PARAMETER :: & ! Parameters for outputs to files "evolu" |
---|
42 | jpinfmx = 100 , & ! maximum number of key variables |
---|
43 | jpchinf = 5 , & ! ??? |
---|
44 | jpchsep = jpchinf + 2 ! ??? |
---|
45 | |
---|
46 | INTEGER :: & |
---|
47 | nfrinf = 4 , & ! number of variables written in one line |
---|
48 | nferme , & ! last time step at which the var. are written on file |
---|
49 | nvinfo , & ! number of total variables |
---|
50 | nbvt , & ! number of time variables |
---|
51 | naveg ! number of step for accumulation before averaging |
---|
52 | |
---|
53 | CHARACTER(len=8) :: & |
---|
54 | fmtinf = '1PE13.5 ' ! format of the output values |
---|
55 | |
---|
56 | CHARACTER(len=30) :: & |
---|
57 | fmtw , & ! formats |
---|
58 | fmtr , & ! ??? |
---|
59 | fmtitr ! ??? |
---|
60 | |
---|
61 | CHARACTER(len=jpchsep), DIMENSION(jpinfmx) :: & |
---|
62 | titvar ! title of key variables |
---|
63 | |
---|
64 | REAL(wp), DIMENSION(jpinfmx) :: & |
---|
65 | vinfom ! temporary working space |
---|
66 | REAL(wp) :: & |
---|
67 | epsi06 = 1.e-06 |
---|
68 | REAL(wp), DIMENSION(jpi,jpj) :: & |
---|
69 | aire ! masked grid cell area |
---|
70 | |
---|
71 | !! * Substitutions |
---|
72 | # include "vectopt_loop_substitute.h90" |
---|
73 | !!---------------------------------------------------------------------- |
---|
74 | !! LIM 2.0, UCL-LOCEAN-IPSL (2005) |
---|
75 | !! $Id$ |
---|
76 | !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt |
---|
77 | !!---------------------------------------------------------------------- |
---|
78 | |
---|
79 | CONTAINS |
---|
80 | |
---|
81 | SUBROUTINE lim_dia |
---|
82 | !!-------------------------------------------------------------------- |
---|
83 | !! *** ROUTINE lim_dia *** |
---|
84 | !! |
---|
85 | !! ** Purpose : Computation and outputs on file ice.evolu |
---|
86 | !! the temporal evolution of some key variables |
---|
87 | !! |
---|
88 | !! History : |
---|
89 | !! 8.0 ! 97-06 (Louvain-La-Neuve) Original code |
---|
90 | !! 8.5 ! 02-09 (C. Ethe , G. Madec ) F90: Free form and module |
---|
91 | !! 9.0 ! 05-06 (M. Vancoppenolle LLN) Rehabilitation and fixing a few bugs |
---|
92 | !! 01-07 (M. Vancoppenolle LLN) added new fields + fram strait |
---|
93 | !! export |
---|
94 | !!------------------------------------------------------------------- |
---|
95 | !! * Local variables |
---|
96 | INTEGER :: jv,ji,jj,jl ! dummy loop indices |
---|
97 | REAL(wp), DIMENSION(jpinfmx) :: & |
---|
98 | vinfor ! temporary working space |
---|
99 | REAL(wp) :: & |
---|
100 | zshift_date , & ! date from the minimum ice extent |
---|
101 | zday, zday_min, & ! current day, day of minimum extent |
---|
102 | zafy, zamy, & ! temporary area of fy and my ice |
---|
103 | zindb |
---|
104 | !!------------------------------------------------------------------- |
---|
105 | |
---|
106 | ! 0) date from the minimum of ice extent |
---|
107 | !--------------------------------------- |
---|
108 | zday_min = 273.0 ! zday_min = date of minimum extent, here September 30th |
---|
109 | zday = FLOAT(numit-nit000) * rdt_ice / ( 86400.0 * FLOAT(nn_fsbc) ) |
---|
110 | IF (zday.GT.zday_min) THEN |
---|
111 | zshift_date = zday - zday_min |
---|
112 | ELSE |
---|
113 | zshift_date = zday - (365.0 - zday_min) |
---|
114 | ENDIF |
---|
115 | |
---|
116 | IF( numit == nstart ) CALL lim_dia_init ! initialisation of ice_evolu file |
---|
117 | |
---|
118 | ! temporal diagnostics |
---|
119 | vinfor(1) = REAL(numit) |
---|
120 | vinfor(2) = nyear |
---|
121 | |
---|
122 | ! put everything to zero |
---|
123 | DO jv = nbvt + 1, nvinfo |
---|
124 | vinfor(jv) = 0.0 |
---|
125 | END DO |
---|
126 | |
---|
127 | !!------------------------------------------------------------------- |
---|
128 | !! 1) Northern hemisphere |
---|
129 | !!------------------------------------------------------------------- |
---|
130 | !! 1.1) Diagnostics independent on age |
---|
131 | !!------------------------------------ |
---|
132 | DO jj = njeq, jpjm1 |
---|
133 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
134 | IF( tms(ji,jj) == 1 ) THEN |
---|
135 | vinfor(3) = vinfor(3) + at_i(ji,jj)*aire(ji,jj) / 1.0e12 !ice area |
---|
136 | IF (at_i(ji,jj).GT.0.15) vinfor(5) = vinfor(5) + aire(ji,jj) / 1.0e12 !ice extent |
---|
137 | vinfor(7) = vinfor(7) + vt_i(ji,jj)*aire(ji,jj) / 1.0e12 !ice volume |
---|
138 | vinfor(9) = vinfor(9) + vt_s(ji,jj)*aire(ji,jj) / 1.0e12 !snow volume |
---|
139 | vinfor(15) = vinfor(15) + ot_i(ji,jj) *vt_i(ji,jj)*aire(ji,jj) / 1.0e12 !mean age |
---|
140 | vinfor(29) = vinfor(29) + smt_i(ji,jj)*vt_i(ji,jj)*aire(ji,jj) / 1.0e12 !mean salinity |
---|
141 | ! the computation of this diagnostic is not reliable |
---|
142 | vinfor(31) = vinfor(31) + vt_i(ji,jj)*( u_ice(ji,jj)*u_ice(ji,jj) + & |
---|
143 | v_ice(ji,jj)*v_ice(ji,jj) )*aire(ji,jj)/1.0e12 |
---|
144 | vinfor(53) = vinfor(53) + emps(ji,jj)*aire(ji,jj) / 1.0e12 !salt flux |
---|
145 | vinfor(55) = vinfor(55) + fsbri(ji,jj)*aire(ji,jj) / 1.0e12 !brine drainage flux |
---|
146 | vinfor(57) = vinfor(57) + fseqv(ji,jj)*aire(ji,jj) / 1.0e12 !equivalent salt flux |
---|
147 | vinfor(59) = vinfor(59) +(sst_m(ji,jj)+rt0)*at_i(ji,jj)*aire(ji,jj) / 1.0e12 !SST |
---|
148 | vinfor(61) = vinfor(61) + sss_m(ji,jj)*at_i(ji,jj)*aire(ji,jj) / 1.0e12 !SSS |
---|
149 | vinfor(65) = vinfor(65) + et_s(ji,jj)/1.0e9*aire(ji,jj) / 1.0e12 ! snow temperature |
---|
150 | vinfor(67) = vinfor(67) + et_i(ji,jj)/1.0e9*aire(ji,jj) / 1.0e12 ! ice heat content |
---|
151 | vinfor(69) = vinfor(69) + v_i(ji,jj,1)*aire(ji,jj) / 1.0e12 !ice volume |
---|
152 | vinfor(71) = vinfor(71) + v_i(ji,jj,2)*aire(ji,jj) / 1.0e12 !ice volume |
---|
153 | vinfor(73) = vinfor(73) + v_i(ji,jj,3)*aire(ji,jj) / 1.0e12 !ice volume |
---|
154 | vinfor(75) = vinfor(75) + v_i(ji,jj,4)*aire(ji,jj) / 1.0e12 !ice volume |
---|
155 | vinfor(77) = vinfor(77) + v_i(ji,jj,5)*aire(ji,jj) / 1.0e12 !ice volume |
---|
156 | vinfor(79) = 0.0 |
---|
157 | vinfor(81) = vinfor(81) + emp(ji,jj)*aire(ji,jj) / 1.0e12 ! mass flux |
---|
158 | ENDIF |
---|
159 | END DO |
---|
160 | END DO |
---|
161 | |
---|
162 | DO jl = ice_cat_bounds(1,1), ice_cat_bounds(1,2) |
---|
163 | DO jj = njeq, jpjm1 |
---|
164 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
165 | IF( tms(ji,jj) == 1 ) THEN |
---|
166 | vinfor(11) = vinfor(11) + v_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 !undef def ice volume |
---|
167 | ENDIF |
---|
168 | END DO |
---|
169 | END DO |
---|
170 | END DO |
---|
171 | |
---|
172 | vinfor(13) = 0.0 |
---|
173 | |
---|
174 | vinfor(15) = vinfor(15) / MAX(vinfor(7),epsi06) ! these have to be divided by total ice volume to have the |
---|
175 | vinfor(29) = vinfor(29) / MAX(vinfor(7),epsi06) ! right value |
---|
176 | vinfor(31) = SQRT( vinfor(31) / MAX( vinfor(7) , epsi06 ) ) |
---|
177 | vinfor(67) = vinfor(67) / MAX(vinfor(7),epsi06) |
---|
178 | |
---|
179 | vinfor(53) = vinfor(53) / MAX(vinfor(5),epsi06) ! these have to be divided by total ice extent to have the |
---|
180 | vinfor(55) = vinfor(55) / MAX(vinfor(5),epsi06) ! right value |
---|
181 | vinfor(57) = vinfor(57) / MAX(vinfor(5),epsi06) ! |
---|
182 | vinfor(79) = vinfor(79) / MAX(vinfor(5),epsi06) ! |
---|
183 | |
---|
184 | zindb = 1.0 - MAX(0.0,SIGN(1.0,-vinfor(3))) ! |
---|
185 | vinfor(59) = zindb*vinfor(59) / MAX(vinfor(3),epsi06) ! divide by ice area |
---|
186 | vinfor(61) = zindb*vinfor(61) / MAX(vinfor(3),epsi06) ! |
---|
187 | |
---|
188 | zindb = 1.0 - MAX(0.0,SIGN(1.0,-vinfor(9))) ! |
---|
189 | vinfor(65) = zindb*vinfor(65) / MAX(vinfor(9),epsi06) ! divide it by snow volume |
---|
190 | |
---|
191 | |
---|
192 | DO jl = 1, jpl |
---|
193 | DO jj = njeq, jpjm1 |
---|
194 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
195 | IF( tms(ji,jj) == 1 ) THEN |
---|
196 | vinfor(33) = vinfor(33) + d_v_i_trp(ji,jj,jl)*aire(ji,jj) / 1.0e12 !ice volume |
---|
197 | vinfor(35) = vinfor(35) + d_v_i_thd(ji,jj,jl)*aire(ji,jj) / 1.0e12 !ice volume |
---|
198 | ENDIF |
---|
199 | END DO |
---|
200 | END DO |
---|
201 | END DO |
---|
202 | |
---|
203 | DO jj = njeq, jpjm1 |
---|
204 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
205 | IF( tms(ji,jj) == 1 ) THEN |
---|
206 | vinfor(37) = vinfor(37) + diag_sni_gr(ji,jj)*aire(ji,jj) / 1.0e12 !th growth rates |
---|
207 | vinfor(39) = vinfor(39) + diag_lat_gr(ji,jj)*aire(ji,jj) / 1.0e12 |
---|
208 | vinfor(41) = vinfor(41) + diag_bot_gr(ji,jj)*aire(ji,jj) / 1.0e12 |
---|
209 | vinfor(43) = vinfor(43) + diag_dyn_gr(ji,jj)*aire(ji,jj) / 1.0e12 |
---|
210 | vinfor(45) = vinfor(45) + dv_dt_thd(ji,jj,5)*aire(ji,jj) / 1.0e12 |
---|
211 | vinfor(47) = vinfor(47) + v_newice(ji,jj) *aire(ji,jj) / 1.0e12 / rdt_ice ! volume acc in OW |
---|
212 | ENDIF |
---|
213 | END DO |
---|
214 | END DO |
---|
215 | |
---|
216 | DO jl = 1, jpl |
---|
217 | DO jj = njeq, jpjm1 |
---|
218 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
219 | IF( tms(ji,jj) == 1 ) THEN |
---|
220 | vinfor(63) = vinfor(63) + t_su(ji,jj,jl)*a_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 |
---|
221 | ENDIF |
---|
222 | END DO |
---|
223 | END DO |
---|
224 | END DO |
---|
225 | vinfor(63) = vinfor(63) / MAX(vinfor(3),epsi06) ! these have to be divided by total ice area |
---|
226 | |
---|
227 | !! 1.2) Diagnostics dependent on age |
---|
228 | !!------------------------------------ |
---|
229 | DO jj = njeq, jpjm1 |
---|
230 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
231 | IF( tms(ji,jj) == 1 ) THEN |
---|
232 | zafy = 0.0 |
---|
233 | zamy = 0.0 |
---|
234 | DO jl = 1, jpl |
---|
235 | IF ((o_i(ji,jj,jl) - zshift_date).LT.0.0) THEN |
---|
236 | vinfor(17) = vinfor(17) + a_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 ! FY ice area |
---|
237 | vinfor(25) = vinfor(25) + v_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 ! FY ice volume |
---|
238 | vinfor(49) = vinfor(49) + sm_i(ji,jj,jl)*v_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 !FY ice salinity |
---|
239 | zafy = zafy + a_i(ji,jj,jl) |
---|
240 | ENDIF |
---|
241 | IF ((o_i(ji,jj,jl) - zshift_date).GT.0.0) THEN |
---|
242 | vinfor(19) = vinfor(19) + a_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 ! MY ice area |
---|
243 | vinfor(27) = vinfor(27) + v_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 ! MY ice volume |
---|
244 | vinfor(51) = vinfor(51) + sm_i(ji,jj,jl)*v_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 !MY ice salinity |
---|
245 | zamy = zamy + a_i(ji,jj,jl) |
---|
246 | ENDIF |
---|
247 | END DO |
---|
248 | IF ((at_i(ji,jj).GT.0.15).AND.(zafy.GT.zamy)) THEN |
---|
249 | vinfor(21) = vinfor(21) + aire(ji,jj) / 1.0e12 ! Seasonal ice extent |
---|
250 | ENDIF |
---|
251 | IF ((at_i(ji,jj).GT.0.15).AND.(zafy.LE.zamy)) THEN |
---|
252 | vinfor(23) = vinfor(23) + aire(ji,jj) / 1.0e12 ! Perennial ice extent |
---|
253 | ENDIF |
---|
254 | ENDIF |
---|
255 | END DO |
---|
256 | END DO |
---|
257 | zindb = 1.0 - MAX(0.0,SIGN(1.0,-vinfor(25))) !=0 if no multiyear ice 1 if yes |
---|
258 | vinfor(49) = zindb*vinfor(49) / MAX(vinfor(25),epsi06) |
---|
259 | zindb = 1.0 - MAX(0.0,SIGN(1.0,-vinfor(27))) !=0 if no multiyear ice 1 if yes |
---|
260 | vinfor(51) = zindb*vinfor(51) / MAX(vinfor(27),epsi06) |
---|
261 | |
---|
262 | !! Fram Strait Export |
---|
263 | !! 83 = area export |
---|
264 | !! 84 = volume export |
---|
265 | !! Fram strait in ORCA2 = 5 points |
---|
266 | !! export = -v_ice*e1t*ddtb*at_i or -v_ice*e1t*ddtb*at_i*h_i |
---|
267 | jj = 136 ! C grid |
---|
268 | vinfor(83) = 0.0 |
---|
269 | vinfor(84) = 0.0 |
---|
270 | DO ji = 134, 138 |
---|
271 | vinfor(83) = vinfor(83) - v_ice(ji,jj) * & |
---|
272 | e1t(ji,jj)*at_i(ji,jj)*rdt_ice / 1.0e12 |
---|
273 | vinfor(84) = vinfor(84) - v_ice(ji,jj) * & |
---|
274 | e1t(ji,jj)*vt_i(ji,jj)*rdt_ice / 1.0e12 |
---|
275 | END DO |
---|
276 | |
---|
277 | !!------------------------------------------------------------------- |
---|
278 | !! 2) Southern hemisphere |
---|
279 | !!------------------------------------------------------------------- |
---|
280 | !! 2.1) Diagnostics independent on age |
---|
281 | !!------------------------------------ |
---|
282 | DO jj = 2, njeqm1 |
---|
283 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
284 | IF( tms(ji,jj) == 1 ) THEN |
---|
285 | vinfor(4) = vinfor(4) + at_i(ji,jj)*aire(ji,jj) / 1.0e12 !ice area |
---|
286 | IF (at_i(ji,jj).GT.0.15) vinfor(6) = vinfor(6) + aire(ji,jj) / 1.0e12 !ice extent |
---|
287 | vinfor(8) = vinfor(8) + vt_i(ji,jj)*aire(ji,jj) / 1.0e12 !ice volume |
---|
288 | vinfor(10) = vinfor(10) + vt_s(ji,jj)*aire(ji,jj) / 1.0e12 !snow volume |
---|
289 | vinfor(16) = vinfor(16) + ot_i(ji,jj)*vt_i(ji,jj)*aire(ji,jj) / 1.0e12 !mean age |
---|
290 | vinfor(30) = vinfor(30) + smt_i(ji,jj)*vt_i(ji,jj)*aire(ji,jj) / 1.0e12 !mean salinity |
---|
291 | ! this diagnostic is not well computed (weighted by vol instead |
---|
292 | ! of area) |
---|
293 | vinfor(32) = vinfor(32) + vt_i(ji,jj)*( u_ice(ji,jj)*u_ice(ji,jj) + & |
---|
294 | v_ice(ji,jj)*v_ice(ji,jj) )*aire(ji,jj)/1.0e12 !ice vel |
---|
295 | vinfor(54) = vinfor(54) + at_i(ji,jj)*emps(ji,jj)*aire(ji,jj) / 1.0e12 ! Total salt flux |
---|
296 | vinfor(56) = vinfor(56) + at_i(ji,jj)*fsbri(ji,jj)*aire(ji,jj) / 1.0e12 ! Brine drainage salt flux |
---|
297 | vinfor(58) = vinfor(58) + at_i(ji,jj)*fseqv(ji,jj)*aire(ji,jj) / 1.0e12 ! Equivalent salt flux |
---|
298 | vinfor(60) = vinfor(60) +(sst_m(ji,jj)+rt0)*at_i(ji,jj)*aire(ji,jj) / 1.0e12 !SST |
---|
299 | vinfor(62) = vinfor(62) + sss_m(ji,jj)*at_i(ji,jj)*aire(ji,jj) / 1.0e12 !SSS |
---|
300 | vinfor(66) = vinfor(66) + et_s(ji,jj)/1.0e9*aire(ji,jj) / 1.0e12 ! snow temperature |
---|
301 | vinfor(68) = vinfor(68) + et_i(ji,jj)/1.0e9*aire(ji,jj) / 1.0e12 ! ice enthalpy |
---|
302 | vinfor(70) = vinfor(70) + v_i(ji,jj,1)*aire(ji,jj) / 1.0e12 !ice volume |
---|
303 | vinfor(72) = vinfor(72) + v_i(ji,jj,2)*aire(ji,jj) / 1.0e12 !ice volume |
---|
304 | vinfor(74) = vinfor(74) + v_i(ji,jj,3)*aire(ji,jj) / 1.0e12 !ice volume |
---|
305 | vinfor(76) = vinfor(76) + v_i(ji,jj,4)*aire(ji,jj) / 1.0e12 !ice volume |
---|
306 | vinfor(78) = vinfor(78) + v_i(ji,jj,5)*aire(ji,jj) / 1.0e12 !ice volume |
---|
307 | vinfor(80) = 0.0 |
---|
308 | vinfor(82) = vinfor(82) + emp(ji,jj)*aire(ji,jj) / 1.0e12 ! mass flux |
---|
309 | ENDIF |
---|
310 | END DO |
---|
311 | END DO |
---|
312 | |
---|
313 | DO jl = ice_cat_bounds(1,1), ice_cat_bounds(1,2) |
---|
314 | DO jj = 2, njeqm1 |
---|
315 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
316 | vinfor(12) = vinfor(12) + v_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 !undef def ice volume |
---|
317 | END DO |
---|
318 | END DO |
---|
319 | END DO |
---|
320 | |
---|
321 | vinfor(14) = 0.0 |
---|
322 | |
---|
323 | zindb = 1.0 - MAX(0.0,SIGN(1.0,-vinfor(8))) |
---|
324 | vinfor(16) = zindb * vinfor(16) / MAX(vinfor(8),epsi06) ! these have to be divided by ice vol |
---|
325 | vinfor(30) = zindb * vinfor(30) / MAX(vinfor(8),epsi06) ! |
---|
326 | vinfor(32) = zindb * SQRT( vinfor(32) / MAX( vinfor(8) , epsi06 ) ) |
---|
327 | vinfor(68) = zindb * vinfor(68) / MAX(vinfor(8),epsi06) ! |
---|
328 | |
---|
329 | zindb = 1.0 - MAX(0.0,SIGN(1.0,-vinfor(6))) |
---|
330 | vinfor(54) = zindb * vinfor(54) / MAX(vinfor(6),epsi06) ! these have to be divided by ice extt |
---|
331 | vinfor(56) = zindb * vinfor(56) / MAX(vinfor(6),epsi06) ! |
---|
332 | vinfor(58) = zindb * vinfor(58) / MAX(vinfor(6),epsi06) ! |
---|
333 | vinfor(80) = zindb * vinfor(80) / MAX(vinfor(6),epsi06) ! |
---|
334 | ! vinfor(84) = vinfor(84) / vinfor(6) ! |
---|
335 | |
---|
336 | zindb = 1.0 - MAX(0.0,SIGN(1.0,-vinfor(4))) ! |
---|
337 | vinfor(60) = zindb*vinfor(60) / ( MAX(vinfor(4), epsi06) ) ! divide by ice area |
---|
338 | vinfor(62) = zindb*vinfor(62) / ( MAX(vinfor(4), epsi06) ) ! |
---|
339 | |
---|
340 | zindb = 1.0 - MAX(0.0,SIGN(1.0,-vinfor(10))) ! |
---|
341 | vinfor(66) = zindb*vinfor(66) / MAX(vinfor(10),epsi06) ! divide it by snow volume |
---|
342 | |
---|
343 | DO jl = 1, jpl |
---|
344 | DO jj = 2, njeqm1 |
---|
345 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
346 | IF( tms(ji,jj) == 1 ) THEN |
---|
347 | vinfor(34) = vinfor(34) + d_v_i_trp(ji,jj,jl)*aire(ji,jj) / 1.0e12 !ice volume |
---|
348 | vinfor(36) = vinfor(36) + d_v_i_thd(ji,jj,jl)*aire(ji,jj) / 1.0e12 !ice volume |
---|
349 | ENDIF |
---|
350 | END DO |
---|
351 | END DO |
---|
352 | END DO |
---|
353 | |
---|
354 | DO jj = 2, njeqm1 |
---|
355 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
356 | IF( tms(ji,jj) == 1 ) THEN |
---|
357 | vinfor(38) = vinfor(38) + diag_sni_gr(ji,jj)*aire(ji,jj) / 1.0e12 !th growth rates |
---|
358 | vinfor(40) = vinfor(40) + diag_lat_gr(ji,jj)*aire(ji,jj) / 1.0e12 |
---|
359 | vinfor(42) = vinfor(42) + diag_bot_gr(ji,jj)*aire(ji,jj) / 1.0e12 |
---|
360 | vinfor(44) = vinfor(44) + diag_dyn_gr(ji,jj)*aire(ji,jj) / 1.0e12 |
---|
361 | vinfor(46) = vinfor(46) + dv_dt_thd(ji,jj,5)*aire(ji,jj) / 1.0e12 |
---|
362 | vinfor(48) = vinfor(48) + v_newice(ji,jj) *aire(ji,jj) / 1.0e12 / rdt_ice ! volume acc in OW |
---|
363 | ENDIF |
---|
364 | END DO |
---|
365 | END DO |
---|
366 | |
---|
367 | |
---|
368 | DO jl = 1, jpl |
---|
369 | DO jj = 2, njeqm1 |
---|
370 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
371 | IF( tms(ji,jj) == 1 ) THEN |
---|
372 | vinfor(64) = vinfor(64) + t_su(ji,jj,jl)*a_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 |
---|
373 | ENDIF |
---|
374 | END DO |
---|
375 | END DO |
---|
376 | END DO |
---|
377 | zindb = 1.0 - MAX(0.0,SIGN(1.0,-vinfor(4))) ! |
---|
378 | vinfor(64) = zindb * vinfor(64) / MAX(vinfor(4),epsi06) ! divide by ice extt |
---|
379 | !! 2.2) Diagnostics dependent on age |
---|
380 | !!------------------------------------ |
---|
381 | DO jj = 2, njeqm1 |
---|
382 | DO ji = fs_2, fs_jpim1 ! vector opt. |
---|
383 | IF( tms(ji,jj) == 1 ) THEN |
---|
384 | zafy = 0.0 |
---|
385 | zamy = 0.0 |
---|
386 | DO jl = 1, jpl |
---|
387 | IF ((o_i(ji,jj,jl) - zshift_date).LT.0.0) THEN |
---|
388 | vinfor(18) = vinfor(18) + a_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 ! FY ice area |
---|
389 | vinfor(26) = vinfor(26) + v_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 ! FY ice volume |
---|
390 | zafy = zafy + a_i(ji,jj,jl) |
---|
391 | vinfor(50) = vinfor(50) + sm_i(ji,jj,jl)*v_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 !FY ice salinity |
---|
392 | ENDIF |
---|
393 | IF ((o_i(ji,jj,jl) - zshift_date).GT.0.0) THEN |
---|
394 | vinfor(20) = vinfor(20) + a_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 ! MY ice area |
---|
395 | vinfor(28) = vinfor(28) + v_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 |
---|
396 | vinfor(52) = vinfor(52) + sm_i(ji,jj,jl)*v_i(ji,jj,jl)*aire(ji,jj) / 1.0e12 !FY ice salinity |
---|
397 | zamy = zamy + a_i(ji,jj,jl) |
---|
398 | ENDIF |
---|
399 | END DO ! jl |
---|
400 | IF ((at_i(ji,jj).GT.0.15).AND.(zafy.GT.zamy)) THEN |
---|
401 | vinfor(22) = vinfor(22) + aire(ji,jj) / 1.0e12 ! Seasonal ice extent |
---|
402 | ENDIF |
---|
403 | IF ((at_i(ji,jj).GT.0.15).AND.(zafy.LE.zamy)) THEN |
---|
404 | vinfor(24) = vinfor(24) + aire(ji,jj) / 1.0e12 ! Perennial ice extent |
---|
405 | ENDIF |
---|
406 | ENDIF ! tms |
---|
407 | END DO ! jj |
---|
408 | END DO ! ji |
---|
409 | zindb = 1.0 - MAX(0.0,SIGN(1.0,-vinfor(26))) !=0 if no multiyear ice 1 if yes |
---|
410 | vinfor(50) = zindb*vinfor(50) / MAX(vinfor(26),epsi06) |
---|
411 | zindb = 1.0 - MAX(0.0,SIGN(1.0,-vinfor(28))) !=0 if no multiyear ice 1 if yes |
---|
412 | vinfor(52) = zindb*vinfor(52) / MAX(vinfor(28),epsi06) |
---|
413 | |
---|
414 | ! Accumulation before averaging |
---|
415 | DO jv = 1, nvinfo |
---|
416 | vinfom(jv) = vinfom(jv) + vinfor(jv) |
---|
417 | END DO |
---|
418 | naveg = naveg + 1 |
---|
419 | |
---|
420 | ! oututs on file ice_evolu |
---|
421 | !MV IF( MOD( numit , ninfo ) == 0 ) THEN |
---|
422 | WRITE(numevo_ice,fmtw) ( titvar(jv), vinfom(jv)/naveg, jv = 1, nvinfo ) |
---|
423 | naveg = 0 |
---|
424 | DO jv = 1, nvinfo |
---|
425 | vinfom(jv)=0.0 |
---|
426 | END DO |
---|
427 | !MV ENDIF |
---|
428 | |
---|
429 | END SUBROUTINE lim_dia |
---|
430 | |
---|
431 | SUBROUTINE lim_dia_init |
---|
432 | !!------------------------------------------------------------------- |
---|
433 | !! *** ROUTINE lim_dia_init *** |
---|
434 | !! |
---|
435 | !! ** Purpose : Preparation of the file ice_evolu for the output of |
---|
436 | !! the temporal evolution of key variables |
---|
437 | !! |
---|
438 | !! ** input : Namelist namicedia |
---|
439 | !! |
---|
440 | !! history : |
---|
441 | !! 8.5 ! 03-08 (C. Ethe) original code |
---|
442 | !! 9.0 ! 08-03 (M. Vancoppenolle) LIM3 |
---|
443 | !!------------------------------------------------------------------- |
---|
444 | NAMELIST/namicedia/fmtinf, nfrinf, ninfo, ntmoy |
---|
445 | |
---|
446 | INTEGER :: jv , & ! dummy loop indice |
---|
447 | & ntot , & |
---|
448 | & ndeb |
---|
449 | |
---|
450 | REAL(wp) :: zxx0, zxx1 ! temporary scalars |
---|
451 | |
---|
452 | CHARACTER(len=jpchinf) :: titinf |
---|
453 | !!------------------------------------------------------------------- |
---|
454 | |
---|
455 | |
---|
456 | ! Read Namelist namicedia |
---|
457 | REWIND ( numnam_ice ) |
---|
458 | READ ( numnam_ice , namicedia ) |
---|
459 | IF(lwp) THEN |
---|
460 | WRITE(numout,*) |
---|
461 | WRITE(numout,*) 'lim_dia_init : ice parameters for ice diagnostics ' |
---|
462 | WRITE(numout,*) '~~~~~~~~~~~~' |
---|
463 | WRITE(numout,*) ' format of the output values fmtinf = ', fmtinf |
---|
464 | WRITE(numout,*) ' number of variables written in one line nfrinf = ', nfrinf |
---|
465 | WRITE(numout,*) ' Instantaneous values of ice evolution or averaging ntmoy = ', ntmoy |
---|
466 | WRITE(numout,*) ' frequency of ouputs on file ice_evolu in case of averaging ninfo = ', ninfo |
---|
467 | ENDIF |
---|
468 | |
---|
469 | ! masked grid cell area |
---|
470 | aire(:,:) = area(:,:) * tms(:,:) |
---|
471 | |
---|
472 | ! Titles of ice key variables : |
---|
473 | titvar(1) = 'NoIt' ! iteration number |
---|
474 | titvar(2) = 'T yr' ! time step in years |
---|
475 | nbvt = 2 ! number of time variables |
---|
476 | |
---|
477 | titvar(3) = 'AI_N' ! sea ice area in the northern Hemisp.(10^12 km2) |
---|
478 | titvar(4) = 'AI_S' ! sea ice area in the southern Hemisp.(10^12 km2) |
---|
479 | titvar(5) = 'EI_N' ! sea ice extent (15%) in the northern Hemisp.(10^12 km2) |
---|
480 | titvar(6) = 'EI_S' ! sea ice extent (15%) in the southern Hemisp.(10^12 km2) |
---|
481 | titvar(7) = 'VI_N' ! sea ice volume in the northern Hemisp.(10^3 km3) |
---|
482 | titvar(8) = 'VI_S' ! sea ice volume in the southern Hemisp.(10^3 km3) |
---|
483 | titvar(9) = 'VS_N' ! snow volume over sea ice in the northern Hemisp.(10^3 km3) |
---|
484 | titvar(10)= 'VS_S' ! snow volume over sea ice in the northern Hemisp.(10^3 km3) |
---|
485 | titvar(11)= 'VuIN' ! undeformed sea ice volume in the northern Hemisp.(10^3 km3) |
---|
486 | titvar(12)= 'VuIS' ! undeformed sea ice volume in the southern Hemisp.(10^3 km3) |
---|
487 | titvar(13)= 'VdIN' ! deformed sea ice volume in the northern Hemisp.(10^3 km3) |
---|
488 | titvar(14)= 'VdIS' ! deformed sea ice volume in the southern Hemisp.(10^3 km3) |
---|
489 | titvar(15)= 'OI_N' ! sea ice mean age in the northern Hemisp.(years) |
---|
490 | titvar(16)= 'OI_S' ! sea ice mean age in the southern Hemisp.(years) |
---|
491 | titvar(17)= 'AFYN' ! total FY ice area northern Hemisp.(10^12 km2) |
---|
492 | titvar(18)= 'AFYS' ! total FY ice area southern Hemisp.(10^12 km2) |
---|
493 | titvar(19)= 'AMYN' ! total MY ice area northern Hemisp.(10^12 km2) |
---|
494 | titvar(20)= 'AMYS' ! total MY ice area southern Hemisp.(10^12 km2) |
---|
495 | titvar(21)= 'EFYN' ! total FY ice extent northern Hemisp.(10^12 km2) (with more 50% FY ice) |
---|
496 | titvar(22)= 'EFYS' ! total FY ice extent southern Hemisp.(10^12 km2) (with more 50% FY ice) |
---|
497 | titvar(23)= 'EMYN' ! total MY ice extent northern Hemisp.(10^12 km2) (with more 50% MY ice) |
---|
498 | titvar(24)= 'EMYS' ! total MY ice extent southern Hemisp.(10^12 km2) (with more 50% MY ice) |
---|
499 | titvar(25)= 'VFYN' ! total undeformed FY ice volume northern Hemisp.(10^3 km3) |
---|
500 | titvar(26)= 'VFYS' ! total undeformed FY ice volume southern Hemisp.(10^3 km3) |
---|
501 | titvar(27)= 'VMYN' ! total undeformed MY ice volume northern Hemisp.(10^3 km3) |
---|
502 | titvar(28)= 'VMYS' ! total undeformed MY ice volume southern Hemisp.(10^3 km3) |
---|
503 | titvar(29)= 'IS_N' ! sea ice mean salinity in the northern hemisphere (ppt) |
---|
504 | titvar(30)= 'IS_S' ! sea ice mean salinity in the southern hemisphere (ppt) |
---|
505 | titvar(31)= 'IVeN' ! sea ice mean velocity in the northern hemisphere (m/s) |
---|
506 | titvar(32)= 'IVeS' ! sea ice mean velocity in the southern hemisphere (m/s) |
---|
507 | titvar(33)= 'DVDN' ! variation of sea ice volume due to dynamics in the northern hemisphere |
---|
508 | titvar(34)= 'DVDS' ! variation of sea ice volume due to dynamics in the southern hemisphere |
---|
509 | titvar(35)= 'DVTN' ! variation of sea ice volume due to thermo in the northern hemisphere |
---|
510 | titvar(36)= 'DVTS' ! variation of sea ice volume due to thermo in the southern hemisphere |
---|
511 | titvar(37)= 'TG1N' ! thermodynamic vertical growth rate in the northern hemisphere, cat 1 |
---|
512 | titvar(38)= 'TG1S' ! thermodynamic vertical growth rate in the souhtern hemisphere, cat 1 |
---|
513 | titvar(39)= 'TG2N' ! thermodynamic vertical growth rate in the northern hemisphere, cat 2 |
---|
514 | titvar(40)= 'TG2S' ! thermodynamic vertical growth rate in the souhtern hemisphere, cat 2 |
---|
515 | titvar(41)= 'TG3N' ! thermodynamic vertical growth rate in the northern hemisphere, cat 3 |
---|
516 | titvar(42)= 'TG3S' ! thermodynamic vertical growth rate in the souhtern hemisphere, cat 3 |
---|
517 | titvar(43)= 'TG4N' ! thermodynamic vertical growth rate in the northern hemisphere, cat 4 |
---|
518 | titvar(44)= 'TG4S' ! thermodynamic vertical growth rate in the souhtern hemisphere, cat 4 |
---|
519 | titvar(45)= 'TG5N' ! thermodynamic vertical growth rate in the northern hemisphere, cat 5 |
---|
520 | titvar(46)= 'TG5S' ! thermodynamic vertical growth rate in the souhtern hemisphere, cat 5 |
---|
521 | titvar(47)= 'LA_N' ! lateral accretion growth rate, northern hemisphere |
---|
522 | titvar(48)= 'LA_S' ! lateral accretion growth rate, southern hemisphere |
---|
523 | titvar(49)= 'SF_N' ! Salinity FY, NH |
---|
524 | titvar(50)= 'SF_S' ! Salinity FY, SH |
---|
525 | titvar(51)= 'SF_N' ! Salinity MY, NH |
---|
526 | titvar(52)= 'SF_S' ! Salinity MY, SH |
---|
527 | titvar(53)= 'Fs_N' ! Total salt flux NH |
---|
528 | titvar(54)= 'Fs_S' ! Total salt flux SH |
---|
529 | titvar(55)= 'FsbN' ! Salt - brine drainage flux NH |
---|
530 | titvar(56)= 'FsbS' ! Salt - brine drainage flux SH |
---|
531 | titvar(57)= 'FseN' ! Salt - Equivalent salt flux NH |
---|
532 | titvar(58)= 'FseS' ! Salt - Equivalent salt flux SH |
---|
533 | titvar(59)= 'SSTN' ! SST, NH |
---|
534 | titvar(60)= 'SSTS' ! SST, SH |
---|
535 | titvar(61)= 'SSSN' ! SSS, NH |
---|
536 | titvar(62)= 'SSSS' ! SSS, SH |
---|
537 | titvar(63)= 'TsuN' ! Tsu, NH |
---|
538 | titvar(64)= 'TsuS' ! Tsu, SH |
---|
539 | titvar(65)= 'TsnN' ! Tsn, NH |
---|
540 | titvar(66)= 'TsnS' ! Tsn, SH |
---|
541 | titvar(67)= 'ei_N' ! ei, NH |
---|
542 | titvar(68)= 'ei_S' ! ei, SH |
---|
543 | titvar(69)= 'vi1N' ! vi1, NH |
---|
544 | titvar(70)= 'vi1S' ! vi1, SH |
---|
545 | titvar(71)= 'vi2N' ! vi2, NH |
---|
546 | titvar(72)= 'vi2S' ! vi2, SH |
---|
547 | titvar(73)= 'vi3N' ! vi3, NH |
---|
548 | titvar(74)= 'vi3S' ! vi3, SH |
---|
549 | titvar(75)= 'vi4N' ! vi4, NH |
---|
550 | titvar(76)= 'vi4S' ! vi4, SH |
---|
551 | titvar(77)= 'vi5N' ! vi5, NH |
---|
552 | titvar(78)= 'vi5S' ! vi5, SH |
---|
553 | titvar(79)= 'vi6N' ! vi6, NH |
---|
554 | titvar(80)= 'vi6S' ! vi6, SH |
---|
555 | titvar(81)= 'fmaN' ! mass flux in the ocean, NH |
---|
556 | titvar(82)= 'fmaS' ! mass flux in the ocean, SH |
---|
557 | titvar(83)= 'AFSE' ! Fram Strait Area export |
---|
558 | titvar(84)= 'VFSE' ! Fram Strait Volume export |
---|
559 | nvinfo = 84 |
---|
560 | |
---|
561 | ! Definition et Ecriture de l'entete : nombre d'enregistrements |
---|
562 | ndeb = ( nstart - 1 ) / ninfo |
---|
563 | IF( nstart == 1 ) ndeb = -1 |
---|
564 | |
---|
565 | nferme = ( nstart - 1 + nitrun) / ninfo |
---|
566 | ntot = nferme - ndeb |
---|
567 | ndeb = ninfo * ( 1 + ndeb ) |
---|
568 | nferme = ninfo * nferme |
---|
569 | |
---|
570 | ! definition of formats |
---|
571 | WRITE( fmtw , '(A,I3,A2,I1,A)' ) '(', nfrinf, '(A', jpchsep, ','//fmtinf//'))' |
---|
572 | WRITE( fmtr , '(A,I3,A,I1,A)' ) '(', nfrinf, '(', jpchsep, 'X,'//fmtinf//'))' |
---|
573 | WRITE( fmtitr, '(A,I3,A,I1,A)' ) '(', nvinfo, 'A', jpchinf, ')' |
---|
574 | |
---|
575 | ! opening "ice_evolu" file |
---|
576 | CALL ctl_opn( numevo_ice, 'ice.evolu', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) |
---|
577 | |
---|
578 | !- ecriture de 2 lignes d''entete : |
---|
579 | WRITE(numevo_ice,1000) fmtr, fmtw, fmtitr, nvinfo, ntot, 0, nfrinf |
---|
580 | zxx0 = 0.001 * REAL(ninfo) |
---|
581 | zxx1 = 0.001 * REAL(ndeb) |
---|
582 | WRITE(numevo_ice,1111) REAL(jpchinf), 0., zxx1, zxx0, 0., 0., 0 |
---|
583 | |
---|
584 | !- ecriture de 2 lignes de titre : |
---|
585 | WRITE(numevo_ice,'(A,I8,A,I8,A,I5)') & |
---|
586 | 'Evolution chronologique - Experience '//cexper & |
---|
587 | //' de', ndeb, ' a', nferme, ' pas', ninfo |
---|
588 | WRITE(numevo_ice,fmtitr) ( titvar(jv), jv = 1, nvinfo ) |
---|
589 | |
---|
590 | |
---|
591 | !--preparation de "titvar" pour l''ecriture parmi les valeurs numeriques : |
---|
592 | DO jv = 2 , nvinfo |
---|
593 | titinf = titvar(jv)(:jpchinf) |
---|
594 | titvar(jv) = ' '//titinf |
---|
595 | END DO |
---|
596 | |
---|
597 | !--Initialisation of the arrays for the accumulation |
---|
598 | DO jv = 1, nvinfo |
---|
599 | vinfom(jv) = 0. |
---|
600 | END DO |
---|
601 | naveg = 0 |
---|
602 | |
---|
603 | 1000 FORMAT( 3(A20),4(1x,I6) ) |
---|
604 | 1111 FORMAT( 3(F7.1,1X,F7.3,1X),I3,A ) |
---|
605 | |
---|
606 | END SUBROUTINE lim_dia_init |
---|
607 | |
---|
608 | #else |
---|
609 | !!---------------------------------------------------------------------- |
---|
610 | !! Default option : NO LIM sea-ice model |
---|
611 | !!---------------------------------------------------------------------- |
---|
612 | CONTAINS |
---|
613 | SUBROUTINE lim_dia ! Empty routine |
---|
614 | END SUBROUTINE lim_dia |
---|
615 | #endif |
---|
616 | |
---|
617 | !!====================================================================== |
---|
618 | END MODULE limdia |
---|