Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

diahsb.F90 @ 6808

Last change on this file since 6808 was 6808, checked in by jamesharle, 8 years ago
merge with trunk@6232 for consistency with SSB code
Property svn:keywords set to `Id`
File size: 20.6 KB

Line
1	MODULE diahsb
2	!!======================================================================
3	!! * MODULE diahsb *
4	!! Ocean diagnostics: Heat, salt and volume budgets
5	!!======================================================================
6	!! History : 3.3 ! 2010-09 (M. Leclair) Original code
7	!! ! 2012-10 (C. Rousset) add iom_put
8	!!----------------------------------------------------------------------
9
10	!!----------------------------------------------------------------------
11	!! dia_hsb : Diagnose the conservation of ocean heat and salt contents, and volume
12	!! dia_hsb_rst : Read or write DIA file in restart file
13	!! dia_hsb_init : Initialization of the conservation diagnostic
14	!!----------------------------------------------------------------------
15	USE oce ! ocean dynamics and tracers
16	USE dom_oce ! ocean space and time domain
17	USE phycst ! physical constants
18	USE sbc_oce ! surface thermohaline fluxes
19	USE sbcrnf ! river runoff
20	USE sbcisf ! ice shelves
21	USE domvvl ! vertical scale factors
22	USE traqsr ! penetrative solar radiation
23	USE trabbc ! bottom boundary condition
24	USE trabbc ! bottom boundary condition
25	USE bdy_par ! (for lk_bdy)
26	USE restart ! ocean restart
27	!
28	USE iom ! I/O manager
29	USE in_out_manager ! I/O manager
30	USE lib_fortran ! glob_sum
31	USE lib_mpp ! distributed memory computing library
32	USE timing ! preformance summary
33	USE wrk_nemo ! work arrays
34
35	IMPLICIT NONE
36	PRIVATE
37
38	PUBLIC dia_hsb ! routine called by step.F90
39	PUBLIC dia_hsb_init ! routine called by nemogcm.F90
40	PUBLIC dia_hsb_rst ! routine called by step.F90
41
42	LOGICAL, PUBLIC :: ln_diahsb !: check the heat and salt budgets
43
44	REAL(wp) :: surf_tot ! ocean surface
45	REAL(wp) :: frc_t, frc_s, frc_v ! global forcing trends
46	REAL(wp) :: frc_wn_t, frc_wn_s ! global forcing trends
47	!
48	REAL(wp), DIMENSION(:,:) , ALLOCATABLE :: surf
49	REAL(wp), DIMENSION(:,:) , ALLOCATABLE :: surf_ini , ssh_ini !
50	REAL(wp), DIMENSION(:,:) , ALLOCATABLE :: ssh_hc_loc_ini, ssh_sc_loc_ini !
51	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: hc_loc_ini, sc_loc_ini, e3t_ini !
52
53	!! * Substitutions
54	# include "vectopt_loop_substitute.h90"
55	!!----------------------------------------------------------------------
56	!! NEMO/OPA 3.3 , NEMO Consortium (2010)
57	!! $Id$
58	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
59	!!----------------------------------------------------------------------
60	CONTAINS
61
62	SUBROUTINE dia_hsb( kt )
63	!!---------------------------------------------------------------------------
64	!! * ROUTINE dia_hsb *
65	!!
66	!! ** Purpose: Compute the ocean global heat content, salt content and volume conservation
67	!!
68	!! ** Method : - Compute the deviation of heat content, salt content and volume
69	!! at the current time step from their values at nit000
70	!! - Compute the contribution of forcing and remove it from these deviations
71	!!
72	!!---------------------------------------------------------------------------
73	INTEGER, INTENT(in) :: kt ! ocean time-step index
74	!
75	INTEGER :: ji, jj, jk ! dummy loop indice
76	REAL(wp) :: zdiff_hc , zdiff_sc ! heat and salt content variations
77	REAL(wp) :: zdiff_hc1 , zdiff_sc1 ! - - - -
78	REAL(wp) :: zdiff_v1 , zdiff_v2 ! volume variation
79	REAL(wp) :: zerr_hc1 , zerr_sc1 ! heat and salt content misfit
80	REAL(wp) :: zvol_tot ! volume
81	REAL(wp) :: z_frc_trd_t , z_frc_trd_s ! - -
82	REAL(wp) :: z_frc_trd_v ! - -
83	REAL(wp) :: z_wn_trd_t , z_wn_trd_s ! - -
84	REAL(wp) :: z_ssh_hc , z_ssh_sc ! - -
85	REAL(wp), DIMENSION(:,:), POINTER :: z2d0, z2d1
86	!!---------------------------------------------------------------------------
87	IF( nn_timing == 1 ) CALL timing_start('dia_hsb')
88	CALL wrk_alloc( jpi,jpj, z2d0, z2d1 )
89	!
90	tsn(:,:,:,1) = tsn(:,:,:,1) * tmask(:,:,:) ; tsb(:,:,:,1) = tsb(:,:,:,1) * tmask(:,:,:) ;
91	tsn(:,:,:,2) = tsn(:,:,:,2) * tmask(:,:,:) ; tsb(:,:,:,2) = tsb(:,:,:,2) * tmask(:,:,:) ;
92	! ------------------------- !
93	! 1 - Trends due to forcing !
94	! ------------------------- !
95	z_frc_trd_v = r1_rau0 * glob_sum( - ( emp(:,:) - rnf(:,:) + fwfisf(:,:) ) * surf(:,:) ) ! volume fluxes
96	z_frc_trd_t = glob_sum( sbc_tsc(:,:,jp_tem) * surf(:,:) ) ! heat fluxes
97	z_frc_trd_s = glob_sum( sbc_tsc(:,:,jp_sal) * surf(:,:) ) ! salt fluxes
98	! Add runoff heat & salt input
99	IF( ln_rnf ) z_frc_trd_t = z_frc_trd_t + glob_sum( rnf_tsc(:,:,jp_tem) * surf(:,:) )
100	IF( ln_rnf_sal) z_frc_trd_s = z_frc_trd_s + glob_sum( rnf_tsc(:,:,jp_sal) * surf(:,:) )
101	! Add ice shelf heat & salt input
102	IF( ln_isf ) z_frc_trd_t = z_frc_trd_t + glob_sum( risf_tsc(:,:,jp_tem) * surf(:,:) )
103	! Add penetrative solar radiation
104	IF( ln_traqsr ) z_frc_trd_t = z_frc_trd_t + r1_rau0_rcp * glob_sum( qsr (:,:) * surf(:,:) )
105	! Add geothermal heat flux
106	IF( ln_trabbc ) z_frc_trd_t = z_frc_trd_t + glob_sum( qgh_trd0(:,:) * surf(:,:) )
107	!
108	IF( ln_linssh ) THEN
109	IF( ln_isfcav ) THEN
110	DO ji=1,jpi
111	DO jj=1,jpj
112	z2d0(ji,jj) = surf(ji,jj) * wn(ji,jj,mikt(ji,jj)) * tsb(ji,jj,mikt(ji,jj),jp_tem)
113	z2d1(ji,jj) = surf(ji,jj) * wn(ji,jj,mikt(ji,jj)) * tsb(ji,jj,mikt(ji,jj),jp_sal)
114	END DO
115	END DO
116	ELSE
117	z2d0(:,:) = surf(:,:) * wn(:,:,1) * tsb(:,:,1,jp_tem)
118	z2d1(:,:) = surf(:,:) * wn(:,:,1) * tsb(:,:,1,jp_sal)
119	END IF
120	z_wn_trd_t = - glob_sum( z2d0 )
121	z_wn_trd_s = - glob_sum( z2d1 )
122	ENDIF
123
124	frc_v = frc_v + z_frc_trd_v * rdt
125	frc_t = frc_t + z_frc_trd_t * rdt
126	frc_s = frc_s + z_frc_trd_s * rdt
127	! ! Advection flux through fixed surface (z=0)
128	IF( ln_linssh ) THEN
129	frc_wn_t = frc_wn_t + z_wn_trd_t * rdt
130	frc_wn_s = frc_wn_s + z_wn_trd_s * rdt
131	ENDIF
132
133	! ------------------------ !
134	! 2 - Content variations !
135	! ------------------------ !
136	! glob_sum_full is needed because you keep the full interior domain to compute the sum (iscpl)
137	zdiff_v2 = 0._wp
138	zdiff_hc = 0._wp
139	zdiff_sc = 0._wp
140
141	! volume variation (calculated with ssh)
142	zdiff_v1 = glob_sum_full( surf(:,:) * sshn(:,:) - surf_ini(:,:) * ssh_ini(:,:) )
143
144	! heat & salt content variation (associated with ssh)
145	IF( ln_linssh ) THEN
146	IF( ln_isfcav ) THEN
147	DO ji = 1, jpi
148	DO jj = 1, jpj
149	z2d0(ji,jj) = surf(ji,jj) * ( tsn(ji,jj,mikt(ji,jj),jp_tem) * sshn(ji,jj) - ssh_hc_loc_ini(ji,jj) )
150	z2d1(ji,jj) = surf(ji,jj) * ( tsn(ji,jj,mikt(ji,jj),jp_sal) * sshn(ji,jj) - ssh_sc_loc_ini(ji,jj) )
151	END DO
152	END DO
153	ELSE
154	z2d0(:,:) = surf(:,:) * ( tsn(:,:,1,jp_tem) * sshn(:,:) - ssh_hc_loc_ini(:,:) )
155	z2d1(:,:) = surf(:,:) * ( tsn(:,:,1,jp_sal) * sshn(:,:) - ssh_sc_loc_ini(:,:) )
156	END IF
157	z_ssh_hc = glob_sum_full( z2d0 )
158	z_ssh_sc = glob_sum_full( z2d1 )
159	ENDIF
160
161	DO jk = 1, jpkm1
162	! volume variation (calculated with scale factors)
163	zdiff_v2 = zdiff_v2 + glob_sum_full( surf(:,:) * tmask(:,:,jk) &
164	& * e3t_n(:,:,jk) - surf_ini(:,:) * e3t_ini(:,:,jk) )
165	! heat content variation
166	zdiff_hc = zdiff_hc + glob_sum_full( surf(:,:) * tmask(:,:,jk) &
167	& * e3t_n(:,:,jk) * tsn(:,:,jk,jp_tem) - surf_ini(:,:) * hc_loc_ini(:,:,jk) )
168	! salt content variation
169	zdiff_sc = zdiff_sc + glob_sum_full( surf (:,:) * tmask(:,:,jk) &
170	* e3t_n(:,:,jk) * tsn(:,:,jk,jp_sal) - surf_ini(:,:) * sc_loc_ini(:,:,jk) )
171	ENDDO
172
173	! Substract forcing from heat content, salt content and volume variations
174	zdiff_v1 = zdiff_v1 - frc_v
175	IF( .NOT.ln_linssh ) zdiff_v2 = zdiff_v2 - frc_v
176	zdiff_hc = zdiff_hc - frc_t
177	zdiff_sc = zdiff_sc - frc_s
178	IF( ln_linssh ) THEN
179	zdiff_hc1 = zdiff_hc + z_ssh_hc
180	zdiff_sc1 = zdiff_sc + z_ssh_sc
181	zerr_hc1 = z_ssh_hc - frc_wn_t
182	zerr_sc1 = z_ssh_sc - frc_wn_s
183	ENDIF
184
185	! ----------------------- !
186	! 3 - Diagnostics writing !
187	! ----------------------- !
188	zvol_tot = 0._wp ! total ocean volume (calculated with scale factors)
189	DO jk = 1, jpkm1
190	zvol_tot = zvol_tot + glob_sum_full( surf(:,:) * tmask(:,:,jk) * e3t_n(:,:,jk) )
191	END DO
192
193	!!gm to be added ?
194	! IF( ln_linssh ) THEN ! fixed volume, add the ssh contribution
195	! zvol_tot = zvol_tot + glob_sum( surf(:,:) * sshn(:,:) )
196	! ENDIF
197	!!gm end
198
199	IF( ln_linssh ) THEN
200	CALL iom_put( 'bgtemper' , zdiff_hc1 / zvol_tot) ! Heat content variation (C)
201	CALL iom_put( 'bgsaline' , zdiff_sc1 / zvol_tot) ! Salt content variation (psu)
202	CALL iom_put( 'bgheatco' , zdiff_hc1 * 1.e-20 * rau0 * rcp ) ! Heat content variation (1.e20 J)
203	CALL iom_put( 'bgsaltco' , zdiff_sc1 * 1.e-9 ) ! Salt content variation (psu*km3)
204	CALL iom_put( 'bgvolssh' , zdiff_v1 * 1.e-9 ) ! volume ssh variation (km3)
205	CALL iom_put( 'bgfrcvol' , frc_v * 1.e-9 ) ! vol - surface forcing (km3)
206	CALL iom_put( 'bgfrctem' , frc_t / zvol_tot ) ! hc - surface forcing (C)
207	CALL iom_put( 'bgfrcsal' , frc_s / zvol_tot ) ! sc - surface forcing (psu)
208	CALL iom_put( 'bgmistem' , zerr_hc1 / zvol_tot ) ! hc - error due to free surface (C)
209	CALL iom_put( 'bgmissal' , zerr_sc1 / zvol_tot ) ! sc - error due to free surface (psu)
210	ELSE
211	CALL iom_put( 'bgtemper' , zdiff_hc / zvol_tot ) ! Temperature variation (C)
212	CALL iom_put( 'bgsaline' , zdiff_sc / zvol_tot ) ! Salinity variation (psu)
213	CALL iom_put( 'bgheatco' , zdiff_hc * 1.e-20 * rau0 * rcp ) ! Heat content variation (1.e20 J)
214	CALL iom_put( 'bgsaltco' , zdiff_sc * 1.e-9 ) ! Salt content variation (psu*km3)
215	CALL iom_put( 'bgvolssh' , zdiff_v1 * 1.e-9 ) ! volume ssh variation (km3)
216	CALL iom_put( 'bgvole3t' , zdiff_v2 * 1.e-9 ) ! volume e3t variation (km3)
217	CALL iom_put( 'bgfrcvol' , frc_v * 1.e-9 ) ! vol - surface forcing (km3)
218	CALL iom_put( 'bgfrctem' , frc_t / zvol_tot ) ! hc - surface forcing (C)
219	CALL iom_put( 'bgfrcsal' , frc_s / zvol_tot ) ! sc - surface forcing (psu)
220	ENDIF
221	!
222	IF( lrst_oce ) CALL dia_hsb_rst( kt, 'WRITE' )
223	!
224	CALL wrk_dealloc( jpi,jpj, z2d0, z2d1 )
225	!
226	IF( nn_timing == 1 ) CALL timing_stop('dia_hsb')
227	!
228	END SUBROUTINE dia_hsb
229
230
231	SUBROUTINE dia_hsb_rst( kt, cdrw )
232	!!---------------------------------------------------------------------
233	!! * ROUTINE limdia_rst *
234	!!
235	!! ** Purpose : Read or write DIA file in restart file
236	!!
237	!! ** Method : use of IOM library
238	!!----------------------------------------------------------------------
239	INTEGER , INTENT(in) :: kt ! ocean time-step
240	CHARACTER(len=*), INTENT(in) :: cdrw ! "READ"/"WRITE" flag
241	!
242	INTEGER :: ji, jj, jk ! dummy loop indices
243	INTEGER :: id1 ! local integers
244	!!----------------------------------------------------------------------
245	!
246	IF( TRIM(cdrw) == 'READ' ) THEN ! Read/initialise
247	IF( ln_rstart ) THEN !* Read the restart file
248	!id1 = iom_varid( numror, 'frc_vol' , ldstop = .FALSE. )
249	!
250	IF(lwp) WRITE(numout,*) '~~~~~~~'
251	IF(lwp) WRITE(numout,*) ' dia_hsb_rst at it= ', kt,' date= ', ndastp
252	IF(lwp) WRITE(numout,*) '~~~~~~~'
253	CALL iom_get( numror, 'frc_v', frc_v )
254	CALL iom_get( numror, 'frc_t', frc_t )
255	CALL iom_get( numror, 'frc_s', frc_s )
256	IF( ln_linssh ) THEN
257	CALL iom_get( numror, 'frc_wn_t', frc_wn_t )
258	CALL iom_get( numror, 'frc_wn_s', frc_wn_s )
259	ENDIF
260	CALL iom_get( numror, jpdom_autoglo, 'surf_ini', surf_ini ) ! ice sheet coupling
261	CALL iom_get( numror, jpdom_autoglo, 'ssh_ini', ssh_ini )
262	CALL iom_get( numror, jpdom_autoglo, 'e3t_ini', e3t_ini )
263	CALL iom_get( numror, jpdom_autoglo, 'hc_loc_ini', hc_loc_ini )
264	CALL iom_get( numror, jpdom_autoglo, 'sc_loc_ini', sc_loc_ini )
265	IF( ln_linssh ) THEN
266	CALL iom_get( numror, jpdom_autoglo, 'ssh_hc_loc_ini', ssh_hc_loc_ini )
267	CALL iom_get( numror, jpdom_autoglo, 'ssh_sc_loc_ini', ssh_sc_loc_ini )
268	ENDIF
269	ELSE
270	IF(lwp) WRITE(numout,*) '~~~~~~~'
271	IF(lwp) WRITE(numout,*) ' dia_hsb at initial state '
272	IF(lwp) WRITE(numout,*) '~~~~~~~'
273	surf_ini(:,:) = e1e2t(:,:) * tmask_i(:,:) ! initial ocean surface
274	ssh_ini(:,:) = sshn(:,:) ! initial ssh
275	DO jk = 1, jpk
276	! if ice sheet/oceqn coupling, need to mask ini variables here (mask could change at the next NEMO instance).
277	e3t_ini (:,:,jk) = e3t_n(:,:,jk) * tmask(:,:,jk) ! initial vertical scale factors
278	hc_loc_ini(:,:,jk) = tsn(:,:,jk,jp_tem) * e3t_n(:,:,jk) * tmask(:,:,jk) ! initial heat content
279	sc_loc_ini(:,:,jk) = tsn(:,:,jk,jp_sal) * e3t_n(:,:,jk) * tmask(:,:,jk) ! initial salt content
280	END DO
281	frc_v = 0._wp ! volume trend due to forcing
282	frc_t = 0._wp ! heat content - - - -
283	frc_s = 0._wp ! salt content - - - -
284	IF( ln_linssh ) THEN
285	IF ( ln_isfcav ) THEN
286	DO ji=1,jpi
287	DO jj=1,jpj
288	ssh_hc_loc_ini(ji,jj) = tsn(ji,jj,mikt(ji,jj),jp_tem) * sshn(ji,jj) ! initial heat content in ssh
289	ssh_sc_loc_ini(ji,jj) = tsn(ji,jj,mikt(ji,jj),jp_sal) * sshn(ji,jj) ! initial salt content in ssh
290	ENDDO
291	ENDDO
292	ELSE
293	ssh_hc_loc_ini(:,:) = tsn(:,:,1,jp_tem) * sshn(:,:) ! initial heat content in ssh
294	ssh_sc_loc_ini(:,:) = tsn(:,:,1,jp_sal) * sshn(:,:) ! initial salt content in ssh
295	END IF
296	frc_wn_t = 0._wp ! initial heat content misfit due to free surface
297	frc_wn_s = 0._wp ! initial salt content misfit due to free surface
298	ENDIF
299	ENDIF
300
301	ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN ! Create restart file
302	! ! -------------------
303	IF(lwp) WRITE(numout,*) '~~~~~~~'
304	IF(lwp) WRITE(numout,*) ' dia_hsb_rst at it= ', kt,' date= ', ndastp
305	IF(lwp) WRITE(numout,*) '~~~~~~~'
306
307	CALL iom_rstput( kt, nitrst, numrow, 'frc_v' , frc_v )
308	CALL iom_rstput( kt, nitrst, numrow, 'frc_t' , frc_t )
309	CALL iom_rstput( kt, nitrst, numrow, 'frc_s' , frc_s )
310	IF( ln_linssh ) THEN
311	CALL iom_rstput( kt, nitrst, numrow, 'frc_wn_t', frc_wn_t )
312	CALL iom_rstput( kt, nitrst, numrow, 'frc_wn_s', frc_wn_s )
313	ENDIF
314	CALL iom_rstput( kt, nitrst, numrow, 'surf_ini', surf_ini ) ! ice sheet coupling
315	CALL iom_rstput( kt, nitrst, numrow, 'ssh_ini', ssh_ini )
316	CALL iom_rstput( kt, nitrst, numrow, 'e3t_ini', e3t_ini )
317	CALL iom_rstput( kt, nitrst, numrow, 'hc_loc_ini', hc_loc_ini )
318	CALL iom_rstput( kt, nitrst, numrow, 'sc_loc_ini', sc_loc_ini )
319	IF( ln_linssh ) THEN
320	CALL iom_rstput( kt, nitrst, numrow, 'ssh_hc_loc_ini', ssh_hc_loc_ini )
321	CALL iom_rstput( kt, nitrst, numrow, 'ssh_sc_loc_ini', ssh_sc_loc_ini )
322	ENDIF
323	!
324	ENDIF
325	!
326	END SUBROUTINE dia_hsb_rst
327
328
329	SUBROUTINE dia_hsb_init
330	!!---------------------------------------------------------------------------
331	!! * ROUTINE dia_hsb *
332	!!
333	!! ** Purpose: Initialization for the heat salt volume budgets
334	!!
335	!! ** Method : Compute initial heat content, salt content and volume
336	!!
337	!! ** Action : - Compute initial heat content, salt content and volume
338	!! - Initialize forcing trends
339	!! - Compute coefficients for conversion
340	!!---------------------------------------------------------------------------
341	INTEGER :: jk ! dummy loop indice
342	INTEGER :: ierror ! local integer
343	INTEGER :: ios
344	!
345	NAMELIST/namhsb/ ln_diahsb
346	!!----------------------------------------------------------------------
347
348	IF(lwp) THEN
349	WRITE(numout,*)
350	WRITE(numout,*) 'dia_hsb_init : check the heat and salt budgets'
351	WRITE(numout,*) '~~~~~~~~ '
352	ENDIF
353
354	REWIND( numnam_ref ) ! Namelist namhsb in reference namelist
355	READ ( numnam_ref, namhsb, IOSTAT = ios, ERR = 901)
356	901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namhsb in reference namelist', lwp )
357
358	REWIND( numnam_cfg ) ! Namelist namhsb in configuration namelist
359	READ ( numnam_cfg, namhsb, IOSTAT = ios, ERR = 902 )
360	902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namhsb in configuration namelist', lwp )
361	IF(lwm) WRITE ( numond, namhsb )
362
363	!
364	IF(lwp) THEN ! Control print
365	WRITE(numout,*)
366	WRITE(numout,*) 'dia_hsb_init : check the heat and salt budgets'
367	WRITE(numout,*) '~~~~~~~~~~~~'
368	WRITE(numout,*) ' Namelist namhsb : set hsb parameters'
369	WRITE(numout,*) ' Switch for hsb diagnostic (T) or not (F) ln_diahsb = ', ln_diahsb
370	WRITE(numout,*)
371	ENDIF
372
373	IF( .NOT. ln_diahsb ) RETURN
374	! IF( .NOT. lk_mpp_rep ) &
375	! CALL ctl_stop (' Your global mpp_sum if performed in single precision - 64 bits -', &
376	! & ' whereas the global sum to be precise must be done in double precision ',&
377	! & ' please add key_mpp_rep')
378
379	! ------------------- !
380	! 1 - Allocate memory !
381	! ------------------- !
382	ALLOCATE( hc_loc_ini(jpi,jpj,jpk), sc_loc_ini(jpi,jpj,jpk), surf_ini(jpi,jpj), &
383	& e3t_ini(jpi,jpj,jpk), surf(jpi,jpj), ssh_ini(jpi,jpj), STAT=ierror )
384	IF( ierror > 0 ) THEN
385	CALL ctl_stop( 'dia_hsb: unable to allocate hc_loc_ini' ) ; RETURN
386	ENDIF
387
388	IF( ln_linssh ) ALLOCATE( ssh_hc_loc_ini(jpi,jpj), ssh_sc_loc_ini(jpi,jpj),STAT=ierror )
389	IF( ierror > 0 ) THEN
390	CALL ctl_stop( 'dia_hsb: unable to allocate hc_loc_ini' ) ; RETURN
391	ENDIF
392
393	! ----------------------------------------------- !
394	! 2 - Time independant variables and file opening !
395	! ----------------------------------------------- !
396	IF(lwp) WRITE(numout,*) "dia_hsb: heat salt volume budgets activated"
397	IF(lwp) WRITE(numout,*) '~~~~~~~'
398	surf(:,:) = e1t(:,:) * e2t(:,:) * tmask_i(:,:) ! masked surface grid cell area
399	surf_tot = glob_sum( surf(:,:) ) ! total ocean surface area
400
401	IF( lk_bdy ) CALL ctl_warn( 'dia_hsb does not take open boundary fluxes into account' )
402	!
403	! ---------------------------------- !
404	! 4 - initial conservation variables !
405	! ---------------------------------- !
406	CALL dia_hsb_rst( nit000, 'READ' ) !* read or initialize all required files
407	!
408	END SUBROUTINE dia_hsb_init
409
410	!!======================================================================
411	END MODULE diahsb

Note: See TracBrowser for help on using the repository browser.

New URL for NEMO forge! http://forge.nemo-ocean.eu

Context Navigation

source: branches/NERC/dev_r5549_BDY_ZEROGRAD/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90 @ 6808

Download in other formats: