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

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.

diawri.F90 in NEMO/branches/2018/dev_r10057_ENHANCE03_ZTILDE/src/OCE/DIA – NEMO

Context Navigation

source: NEMO/branches/2018/dev_r10057_ENHANCE03_ZTILDE/src/OCE/DIA/diawri.F90 @ 10164

Last change on this file since 10164 was 10164, checked in by jchanut, 6 years ago
ztilde update, #2126: Add higher order filtering + alternative scale factor decomposition
Property svn:keywords set to `Id`
File size: 56.0 KB

Rev	Line
[3]	1	MODULE diawri
	2	!!======================================================================
	3	!! * MODULE diawri *
	4	!! Ocean diagnostics : write ocean output files
	5	!!=====================================================================
[2528]	6	!! History : OPA ! 1991-03 (M.-A. Foujols) Original code
	7	!! 4.0 ! 1991-11 (G. Madec)
	8	!! ! 1992-06 (M. Imbard) correction restart file
	9	!! ! 1992-07 (M. Imbard) split into diawri and rstwri
	10	!! ! 1993-03 (M. Imbard) suppress writibm
	11	!! ! 1998-01 (C. Levy) NETCDF format using ioipsl INTERFACE
	12	!! ! 1999-02 (E. Guilyardi) name of netCDF files + variables
	13	!! 8.2 ! 2000-06 (M. Imbard) Original code (diabort.F)
	14	!! NEMO 1.0 ! 2002-06 (A.Bozec, E. Durand) Original code (diainit.F)
	15	!! - ! 2002-09 (G. Madec) F90: Free form and module
	16	!! - ! 2002-12 (G. Madec) merge of diabort and diainit, F90
	17	!! ! 2005-11 (V. Garnier) Surface pressure gradient organization
	18	!! 3.2 ! 2008-11 (B. Lemaire) creation from old diawri
[5836]	19	!! 3.7 ! 2014-01 (G. Madec) remove eddy induced velocity from no-IOM output
	20	!! ! change name of output variables in dia_wri_state
[2528]	21	!!----------------------------------------------------------------------
[3]	22
	23	!!----------------------------------------------------------------------
[2528]	24	!! dia_wri : create the standart output files
	25	!! dia_wri_state : create an output NetCDF file for a single instantaeous ocean state and forcing fields
	26	!!----------------------------------------------------------------------
[9019]	27	USE oce ! ocean dynamics and tracers
	28	USE dom_oce ! ocean space and time domain
	29	USE phycst ! physical constants
	30	USE dianam ! build name of file (routine)
	31	USE diahth ! thermocline diagnostics
	32	USE dynadv , ONLY: ln_dynadv_vec
	33	USE icb_oce ! Icebergs
	34	USE icbdia ! Iceberg budgets
	35	USE ldftra ! lateral physics: eddy diffusivity coef.
	36	USE ldfdyn ! lateral physics: eddy viscosity coef.
	37	USE sbc_oce ! Surface boundary condition: ocean fields
	38	USE sbc_ice ! Surface boundary condition: ice fields
	39	USE sbcssr ! restoring term toward SST/SSS climatology
	40	USE sbcwave ! wave parameters
	41	USE wet_dry ! wetting and drying
	42	USE zdf_oce ! ocean vertical physics
	43	USE zdfdrg ! ocean vertical physics: top/bottom friction
	44	USE zdfmxl ! mixed layer
[6140]	45	!
[9019]	46	USE lbclnk ! ocean lateral boundary conditions (or mpp link)
	47	USE in_out_manager ! I/O manager
	48	USE diatmb ! Top,middle,bottom output
	49	USE dia25h ! 25h Mean output
	50	USE iom !
	51	USE ioipsl !
[5463]	52
[9570]	53	#if defined key_si3
[9019]	54	USE icewri
[1482]	55	#endif
[2715]	56	USE lib_mpp ! MPP library
[3294]	57	USE timing ! preformance summary
[6140]	58	USE diurnal_bulk ! diurnal warm layer
	59	USE cool_skin ! Cool skin
[10164]	60	USE domvvl
[2528]	61
[3]	62	IMPLICIT NONE
	63	PRIVATE
	64
[2528]	65	PUBLIC dia_wri ! routines called by step.F90
	66	PUBLIC dia_wri_state
[2715]	67	PUBLIC dia_wri_alloc ! Called by nemogcm module
[3]	68
[2528]	69	INTEGER :: nid_T, nz_T, nh_T, ndim_T, ndim_hT ! grid_T file
[3609]	70	INTEGER :: nb_T , ndim_bT ! grid_T file
[2528]	71	INTEGER :: nid_U, nz_U, nh_U, ndim_U, ndim_hU ! grid_U file
	72	INTEGER :: nid_V, nz_V, nh_V, ndim_V, ndim_hV ! grid_V file
	73	INTEGER :: nid_W, nz_W, nh_W ! grid_W file
	74	INTEGER :: ndex(1) ! ???
[2715]	75	INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) :: ndex_hT, ndex_hU, ndex_hV
	76	INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) :: ndex_T, ndex_U, ndex_V
[3609]	77	INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) :: ndex_bT
[3]	78
	79	!! * Substitutions
[1756]	80	# include "vectopt_loop_substitute.h90"
[3]	81	!!----------------------------------------------------------------------
[9598]	82	!! NEMO/OCE 4.0 , NEMO Consortium (2018)
[5217]	83	!! $Id$
[10126]	84	!! Software governed by the CeCILL license (see ./LICENSE)
[3]	85	!!----------------------------------------------------------------------
	86	CONTAINS
	87
[6140]	88	#if defined key_iomput
[3]	89	!!----------------------------------------------------------------------
[2528]	90	!! 'key_iomput' use IOM library
	91	!!----------------------------------------------------------------------
[9652]	92	INTEGER FUNCTION dia_wri_alloc()
	93	!
	94	dia_wri_alloc = 0
	95	!
	96	END FUNCTION dia_wri_alloc
[2715]	97
[9652]	98
[1561]	99	SUBROUTINE dia_wri( kt )
[1482]	100	!!---------------------------------------------------------------------
	101	!! * ROUTINE dia_wri *
	102	!!
	103	!! ** Purpose : Standard output of opa: dynamics and tracer fields
	104	!! NETCDF format is used by default
	105	!!
	106	!! ** Method : use iom_put
	107	!!----------------------------------------------------------------------
	108	INTEGER, INTENT( in ) :: kt ! ocean time-step index
[1756]	109	!!
[9019]	110	INTEGER :: ji, jj, jk ! dummy loop indices
	111	INTEGER :: ikbot ! local integer
	112	REAL(wp):: zztmp , zztmpx ! local scalar
	113	REAL(wp):: zztmp2, zztmpy ! - -
	114	REAL(wp), DIMENSION(jpi,jpj) :: z2d ! 2D workspace
	115	REAL(wp), DIMENSION(jpi,jpj,jpk) :: z3d ! 3D workspace
[1482]	116	!!----------------------------------------------------------------------
	117	!
[9124]	118	IF( ln_timing ) CALL timing_start('dia_wri')
[3294]	119	!
[1482]	120	! Output the initial state and forcings
	121	IF( ninist == 1 ) THEN
	122	CALL dia_wri_state( 'output.init', kt )
	123	ninist = 0
	124	ENDIF
[3]	125
[6351]	126	! Output of initial vertical scale factor
[10164]	127	CALL iom_put("e3t_0", e3t_0(:,:,:))
[10116]	128	CALL iom_put("e3u_0", e3u_0(:,:,:) )
	129	CALL iom_put("e3v_0", e3v_0(:,:,:) )
[6351]	130	!
[6387]	131	CALL iom_put( "e3t" , e3t_n(:,:,:) )
	132	CALL iom_put( "e3u" , e3u_n(:,:,:) )
	133	CALL iom_put( "e3v" , e3v_n(:,:,:) )
	134	CALL iom_put( "e3w" , e3w_n(:,:,:) )
[10164]	135	!
	136	IF ( (ln_vvl_ztilde).OR.(ln_vvl_layer) ) THEN
	137	IF( iom_use("e3t_star") ) THEN ! Barotropic cell thickness anomaly
	138	z3d(:,:,:) = (e3t_n(:,:,:)-tilde_e3t_n(:,:,:)-e3t_0(:,:,:))*tmask(:,:,:)
	139	CALL iom_put( "e3t_star" , z3d(:,:,:) )
	140	ENDIF
	141	IF( iom_use("e3t_tilde") ) THEN ! Baroclinic cell thickness anomaly
	142	CALL iom_put( "e3t_tilde" , tilde_e3t_n(:,:,:) )
	143	ENDIF
	144	ENDIF
	145
[6351]	146	IF( iom_use("e3tdef") ) &
[6387]	147	CALL iom_put( "e3tdef" , ( ( e3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 )
[5461]	148
[9023]	149	IF( ll_wd ) THEN
	150	CALL iom_put( "ssh" , (sshn+ssh_ref)*tmask(:,:,1) ) ! sea surface height (brought back to the reference used for wetting and drying)
	151	ELSE
	152	CALL iom_put( "ssh" , sshn ) ! sea surface height
	153	ENDIF
	154
[7646]	155	IF( iom_use("wetdep") ) & ! wet depth
[9023]	156	CALL iom_put( "wetdep" , ht_0(:,:) + sshn(:,:) )
[5107]	157
	158	CALL iom_put( "toce", tsn(:,:,:,jp_tem) ) ! 3D temperature
	159	CALL iom_put( "sst", tsn(:,:,1,jp_tem) ) ! surface temperature
	160	IF ( iom_use("sbt") ) THEN
[4990]	161	DO jj = 1, jpj
	162	DO ji = 1, jpi
[9019]	163	ikbot = mbkt(ji,jj)
	164	z2d(ji,jj) = tsn(ji,jj,ikbot,jp_tem)
[4990]	165	END DO
[5107]	166	END DO
	167	CALL iom_put( "sbt", z2d ) ! bottom temperature
	168	ENDIF
	169
	170	CALL iom_put( "soce", tsn(:,:,:,jp_sal) ) ! 3D salinity
	171	CALL iom_put( "sss", tsn(:,:,1,jp_sal) ) ! surface salinity
	172	IF ( iom_use("sbs") ) THEN
[4990]	173	DO jj = 1, jpj
	174	DO ji = 1, jpi
[9019]	175	ikbot = mbkt(ji,jj)
	176	z2d(ji,jj) = tsn(ji,jj,ikbot,jp_sal)
[4990]	177	END DO
[5107]	178	END DO
	179	CALL iom_put( "sbs", z2d ) ! bottom salinity
	180	ENDIF
[5463]	181
	182	IF ( iom_use("taubot") ) THEN ! bottom stress
[9019]	183	zztmp = rau0 * 0.25
[7753]	184	z2d(:,:) = 0._wp
[5463]	185	DO jj = 2, jpjm1
	186	DO ji = fs_2, fs_jpim1 ! vector opt.
[9019]	187	zztmp2 = ( ( rCdU_bot(ji+1,jj)+rCdU_bot(ji ,jj) ) * un(ji ,jj,mbku(ji ,jj)) )**2 &
	188	& + ( ( rCdU_bot(ji ,jj)+rCdU_bot(ji-1,jj) ) * un(ji-1,jj,mbku(ji-1,jj)) )**2 &
	189	& + ( ( rCdU_bot(ji,jj+1)+rCdU_bot(ji,jj ) ) * vn(ji,jj ,mbkv(ji,jj )) )**2 &
	190	& + ( ( rCdU_bot(ji,jj )+rCdU_bot(ji,jj-1) ) * vn(ji,jj-1,mbkv(ji,jj-1)) )**2
	191	z2d(ji,jj) = zztmp * SQRT( zztmp2 ) * tmask(ji,jj,1)
[5463]	192	!
[9019]	193	END DO
	194	END DO
[5463]	195	CALL lbc_lnk( z2d, 'T', 1. )
	196	CALL iom_put( "taubot", z2d )
	197	ENDIF
[5107]	198
[9019]	199	CALL iom_put( "uoce", un(:,:,:) ) ! 3D i-current
	200	CALL iom_put( "ssu", un(:,:,1) ) ! surface i-current
[5107]	201	IF ( iom_use("sbu") ) THEN
[4990]	202	DO jj = 1, jpj
	203	DO ji = 1, jpi
[9019]	204	ikbot = mbku(ji,jj)
	205	z2d(ji,jj) = un(ji,jj,ikbot)
[4990]	206	END DO
[5107]	207	END DO
	208	CALL iom_put( "sbu", z2d ) ! bottom i-current
	209	ENDIF
	210
[9019]	211	CALL iom_put( "voce", vn(:,:,:) ) ! 3D j-current
	212	CALL iom_put( "ssv", vn(:,:,1) ) ! surface j-current
[5107]	213	IF ( iom_use("sbv") ) THEN
[4990]	214	DO jj = 1, jpj
	215	DO ji = 1, jpi
[9019]	216	ikbot = mbkv(ji,jj)
	217	z2d(ji,jj) = vn(ji,jj,ikbot)
[4990]	218	END DO
[5107]	219	END DO
	220	CALL iom_put( "sbv", z2d ) ! bottom j-current
[4990]	221	ENDIF
[1482]	222
[5461]	223	CALL iom_put( "woce", wn ) ! vertical velocity
	224	IF( iom_use('w_masstr') .OR. iom_use('w_masstr2') ) THEN ! vertical mass transport & its square value
	225	! Caution: in the VVL case, it only correponds to the baroclinic mass transport.
[7753]	226	z2d(:,:) = rau0 * e1e2t(:,:)
[5461]	227	DO jk = 1, jpk
[7753]	228	z3d(:,:,jk) = wn(:,:,jk) * z2d(:,:)
[5461]	229	END DO
	230	CALL iom_put( "w_masstr" , z3d )
	231	IF( iom_use('w_masstr2') ) CALL iom_put( "w_masstr2", z3d(:,:,:) * z3d(:,:,:) )
	232	ENDIF
	233
[9019]	234	CALL iom_put( "avt" , avt ) ! T vert. eddy diff. coef.
	235	CALL iom_put( "avs" , avs ) ! S vert. eddy diff. coef.
	236	CALL iom_put( "avm" , avm ) ! T vert. eddy visc. coef.
[5107]	237
[9019]	238	IF( iom_use('logavt') ) CALL iom_put( "logavt", LOG( MAX( 1.e-20_wp, avt(:,:,:) ) ) )
	239	IF( iom_use('logavs') ) CALL iom_put( "logavs", LOG( MAX( 1.e-20_wp, avs(:,:,:) ) ) )
[6351]	240
[5107]	241	IF ( iom_use("sstgrad") .OR. iom_use("sstgrad2") ) THEN
[4990]	242	DO jj = 2, jpjm1 ! sst gradient
	243	DO ji = fs_2, fs_jpim1 ! vector opt.
[5836]	244	zztmp = tsn(ji,jj,1,jp_tem)
	245	zztmpx = ( tsn(ji+1,jj,1,jp_tem) - zztmp ) * r1_e1u(ji,jj) + ( zztmp - tsn(ji-1,jj ,1,jp_tem) ) * r1_e1u(ji-1,jj)
	246	zztmpy = ( tsn(ji,jj+1,1,jp_tem) - zztmp ) * r1_e2v(ji,jj) + ( zztmp - tsn(ji ,jj-1,1,jp_tem) ) * r1_e2v(ji,jj-1)
[4990]	247	z2d(ji,jj) = 0.25 * ( zztmpx * zztmpx + zztmpy * zztmpy ) &
	248	& * umask(ji,jj,1) * umask(ji-1,jj,1) * vmask(ji,jj,1) * umask(ji,jj-1,1)
	249	END DO
[1756]	250	END DO
[4990]	251	CALL lbc_lnk( z2d, 'T', 1. )
[9019]	252	CALL iom_put( "sstgrad2", z2d ) ! square of module of sst gradient
[7753]	253	z2d(:,:) = SQRT( z2d(:,:) )
[9019]	254	CALL iom_put( "sstgrad" , z2d ) ! module of sst gradient
[4990]	255	ENDIF
	256
[9019]	257	! heat and salt contents
[4990]	258	IF( iom_use("heatc") ) THEN
[7753]	259	z2d(:,:) = 0._wp
[4990]	260	DO jk = 1, jpkm1
[5107]	261	DO jj = 1, jpj
	262	DO ji = 1, jpi
[6140]	263	z2d(ji,jj) = z2d(ji,jj) + e3t_n(ji,jj,jk) * tsn(ji,jj,jk,jp_tem) * tmask(ji,jj,jk)
[4990]	264	END DO
[4761]	265	END DO
	266	END DO
[9019]	267	CALL iom_put( "heatc", rau0_rcp * z2d ) ! vertically integrated heat content (J/m2)
[4990]	268	ENDIF
	269
	270	IF( iom_use("saltc") ) THEN
[7753]	271	z2d(:,:) = 0._wp
[4990]	272	DO jk = 1, jpkm1
[5107]	273	DO jj = 1, jpj
	274	DO ji = 1, jpi
[6140]	275	z2d(ji,jj) = z2d(ji,jj) + e3t_n(ji,jj,jk) * tsn(ji,jj,jk,jp_sal) * tmask(ji,jj,jk)
[4990]	276	END DO
	277	END DO
	278	END DO
[9019]	279	CALL iom_put( "saltc", rau0 * z2d ) ! vertically integrated salt content (PSU*kg/m2)
[4990]	280	ENDIF
[4840]	281	!
[4990]	282	IF ( iom_use("eken") ) THEN
[9399]	283	z3d(:,:,jpk) = 0._wp
[4990]	284	DO jk = 1, jpkm1
	285	DO jj = 2, jpjm1
	286	DO ji = fs_2, fs_jpim1 ! vector opt.
[9019]	287	zztmp = 0.25_wp * r1_e1e2t(ji,jj) / e3t_n(ji,jj,jk)
	288	z3d(ji,jj,jk) = zztmp * ( un(ji-1,jj,jk)*2 e2u(ji-1,jj) * e3u_n(ji-1,jj,jk) &
	289	& + un(ji ,jj,jk)*2 e2u(ji ,jj) * e3u_n(ji ,jj,jk) &
	290	& + vn(ji,jj-1,jk)*2 e1v(ji,jj-1) * e3v_n(ji,jj-1,jk) &
	291	& + vn(ji,jj ,jk)*2 e1v(ji,jj ) * e3v_n(ji,jj ,jk) )
	292	END DO
	293	END DO
	294	END DO
	295	CALL lbc_lnk( z3d, 'T', 1. )
	296	CALL iom_put( "eken", z3d ) ! kinetic energy
[4990]	297	ENDIF
[6351]	298	!
	299	CALL iom_put( "hdiv", hdivn ) ! Horizontal divergence
	300	!
[7646]	301	IF( iom_use("u_masstr") .OR. iom_use("u_masstr_vint") .OR. iom_use("u_heattr") .OR. iom_use("u_salttr") ) THEN
[7753]	302	z3d(:,:,jpk) = 0.e0
	303	z2d(:,:) = 0.e0
[1756]	304	DO jk = 1, jpkm1
[7753]	305	z3d(:,:,jk) = rau0 * un(:,:,jk) * e2u(:,:) * e3u_n(:,:,jk) * umask(:,:,jk)
	306	z2d(:,:) = z2d(:,:) + z3d(:,:,jk)
[1756]	307	END DO
[9019]	308	CALL iom_put( "u_masstr" , z3d ) ! mass transport in i-direction
	309	CALL iom_put( "u_masstr_vint", z2d ) ! mass transport in i-direction vertical sum
[4990]	310	ENDIF
	311
	312	IF( iom_use("u_heattr") ) THEN
[9019]	313	z2d(:,:) = 0._wp
[4990]	314	DO jk = 1, jpkm1
	315	DO jj = 2, jpjm1
	316	DO ji = fs_2, fs_jpim1 ! vector opt.
	317	z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( tsn(ji,jj,jk,jp_tem) + tsn(ji+1,jj,jk,jp_tem) )
	318	END DO
	319	END DO
	320	END DO
	321	CALL lbc_lnk( z2d, 'U', -1. )
[9019]	322	CALL iom_put( "u_heattr", 0.5rcp z2d ) ! heat transport in i-direction
[4990]	323	ENDIF
[4761]	324
[4990]	325	IF( iom_use("u_salttr") ) THEN
[7753]	326	z2d(:,:) = 0.e0
[1756]	327	DO jk = 1, jpkm1
	328	DO jj = 2, jpjm1
	329	DO ji = fs_2, fs_jpim1 ! vector opt.
[4990]	330	z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( tsn(ji,jj,jk,jp_sal) + tsn(ji+1,jj,jk,jp_sal) )
[1756]	331	END DO
	332	END DO
	333	END DO
	334	CALL lbc_lnk( z2d, 'U', -1. )
[9019]	335	CALL iom_put( "u_salttr", 0.5 * z2d ) ! heat transport in i-direction
[4990]	336	ENDIF
[4761]	337
[4990]	338
	339	IF( iom_use("v_masstr") .OR. iom_use("v_heattr") .OR. iom_use("v_salttr") ) THEN
[7753]	340	z3d(:,:,jpk) = 0.e0
[1756]	341	DO jk = 1, jpkm1
[7753]	342	z3d(:,:,jk) = rau0 * vn(:,:,jk) * e1v(:,:) * e3v_n(:,:,jk) * vmask(:,:,jk)
[1756]	343	END DO
[9019]	344	CALL iom_put( "v_masstr", z3d ) ! mass transport in j-direction
[4990]	345	ENDIF
	346
	347	IF( iom_use("v_heattr") ) THEN
[7753]	348	z2d(:,:) = 0.e0
[4990]	349	DO jk = 1, jpkm1
	350	DO jj = 2, jpjm1
	351	DO ji = fs_2, fs_jpim1 ! vector opt.
	352	z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( tsn(ji,jj,jk,jp_tem) + tsn(ji,jj+1,jk,jp_tem) )
	353	END DO
	354	END DO
	355	END DO
	356	CALL lbc_lnk( z2d, 'V', -1. )
[9019]	357	CALL iom_put( "v_heattr", 0.5rcp z2d ) ! heat transport in j-direction
[4990]	358	ENDIF
[4761]	359
[4990]	360	IF( iom_use("v_salttr") ) THEN
[9019]	361	z2d(:,:) = 0._wp
[1756]	362	DO jk = 1, jpkm1
	363	DO jj = 2, jpjm1
	364	DO ji = fs_2, fs_jpim1 ! vector opt.
[4990]	365	z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( tsn(ji,jj,jk,jp_sal) + tsn(ji,jj+1,jk,jp_sal) )
[1756]	366	END DO
	367	END DO
	368	END DO
	369	CALL lbc_lnk( z2d, 'V', -1. )
[9019]	370	CALL iom_put( "v_salttr", 0.5 * z2d ) ! heat transport in j-direction
[1756]	371	ENDIF
[7646]	372
	373	IF( iom_use("tosmint") ) THEN
[9019]	374	z2d(:,:) = 0._wp
[7646]	375	DO jk = 1, jpkm1
	376	DO jj = 2, jpjm1
	377	DO ji = fs_2, fs_jpim1 ! vector opt.
[9019]	378	z2d(ji,jj) = z2d(ji,jj) + e3t_n(ji,jj,jk) * tsn(ji,jj,jk,jp_tem)
[7646]	379	END DO
	380	END DO
	381	END DO
	382	CALL lbc_lnk( z2d, 'T', -1. )
[9019]	383	CALL iom_put( "tosmint", rau0 * z2d ) ! Vertical integral of temperature
[7646]	384	ENDIF
	385	IF( iom_use("somint") ) THEN
[7753]	386	z2d(:,:)=0._wp
[7646]	387	DO jk = 1, jpkm1
	388	DO jj = 2, jpjm1
	389	DO ji = fs_2, fs_jpim1 ! vector opt.
[9019]	390	z2d(ji,jj) = z2d(ji,jj) + e3t_n(ji,jj,jk) * tsn(ji,jj,jk,jp_sal)
[7646]	391	END DO
	392	END DO
	393	END DO
	394	CALL lbc_lnk( z2d, 'T', -1. )
[9019]	395	CALL iom_put( "somint", rau0 * z2d ) ! Vertical integral of salinity
[7646]	396	ENDIF
	397
[9019]	398	CALL iom_put( "bn2", rn2 ) ! Brunt-Vaisala buoyancy frequency (N^2)
[2528]	399	!
[6140]	400
[9019]	401	IF (ln_diatmb) CALL dia_tmb ! tmb values
	402
	403	IF (ln_dia25h) CALL dia_25h( kt ) ! 25h averaging
[6140]	404
[9124]	405	IF( ln_timing ) CALL timing_stop('dia_wri')
[3294]	406	!
[1482]	407	END SUBROUTINE dia_wri
	408
	409	#else
[2528]	410	!!----------------------------------------------------------------------
	411	!! Default option use IOIPSL library
	412	!!----------------------------------------------------------------------
	413
[9652]	414	INTEGER FUNCTION dia_wri_alloc()
	415	!!----------------------------------------------------------------------
	416	INTEGER, DIMENSION(2) :: ierr
	417	!!----------------------------------------------------------------------
	418	ierr = 0
	419	ALLOCATE( ndex_hT(jpijpj) , ndex_T(jpijpj*jpk) , &
	420	& ndex_hU(jpijpj) , ndex_U(jpijpj*jpk) , &
	421	& ndex_hV(jpijpj) , ndex_V(jpijpj*jpk) , STAT=ierr(1) )
	422	!
	423	dia_wri_alloc = MAXVAL(ierr)
	424	IF( lk_mpp ) CALL mpp_sum( dia_wri_alloc )
	425	!
	426	END FUNCTION dia_wri_alloc
	427
	428
[1561]	429	SUBROUTINE dia_wri( kt )
[3]	430	!!---------------------------------------------------------------------
	431	!! * ROUTINE dia_wri *
	432	!!
	433	!! ** Purpose : Standard output of opa: dynamics and tracer fields
	434	!! NETCDF format is used by default
	435	!!
	436	!! ** Method : At the beginning of the first time step (nit000),
	437	!! define all the NETCDF files and fields
	438	!! At each time step call histdef to compute the mean if ncessary
	439	!! Each nwrite time step, output the instantaneous or mean fields
	440	!!----------------------------------------------------------------------
[5836]	441	INTEGER, INTENT( in ) :: kt ! ocean time-step index
	442	!
[2528]	443	LOGICAL :: ll_print = .FALSE. ! =T print and flush numout
	444	CHARACTER (len=40) :: clhstnam, clop, clmx ! local names
	445	INTEGER :: inum = 11 ! temporary logical unit
[3294]	446	INTEGER :: ji, jj, jk ! dummy loop indices
	447	INTEGER :: ierr ! error code return from allocation
[2528]	448	INTEGER :: iimi, iima, ipk, it, itmod, ijmi, ijma ! local integers
[3609]	449	INTEGER :: jn, ierror ! local integers
[6140]	450	REAL(wp) :: zsto, zout, zmax, zjulian ! local scalars
[5836]	451	!
[9019]	452	REAL(wp), DIMENSION(jpi,jpj) :: zw2d ! 2D workspace
	453	REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d ! 3D workspace
[3]	454	!!----------------------------------------------------------------------
[3294]	455	!
[9124]	456	IF( ln_timing ) CALL timing_start('dia_wri')
[1482]	457	!
[9019]	458	IF( ninist == 1 ) THEN !== Output the initial state and forcings ==!
[1482]	459	CALL dia_wri_state( 'output.init', kt )
	460	ninist = 0
	461	ENDIF
	462	!
[3]	463	! 0. Initialisation
	464	! -----------------
[632]	465
[9019]	466	ll_print = .FALSE. ! local variable for debugging
[3]	467	ll_print = ll_print .AND. lwp
	468
	469	! Define frequency of output and means
[5566]	470	clop = "x" ! no use of the mask value (require less cpu time and otherwise the model crashes)
[3]	471	#if defined key_diainstant
[6140]	472	zsto = nwrite * rdt
[1312]	473	clop = "inst("//TRIM(clop)//")"
[3]	474	#else
[6140]	475	zsto=rdt
[1312]	476	clop = "ave("//TRIM(clop)//")"
[3]	477	#endif
[6140]	478	zout = nwrite * rdt
	479	zmax = ( nitend - nit000 + 1 ) * rdt
[3]	480
	481	! Define indices of the horizontal output zoom and vertical limit storage
	482	iimi = 1 ; iima = jpi
	483	ijmi = 1 ; ijma = jpj
	484	ipk = jpk
	485
	486	! define time axis
[1334]	487	it = kt
	488	itmod = kt - nit000 + 1
[3]	489
	490
	491	! 1. Define NETCDF files and fields at beginning of first time step
	492	! -----------------------------------------------------------------
	493
	494	IF( kt == nit000 ) THEN
	495
	496	! Define the NETCDF files (one per grid)
[632]	497
[3]	498	! Compute julian date from starting date of the run
[1309]	499	CALL ymds2ju( nyear, nmonth, nday, rdt, zjulian )
	500	zjulian = zjulian - adatrj ! set calendar origin to the beginning of the experiment
[3]	501	IF(lwp)WRITE(numout,*)
	502	IF(lwp)WRITE(numout,*) 'Date 0 used :', nit000, ' YEAR ', nyear, &
	503	& ' MONTH ', nmonth, ' DAY ', nday, 'Julian day : ', zjulian
	504	IF(lwp)WRITE(numout,*) ' indexes of zoom = ', iimi, iima, ijmi, ijma, &
	505	' limit storage in depth = ', ipk
	506
	507	! WRITE root name in date.file for use by postpro
[1581]	508	IF(lwp) THEN
[895]	509	CALL dia_nam( clhstnam, nwrite,' ' )
[1581]	510	CALL ctl_opn( inum, 'date.file', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea )
[895]	511	WRITE(inum,*) clhstnam
	512	CLOSE(inum)
	513	ENDIF
[632]	514
[3]	515	! Define the T grid FILE ( nid_T )
[632]	516
[3]	517	CALL dia_nam( clhstnam, nwrite, 'grid_T' )
	518	IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam ! filename
	519	CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit
	520	& iimi, iima-iimi+1, ijmi, ijma-ijmi+1, &
[6140]	521	& nit000-1, zjulian, rdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set )
[3]	522	CALL histvert( nid_T, "deptht", "Vertical T levels", & ! Vertical grid: gdept
[4292]	523	& "m", ipk, gdept_1d, nz_T, "down" )
[3]	524	! ! Index of ocean points
	525	CALL wheneq( jpijpjipk, tmask, 1, 1., ndex_T , ndim_T ) ! volume
	526	CALL wheneq( jpi*jpj , tmask, 1, 1., ndex_hT, ndim_hT ) ! surface
[3609]	527	!
	528	IF( ln_icebergs ) THEN
	529	!
	530	!! allocation cant go in dia_wri_alloc because ln_icebergs is only set after
	531	!! that routine is called from nemogcm, so do it here immediately before its needed
	532	ALLOCATE( ndex_bT(jpijpjnclasses), STAT=ierror )
	533	IF( lk_mpp ) CALL mpp_sum( ierror )
	534	IF( ierror /= 0 ) THEN
	535	CALL ctl_stop('dia_wri: failed to allocate iceberg diagnostic array')
	536	RETURN
	537	ENDIF
	538	!
	539	!! iceberg vertical coordinate is class number
	540	CALL histvert( nid_T, "class", "Iceberg class", & ! Vertical grid: class
	541	& "number", nclasses, class_num, nb_T )
	542	!
	543	!! each class just needs the surface index pattern
	544	ndim_bT = 3
	545	DO jn = 1,nclasses
	546	ndex_bT((jn-1)jpijpj+1:jnjpijpj) = ndex_hT(1:jpi*jpj)
	547	ENDDO
	548	!
	549	ENDIF
[3]	550
	551	! Define the U grid FILE ( nid_U )
	552
	553	CALL dia_nam( clhstnam, nwrite, 'grid_U' )
	554	IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam ! filename
	555	CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu, & ! Horizontal grid: glamu and gphiu
	556	& iimi, iima-iimi+1, ijmi, ijma-ijmi+1, &
[6140]	557	& nit000-1, zjulian, rdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set )
[3]	558	CALL histvert( nid_U, "depthu", "Vertical U levels", & ! Vertical grid: gdept
[4292]	559	& "m", ipk, gdept_1d, nz_U, "down" )
[3]	560	! ! Index of ocean points
	561	CALL wheneq( jpijpjipk, umask, 1, 1., ndex_U , ndim_U ) ! volume
	562	CALL wheneq( jpi*jpj , umask, 1, 1., ndex_hU, ndim_hU ) ! surface
	563
	564	! Define the V grid FILE ( nid_V )
	565
	566	CALL dia_nam( clhstnam, nwrite, 'grid_V' ) ! filename
	567	IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam
	568	CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv, & ! Horizontal grid: glamv and gphiv
	569	& iimi, iima-iimi+1, ijmi, ijma-ijmi+1, &
[6140]	570	& nit000-1, zjulian, rdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set )
[3]	571	CALL histvert( nid_V, "depthv", "Vertical V levels", & ! Vertical grid : gdept
[4292]	572	& "m", ipk, gdept_1d, nz_V, "down" )
[3]	573	! ! Index of ocean points
	574	CALL wheneq( jpijpjipk, vmask, 1, 1., ndex_V , ndim_V ) ! volume
	575	CALL wheneq( jpi*jpj , vmask, 1, 1., ndex_hV, ndim_hV ) ! surface
	576
	577	! Define the W grid FILE ( nid_W )
	578
	579	CALL dia_nam( clhstnam, nwrite, 'grid_W' ) ! filename
	580	IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam
	581	CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit
	582	& iimi, iima-iimi+1, ijmi, ijma-ijmi+1, &
[6140]	583	& nit000-1, zjulian, rdt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set )
[3]	584	CALL histvert( nid_W, "depthw", "Vertical W levels", & ! Vertical grid: gdepw
[4292]	585	& "m", ipk, gdepw_1d, nz_W, "down" )
[3]	586
[632]	587
[3]	588	! Declare all the output fields as NETCDF variables
	589
	590	! !!! nid_T : 3D
	591	CALL histdef( nid_T, "votemper", "Temperature" , "C" , & ! tn
	592	& jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout )
	593	CALL histdef( nid_T, "vosaline", "Salinity" , "PSU" , & ! sn
	594	& jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout )
[6140]	595	IF( .NOT.ln_linssh ) THEN
[4292]	596	CALL histdef( nid_T, "vovvle3t", "Level thickness" , "m" ,& ! e3t_n
	597	& jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout )
	598	CALL histdef( nid_T, "vovvldep", "T point depth" , "m" ,& ! e3t_n
	599	& jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout )
	600	CALL histdef( nid_T, "vovvldef", "Squared level deformation" , "%^2" ,& ! e3t_n
	601	& jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout )
	602	ENDIF
[3]	603	! !!! nid_T : 2D
	604	CALL histdef( nid_T, "sosstsst", "Sea Surface temperature" , "C" , & ! sst
	605	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	606	CALL histdef( nid_T, "sosaline", "Sea Surface Salinity" , "PSU" , & ! sss
	607	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[359]	608	CALL histdef( nid_T, "sossheig", "Sea Surface Height" , "m" , & ! ssh
[3]	609	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[2528]	610	CALL histdef( nid_T, "sowaflup", "Net Upward Water Flux" , "Kg/m2/s", & ! (emp-rnf)
[3]	611	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[4570]	612	CALL histdef( nid_T, "sorunoff", "River runoffs" , "Kg/m2/s", & ! runoffs
	613	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[3625]	614	CALL histdef( nid_T, "sosfldow", "downward salt flux" , "PSU/m2/s", & ! sfx
[3]	615	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[6140]	616	IF( ln_linssh ) THEN
[4292]	617	CALL histdef( nid_T, "sosst_cd", "Concentration/Dilution term on temperature" & ! emp * tsn(:,:,1,jp_tem)
[3625]	618	& , "KgC/m2/s", & ! sosst_cd
[4292]	619	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	620	CALL histdef( nid_T, "sosss_cd", "Concentration/Dilution term on salinity" & ! emp * tsn(:,:,1,jp_sal)
[3625]	621	& , "KgPSU/m2/s",& ! sosss_cd
[4292]	622	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	623	ENDIF
[888]	624	CALL histdef( nid_T, "sohefldo", "Net Downward Heat Flux" , "W/m2" , & ! qns + qsr
[3]	625	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	626	CALL histdef( nid_T, "soshfldo", "Shortwave Radiation" , "W/m2" , & ! qsr
	627	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[1585]	628	CALL histdef( nid_T, "somixhgt", "Turbocline Depth" , "m" , & ! hmld
	629	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[3]	630	CALL histdef( nid_T, "somxl010", "Mixed Layer Depth 0.01" , "m" , & ! hmlp
	631	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[1037]	632	CALL histdef( nid_T, "soicecov", "Ice fraction" , "[0,1]" , & ! fr_i
[3]	633	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[1649]	634	CALL histdef( nid_T, "sowindsp", "wind speed at 10m" , "m/s" , & ! wndm
	635	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[3609]	636	!
	637	IF( ln_icebergs ) THEN
	638	CALL histdef( nid_T, "calving" , "calving mass input" , "kg/s" , &
	639	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	640	CALL histdef( nid_T, "calving_heat" , "calving heat flux" , "XXXX" , &
	641	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	642	CALL histdef( nid_T, "berg_floating_melt" , "Melt rate of icebergs + bits" , "kg/m2/s", &
	643	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	644	CALL histdef( nid_T, "berg_stored_ice" , "Accumulated ice mass by class" , "kg" , &
	645	& jpi, jpj, nh_T, nclasses , 1, nclasses , nb_T , 32, clop, zsto, zout )
	646	IF( ln_bergdia ) THEN
	647	CALL histdef( nid_T, "berg_melt" , "Melt rate of icebergs" , "kg/m2/s", &
	648	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	649	CALL histdef( nid_T, "berg_buoy_melt" , "Buoyancy component of iceberg melt rate" , "kg/m2/s", &
	650	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	651	CALL histdef( nid_T, "berg_eros_melt" , "Erosion component of iceberg melt rate" , "kg/m2/s", &
	652	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	653	CALL histdef( nid_T, "berg_conv_melt" , "Convective component of iceberg melt rate", "kg/m2/s", &
	654	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	655	CALL histdef( nid_T, "berg_virtual_area" , "Virtual coverage by icebergs" , "m2" , &
	656	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	657	CALL histdef( nid_T, "bits_src" , "Mass source of bergy bits" , "kg/m2/s", &
	658	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	659	CALL histdef( nid_T, "bits_melt" , "Melt rate of bergy bits" , "kg/m2/s", &
	660	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	661	CALL histdef( nid_T, "bits_mass" , "Bergy bit density field" , "kg/m2" , &
	662	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	663	CALL histdef( nid_T, "berg_mass" , "Iceberg density field" , "kg/m2" , &
	664	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	665	CALL histdef( nid_T, "berg_real_calving" , "Calving into iceberg class" , "kg/s" , &
	666	& jpi, jpj, nh_T, nclasses , 1, nclasses , nb_T , 32, clop, zsto, zout )
	667	ENDIF
	668	ENDIF
	669
[5407]	670	IF( .NOT. ln_cpl ) THEN
[4990]	671	CALL histdef( nid_T, "sohefldp", "Surface Heat Flux: Damping" , "W/m2" , & ! qrp
	672	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	673	CALL histdef( nid_T, "sowafldp", "Surface Water Flux: Damping" , "Kg/m2/s", & ! erp
	674	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	675	CALL histdef( nid_T, "sosafldp", "Surface salt flux: damping" , "Kg/m2/s", & ! erp * sn
	676	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	677	ENDIF
[3]	678
[5407]	679	IF( ln_cpl .AND. nn_ice <= 1 ) THEN
[4990]	680	CALL histdef( nid_T, "sohefldp", "Surface Heat Flux: Damping" , "W/m2" , & ! qrp
	681	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	682	CALL histdef( nid_T, "sowafldp", "Surface Water Flux: Damping" , "Kg/m2/s", & ! erp
	683	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	684	CALL histdef( nid_T, "sosafldp", "Surface salt flux: Damping" , "Kg/m2/s", & ! erp * sn
	685	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	686	ENDIF
	687
[3]	688	clmx ="l_max(only(x))" ! max index on a period
[5836]	689	! CALL histdef( nid_T, "sobowlin", "Bowl Index" , "W-point", & ! bowl INDEX
	690	! & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clmx, zsto, zout )
[3]	691	#if defined key_diahth
	692	CALL histdef( nid_T, "sothedep", "Thermocline Depth" , "m" , & ! hth
	693	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	694	CALL histdef( nid_T, "so20chgt", "Depth of 20C isotherm" , "m" , & ! hd20
	695	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	696	CALL histdef( nid_T, "so28chgt", "Depth of 28C isotherm" , "m" , & ! hd28
	697	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[7646]	698	CALL histdef( nid_T, "sohtc300", "Heat content 300 m" , "J/m2" , & ! htc3
[3]	699	& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	700	#endif
	701
[2528]	702	CALL histend( nid_T, snc4chunks=snc4set )
[3]	703
	704	! !!! nid_U : 3D
	705	CALL histdef( nid_U, "vozocrtx", "Zonal Current" , "m/s" , & ! un
	706	& jpi, jpj, nh_U, ipk, 1, ipk, nz_U, 32, clop, zsto, zout )
[7646]	707	IF( ln_wave .AND. ln_sdw) THEN
	708	CALL histdef( nid_U, "sdzocrtx", "Stokes Drift Zonal Current" , "m/s" , & ! usd
	709	& jpi, jpj, nh_U, ipk, 1, ipk, nz_U, 32, clop, zsto, zout )
	710	ENDIF
[3]	711	! !!! nid_U : 2D
[888]	712	CALL histdef( nid_U, "sozotaux", "Wind Stress along i-axis" , "N/m2" , & ! utau
[3]	713	& jpi, jpj, nh_U, 1 , 1, 1 , - 99, 32, clop, zsto, zout )
	714
[2528]	715	CALL histend( nid_U, snc4chunks=snc4set )
[3]	716
	717	! !!! nid_V : 3D
	718	CALL histdef( nid_V, "vomecrty", "Meridional Current" , "m/s" , & ! vn
	719	& jpi, jpj, nh_V, ipk, 1, ipk, nz_V, 32, clop, zsto, zout )
[7646]	720	IF( ln_wave .AND. ln_sdw) THEN
	721	CALL histdef( nid_V, "sdmecrty", "Stokes Drift Meridional Current" , "m/s" , & ! vsd
	722	& jpi, jpj, nh_V, ipk, 1, ipk, nz_V, 32, clop, zsto, zout )
	723	ENDIF
[3]	724	! !!! nid_V : 2D
[888]	725	CALL histdef( nid_V, "sometauy", "Wind Stress along j-axis" , "N/m2" , & ! vtau
[3]	726	& jpi, jpj, nh_V, 1 , 1, 1 , - 99, 32, clop, zsto, zout )
	727
[2528]	728	CALL histend( nid_V, snc4chunks=snc4set )
[3]	729
	730	! !!! nid_W : 3D
	731	CALL histdef( nid_W, "vovecrtz", "Vertical Velocity" , "m/s" , & ! wn
	732	& jpi, jpj, nh_W, ipk, 1, ipk, nz_W, 32, clop, zsto, zout )
	733	CALL histdef( nid_W, "votkeavt", "Vertical Eddy Diffusivity" , "m2/s" , & ! avt
	734	& jpi, jpj, nh_W, ipk, 1, ipk, nz_W, 32, clop, zsto, zout )
[9019]	735	CALL histdef( nid_W, "votkeavm", "Vertical Eddy Viscosity" , "m2/s" , & ! avm
[255]	736	& jpi, jpj, nh_W, ipk, 1, ipk, nz_W, 32, clop, zsto, zout )
	737
[9019]	738	IF( ln_zdfddm ) THEN
[3]	739	CALL histdef( nid_W,"voddmavs","Salt Vertical Eddy Diffusivity" , "m2/s" , & ! avs
	740	& jpi, jpj, nh_W, ipk, 1, ipk, nz_W, 32, clop, zsto, zout )
	741	ENDIF
[7646]	742
	743	IF( ln_wave .AND. ln_sdw) THEN
	744	CALL histdef( nid_W, "sdvecrtz", "Stokes Drift Vertical Current" , "m/s" , & ! wsd
	745	& jpi, jpj, nh_W, ipk, 1, ipk, nz_W, 32, clop, zsto, zout )
	746	ENDIF
[3]	747	! !!! nid_W : 2D
[2528]	748	CALL histend( nid_W, snc4chunks=snc4set )
[3]	749
	750	IF(lwp) WRITE(numout,*)
	751	IF(lwp) WRITE(numout,*) 'End of NetCDF Initialization'
	752	IF(ll_print) CALL FLUSH(numout )
	753
	754	ENDIF
	755
	756	! 2. Start writing data
	757	! ---------------------
	758
[4292]	759	! ndex(1) est utilise ssi l'avant dernier argument est different de
[3]	760	! la taille du tableau en sortie. Dans ce cas , l'avant dernier argument
	761	! donne le nombre d'elements, et ndex la liste des indices a sortir
	762
[1334]	763	IF( lwp .AND. MOD( itmod, nwrite ) == 0 ) THEN
[3]	764	WRITE(numout,*) 'dia_wri : write model outputs in NetCDF files at ', kt, 'time-step'
	765	WRITE(numout,*) '~~~~~~ '
	766	ENDIF
	767
[6140]	768	IF( .NOT.ln_linssh ) THEN
	769	CALL histwrite( nid_T, "votemper", it, tsn(:,:,:,jp_tem) * e3t_n(:,:,:) , ndim_T , ndex_T ) ! heat content
	770	CALL histwrite( nid_T, "vosaline", it, tsn(:,:,:,jp_sal) * e3t_n(:,:,:) , ndim_T , ndex_T ) ! salt content
	771	CALL histwrite( nid_T, "sosstsst", it, tsn(:,:,1,jp_tem) * e3t_n(:,:,1) , ndim_hT, ndex_hT ) ! sea surface heat content
	772	CALL histwrite( nid_T, "sosaline", it, tsn(:,:,1,jp_sal) * e3t_n(:,:,1) , ndim_hT, ndex_hT ) ! sea surface salinity content
[4292]	773	ELSE
	774	CALL histwrite( nid_T, "votemper", it, tsn(:,:,:,jp_tem) , ndim_T , ndex_T ) ! temperature
	775	CALL histwrite( nid_T, "vosaline", it, tsn(:,:,:,jp_sal) , ndim_T , ndex_T ) ! salinity
	776	CALL histwrite( nid_T, "sosstsst", it, tsn(:,:,1,jp_tem) , ndim_hT, ndex_hT ) ! sea surface temperature
	777	CALL histwrite( nid_T, "sosaline", it, tsn(:,:,1,jp_sal) , ndim_hT, ndex_hT ) ! sea surface salinity
	778	ENDIF
[6140]	779	IF( .NOT.ln_linssh ) THEN
[7753]	780	zw3d(:,:,:) = ( ( e3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2
[6140]	781	CALL histwrite( nid_T, "vovvle3t", it, e3t_n (:,:,:) , ndim_T , ndex_T ) ! level thickness
	782	CALL histwrite( nid_T, "vovvldep", it, gdept_n(:,:,:) , ndim_T , ndex_T ) ! t-point depth
[4292]	783	CALL histwrite( nid_T, "vovvldef", it, zw3d , ndim_T , ndex_T ) ! level thickness deformation
	784	ENDIF
[3]	785	CALL histwrite( nid_T, "sossheig", it, sshn , ndim_hT, ndex_hT ) ! sea surface height
[2528]	786	CALL histwrite( nid_T, "sowaflup", it, ( emp-rnf ) , ndim_hT, ndex_hT ) ! upward water flux
[4570]	787	CALL histwrite( nid_T, "sorunoff", it, rnf , ndim_hT, ndex_hT ) ! river runoffs
[3625]	788	CALL histwrite( nid_T, "sosfldow", it, sfx , ndim_hT, ndex_hT ) ! downward salt flux
	789	! (includes virtual salt flux beneath ice
	790	! in linear free surface case)
[6140]	791	IF( ln_linssh ) THEN
[7753]	792	zw2d(:,:) = emp (:,:) * tsn(:,:,1,jp_tem)
[4292]	793	CALL histwrite( nid_T, "sosst_cd", it, zw2d, ndim_hT, ndex_hT ) ! c/d term on sst
[7753]	794	zw2d(:,:) = emp (:,:) * tsn(:,:,1,jp_sal)
[4292]	795	CALL histwrite( nid_T, "sosss_cd", it, zw2d, ndim_hT, ndex_hT ) ! c/d term on sss
	796	ENDIF
[888]	797	CALL histwrite( nid_T, "sohefldo", it, qns + qsr , ndim_hT, ndex_hT ) ! total heat flux
[3]	798	CALL histwrite( nid_T, "soshfldo", it, qsr , ndim_hT, ndex_hT ) ! solar heat flux
[1585]	799	CALL histwrite( nid_T, "somixhgt", it, hmld , ndim_hT, ndex_hT ) ! turbocline depth
[3]	800	CALL histwrite( nid_T, "somxl010", it, hmlp , ndim_hT, ndex_hT ) ! mixed layer depth
[1037]	801	CALL histwrite( nid_T, "soicecov", it, fr_i , ndim_hT, ndex_hT ) ! ice fraction
[1649]	802	CALL histwrite( nid_T, "sowindsp", it, wndm , ndim_hT, ndex_hT ) ! wind speed
[3609]	803	!
	804	IF( ln_icebergs ) THEN
	805	!
	806	CALL histwrite( nid_T, "calving" , it, berg_grid%calving , ndim_hT, ndex_hT )
	807	CALL histwrite( nid_T, "calving_heat" , it, berg_grid%calving_hflx , ndim_hT, ndex_hT )
	808	CALL histwrite( nid_T, "berg_floating_melt" , it, berg_grid%floating_melt, ndim_hT, ndex_hT )
	809	!
	810	CALL histwrite( nid_T, "berg_stored_ice" , it, berg_grid%stored_ice , ndim_bT, ndex_bT )
	811	!
	812	IF( ln_bergdia ) THEN
	813	CALL histwrite( nid_T, "berg_melt" , it, berg_melt , ndim_hT, ndex_hT )
	814	CALL histwrite( nid_T, "berg_buoy_melt" , it, buoy_melt , ndim_hT, ndex_hT )
	815	CALL histwrite( nid_T, "berg_eros_melt" , it, eros_melt , ndim_hT, ndex_hT )
	816	CALL histwrite( nid_T, "berg_conv_melt" , it, conv_melt , ndim_hT, ndex_hT )
	817	CALL histwrite( nid_T, "berg_virtual_area" , it, virtual_area , ndim_hT, ndex_hT )
	818	CALL histwrite( nid_T, "bits_src" , it, bits_src , ndim_hT, ndex_hT )
	819	CALL histwrite( nid_T, "bits_melt" , it, bits_melt , ndim_hT, ndex_hT )
	820	CALL histwrite( nid_T, "bits_mass" , it, bits_mass , ndim_hT, ndex_hT )
	821	CALL histwrite( nid_T, "berg_mass" , it, berg_mass , ndim_hT, ndex_hT )
	822	!
	823	CALL histwrite( nid_T, "berg_real_calving" , it, real_calving , ndim_bT, ndex_bT )
	824	ENDIF
	825	ENDIF
	826
[5407]	827	IF( .NOT. ln_cpl ) THEN
[4990]	828	CALL histwrite( nid_T, "sohefldp", it, qrp , ndim_hT, ndex_hT ) ! heat flux damping
	829	CALL histwrite( nid_T, "sowafldp", it, erp , ndim_hT, ndex_hT ) ! freshwater flux damping
[7753]	830	IF( ln_ssr ) zw2d(:,:) = erp(:,:) * tsn(:,:,1,jp_sal) * tmask(:,:,1)
[4990]	831	CALL histwrite( nid_T, "sosafldp", it, zw2d , ndim_hT, ndex_hT ) ! salt flux damping
	832	ENDIF
[5407]	833	IF( ln_cpl .AND. nn_ice <= 1 ) THEN
[4990]	834	CALL histwrite( nid_T, "sohefldp", it, qrp , ndim_hT, ndex_hT ) ! heat flux damping
	835	CALL histwrite( nid_T, "sowafldp", it, erp , ndim_hT, ndex_hT ) ! freshwater flux damping
[7753]	836	IF( ln_ssr ) zw2d(:,:) = erp(:,:) * tsn(:,:,1,jp_sal) * tmask(:,:,1)
[4990]	837	CALL histwrite( nid_T, "sosafldp", it, zw2d , ndim_hT, ndex_hT ) ! salt flux damping
	838	ENDIF
	839	! zw2d(:,:) = FLOAT( nmln(:,:) ) * tmask(:,:,1)
	840	! CALL histwrite( nid_T, "sobowlin", it, zw2d , ndim_hT, ndex_hT ) ! ???
[3]	841
	842	#if defined key_diahth
	843	CALL histwrite( nid_T, "sothedep", it, hth , ndim_hT, ndex_hT ) ! depth of the thermocline
	844	CALL histwrite( nid_T, "so20chgt", it, hd20 , ndim_hT, ndex_hT ) ! depth of the 20 isotherm
	845	CALL histwrite( nid_T, "so28chgt", it, hd28 , ndim_hT, ndex_hT ) ! depth of the 28 isotherm
	846	CALL histwrite( nid_T, "sohtc300", it, htc3 , ndim_hT, ndex_hT ) ! first 300m heaat content
	847	#endif
[888]	848
[3]	849	CALL histwrite( nid_U, "vozocrtx", it, un , ndim_U , ndex_U ) ! i-current
[888]	850	CALL histwrite( nid_U, "sozotaux", it, utau , ndim_hU, ndex_hU ) ! i-wind stress
[3]	851
	852	CALL histwrite( nid_V, "vomecrty", it, vn , ndim_V , ndex_V ) ! j-current
[888]	853	CALL histwrite( nid_V, "sometauy", it, vtau , ndim_hV, ndex_hV ) ! j-wind stress
[3]	854
	855	CALL histwrite( nid_W, "vovecrtz", it, wn , ndim_T, ndex_T ) ! vert. current
	856	CALL histwrite( nid_W, "votkeavt", it, avt , ndim_T, ndex_T ) ! T vert. eddy diff. coef.
[9019]	857	CALL histwrite( nid_W, "votkeavm", it, avm , ndim_T, ndex_T ) ! T vert. eddy visc. coef.
	858	IF( ln_zdfddm ) THEN
	859	CALL histwrite( nid_W, "voddmavs", it, avs , ndim_T, ndex_T ) ! S vert. eddy diff. coef.
[3]	860	ENDIF
	861
[7646]	862	IF( ln_wave .AND. ln_sdw ) THEN
[9019]	863	CALL histwrite( nid_U, "sdzocrtx", it, usd , ndim_U , ndex_U ) ! i-StokesDrift-current
	864	CALL histwrite( nid_V, "sdmecrty", it, vsd , ndim_V , ndex_V ) ! j-StokesDrift-current
	865	CALL histwrite( nid_W, "sdvecrtz", it, wsd , ndim_T , ndex_T ) ! StokesDrift vert. current
[7646]	866	ENDIF
	867
[1318]	868	! 3. Close all files
	869	! ---------------------------------------
[1561]	870	IF( kt == nitend ) THEN
[3]	871	CALL histclo( nid_T )
	872	CALL histclo( nid_U )
	873	CALL histclo( nid_V )
	874	CALL histclo( nid_W )
	875	ENDIF
[2528]	876	!
[9124]	877	IF( ln_timing ) CALL timing_stop('dia_wri')
[3294]	878	!
[3]	879	END SUBROUTINE dia_wri
[1567]	880	#endif
	881
[1334]	882	SUBROUTINE dia_wri_state( cdfile_name, kt )
[3]	883	!!---------------------------------------------------------------------
	884	!! * ROUTINE dia_wri_state *
	885	!!
	886	!! ** Purpose : create a NetCDF file named cdfile_name which contains
	887	!! the instantaneous ocean state and forcing fields.
	888	!! Used to find errors in the initial state or save the last
	889	!! ocean state in case of abnormal end of a simulation
	890	!!
	891	!! ** Method : NetCDF files using ioipsl
	892	!! File 'output.init.nc' is created if ninist = 1 (namelist)
	893	!! File 'output.abort.nc' is created in case of abnormal job end
	894	!!----------------------------------------------------------------------
[1334]	895	CHARACTER (len=* ), INTENT( in ) :: cdfile_name ! name of the file created
	896	INTEGER , INTENT( in ) :: kt ! ocean time-step index
[2528]	897	!!
[648]	898	CHARACTER (len=32) :: clname
[3]	899	CHARACTER (len=40) :: clop
[2528]	900	INTEGER :: id_i , nz_i, nh_i
	901	INTEGER, DIMENSION(1) :: idex ! local workspace
[6140]	902	REAL(wp) :: zsto, zout, zmax, zjulian
[3]	903	!!----------------------------------------------------------------------
[3294]	904	!
[3]	905	! 0. Initialisation
	906	! -----------------
[632]	907
[648]	908	! Define name, frequency of output and means
	909	clname = cdfile_name
[1792]	910	IF( .NOT. Agrif_Root() ) clname = TRIM(Agrif_CFixed())//'_'//TRIM(clname)
[3]	911	zsto = rdt
	912	clop = "inst(x)" ! no use of the mask value (require less cpu time)
	913	zout = rdt
[6140]	914	zmax = ( nitend - nit000 + 1 ) * rdt
[3]	915
[648]	916	IF(lwp) WRITE(numout,*)
	917	IF(lwp) WRITE(numout,*) 'dia_wri_state : single instantaneous ocean state'
	918	IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~ and forcing fields file created '
	919	IF(lwp) WRITE(numout,*) ' and named :', clname, '.nc'
	920
	921
[3]	922	! 1. Define NETCDF files and fields at beginning of first time step
	923	! -----------------------------------------------------------------
	924
	925	! Compute julian date from starting date of the run
[1310]	926	CALL ymds2ju( nyear, nmonth, nday, rdt, zjulian ) ! time axis
	927	zjulian = zjulian - adatrj ! set calendar origin to the beginning of the experiment
[648]	928	CALL histbeg( clname, jpi, glamt, jpj, gphit, &
[6140]	929	1, jpi, 1, jpj, nit000-1, zjulian, rdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit
[3]	930	CALL histvert( id_i, "deptht", "Vertical T levels", & ! Vertical grid : gdept
[4292]	931	"m", jpk, gdept_1d, nz_i, "down")
[3]	932
	933	! Declare all the output fields as NetCDF variables
	934
	935	CALL histdef( id_i, "vosaline", "Salinity" , "PSU" , & ! salinity
	936	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
	937	CALL histdef( id_i, "votemper", "Temperature" , "C" , & ! temperature
	938	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
[359]	939	CALL histdef( id_i, "sossheig", "Sea Surface Height" , "m" , & ! ssh
	940	& jpi, jpj, nh_i, 1 , 1, 1 , nz_i, 32, clop, zsto, zout )
[3]	941	CALL histdef( id_i, "vozocrtx", "Zonal Current" , "m/s" , & ! zonal current
	942	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
	943	CALL histdef( id_i, "vomecrty", "Meridional Current" , "m/s" , & ! meridonal current
	944	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
	945	CALL histdef( id_i, "vovecrtz", "Vertical Velocity" , "m/s" , & ! vertical current
	946	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
[6140]	947	!
[7646]	948	IF( ALLOCATED(ahtu) ) THEN
	949	CALL histdef( id_i, "ahtu" , "u-eddy diffusivity" , "m2/s" , & ! zonal current
	950	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
	951	CALL histdef( id_i, "ahtv" , "v-eddy diffusivity" , "m2/s" , & ! meridonal current
	952	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
	953	ENDIF
	954	IF( ALLOCATED(ahmt) ) THEN
	955	CALL histdef( id_i, "ahmt" , "t-eddy viscosity" , "m2/s" , & ! zonal current
	956	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
	957	CALL histdef( id_i, "ahmf" , "f-eddy viscosity" , "m2/s" , & ! meridonal current
	958	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
	959	ENDIF
[6140]	960	!
[3]	961	CALL histdef( id_i, "sowaflup", "Net Upward Water Flux" , "Kg/m2/S", & ! net freshwater
	962	& jpi, jpj, nh_i, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	963	CALL histdef( id_i, "sohefldo", "Net Downward Heat Flux", "W/m2" , & ! net heat flux
	964	& jpi, jpj, nh_i, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	965	CALL histdef( id_i, "soshfldo", "Shortwave Radiation" , "W/m2" , & ! solar flux
	966	& jpi, jpj, nh_i, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[1037]	967	CALL histdef( id_i, "soicecov", "Ice fraction" , "[0,1]" , & ! fr_i
[3]	968	& jpi, jpj, nh_i, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	969	CALL histdef( id_i, "sozotaux", "Zonal Wind Stress" , "N/m2" , & ! i-wind stress
	970	& jpi, jpj, nh_i, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
	971	CALL histdef( id_i, "sometauy", "Meridional Wind Stress", "N/m2" , & ! j-wind stress
	972	& jpi, jpj, nh_i, 1 , 1, 1 , -99 , 32, clop, zsto, zout )
[6140]	973	IF( .NOT.ln_linssh ) THEN
	974	CALL histdef( id_i, "vovvldep", "T point depth" , "m" , & ! t-point depth
[4292]	975	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
[6140]	976	CALL histdef( id_i, "vovvle3t", "T point thickness" , "m" , & ! t-point depth
	977	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
[5836]	978	ENDIF
[7646]	979	!
	980	IF( ln_wave .AND. ln_sdw ) THEN
	981	CALL histdef( id_i, "sdzocrtx", "Stokes Drift Zonal" , "m/s" , & ! StokesDrift zonal current
	982	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
	983	CALL histdef( id_i, "sdmecrty", "Stokes Drift Merid" , "m/s" , & ! StokesDrift meridonal current
	984	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
	985	CALL histdef( id_i, "sdvecrtz", "Stokes Drift Vert" , "m/s" , & ! StokesDrift vertical current
	986	& jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )
	987	ENDIF
[3]	988
[9570]	989	#if defined key_si3
[9019]	990	IF( nn_ice == 2 ) THEN ! condition needed in case agrif + ice-model but no-ice in child grid
	991	CALL ice_wri_state( kt, id_i, nh_i )
	992	ENDIF
[1482]	993	#else
[2528]	994	CALL histend( id_i, snc4chunks=snc4set )
[1482]	995	#endif
[3]	996
	997	! 2. Start writing data
	998	! ---------------------
	999	! idex(1) est utilise ssi l'avant dernier argument est diffferent de
	1000	! la taille du tableau en sortie. Dans ce cas , l'avant dernier argument
	1001	! donne le nombre d'elements, et idex la liste des indices a sortir
	1002	idex(1) = 1 ! init to avoid compil warning
[632]	1003
[3]	1004	! Write all fields on T grid
[3294]	1005	CALL histwrite( id_i, "votemper", kt, tsn(:,:,:,jp_tem), jpijpjjpk, idex ) ! now temperature
	1006	CALL histwrite( id_i, "vosaline", kt, tsn(:,:,:,jp_sal), jpijpjjpk, idex ) ! now salinity
	1007	CALL histwrite( id_i, "sossheig", kt, sshn , jpi*jpj , idex ) ! sea surface height
	1008	CALL histwrite( id_i, "vozocrtx", kt, un , jpijpjjpk, idex ) ! now i-velocity
	1009	CALL histwrite( id_i, "vomecrty", kt, vn , jpijpjjpk, idex ) ! now j-velocity
	1010	CALL histwrite( id_i, "vovecrtz", kt, wn , jpijpjjpk, idex ) ! now k-velocity
[6140]	1011	!
[7646]	1012	IF( ALLOCATED(ahtu) ) THEN
	1013	CALL histwrite( id_i, "ahtu" , kt, ahtu , jpijpjjpk, idex ) ! aht at u-point
	1014	CALL histwrite( id_i, "ahtv" , kt, ahtv , jpijpjjpk, idex ) ! - at v-point
	1015	ENDIF
	1016	IF( ALLOCATED(ahmt) ) THEN
	1017	CALL histwrite( id_i, "ahmt" , kt, ahmt , jpijpjjpk, idex ) ! ahm at t-point
	1018	CALL histwrite( id_i, "ahmf" , kt, ahmf , jpijpjjpk, idex ) ! - at f-point
	1019	ENDIF
[6140]	1020	!
[7646]	1021	CALL histwrite( id_i, "sowaflup", kt, emp - rnf , jpi*jpj , idex ) ! freshwater budget
[3294]	1022	CALL histwrite( id_i, "sohefldo", kt, qsr + qns , jpi*jpj , idex ) ! total heat flux
	1023	CALL histwrite( id_i, "soshfldo", kt, qsr , jpi*jpj , idex ) ! solar heat flux
	1024	CALL histwrite( id_i, "soicecov", kt, fr_i , jpi*jpj , idex ) ! ice fraction
	1025	CALL histwrite( id_i, "sozotaux", kt, utau , jpi*jpj , idex ) ! i-wind stress
	1026	CALL histwrite( id_i, "sometauy", kt, vtau , jpi*jpj , idex ) ! j-wind stress
[3]	1027
[6140]	1028	IF( .NOT.ln_linssh ) THEN
	1029	CALL histwrite( id_i, "vovvldep", kt, gdept_n(:,:,:), jpijpjjpk, idex )! T-cell depth
	1030	CALL histwrite( id_i, "vovvle3t", kt, e3t_n (:,:,:) , jpijpjjpk, idex )! T-cell thickness
	1031	END IF
[7646]	1032
	1033	IF( ln_wave .AND. ln_sdw ) THEN
	1034	CALL histwrite( id_i, "sdzocrtx", kt, usd , jpijpjjpk, idex) ! now StokesDrift i-velocity
	1035	CALL histwrite( id_i, "sdmecrty", kt, vsd , jpijpjjpk, idex) ! now StokesDrift j-velocity
	1036	CALL histwrite( id_i, "sdvecrtz", kt, wsd , jpijpjjpk, idex) ! now StokesDrift k-velocity
	1037	ENDIF
	1038
[3]	1039	! 3. Close the file
	1040	! -----------------
	1041	CALL histclo( id_i )
[6140]	1042	#if ! defined key_iomput
[1561]	1043	IF( ninist /= 1 ) THEN
	1044	CALL histclo( nid_T )
	1045	CALL histclo( nid_U )
	1046	CALL histclo( nid_V )
	1047	CALL histclo( nid_W )
	1048	ENDIF
	1049	#endif
[3294]	1050	!
[3]	1051	END SUBROUTINE dia_wri_state
[6140]	1052
[3]	1053	!!======================================================================
	1054	END MODULE diawri

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

Download in other formats: