Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

diahsb.F90 @ 7924

Last change on this file since 7924 was 7924, checked in by andmirek, 7 years ago
first commit with XIOS restart read functionality
File size: 20.2 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	USE iom_def, ONLY : lxios_read
35
36	IMPLICIT NONE
37	PRIVATE
38
39	PUBLIC dia_hsb ! routine called by step.F90
40	PUBLIC dia_hsb_init ! routine called by nemogcm.F90
41	PUBLIC dia_hsb_rst ! routine called by step.F90
42
43	LOGICAL, PUBLIC :: ln_diahsb !: check the heat and salt budgets
44
45	REAL(wp) :: surf_tot ! ocean surface
46	REAL(wp) :: frc_t, frc_s, frc_v ! global forcing trends
47	REAL(wp) :: frc_wn_t, frc_wn_s ! global forcing trends
48	!
49	REAL(wp), DIMENSION(:,:) , ALLOCATABLE :: surf , 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 "domzgr_substitute.h90"
55	# include "vectopt_loop_substitute.h90"
56	!!----------------------------------------------------------------------
57	!! NEMO/OPA 3.3 , NEMO Consortium (2010)
58	!! $Id$
59	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
60	!!----------------------------------------------------------------------
61	CONTAINS
62
63	SUBROUTINE dia_hsb( kt )
64	!!---------------------------------------------------------------------------
65	!! * ROUTINE dia_hsb *
66	!!
67	!! ** Purpose: Compute the ocean global heat content, salt content and volume conservation
68	!!
69	!! ** Method : - Compute the deviation of heat content, salt content and volume
70	!! at the current time step from their values at nit000
71	!! - Compute the contribution of forcing and remove it from these deviations
72	!!
73	!!---------------------------------------------------------------------------
74	INTEGER, INTENT(in) :: kt ! ocean time-step index
75	!
76	INTEGER :: ji, jj, jk ! dummy loop indice
77	REAL(wp) :: zdiff_hc , zdiff_sc ! heat and salt content variations
78	REAL(wp) :: zdiff_hc1 , zdiff_sc1 ! - - - -
79	REAL(wp) :: zdiff_v1 , zdiff_v2 ! volume variation
80	REAL(wp) :: zerr_hc1 , zerr_sc1 ! heat and salt content misfit
81	REAL(wp) :: zvol_tot ! volume
82	REAL(wp) :: z_frc_trd_t , z_frc_trd_s ! - -
83	REAL(wp) :: z_frc_trd_v ! - -
84	REAL(wp) :: z_wn_trd_t , z_wn_trd_s ! - -
85	REAL(wp) :: z_ssh_hc , z_ssh_sc ! - -
86	REAL(wp), DIMENSION(:,:), POINTER :: z2d0, z2d1
87	!!---------------------------------------------------------------------------
88	IF( nn_timing == 1 ) CALL timing_start('dia_hsb')
89	CALL wrk_alloc( jpi,jpj, z2d0, z2d1 )
90	!
91	tsn(:,:,:,1) = tsn(:,:,:,1) * tmask(:,:,:) ; tsb(:,:,:,1) = tsb(:,:,:,1) * tmask(:,:,:) ;
92	tsn(:,:,:,2) = tsn(:,:,:,2) * tmask(:,:,:) ; tsb(:,:,:,2) = tsb(:,:,:,2) * tmask(:,:,:) ;
93	! ------------------------- !
94	! 1 - Trends due to forcing !
95	! ------------------------- !
96	z_frc_trd_v = r1_rau0 * glob_sum( - ( emp(:,:) - rnf(:,:) + fwfisf(:,:) ) * surf(:,:) ) ! volume fluxes
97	z_frc_trd_t = glob_sum( sbc_tsc(:,:,jp_tem) * surf(:,:) ) ! heat fluxes
98	z_frc_trd_s = glob_sum( sbc_tsc(:,:,jp_sal) * surf(:,:) ) ! salt fluxes
99	! Add runoff heat & salt input
100	IF( ln_rnf ) z_frc_trd_t = z_frc_trd_t + glob_sum( rnf_tsc(:,:,jp_tem) * surf(:,:) )
101	IF( ln_rnf_sal) z_frc_trd_s = z_frc_trd_s + glob_sum( rnf_tsc(:,:,jp_sal) * surf(:,:) )
102	! Add ice shelf heat & salt input
103	IF( nn_isf .GE. 1 ) THEN
104	z_frc_trd_t = z_frc_trd_t + glob_sum( risf_tsc(:,:,jp_tem) * surf(:,:) )
105	z_frc_trd_s = z_frc_trd_s + glob_sum( risf_tsc(:,:,jp_sal) * surf(:,:) )
106	ENDIF
107
108	! Add penetrative solar radiation
109	IF( ln_traqsr ) z_frc_trd_t = z_frc_trd_t + r1_rau0_rcp * glob_sum( qsr (:,:) * surf(:,:) )
110	! Add geothermal heat flux
111	IF( ln_trabbc ) z_frc_trd_t = z_frc_trd_t + glob_sum( qgh_trd0(:,:) * surf(:,:) )
112	!
113	IF( .NOT. lk_vvl ) THEN
114	IF ( ln_isfcav ) THEN
115	DO ji=1,jpi
116	DO jj=1,jpj
117	z2d0(ji,jj) = surf(ji,jj) * wn(ji,jj,mikt(ji,jj)) * tsb(ji,jj,mikt(ji,jj),jp_tem)
118	z2d1(ji,jj) = surf(ji,jj) * wn(ji,jj,mikt(ji,jj)) * tsb(ji,jj,mikt(ji,jj),jp_sal)
119	ENDDO
120	ENDDO
121	ELSE
122	z2d0(:,:) = surf(:,:) * wn(:,:,1) * tsb(:,:,1,jp_tem)
123	z2d1(:,:) = surf(:,:) * wn(:,:,1) * tsb(:,:,1,jp_sal)
124	END IF
125	z_wn_trd_t = - glob_sum( z2d0 )
126	z_wn_trd_s = - glob_sum( z2d1 )
127	ENDIF
128
129	frc_v = frc_v + z_frc_trd_v * rdt
130	frc_t = frc_t + z_frc_trd_t * rdt
131	frc_s = frc_s + z_frc_trd_s * rdt
132	! ! Advection flux through fixed surface (z=0)
133	IF( .NOT. lk_vvl ) THEN
134	frc_wn_t = frc_wn_t + z_wn_trd_t * rdt
135	frc_wn_s = frc_wn_s + z_wn_trd_s * rdt
136	ENDIF
137
138	! ------------------------ !
139	! 2 - Content variations !
140	! ------------------------ !
141	zdiff_v2 = 0._wp
142	zdiff_hc = 0._wp
143	zdiff_sc = 0._wp
144
145	! volume variation (calculated with ssh)
146	zdiff_v1 = glob_sum( surf(:,:) * ( sshn(:,:) - ssh_ini(:,:) ) )
147
148	! heat & salt content variation (associated with ssh)
149	IF( .NOT. lk_vvl ) THEN
150	IF ( ln_isfcav ) THEN
151	DO ji = 1, jpi
152	DO jj = 1, jpj
153	z2d0(ji,jj) = surf(ji,jj) * ( tsn(ji,jj,mikt(ji,jj),jp_tem) * sshn(ji,jj) - ssh_hc_loc_ini(ji,jj) )
154	z2d1(ji,jj) = surf(ji,jj) * ( tsn(ji,jj,mikt(ji,jj),jp_sal) * sshn(ji,jj) - ssh_sc_loc_ini(ji,jj) )
155	END DO
156	END DO
157	ELSE
158	z2d0(:,:) = surf(:,:) * ( tsn(:,:,1,jp_tem) * sshn(:,:) - ssh_hc_loc_ini(:,:) )
159	z2d1(:,:) = surf(:,:) * ( tsn(:,:,1,jp_sal) * sshn(:,:) - ssh_sc_loc_ini(:,:) )
160	END IF
161	z_ssh_hc = glob_sum( z2d0 )
162	z_ssh_sc = glob_sum( z2d1 )
163	ENDIF
164
165	DO jk = 1, jpkm1
166	! volume variation (calculated with scale factors)
167	zdiff_v2 = zdiff_v2 + glob_sum( surf(:,:) * tmask(:,:,jk) &
168	& * ( fse3t_n(:,:,jk) - e3t_ini(:,:,jk) ) )
169	! heat content variation
170	zdiff_hc = zdiff_hc + glob_sum( surf(:,:) * tmask(:,:,jk) &
171	& * ( fse3t_n(:,:,jk) * tsn(:,:,jk,jp_tem) - hc_loc_ini(:,:,jk) ) )
172	! salt content variation
173	zdiff_sc = zdiff_sc + glob_sum( surf(:,:) * tmask(:,:,jk) &
174	& * ( fse3t_n(:,:,jk) * tsn(:,:,jk,jp_sal) - sc_loc_ini(:,:,jk) ) )
175	ENDDO
176
177	! Substract forcing from heat content, salt content and volume variations
178	zdiff_v1 = zdiff_v1 - frc_v
179	IF( lk_vvl ) zdiff_v2 = zdiff_v2 - frc_v
180	zdiff_hc = zdiff_hc - frc_t
181	zdiff_sc = zdiff_sc - frc_s
182	IF( .NOT. lk_vvl ) THEN
183	zdiff_hc1 = zdiff_hc + z_ssh_hc
184	zdiff_sc1 = zdiff_sc + z_ssh_sc
185	zerr_hc1 = z_ssh_hc - frc_wn_t
186	zerr_sc1 = z_ssh_sc - frc_wn_s
187	ENDIF
188
189	! ----------------------- !
190	! 3 - Diagnostics writing !
191	! ----------------------- !
192	zvol_tot = 0._wp ! total ocean volume (calculated with scale factors)
193	DO jk = 1, jpkm1
194	zvol_tot = zvol_tot + glob_sum( surf(:,:) * tmask(:,:,jk) * fse3t_n(:,:,jk) )
195	END DO
196
197	!!gm to be added ?
198	! IF( .NOT. lk_vvl ) THEN ! fixed volume, add the ssh contribution
199	! zvol_tot = zvol_tot + glob_sum( surf(:,:) * sshn(:,:) )
200	! ENDIF
201	!!gm end
202
203	IF( lk_vvl ) THEN
204	CALL iom_put( 'bgtemper' , zdiff_hc / zvol_tot ) ! Temperature variation (C)
205	CALL iom_put( 'bgsaline' , zdiff_sc / zvol_tot ) ! Salinity variation (psu)
206	CALL iom_put( 'bgheatco' , zdiff_hc * 1.e-20 * rau0 * rcp ) ! Heat content variation (1.e20 J)
207	CALL iom_put( 'bgsaltco' , zdiff_sc * 1.e-9 ) ! Salt content variation (psu*km3)
208	CALL iom_put( 'bgvolssh' , zdiff_v1 * 1.e-9 ) ! volume ssh variation (km3)
209	CALL iom_put( 'bgvole3t' , zdiff_v2 * 1.e-9 ) ! volume e3t variation (km3)
210	CALL iom_put( 'bgfrcvol' , frc_v * 1.e-9 ) ! vol - surface forcing (km3)
211	CALL iom_put( 'bgfrctem' , frc_t / zvol_tot ) ! hc - surface forcing (C)
212	CALL iom_put( 'bgfrcsal' , frc_s / zvol_tot ) ! sc - surface forcing (psu)
213	ELSE
214	CALL iom_put( 'bgtemper' , zdiff_hc1 / zvol_tot) ! Heat content variation (C)
215	CALL iom_put( 'bgsaline' , zdiff_sc1 / zvol_tot) ! Salt content variation (psu)
216	CALL iom_put( 'bgheatco' , zdiff_hc1 * 1.e-20 * rau0 * rcp ) ! Heat content variation (1.e20 J)
217	CALL iom_put( 'bgsaltco' , zdiff_sc1 * 1.e-9 ) ! Salt content variation (psu*km3)
218	CALL iom_put( 'bgvolssh' , zdiff_v1 * 1.e-9 ) ! volume ssh variation (km3)
219	CALL iom_put( 'bgfrcvol' , frc_v * 1.e-9 ) ! vol - surface forcing (km3)
220	CALL iom_put( 'bgfrctem' , frc_t / zvol_tot ) ! hc - surface forcing (C)
221	CALL iom_put( 'bgfrcsal' , frc_s / zvol_tot ) ! sc - surface forcing (psu)
222	CALL iom_put( 'bgmistem' , zerr_hc1 / zvol_tot ) ! hc - error due to free surface (C)
223	CALL iom_put( 'bgmissal' , zerr_sc1 / zvol_tot ) ! sc - error due to free surface (psu)
224	ENDIF
225	!
226	IF( lrst_oce ) CALL dia_hsb_rst( kt, 'WRITE' )
227
228	CALL wrk_dealloc( jpi,jpj, z2d0, z2d1 )
229
230	IF( nn_timing == 1 ) CALL timing_stop('dia_hsb')
231	!
232	END SUBROUTINE dia_hsb
233
234
235	SUBROUTINE dia_hsb_rst( kt, cdrw )
236	!!---------------------------------------------------------------------
237	!! * ROUTINE limdia_rst *
238	!!
239	!! ** Purpose : Read or write DIA file in restart file
240	!!
241	!! ** Method : use of IOM library
242	!!----------------------------------------------------------------------
243	INTEGER , INTENT(in) :: kt ! ocean time-step
244	CHARACTER(len=*), INTENT(in) :: cdrw ! "READ"/"WRITE" flag
245	!
246	INTEGER :: ji, jj, jk ! dummy loop indices
247	INTEGER :: id1 ! local integers
248	!!----------------------------------------------------------------------
249	!
250	IF( TRIM(cdrw) == 'READ' ) THEN ! Read/initialise
251	IF( ln_rstart ) THEN !* Read the restart file
252	!id1 = iom_varid( numror, 'frc_vol' , ldstop = .FALSE. )
253	!
254	IF(lwp) WRITE(numout,*) '~~~~~~~'
255	IF(lwp) WRITE(numout,*) ' dia_hsb_rst at it= ', kt,' date= ', ndastp
256	IF(lwp) WRITE(numout,*) '~~~~~~~'
257	CALL iom_get( numror, 'frc_v', frc_v, lrxios = lxios_read )
258	CALL iom_get( numror, 'frc_t', frc_t, lrxios = lxios_read )
259	CALL iom_get( numror, 'frc_s', frc_s, lrxios = lxios_read )
260	IF( .NOT. lk_vvl ) THEN
261	CALL iom_get( numror, 'frc_wn_t', frc_wn_t, lrxios = lxios_read )
262	CALL iom_get( numror, 'frc_wn_s', frc_wn_s, lrxios = lxios_read )
263	ENDIF
264	CALL iom_get( numror, jpdom_autoglo, 'ssh_ini', ssh_ini, lrxios = lxios_read )
265	CALL iom_get( numror, jpdom_autoglo, 'e3t_ini', e3t_ini, lrxios = lxios_read )
266	CALL iom_get( numror, jpdom_autoglo, 'hc_loc_ini', hc_loc_ini, lrxios = lxios_read )
267	CALL iom_get( numror, jpdom_autoglo, 'sc_loc_ini', sc_loc_ini, lrxios = lxios_read )
268	IF( .NOT. lk_vvl ) THEN
269	CALL iom_get( numror, jpdom_autoglo, 'ssh_hc_loc_ini', ssh_hc_loc_ini, lrxios = lxios_read )
270	CALL iom_get( numror, jpdom_autoglo, 'ssh_sc_loc_ini', ssh_sc_loc_ini, lrxios = lxios_read )
271	ENDIF
272	ELSE
273	IF(lwp) WRITE(numout,*) '~~~~~~~'
274	IF(lwp) WRITE(numout,*) ' dia_hsb at initial state '
275	IF(lwp) WRITE(numout,*) '~~~~~~~'
276	ssh_ini(:,:) = sshn(:,:) ! initial ssh
277	DO jk = 1, jpk
278	e3t_ini (:,:,jk) = fse3t_n(:,:,jk) ! initial vertical scale factors
279	hc_loc_ini(:,:,jk) = tsn(:,:,jk,jp_tem) * fse3t_n(:,:,jk) ! initial heat content
280	sc_loc_ini(:,:,jk) = tsn(:,:,jk,jp_sal) * fse3t_n(:,:,jk) ! initial salt content
281	END DO
282	frc_v = 0._wp ! volume trend due to forcing
283	frc_t = 0._wp ! heat content - - - -
284	frc_s = 0._wp ! salt content - - - -
285	IF( .NOT. lk_vvl ) THEN
286	IF ( ln_isfcav ) THEN
287	DO ji=1,jpi
288	DO jj=1,jpj
289	ssh_hc_loc_ini(ji,jj) = tsn(ji,jj,mikt(ji,jj),jp_tem) * sshn(ji,jj) ! initial heat content in ssh
290	ssh_sc_loc_ini(ji,jj) = tsn(ji,jj,mikt(ji,jj),jp_sal) * sshn(ji,jj) ! initial salt content in ssh
291	ENDDO
292	ENDDO
293	ELSE
294	ssh_hc_loc_ini(:,:) = tsn(:,:,1,jp_tem) * sshn(:,:) ! initial heat content in ssh
295	ssh_sc_loc_ini(:,:) = tsn(:,:,1,jp_sal) * sshn(:,:) ! initial salt content in ssh
296	END IF
297	frc_wn_t = 0._wp ! initial heat content misfit due to free surface
298	frc_wn_s = 0._wp ! initial salt content misfit due to free surface
299	ENDIF
300	ENDIF
301
302	ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN ! Create restart file
303	! ! -------------------
304	IF(lwp) WRITE(numout,*) '~~~~~~~'
305	IF(lwp) WRITE(numout,*) ' dia_hsb_rst at it= ', kt,' date= ', ndastp
306	IF(lwp) WRITE(numout,*) '~~~~~~~'
307
308	CALL iom_rstput( kt, nitrst, numrow, 'frc_v' , frc_v )
309	CALL iom_rstput( kt, nitrst, numrow, 'frc_t' , frc_t )
310	CALL iom_rstput( kt, nitrst, numrow, 'frc_s' , frc_s )
311	IF( .NOT. lk_vvl ) THEN
312	CALL iom_rstput( kt, nitrst, numrow, 'frc_wn_t', frc_wn_t )
313	CALL iom_rstput( kt, nitrst, numrow, 'frc_wn_s', frc_wn_s )
314	ENDIF
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( .NOT. lk_vvl ) 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), &
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(.NOT. lk_vvl ) 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/UKMO/dev_r5518_GO6_package_XIOS_read/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90 @ 7924

Download in other formats: