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.

limwri.F90 in branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3 – NEMO

Context Navigation

source: branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90 @ 8291

Last change on this file since 8291 was 8291, checked in by clem, 7 years ago
correct last commit changes. This is the start version for the modifications of the ice model. All sette tests passed except restartability for agrif and sas (though sas_biper is restartable). Configurations CREG, SPITZ, FER & SASBIPER are running.
Property svn:keywords set to `Id`
File size: 37.0 KB

Rev	Line
[825]	1	MODULE limwri
	2	!!======================================================================
	3	!! * MODULE limwri *
	4	!! Ice diagnostics : write ice output files
	5	!!======================================================================
[888]	6	#if defined key_lim3
[825]	7	!!----------------------------------------------------------------------
[2528]	8	!! 'key_lim3' LIM3 sea-ice model
[825]	9	!!----------------------------------------------------------------------
	10	!! lim_wri : write of the diagnostics variables in ouput file
[4161]	11	!! lim_wri_state : write for initial state or/and abandon
[825]	12	!!----------------------------------------------------------------------
	13	USE ioipsl
[888]	14	USE dianam ! build name of file (routine)
[825]	15	USE phycst
	16	USE dom_oce
[888]	17	USE sbc_oce ! Surface boundary condition: ocean fields
	18	USE sbc_ice ! Surface boundary condition: ice fields
[825]	19	USE ice
[2715]	20	USE limvar
	21	USE in_out_manager
[825]	22	USE lbclnk
[2715]	23	USE lib_mpp ! MPP library
[3294]	24	USE wrk_nemo ! work arrays
[4161]	25	USE iom
	26	USE timing ! Timing
	27	USE lib_fortran ! Fortran utilities
[825]	28
	29	IMPLICIT NONE
	30	PRIVATE
	31
	32	PUBLIC lim_wri ! routine called by lim_step.F90
[4161]	33	PUBLIC lim_wri_state ! called by dia_wri_state
[825]	34
[2528]	35	!!----------------------------------------------------------------------
[4161]	36	!! NEMO/LIM3 4.0 , UCL - NEMO Consortium (2011)
[2528]	37	!! $Id$
[2715]	38	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
[2528]	39	!!----------------------------------------------------------------------
	40	CONTAINS
[825]	41
	42	SUBROUTINE lim_wri( kindic )
	43	!!-------------------------------------------------------------------
	44	!! This routine computes the average of some variables and write it
	45	!! on the ouput files.
	46	!! ATTENTION cette routine n'est valable que si le pas de temps est
	47	!! egale a une fraction entiere de 1 jours.
	48	!! Diff 1-D 3-D : suppress common also included in etat
	49	!! suppress cmoymo 11-18
	50	!! modif : 03/06/98
	51	!!-------------------------------------------------------------------
[2715]	52	INTEGER, INTENT(in) :: kindic ! if kindic < 0 there has been an error somewhere
	53	!
[4990]	54	INTEGER :: ji, jj, jk, jl ! dummy loop indices
[8266]	55	REAL(wp) :: z2da, z2db, ztmp, zrho1, zrho2, zmiss_val
[8291]	56	REAL(wp) :: zs12, zshear
[8266]	57	REAL(wp), POINTER, DIMENSION(:,:,:) :: zswi2, zmiss2
	58	REAL(wp), POINTER, DIMENSION(:,:) :: z2d, zswi, zmiss ! 2D workspace
	59	REAL(wp), POINTER, DIMENSION(:,:) :: zfb ! ice freeboard
[8239]	60	REAL(wp), POINTER, DIMENSION(:,:) :: zamask, zamask15 ! 15% concentration mask
[8291]	61	REAL(wp), POINTER, DIMENSION(:,:) :: zsig1, zsig2, zsig3
[8239]	62
	63	! Global ice diagnostics (SIMIP)
	64	REAL(wp) :: zdiag_area_nh, & ! area, extent, volume
	65	& zdiag_extt_nh, &
	66	& zdiag_area_sh, &
	67	& zdiag_extt_sh, &
	68	& zdiag_volu_nh, &
	69	& zdiag_volu_sh
	70
[2715]	71	!!-------------------------------------------------------------------
[825]	72
[4161]	73	IF( nn_timing == 1 ) CALL timing_start('limwri')
	74
[8291]	75	CALL wrk_alloc( jpi,jpj, jpl, zswi2, zmiss2 )
	76	CALL wrk_alloc( jpi,jpj , z2d, zswi, zmiss )
	77	CALL wrk_alloc( jpi,jpj , zfb, zamask, zamask15 )
	78	CALL wrk_alloc( jpi,jpj , zsig1, zsig2, zsig3 )
[825]	79
[8266]	80	!----------------------------------------
	81	! Brine volume, switches, missing values
	82	!----------------------------------------
	83
[7646]	84	! brine volume
	85	CALL lim_var_bv
[921]	86
[7646]	87	! tresholds for outputs
	88	DO jj = 1, jpj
[4688]	89	DO ji = 1, jpi
[8239]	90	zswi(ji,jj) = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - epsi06 ) ) ! 1 if ice, 0 if no ice
	91	zamask(ji,jj) = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - 0.05 ) ) ! 1 if 5% ice, 0 if less - required to mask thickness and snow depth
	92	zamask15(ji,jj) = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - 0.15 ) ) ! 1 if 15% ice, 0 if less
[4688]	93	END DO
	94	END DO
[7646]	95	DO jl = 1, jpl
	96	DO jj = 1, jpj
[4688]	97	DO ji = 1, jpi
[7646]	98	zswi2(ji,jj,jl) = MAX( 0._wp , SIGN( 1._wp , a_i(ji,jj,jl) - epsi06 ) )
[4688]	99	END DO
	100	END DO
[7646]	101	END DO
[8266]	102
[8289]	103	zmiss_val = 1.0e20
[8266]	104	zmiss(:,:) = zmiss_val * ( 1. - zswi(:,:) )
	105	zmiss2(:,:,:) = zmiss_val * ( 1. - zswi2(:,:,:) )
	106
	107	!----------------------------------------
	108	! Standard outputs
	109	!----------------------------------------
	110	! fluxes
[7646]	111	! pfrld is the lead fraction at the previous time step (actually between TRP and THD)
	112	IF( iom_use('qsr_oce') ) CALL iom_put( "qsr_oce" , qsr_oce(:,:) * pfrld(:,:) ) ! solar flux at ocean surface
	113	IF( iom_use('qns_oce') ) CALL iom_put( "qns_oce" , qns_oce(:,:) * pfrld(:,:) + qemp_oce(:,:) ) ! non-solar flux at ocean surface
	114	IF( iom_use('qsr_ice') ) CALL iom_put( "qsr_ice" , SUM( qsr_ice(:,:,:) * a_i_b(:,:,:), dim=3 ) ) ! solar flux at ice surface
	115	IF( iom_use('qns_ice') ) CALL iom_put( "qns_ice" , SUM( qns_ice(:,:,:) * a_i_b(:,:,:), dim=3 ) + qemp_ice(:,:) ) ! non-solar flux at ice surface
	116	IF( iom_use('qtr_ice') ) CALL iom_put( "qtr_ice" , SUM( ftr_ice(:,:,:) * a_i_b(:,:,:), dim=3 ) ) ! solar flux transmitted thru ice
	117	IF( iom_use('qt_oce' ) ) CALL iom_put( "qt_oce" , ( qsr_oce(:,:) + qns_oce(:,:) ) * pfrld(:,:) + qemp_oce(:,:) )
	118	IF( iom_use('qt_ice' ) ) CALL iom_put( "qt_ice" , SUM( ( qns_ice(:,:,:) + qsr_ice(:,:,:) ) &
	119	& * a_i_b(:,:,:), dim=3 ) + qemp_ice(:,:) )
	120	IF( iom_use('qemp_oce') ) CALL iom_put( "qemp_oce" , qemp_oce(:,:) )
	121	IF( iom_use('qemp_ice') ) CALL iom_put( "qemp_ice" , qemp_ice(:,:) )
	122	IF( iom_use('emp_oce' ) ) CALL iom_put( "emp_oce" , emp_oce(:,:) ) ! emp over ocean (taking into account the snow blown away from the ice)
	123	IF( iom_use('emp_ice' ) ) CALL iom_put( "emp_ice" , emp_ice(:,:) ) ! emp over ice (taking into account the snow blown away from the ice)
[921]	124
[7646]	125	! velocity
[8266]	126	IF ( iom_use('uice_ipa') ) CALL iom_put( "uice_ipa" , u_ice ) ! ice velocity u component
	127	IF ( iom_use('vice_ipa') ) CALL iom_put( "vice_ipa" , v_ice ) ! ice velocity v component
	128
	129	IF ( ( iom_use( "icevel" ) ) .OR. ( iom_use( "icevel_mv" ) ) ) THEN
[4688]	130	DO jj = 2 , jpjm1
	131	DO ji = 2 , jpim1
[7646]	132	z2da = ( u_ice(ji,jj) * umask(ji,jj,1) + u_ice(ji-1,jj) * umask(ji-1,jj,1) ) * 0.5_wp
	133	z2db = ( v_ice(ji,jj) * vmask(ji,jj,1) + v_ice(ji,jj-1) * vmask(ji,jj-1,1) ) * 0.5_wp
	134	z2d(ji,jj) = SQRT( z2da * z2da + z2db * z2db )
[4688]	135	END DO
	136	END DO
[7646]	137	CALL lbc_lnk( z2d, 'T', 1. )
[8266]	138	IF ( iom_use( "icevel" ) ) CALL iom_put( "icevel" , z2d ) ! ice velocity module
	139	IF ( iom_use( "icevel_mv" ) ) CALL iom_put( "icevel_mv" , z2d(:,:) * zswi(:,:) + zmiss(:,:) ) ! ice velocity module (missing value)
[4688]	140	ENDIF
[825]	141
[7646]	142	IF ( iom_use( "tau_icebfr" ) ) CALL iom_put( "tau_icebfr" , tau_icebfr ) ! ice friction with ocean bottom (landfast ice)
[4688]	143	!
[8239]	144	IF ( iom_use( "miceage" ) ) CALL iom_put( "miceage" , om_i * zswi * zamask15 ) ! mean ice age
[7646]	145	IF ( iom_use( "micet" ) ) CALL iom_put( "micet" , ( tm_i - rt0 ) * zswi ) ! ice mean temperature
	146	IF ( iom_use( "icest" ) ) CALL iom_put( "icest" , ( tm_su - rt0 ) * zswi ) ! ice surface temperature
	147	IF ( iom_use( "icecolf" ) ) CALL iom_put( "icecolf" , hicol ) ! frazil ice collection thickness
	148	!
[4688]	149	CALL iom_put( "isst" , sst_m ) ! sea surface temperature
	150	CALL iom_put( "isss" , sss_m ) ! sea surface salinity
[7646]	151	CALL iom_put( "iceconc" , at_i * zswi ) ! ice concentration
	152	CALL iom_put( "icevolu" , vt_i * zswi ) ! ice volume = mean ice thickness over the cell
	153	CALL iom_put( "icehc" , et_i * zswi ) ! ice total heat content
	154	CALL iom_put( "isnowhc" , et_s * zswi ) ! snow total heat content
	155	CALL iom_put( "ibrinv" , bvm_i * zswi * 100. ) ! brine volume
[5407]	156	CALL iom_put( "snowpre" , sprecip * 86400. ) ! snow precipitation
[7646]	157	CALL iom_put( "micesalt" , smt_i * zswi ) ! mean ice salinity
[8272]	158	CALL iom_put( "snowvol" , vt_s * zswi ) ! snow volume
[4161]	159
[8272]	160	CALL iom_put( "idive" , divu_i(:,:) * zswi(:,:) + zmiss(:,:) ) ! divergence
	161	CALL iom_put( "ishear" , shear_i(:,:) * zswi(:,:) + zmiss(:,:) ) ! shear
[4688]	162
	163	CALL iom_put( "icetrp" , diag_trp_vi * rday ) ! ice volume transport
	164	CALL iom_put( "snwtrp" , diag_trp_vs * rday ) ! snw volume transport
[5123]	165	CALL iom_put( "saltrp" , diag_trp_smv * rday * rhoic ) ! salt content transport
[4688]	166	CALL iom_put( "deitrp" , diag_trp_ei ) ! advected ice enthalpy (W/m2)
	167	CALL iom_put( "destrp" , diag_trp_es ) ! advected snw enthalpy (W/m2)
[921]	168
[6418]	169	CALL iom_put( "sfxbog" , sfx_bog * rday ) ! salt flux from bottom growth
[7646]	170	CALL iom_put( "sfxbom" , sfx_bom * rday ) ! salt flux from bottom melting
	171	CALL iom_put( "sfxsum" , sfx_sum * rday ) ! salt flux from surface melting
	172	CALL iom_put( "sfxlam" , sfx_lam * rday ) ! salt flux from lateral melting
[6418]	173	CALL iom_put( "sfxsni" , sfx_sni * rday ) ! salt flux from snow ice formation
	174	CALL iom_put( "sfxopw" , sfx_opw * rday ) ! salt flux from open water formation
[4688]	175	CALL iom_put( "sfxdyn" , sfx_dyn * rday ) ! salt flux from ridging rafting
[7646]	176	CALL iom_put( "sfxres" , sfx_res * rday ) ! salt flux from limupdate (resultant)
[4688]	177	CALL iom_put( "sfxbri" , sfx_bri * rday ) ! salt flux from brines
[6418]	178	CALL iom_put( "sfxsub" , sfx_sub * rday ) ! salt flux from sublimation
[4688]	179	CALL iom_put( "sfx" , sfx * rday ) ! total salt flux
[825]	180
[4990]	181	ztmp = rday / rhoic
[5123]	182	CALL iom_put( "vfxres" , wfx_res * ztmp ) ! daily prod./melting due to limupdate
	183	CALL iom_put( "vfxopw" , wfx_opw * ztmp ) ! daily lateral thermodynamic ice production
	184	CALL iom_put( "vfxsni" , wfx_sni * ztmp ) ! daily snowice ice production
	185	CALL iom_put( "vfxbog" , wfx_bog * ztmp ) ! daily bottom thermodynamic ice production
	186	CALL iom_put( "vfxdyn" , wfx_dyn * ztmp ) ! daily dynamic ice production (rid/raft)
	187	CALL iom_put( "vfxsum" , wfx_sum * ztmp ) ! surface melt
	188	CALL iom_put( "vfxbom" , wfx_bom * ztmp ) ! bottom melt
[7646]	189	CALL iom_put( "vfxlam" , wfx_lam * ztmp ) ! lateral melt
[5123]	190	CALL iom_put( "vfxice" , wfx_ice * ztmp ) ! total ice growth/melt
[7646]	191
[8233]	192	IF ( ln_pnd ) &
	193	CALL iom_put( "vfxpnd" , wfx_pnd * ztmp ) ! melt pond water flux
	194
[7646]	195	IF ( iom_use( "vfxthin" ) ) THEN ! ice production for open water + thin ice (<20cm) => comparable to observations
[7753]	196	WHERE( htm_i(:,:) < 0.2 .AND. htm_i(:,:) > 0. ) ; z2d = wfx_bog
	197	ELSEWHERE ; z2d = 0._wp
	198	END WHERE
[7646]	199	CALL iom_put( "vfxthin", ( wfx_opw + z2d ) * ztmp )
	200	ENDIF
	201
	202	ztmp = rday / rhosn
	203	CALL iom_put( "vfxspr" , wfx_spr * ztmp ) ! precip (snow)
[5123]	204	CALL iom_put( "vfxsnw" , wfx_snw * ztmp ) ! total snw growth/melt
[7646]	205	CALL iom_put( "vfxsub" , wfx_sub * ztmp ) ! sublimation (snow/ice)
[8239]	206	CALL iom_put( "vfxsub_err" , wfx_err_sub * ztmp ) ! "excess" of sublimation sent to ocean
	207
	208	CALL iom_put( "afxtot" , afx_tot ) ! concentration tendency (total)
	209	CALL iom_put( "afxdyn" , afx_dyn ) ! concentration tendency (dynamics)
	210	CALL iom_put( "afxthd" , afx_thd ) ! concentration tendency (thermo)
[825]	211
[5123]	212	CALL iom_put ('hfxthd' , hfx_thd(:,:) ) !
	213	CALL iom_put ('hfxdyn' , hfx_dyn(:,:) ) !
	214	CALL iom_put ('hfxres' , hfx_res(:,:) ) !
	215	CALL iom_put ('hfxout' , hfx_out(:,:) ) !
	216	CALL iom_put ('hfxin' , hfx_in(:,:) ) !
	217	CALL iom_put ('hfxsnw' , hfx_snw(:,:) ) !
	218	CALL iom_put ('hfxsub' , hfx_sub(:,:) ) !
	219	CALL iom_put ('hfxerr' , hfx_err(:,:) ) !
	220	CALL iom_put ('hfxerr_rem' , hfx_err_rem(:,:) ) !
[4688]	221
[5123]	222	CALL iom_put ('hfxsum' , hfx_sum(:,:) ) !
	223	CALL iom_put ('hfxbom' , hfx_bom(:,:) ) !
	224	CALL iom_put ('hfxbog' , hfx_bog(:,:) ) !
	225	CALL iom_put ('hfxdif' , hfx_dif(:,:) ) !
	226	CALL iom_put ('hfxopw' , hfx_opw(:,:) ) !
[7646]	227	CALL iom_put ('hfxtur' , fhtur(:,:) * at_i_b(:,:) ) ! turbulent heat flux at ice base
[5167]	228	CALL iom_put ('hfxdhc' , diag_heat(:,:) ) ! Heat content variation in snow and ice
[5123]	229	CALL iom_put ('hfxspr' , hfx_spr(:,:) ) ! Heat content of snow precip
[7646]	230
[8291]	231	! specific outputs for EVP rheology
	232	IF( iom_use( "isig1" ) .OR. iom_use( "isig2" ) .OR. iom_use( "isig3" ) ) THEN
	233	zsig1(:,:) = 0._wp; zsig2(:,:) = 0._wp; zsig3(:,:) = 0._wp;
	234	DO jj = 2, jpjm1
	235	DO ji = 2, jpim1
	236	zs12 = ( zswi(ji-1,jj) * stress12_i(ji-1,jj) + zswi(ji ,jj-1) * stress12_i(ji ,jj-1) + & ! stress12_i at T-point
	237	& zswi(ji ,jj) * stress12_i(ji ,jj) + zswi(ji-1,jj-1) * stress12_i(ji-1,jj-1) ) &
	238	& / MAX( 1._wp, zswi(ji-1,jj) + zswi(ji,jj-1) + zswi(ji,jj) + zswi(ji-1,jj-1) )
	239
	240	zshear = SQRT( stress2_i(ji,jj) * stress2_i(ji,jj) + 4._wp * zs12 * zs12 ) ! shear stress
	241
	242	z2da = zswi(ji,jj) / MAX( 1._wp, strength(ji,jj) )
	243
	244	!! zsig1(ji,jj) = 0.5_wp * z2da * ( stress1_i(ji,jj) + zshear ) ! principal stress (y-direction, see Hunke & Dukowicz 2002)
	245	!! zsig2(ji,jj) = 0.5_wp * z2da * ( stress1_i(ji,jj) - zshear ) ! principal stress (x-direction, see Hunke & Dukowicz 2002)
	246	!! zsig3(ji,jj) = z2da*2 ( ( stress1_i(ji,jj) + strength(ji,jj) )*2 + ( rn_ecc zshear )**2 ) ! quadratic relation linking compressive stress to shear stress
	247	!! ! (scheme converges if this value is ~1, see Bouillon et al 2009 (eq. 11))
	248	zsig1(ji,jj) = 0.5_wp * z2da * ( stress1_i(ji,jj) ) ! compressive stress, see Bouillon et al. 2015
	249	zsig2(ji,jj) = 0.5_wp * z2da * ( zshear ) ! shear stress
	250	zsig3(ji,jj) = z2da*2 ( ( stress1_i(ji,jj) + strength(ji,jj) )*2 + ( rn_ecc zshear )**2 )
	251	END DO
	252	END DO
	253	CALL lbc_lnk_multi( zsig1, 'T', 1., zsig2, 'T', 1., zsig3, 'T', 1. )
	254	CALL iom_put( "isig1" , zsig1 )
	255	CALL iom_put( "isig2" , zsig2 )
	256	CALL iom_put( "isig3" , zsig3 )
	257	ENDIF
	258
[8233]	259	! MV MP 2016
	260	IF ( ln_pnd ) THEN
	261	CALL iom_put( "iceamp" , at_ip * zswi ) ! melt pond total fraction
	262	CALL iom_put( "icevmp" , vt_ip * zswi ) ! melt pond total volume per unit area
	263	ENDIF
	264	! END MV MP 2016
	265
[8239]	266	!----------------------------------
	267	! Output category-dependent fields
	268	!----------------------------------
[7646]	269	IF ( iom_use( "iceconc_cat" ) ) CALL iom_put( "iceconc_cat" , a_i * zswi2 ) ! area for categories
	270	IF ( iom_use( "icethic_cat" ) ) CALL iom_put( "icethic_cat" , ht_i * zswi2 ) ! thickness for categories
	271	IF ( iom_use( "snowthic_cat" ) ) CALL iom_put( "snowthic_cat" , ht_s * zswi2 ) ! snow depth for categories
	272	IF ( iom_use( "salinity_cat" ) ) CALL iom_put( "salinity_cat" , sm_i * zswi2 ) ! salinity for categories
[5517]	273	! ice temperature
[7646]	274	IF ( iom_use( "icetemp_cat" ) ) CALL iom_put( "icetemp_cat", ( SUM( t_i(:,:,:,:), dim=3 ) * r1_nlay_i - rt0 ) * zswi2 )
[5517]	275	! snow temperature
[7646]	276	IF ( iom_use( "snwtemp_cat" ) ) CALL iom_put( "snwtemp_cat", ( SUM( t_s(:,:,:,:), dim=3 ) * r1_nlay_s - rt0 ) * zswi2 )
	277	! ice age
[8239]	278	IF ( iom_use( "iceage_cat" ) ) CALL iom_put( "iceage_cat" , o_i * zswi2 )
[7646]	279	! brine volume
	280	IF ( iom_use( "brinevol_cat" ) ) CALL iom_put( "brinevol_cat", bv_i * 100. * zswi2 )
[5517]	281
[8233]	282	! MV MP 2016
	283	IF ( ln_pnd ) THEN
	284	IF ( iom_use( "iceamp_cat" ) ) CALL iom_put( "iceamp_cat" , a_ip * zswi2 ) ! melt pond frac for categories
	285	IF ( iom_use( "icevmp_cat" ) ) CALL iom_put( "icevmp_cat" , v_ip * zswi2 ) ! melt pond frac for categories
	286	IF ( iom_use( "icehmp_cat" ) ) CALL iom_put( "icehmp_cat" , h_ip * zswi2 ) ! melt pond frac for categories
	287	IF ( iom_use( "iceafp_cat" ) ) CALL iom_put( "iceafp_cat" , a_ip_frac * zswi2 ) ! melt pond frac for categories
	288	ENDIF
	289	! END MV MP 2016
	290
[8239]	291	!--------------------------------
	292	! Add-ons for SIMIP
	293	!--------------------------------
	294	zrho1 = ( rau0 - rhoic ) / rau0; zrho2 = rhosn / rau0
	295
[8272]	296	IF ( iom_use( "icepres" ) ) CALL iom_put( "icepres" , zswi(:,:) ) ! Ice presence (1 or 0)
	297	IF ( iom_use( "icemass" ) ) CALL iom_put( "icemass" , rhoic * vt_i(:,:) * zswi(:,:) ) ! Ice mass per cell area
	298	IF ( iom_use( "icethic" ) ) CALL iom_put( "icethic" , htm_i(:,:) * zamask(:,:) + ( 1. - zamask(:,:) ) * zmiss_val ) ! Ice thickness
	299	IF ( iom_use( "snomass" ) ) CALL iom_put( "snomass" , rhosn * vt_s(:,:) * zswi(:,:) + zmiss(:,:) ) ! Snow mass per cell area
	300	IF ( iom_use( "snothic" ) ) CALL iom_put( "snothic" , htm_s(:,:) * zamask(:,:) + ( 1. - zamask(:,:) ) * zmiss_val ) ! Snow thickness
[8266]	301
[8272]	302	IF ( iom_use( "iceconc_cat_mv" ) ) CALL iom_put( "iceconc_cat_mv" , a_i(:,:,:) * zswi2(:,:,:) + zmiss2(:,:,:) ) ! Area for categories
	303	IF ( iom_use( "icethic_cat_mv" ) ) CALL iom_put( "icethic_cat_mv" , ht_i(:,:,:) * zswi2(:,:,:) + zmiss2(:,:,:) ) ! Thickness for categories
	304	IF ( iom_use( "snowthic_cat_mv" ) ) CALL iom_put( "snowthic_cat_mv", ht_s(:,:,:) * zswi2(:,:,:) + zmiss2(:,:,:) ) ! Snow depth for categories
[8266]	305
[8272]	306	IF ( iom_use( "icestK" ) ) CALL iom_put( "icestK" , tm_su(:,:) * zswi(:,:) + zmiss(:,:) ) ! Ice surface temperature
	307	IF ( iom_use( "icesntK" ) ) CALL iom_put( "icesntK" , tm_si(:,:) * zswi(:,:) + zmiss(:,:) ) ! Snow-ice interface temperature
	308	IF ( iom_use( "icebotK" ) ) CALL iom_put( "icebotK" , t_bo(:,:) * zswi(:,:) + zmiss(:,:) ) ! Ice bottom temperature
	309	IF ( iom_use( "iceage" ) ) CALL iom_put( "iceage" , om_i(:,:) * zamask15(:,:) + ( 1. - zamask15(:,:) ) * zmiss_val ) ! Ice age
	310	IF ( iom_use( "icesmass" ) ) CALL iom_put( "icesmass" , SUM( smv_i, DIM = 3 ) * rhoic * 1.0e-3 * zswi(:,:) ) ! Mass of salt in sea ice per cell area
	311	IF ( iom_use( "icesal" ) ) CALL iom_put( "icesal" , smt_i(:,:) * zswi(:,:) + zmiss(:,:) ) ! Ice salinity
[8266]	312
[8239]	313	IF ( iom_use( "icefb" ) ) THEN
[8266]	314	zfb(:,:) = ( zrho1 * htm_i(:,:) - zrho2 * htm_s(:,:) )
[8239]	315	WHERE( zfb < 0._wp ) ; zfb = 0._wp ; END WHERE
[8272]	316	CALL iom_put( "icefb" , zfb(:,:) * zswi(:,:) + zmiss(:,:) ) ! Ice freeboard
[8239]	317	ENDIF
[8266]	318
[8272]	319	IF ( iom_use( "isnhcneg" ) ) CALL iom_put( "isnhcneg" , - et_s(:,:) * zswi(:,:) + zmiss(:,:) ) ! Snow total heat content
[8266]	320
	321	IF ( iom_use( "dmithd" ) ) CALL iom_put( "dmithd" , - wfx_bog - wfx_bom - wfx_sum & ! Sea-ice mass change from thermodynamics
[8239]	322	& - wfx_sni - wfx_opw - wfx_res )
[8266]	323	IF ( iom_use( "dmidyn" ) ) CALL iom_put( "dmidyn" , diag_dmi_dyn ) ! Sea-ice mass change from dynamics
	324	IF ( iom_use( "dmiopw" ) ) CALL iom_put( "dmiopw" , - wfx_opw ) ! Sea-ice mass change through growth in open water
	325	IF ( iom_use( "dmibog" ) ) CALL iom_put( "dmibog" , - wfx_bog ) ! Sea-ice mass change through basal growth
	326	IF ( iom_use( "dmisni" ) ) CALL iom_put( "dmisni" , - wfx_sni ) ! Sea-ice mass change through snow-to-ice conversion
	327	IF ( iom_use( "dmisum" ) ) CALL iom_put( "dmisum" , - wfx_sum ) ! Sea-ice mass change through surface melting
	328	IF ( iom_use( "dmibom" ) ) CALL iom_put( "dmibom" , - wfx_bom ) ! Sea-ice mass change through bottom melting
[8239]	329
[8266]	330	IF ( iom_use( "dmtsub" ) ) CALL iom_put( "dmtsub" , - wfx_sub ) ! Sea-ice mass change through evaporation and sublimation
	331	IF ( iom_use( "dmssub" ) ) CALL iom_put( "dmssub" , - wfx_snw_sub ) ! Snow mass change through sublimation
	332	IF ( iom_use( "dmisub" ) ) CALL iom_put( "dmisub" , - wfx_ice_sub ) ! Sea-ice mass change through sublimation
[8239]	333
[8266]	334	IF ( iom_use( "dmsspr" ) ) CALL iom_put( "dmsspr" , - wfx_spr ) ! Snow mass change through snow fall
	335	IF ( iom_use( "dmsssi" ) ) CALL iom_put( "dmsssi" , wfx_sni*rhosn/rhoic ) ! Snow mass change through snow-to-ice conversion
[8239]	336
[8266]	337	IF ( iom_use( "dmsmel" ) ) CALL iom_put( "dmsmel" , - wfx_snw_sum ) ! Snow mass change through melt
	338	IF ( iom_use( "dmsdyn" ) ) CALL iom_put( "dmsdyn" , diag_dms_dyn ) ! Snow mass change through dynamics
[8239]	339
[8266]	340	IF ( iom_use( "hfxsenso" ) ) CALL iom_put( "hfxsenso" , -fhtur(:,:) * zswi(:,:) + zmiss(:,:) ) ! Sensible oceanic heat flux
	341	IF ( iom_use( "hfxconbo" ) ) CALL iom_put( "hfxconbo" , diag_fc_bo * zswi(:,:) + zmiss(:,:) ) ! Bottom conduction flux
	342	IF ( iom_use( "hfxconsu" ) ) CALL iom_put( "hfxconsu" , diag_fc_su * zswi(:,:) + zmiss(:,:) ) ! Surface conduction flux
[8239]	343
[8266]	344	IF ( iom_use( "wfxtot" ) ) CALL iom_put( "wfxtot" , wfx_ice(:,:) * zswi(:,:) + zmiss(:,:) ) ! Total freshwater flux from sea ice
	345	IF ( iom_use( "wfxsum" ) ) CALL iom_put( "wfxsum" , wfx_sum(:,:) * zswi(:,:) + zmiss(:,:) ) ! Freshwater flux from sea-ice surface
	346	IF ( iom_use( "sfx_mv" ) ) CALL iom_put( "sfx_mv" , sfx(:,:) * 0.001 * zswi(:,:) + zmiss(:,:) ) ! Total salt flux
[8239]	347
[8266]	348	IF ( iom_use( "uice_mv" ) ) CALL iom_put( "uice_mv" , u_ice(:,:) * zswi(:,:) + zmiss(:,:) ) ! ice velocity u component
	349	IF ( iom_use( "vice_mv" ) ) CALL iom_put( "vice_mv" , v_ice(:,:) * zswi(:,:) + zmiss(:,:) ) ! ice velocity v component
	350
	351	IF ( iom_use( "xmtrpice" ) ) CALL iom_put( "xmtrpice" , diag_xmtrp_ice(:,:) ) ! X-component of sea-ice mass transport (kg/s)
	352	IF ( iom_use( "ymtrpice" ) ) CALL iom_put( "ymtrpice" , diag_ymtrp_ice(:,:) ) ! Y-component of sea-ice mass transport
[8239]	353
[8266]	354	IF ( iom_use( "xmtrpsnw" ) ) CALL iom_put( "xmtrpsnw" , diag_xmtrp_snw(:,:) ) ! X-component of snow mass transport (kg/s)
	355	IF ( iom_use( "ymtrpsnw" ) ) CALL iom_put( "ymtrpsnw" , diag_ymtrp_snw(:,:) ) ! Y-component of snow mass transport
[8239]	356
[8266]	357	IF ( iom_use( "xatrp" ) ) CALL iom_put( "xatrp" , diag_xatrp(:,:) ) ! X-component of ice area transport
	358	IF ( iom_use( "yatrp" ) ) CALL iom_put( "yatrp" , diag_yatrp(:,:) ) ! Y-component of ice area transport
[8239]	359
[8266]	360	IF ( iom_use( "utau_ice" ) ) CALL iom_put( "utau_ice" , utau_ice(:,:) * zswi(:,:) + zmiss(:,:) ) ! Wind stress term in force balance (x)
	361	IF ( iom_use( "vtau_ice" ) ) CALL iom_put( "vtau_ice" , vtau_ice(:,:) * zswi(:,:) + zmiss(:,:) ) ! Wind stress term in force balance (y)
[8239]	362
[8266]	363	IF ( iom_use( "utau_oi" ) ) CALL iom_put( "utau_oi" , diag_utau_oi(:,:) * zswi(:,:) + zmiss(:,:) ) ! Ocean stress term in force balance (x)
	364	IF ( iom_use( "vtau_oi" ) ) CALL iom_put( "vtau_oi" , diag_vtau_oi(:,:) * zswi(:,:) + zmiss(:,:) ) ! Ocean stress term in force balance (y)
[8239]	365
[8266]	366	IF ( iom_use( "icestr" ) ) CALL iom_put( "icestr" , strength(:,:) * zswi(:,:) + zmiss(:,:) ) ! Ice strength
[8239]	367
[8266]	368	IF ( iom_use( "dssh_dx" ) ) CALL iom_put( "dssh_dx" , diag_dssh_dx(:,:) * zswi(:,:) + zmiss(:,:) ) ! Sea-surface tilt term in force balance (x)
	369	IF ( iom_use( "dssh_dy" ) ) CALL iom_put( "dssh_dy" , diag_dssh_dy(:,:) * zswi(:,:) + zmiss(:,:) ) ! Sea-surface tilt term in force balance (y)
[8239]	370
[8266]	371	IF ( iom_use( "corstrx" ) ) CALL iom_put( "corstrx" , diag_corstrx(:,:) * zswi(:,:) + zmiss(:,:) ) ! Coriolis force term in force balance (x)
	372	IF ( iom_use( "corstry" ) ) CALL iom_put( "corstry" , diag_corstry(:,:) * zswi(:,:) + zmiss(:,:) ) ! Coriolis force term in force balance (y)
[8239]	373
[8266]	374	IF ( iom_use( "intstrx" ) ) CALL iom_put( "intstrx" , diag_intstrx(:,:) * zswi(:,:) + zmiss(:,:) ) ! Internal force term in force balance (x)
	375	IF ( iom_use( "intstry" ) ) CALL iom_put( "intstry" , diag_intstry(:,:) * zswi(:,:) + zmiss(:,:) ) ! Internal force term in force balance (y)
	376
[8272]	377	IF ( iom_use( "normstr" ) ) CALL iom_put( "normstr" , diag_sig1(:,:) * zswi(:,:) + zmiss(:,:) ) ! Normal stress
	378	IF ( iom_use( "sheastr" ) ) CALL iom_put( "sheastr" , diag_sig2(:,:) * zswi(:,:) + zmiss(:,:) ) ! Shear stress
[8266]	379
[8239]	380	!--------------------------------
	381	! Global ice diagnostics (SIMIP)
	382	!--------------------------------
	383
	384	IF ( iom_use( "NH_icearea" ) .OR. iom_use( "NH_icevolu" ) .OR. iom_use( "NH_iceextt" ) ) THEN ! NH integrated diagnostics
	385
[8240]	386	WHERE( ff_t > 0._wp ); zswi(:,:) = 1.0e-12
[8239]	387	ELSEWHERE ; zswi(:,:) = 0.
	388	END WHERE
	389
[8240]	390	zdiag_area_nh = glob_sum( at_i(:,:) * zswi(:,:) * e1e2t(:,:) )
	391	zdiag_volu_nh = glob_sum( vt_i(:,:) * zswi(:,:) * e1e2t(:,:) )
[8239]	392
[8240]	393	WHERE( ff_t > 0._wp .AND. at_i > 0.15 ); zswi(:,:) = 1.0e-12
[8239]	394	ELSEWHERE ; zswi(:,:) = 0.
	395	END WHERE
	396
[8240]	397	zdiag_extt_nh = glob_sum( zswi(:,:) * e1e2t(:,:) )
[8239]	398
	399	IF ( iom_use( "NH_icearea" ) ) CALL iom_put( "NH_icearea" , zdiag_area_nh )
	400	IF ( iom_use( "NH_icevolu" ) ) CALL iom_put( "NH_icevolu" , zdiag_volu_nh )
	401	IF ( iom_use( "NH_iceextt" ) ) CALL iom_put( "NH_iceextt" , zdiag_extt_nh )
	402
	403	ENDIF
	404
	405	IF ( iom_use( "SH_icearea" ) .OR. iom_use( "SH_icevolu" ) .OR. iom_use( "SH_iceextt" ) ) THEN ! SH integrated diagnostics
	406
[8240]	407	WHERE( ff_t < 0._wp ); zswi(:,:) = 1.0e-12;
[8239]	408	ELSEWHERE ; zswi(:,:) = 0.
	409	END WHERE
	410
[8240]	411	zdiag_area_sh = glob_sum( at_i(:,:) * zswi(:,:) * e1e2t(:,:) )
	412	zdiag_volu_sh = glob_sum( vt_i(:,:) * zswi(:,:) * e1e2t(:,:) )
[8239]	413
[8240]	414	WHERE( ff_t < 0._wp .AND. at_i > 0.15 ); zswi(:,:) = 1.0e-12
[8239]	415	ELSEWHERE ; zswi(:,:) = 0.
	416	END WHERE
	417
[8240]	418	zdiag_extt_sh = glob_sum( zswi(:,:) * e1e2t(:,:) )
[8239]	419
	420	IF ( iom_use( "SH_icearea" ) ) CALL iom_put( "SH_icearea", zdiag_area_sh )
	421	IF ( iom_use( "SH_icevolu" ) ) CALL iom_put( "SH_icevolu", zdiag_volu_sh )
	422	IF ( iom_use( "SH_iceextt" ) ) CALL iom_put( "SH_iceextt", zdiag_extt_sh )
	423
	424	ENDIF
	425
[4688]	426	! ! Create an output files (output.lim.abort.nc) if S < 0 or u > 20 m/s
	427	! IF( kindic < 0 ) CALL lim_wri_state( 'output.abort' )
	428	! not yet implemented
	429
[8266]	430	CALL wrk_dealloc( jpi, jpj, jpl, zswi2, zmiss2 )
	431	CALL wrk_dealloc( jpi, jpj , z2d, zswi, zmiss )
[8239]	432	CALL wrk_dealloc( jpi, jpj , zfb, zamask, zamask15 )
[8291]	433	CALL wrk_dealloc( jpi, jpj , zsig1, zsig2, zsig3 )
[4161]	434
	435	IF( nn_timing == 1 ) CALL timing_stop('limwri')
[2777]	436
[825]	437	END SUBROUTINE lim_wri
	438
[4161]	439
	440	SUBROUTINE lim_wri_state( kt, kid, kh_i )
	441	!!---------------------------------------------------------------------
	442	!! * ROUTINE lim_wri_state *
	443	!!
	444	!! ** Purpose : create a NetCDF file named cdfile_name which contains
	445	!! the instantaneous ice state and forcing fields for ice model
	446	!! Used to find errors in the initial state or save the last
	447	!! ocean state in case of abnormal end of a simulation
	448	!!
	449	!! History :
[7646]	450	!! 4.0 ! 2013-06 (C. Rousset)
[4161]	451	!!----------------------------------------------------------------------
[7646]	452	INTEGER, INTENT( in ) :: kt ! ocean time-step index)
	453	INTEGER, INTENT( in ) :: kid , kh_i
	454	INTEGER :: nz_i, jl
	455	REAL(wp), DIMENSION(jpl) :: jcat
[4161]	456	!!----------------------------------------------------------------------
[7646]	457	DO jl = 1, jpl
	458	jcat(jl) = REAL(jl)
	459	ENDDO
	460
	461	CALL histvert( kid, "ncatice", "Ice Categories","", jpl, jcat, nz_i, "up")
[825]	462
[7646]	463	CALL histdef( kid, "sithic", "Ice thickness" , "m" , &
[4765]	464	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	465	CALL histdef( kid, "siconc", "Ice concentration" , "%" , &
[4765]	466	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	467	CALL histdef( kid, "sitemp", "Ice temperature" , "C" , &
[4765]	468	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	469	CALL histdef( kid, "sivelu", "i-Ice speed " , "m/s" , &
[4765]	470	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	471	CALL histdef( kid, "sivelv", "j-Ice speed " , "m/s" , &
[4765]	472	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	473	CALL histdef( kid, "sistru", "i-Wind stress over ice " , "Pa" , &
[4765]	474	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	475	CALL histdef( kid, "sistrv", "j-Wind stress over ice " , "Pa" , &
[4765]	476	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	477	CALL histdef( kid, "sisflx", "Solar flux over ocean" , "w/m2" , &
[4765]	478	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	479	CALL histdef( kid, "sinflx", "Non-solar flux over ocean" , "w/m2" , &
[4765]	480	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
	481	CALL histdef( kid, "isnowpre", "Snow precipitation" , "kg/m2/s", &
	482	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	483	CALL histdef( kid, "sisali", "Ice salinity" , "PSU" , &
[4765]	484	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	485	CALL histdef( kid, "sivolu", "Ice volume" , "m" , &
[4765]	486	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	487	CALL histdef( kid, "sidive", "Ice divergence" , "10-8s-1", &
[4765]	488	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	489
[8233]	490	! MV MP 2016
	491	IF ( ln_pnd ) THEN
	492	CALL histdef( kid, "si_amp", "Melt pond fraction" , "%" , &
	493	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
	494	CALL histdef( kid, "si_vmp", "Melt pond volume" , "m" , &
	495	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
	496	ENDIF
	497	! END MV MP 2016
	498
[7646]	499	CALL histdef( kid, "vfxbog", "Ice bottom production" , "m/s" , &
[4765]	500	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	501	CALL histdef( kid, "vfxdyn", "Ice dynamic production" , "m/s" , &
[4765]	502	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	503	CALL histdef( kid, "vfxopw", "Ice open water prod" , "m/s" , &
[4765]	504	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	505	CALL histdef( kid, "vfxsni", "Snow ice production " , "m/s" , &
[4765]	506	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	507	CALL histdef( kid, "vfxres", "Ice prod from limupdate" , "m/s" , &
[4765]	508	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	509	CALL histdef( kid, "vfxbom", "Ice bottom melt" , "m/s" , &
[4765]	510	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[7646]	511	CALL histdef( kid, "vfxsum", "Ice surface melt" , "m/s" , &
[4765]	512	& jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
[4161]	513
[7646]	514	CALL histdef( kid, "sithicat", "Ice thickness" , "m" , &
	515	& jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt )
	516	CALL histdef( kid, "siconcat", "Ice concentration" , "%" , &
	517	& jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt )
	518	CALL histdef( kid, "sisalcat", "Ice salinity" , "" , &
	519	& jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt )
	520	CALL histdef( kid, "sitemcat", "Ice temperature" , "C" , &
	521	& jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt )
	522	CALL histdef( kid, "snthicat", "Snw thickness" , "m" , &
	523	& jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt )
	524	CALL histdef( kid, "sntemcat", "Snw temperature" , "C" , &
	525	& jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt )
	526
[4161]	527	CALL histend( kid, snc4set ) ! end of the file definition
	528
[7646]	529	CALL histwrite( kid, "sithic", kt, htm_i , jpi*jpj, (/1/) )
	530	CALL histwrite( kid, "siconc", kt, at_i , jpi*jpj, (/1/) )
	531	CALL histwrite( kid, "sitemp", kt, tm_i - rt0 , jpi*jpj, (/1/) )
	532	CALL histwrite( kid, "sivelu", kt, u_ice , jpi*jpj, (/1/) )
	533	CALL histwrite( kid, "sivelv", kt, v_ice , jpi*jpj, (/1/) )
	534	CALL histwrite( kid, "sistru", kt, utau_ice , jpi*jpj, (/1/) )
	535	CALL histwrite( kid, "sistrv", kt, vtau_ice , jpi*jpj, (/1/) )
	536	CALL histwrite( kid, "sisflx", kt, qsr , jpi*jpj, (/1/) )
	537	CALL histwrite( kid, "sinflx", kt, qns , jpi*jpj, (/1/) )
[4161]	538	CALL histwrite( kid, "isnowpre", kt, sprecip , jpi*jpj, (/1/) )
[7646]	539	CALL histwrite( kid, "sisali", kt, smt_i , jpi*jpj, (/1/) )
	540	CALL histwrite( kid, "sivolu", kt, vt_i , jpi*jpj, (/1/) )
	541	CALL histwrite( kid, "sidive", kt, divu_i1.0e8 , jpijpj, (/1/) )
[4161]	542
[8233]	543	! MV MP 2016
	544	IF ( ln_pnd ) THEN
	545	CALL histwrite( kid, "si_amp", kt, at_ip , jpi*jpj, (/1/) )
	546	CALL histwrite( kid, "si_vmp", kt, vt_ip , jpi*jpj, (/1/) )
	547	ENDIF
	548	! END MV MP 2016
	549
[7646]	550	CALL histwrite( kid, "vfxbog", kt, wfx_bog , jpi*jpj, (/1/) )
	551	CALL histwrite( kid, "vfxdyn", kt, wfx_dyn , jpi*jpj, (/1/) )
	552	CALL histwrite( kid, "vfxopw", kt, wfx_opw , jpi*jpj, (/1/) )
	553	CALL histwrite( kid, "vfxsni", kt, wfx_sni , jpi*jpj, (/1/) )
	554	CALL histwrite( kid, "vfxres", kt, wfx_res , jpi*jpj, (/1/) )
	555	CALL histwrite( kid, "vfxbom", kt, wfx_bom , jpi*jpj, (/1/) )
	556	CALL histwrite( kid, "vfxsum", kt, wfx_sum , jpi*jpj, (/1/) )
[8233]	557	IF ( ln_pnd ) &
	558	CALL histwrite( kid, "vfxpnd", kt, wfx_pnd , jpi*jpj, (/1/) )
[4333]	559
[7646]	560	CALL histwrite( kid, "sithicat", kt, ht_i , jpijpjjpl, (/1/) )
	561	CALL histwrite( kid, "siconcat", kt, a_i , jpijpjjpl, (/1/) )
	562	CALL histwrite( kid, "sisalcat", kt, sm_i , jpijpjjpl, (/1/) )
	563	CALL histwrite( kid, "sitemcat", kt, tm_i - rt0 , jpijpjjpl, (/1/) )
	564	CALL histwrite( kid, "snthicat", kt, ht_s , jpijpjjpl, (/1/) )
	565	CALL histwrite( kid, "sntemcat", kt, tm_su - rt0 , jpijpjjpl, (/1/) )
	566
[4688]	567	! Close the file
	568	! -----------------
	569	!CALL histclo( kid )
[4161]	570
	571	END SUBROUTINE lim_wri_state
	572
[825]	573	#else
	574	!!----------------------------------------------------------------------
	575	!! Default option : Empty module NO LIM sea-ice model
	576	!!----------------------------------------------------------------------
	577	CONTAINS
	578	SUBROUTINE lim_wri ! Empty routine
	579	END SUBROUTINE lim_wri
	580	#endif
	581
	582	!!======================================================================
	583	END MODULE limwri

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

Download in other formats: