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 trunk/NEMO/LIM_SRC_3 – NEMO

Context Navigation

source: trunk/NEMO/LIM_SRC_3/limwri.F90 @ 1530

Last change on this file since 1530 was 1469, checked in by smasson, 15 years ago
[uv]taui_ice renamed [uv]tau_ice, see ticket:452
Property svn:keywords set to `Id`
File size: 20.0 KB

Rev	Line
[825]	1	MODULE limwri
[834]	2	!!----------------------------------------------------------------------
	3	!! 'key_lim3' LIM3 sea-ice model
	4	!!----------------------------------------------------------------------
[825]	5	!!======================================================================
	6	!! * MODULE limwri *
	7	!! Ice diagnostics : write ice output files
	8	!!======================================================================
[888]	9	#if defined key_lim3
[825]	10	!!----------------------------------------------------------------------
[1156]	11	!! LIM 3.0, UCL-LOCEAN-IPSL (2008)
	12	!! $Id$
[825]	13	!! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
	14	!!----------------------------------------------------------------------
	15	!! lim_wri : write of the diagnostics variables in ouput file
	16	!! lim_wri_init : initialization and namelist read
	17	!!----------------------------------------------------------------------
	18	!! * Modules used
	19	USE ioipsl
[888]	20	USE dianam ! build name of file (routine)
[825]	21	USE phycst
	22	USE dom_oce
	23	USE daymod
	24	USE in_out_manager
[888]	25	USE sbc_oce ! Surface boundary condition: ocean fields
	26	USE sbc_ice ! Surface boundary condition: ice fields
[825]	27	USE dom_ice
	28	USE ice
	29	USE iceini
	30	USE lbclnk
	31	USE par_ice
	32	USE limvar
	33
	34	IMPLICIT NONE
	35	PRIVATE
	36
	37	!! * Accessibility
	38	PUBLIC lim_wri ! routine called by lim_step.F90
	39
	40	!! * Module variables
	41	INTEGER, PARAMETER :: & !:
	42	jpnoumax = 40 !: maximum number of variable for ice output
	43	INTEGER :: &
	44	noumef , & ! number of fields
	45	noumefa , & ! number of additional fields
	46	add_diag_swi , & ! additional diagnostics
	47	nz ! dimension for the itd field
	48
	49	REAL(wp) , DIMENSION(jpnoumax) :: &
	50	cmulti , & ! multiplicative constant
	51	cadd , & ! additive constant
	52	cmultia , & ! multiplicative constant
	53	cadda ! additive constant
	54	CHARACTER(len = 35), DIMENSION(jpnoumax) :: &
	55	titn, titna ! title of the field
	56	CHARACTER(len = 8 ), DIMENSION(jpnoumax) :: &
	57	nam, nama ! name of the field
	58	CHARACTER(len = 8 ), DIMENSION(jpnoumax) :: &
	59	uni, unia ! unit of the field
	60	INTEGER , DIMENSION(jpnoumax) :: &
	61	nc, nca ! switch for saving field ( = 1 ) or not ( = 0 )
	62
	63	REAL(wp) :: & ! constant values
	64	epsi16 = 1e-16 , &
	65	zzero = 0.e0 , &
	66	zone = 1.e0
	67
	68	CONTAINS
	69	#if defined key_dimgout
	70
	71	# include "limwri_dimg.h90"
	72
	73	#else
	74
	75	SUBROUTINE lim_wri( kindic )
	76	!!-------------------------------------------------------------------
	77	!! This routine computes the average of some variables and write it
	78	!! on the ouput files.
	79	!! ATTENTION cette routine n'est valable que si le pas de temps est
	80	!! egale a une fraction entiere de 1 jours.
	81	!! Diff 1-D 3-D : suppress common also included in etat
	82	!! suppress cmoymo 11-18
	83	!! modif : 03/06/98
	84	!!-------------------------------------------------------------------
	85	INTEGER, INTENT(in) :: &
[921]	86	kindic ! if kindic < 0 there has been an error somewhere
[825]	87
	88	!! * Local variables
	89	REAL(wp),DIMENSION(1) :: zdept
[921]	90
[825]	91	REAL(wp) :: &
[1310]	92	zsto, zjulian,zout, &
[825]	93	zindh,zinda,zindb
	94	REAL(wp), DIMENSION(jpi,jpj,jpnoumax) :: &
	95	zcmo, &
	96	zcmoa ! additional fields
[921]	97
[825]	98	REAL(wp), DIMENSION(jpi,jpj) :: &
[834]	99	zfield
[825]	100
	101	REAL(wp), DIMENSION(jpi,jpj,jpl) :: &
	102	zmaskitd, zoi, zei
	103
	104	INTEGER :: ji, jj, jk, jl, jf, ipl ! dummy loop indices
	105
	106	CHARACTER(len = 40) :: &
	107	clhstnam, clop, &
[834]	108	clhstnama
[825]	109
	110	INTEGER , SAVE :: &
	111	nice, nhorid, ndim, niter, ndepid
	112	INTEGER , SAVE :: &
[1334]	113	nicea, nhorida, ndimitd
[825]	114	INTEGER , DIMENSION( jpij ) , SAVE :: &
	115	ndex51
	116	INTEGER , DIMENSION( jpij*jpl ) , SAVE :: &
	117	ndexitd
	118	!!-------------------------------------------------------------------
[921]	119
[825]	120	ipl = jpl
	121
	122	IF ( numit == nstart ) THEN
	123
	124	CALL lim_wri_init
[921]	125
[1112]	126	IF(lwp) WRITE(numout,*) ' lim_wri, first time step '
	127	IF(lwp) WRITE(numout,*) ' add_diag_swi ', add_diag_swi
[825]	128
	129	!--------------------
	130	! 1) Initialization
	131	!--------------------
	132
	133	!-------------
	134	! Normal file
	135	!-------------
[921]	136
[825]	137	zsto = rdt_ice
[1312]	138	IF( ln_mskland ) THEN ; clop = "ave(only(x))" ! put 1.e+20 on land (very expensive!!)
	139	ELSE ; clop = "ave(x)" ! no use of the mask value (require less cpu time)
	140	ENDIF
[888]	141	zout = nwrite * rdt_ice / nn_fsbc
[1339]	142	niter = ( nit000 - 1 ) / nn_fsbc
[825]	143	zdept(1) = 0.
[869]	144
[1310]	145	CALL ymds2ju ( nyear, nmonth, nday, rdt, zjulian )
	146	zjulian = zjulian - adatrj ! set calendar origin to the beginning of the experiment
[825]	147	CALL dia_nam ( clhstnam, nwrite, 'icemod' )
[1339]	148	CALL histbeg ( clhstnam, jpi, glamt, jpj, gphit, 1, jpi, 1, jpj, niter, zjulian, rdt_ice, &
[1334]	149	& nhorid, nice, domain_id=nidom )
	150	CALL histvert( nice, "deptht", "Vertical T levels", "m", 1, zdept, ndepid, "down")
[825]	151	CALL wheneq ( jpij , tmask(:,:,1), 1, 1., ndex51, ndim)
[921]	152
[825]	153	DO jf = 1 , noumef
[1112]	154	IF(lwp) WRITE(numout,*) 'jf', jf
[825]	155	IF ( nc(jf) == 1 ) THEN
	156	CALL histdef( nice, nam(jf), titn(jf), uni(jf), jpi, jpj &
	157	, nhorid, 1, 1, 1, -99, 32, clop, zsto, zout )
[1112]	158	IF(lwp) WRITE(numout,*) 'nice, nam(jf), titn(jf), uni(jf), nhorid, clop, zsto, zout'
	159	IF(lwp) WRITE(numout,*) nice, nam(jf), titn(jf), uni(jf), nhorid, clop, zsto, zout
[825]	160	ENDIF
	161	END DO
	162
	163	CALL histend(nice)
[921]	164
[825]	165	!-----------------
	166	! ITD file output
	167	!-----------------
	168	zsto = rdt_ice
	169	clop = "ave(x)"
[888]	170	zout = nwrite * rdt_ice / nn_fsbc
[825]	171	zdept(1) = 0.
	172
	173	CALL dia_nam ( clhstnama, nwrite, 'icemoa' )
	174	CALL histbeg ( clhstnama, jpi, glamt, jpj, gphit, &
[1339]	175	1, jpi, 1, jpj, & ! zoom
	176	niter, zjulian, rdt_ice, & ! time
	177	nhorida, & ! ? linked with horizontal ...
[921]	178	nicea , domain_id=nidom) ! file
[825]	179	CALL histvert( nicea, "icethi", "L levels", &
[921]	180	"m", ipl , hi_mean , nz )
[825]	181	DO jl = 1, jpl
	182	zmaskitd(:,:,jl) = tmask(:,:,1)
	183	END DO
	184	CALL wheneq ( jpij , tmask(:,:,1), 1, 1., ndex51, ndim)
	185	CALL wheneq( jpijpjjpl, zmaskitd, 1, 1., ndexitd, ndimitd )
	186	CALL histdef( nicea, "iice_itd", "Ice area in categories" , "-" , &
[921]	187	jpi, jpj, nhorida, jpl, 1, jpl, nz, 15, clop, zsto, zout )
[825]	188	CALL histdef( nicea, "iice_hid", "Ice thickness in categories" , "m" , &
[921]	189	jpi, jpj, nhorida, jpl, 1, jpl, nz, 15, clop, zsto, zout )
[825]	190	CALL histdef( nicea, "iice_hsd", "Snow depth in in categories" , "m" , &
[921]	191	jpi, jpj, nhorida, jpl, 1, jpl, nz, 15, clop, zsto, zout )
[825]	192	CALL histdef( nicea, "iice_std", "Ice salinity distribution" , "ppt" , &
[921]	193	jpi, jpj, nhorida, jpl, 1, jpl, nz, 15, clop, zsto, zout )
[825]	194	CALL histdef( nicea, "iice_otd", "Ice age distribution" , "days", &
[921]	195	jpi, jpj, nhorida, jpl, 1, jpl, nz, 15, clop, zsto, zout )
[825]	196	CALL histdef( nicea, "iice_etd", "Brine volume distr. " , "%" , &
[921]	197	jpi, jpj, nhorida, jpl, 1, jpl, nz, 15, clop, zsto, zout )
[825]	198	CALL histend(nicea)
	199	ENDIF
	200
[921]	201	! !-----------------------------------------------------------------------!
	202	! !--2. Computation of instantaneous values !
	203	! !-----------------------------------------------------------------------!
	204
	205	!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
[1055]	206	IF( ln_nicep ) THEN
[825]	207	WRITE(numout,*)
	208	WRITE(numout,*) 'lim_wri : write ice outputs in NetCDF files at time : ', nyear, nmonth, nday, numit
	209	WRITE(numout,*) '~~~~~~~ '
	210	WRITE(numout,*) ' kindic = ', kindic
	211	ENDIF
[921]	212	!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
[825]	213
	214	!-- calculs des valeurs instantanees
	215	zcmo( 1:jpi, 1:jpj, 1:jpnoumax ) = 0.0
	216	zcmoa( 1:jpi, 1:jpj, 1:jpnoumax ) = 0.0
	217
	218	DO jl = 1, jpl
	219	DO jj = 1, jpj
	220	DO ji = 1, jpi
	221	zindh = MAX( zzero , SIGN( zone , vt_i(ji,jj) * at_i(ji,jj) - 0.10 ) )
	222	zinda = MAX( zzero , SIGN( zone , at_i(ji,jj) - 0.10 ) )
	223	zcmo(ji,jj,17) = zcmo(ji,jj,17) + a_i(ji,jj,jl)*qsr_ice (ji,jj,jl)
[888]	224	zcmo(ji,jj,18) = zcmo(ji,jj,18) + a_i(ji,jj,jl)*qns_ice(ji,jj,jl)
[825]	225	zcmo(ji,jj,27) = zcmo(ji,jj,27) + t_su(ji,jj,jl)a_i(ji,jj,jl)/MAX(at_i(ji,jj),epsi16)zinda
	226	END DO
	227	END DO
	228	END DO
	229
	230	CALL lim_var_bv
[921]	231
[825]	232	DO jj = 2 , jpjm1
	233	DO ji = 2 , jpim1
	234	zindh = MAX( zzero , SIGN( zone , vt_i(ji,jj) * at_i(ji,jj) - 0.10 ) )
	235	zinda = MAX( zzero , SIGN( zone , at_i(ji,jj) - 0.10 ) )
	236	zindb = zindh * zinda
	237
	238	zcmo(ji,jj,1) = at_i(ji,jj)
	239	zcmo(ji,jj,2) = vt_i(ji,jj)/MAX(at_i(ji,jj),epsi16)*zinda
	240	zcmo(ji,jj,3) = vt_s(ji,jj)/MAX(at_i(ji,jj),epsi16)*zinda
	241	zcmo(ji,jj,4) = diag_bot_gr(ji,jj) * &
[921]	242	86400.0 * zinda !Bottom thermodynamic ice production
[825]	243	zcmo(ji,jj,5) = diag_dyn_gr(ji,jj) * &
[921]	244	86400.0 * zinda !Dynamic ice production (rid/raft)
[825]	245	zcmo(ji,jj,22) = diag_lat_gr(ji,jj) * &
[921]	246	86400.0 * zinda !Lateral thermodynamic ice production
[825]	247	zcmo(ji,jj,23) = diag_sni_gr(ji,jj) * &
[921]	248	86400.0 * zinda !Snow ice production ice production
[825]	249	zcmo(ji,jj,24) = tm_i(ji,jj) - rtt
	250
	251	zcmo(ji,jj,6) = fbif (ji,jj)
	252	zcmo(ji,jj,7) = zindb * ( u_ice(ji,jj ) * tmu(ji,jj) &
[921]	253	& + u_ice(ji-1,jj) * tmu(ji-1,jj) ) &
	254	& / 2.0
[825]	255	zcmo(ji,jj,8) = zindb * ( v_ice(ji,jj ) * tmv(ji,jj) &
[921]	256	& + v_ice(ji,jj-1) * tmv(ji,jj-1) ) &
	257	& / 2.0
[888]	258	zcmo(ji,jj,9) = sst_m(ji,jj)
	259	zcmo(ji,jj,10) = sss_m(ji,jj)
[825]	260
[888]	261	zcmo(ji,jj,11) = qns(ji,jj) + qsr(ji,jj)
	262	zcmo(ji,jj,12) = qsr(ji,jj)
	263	zcmo(ji,jj,13) = qns(ji,jj)
[825]	264	zcmo(ji,jj,14) = fhbri(ji,jj)
[1469]	265	zcmo(ji,jj,15) = utau_ice(ji,jj)
	266	zcmo(ji,jj,16) = vtau_ice(ji,jj)
[888]	267	zcmo(ji,jj,17) = zcmo(ji,jj,17) + (1.0-at_i(ji,jj))*qsr(ji,jj)
	268	zcmo(ji,jj,18) = zcmo(ji,jj,18) + (1.0-at_i(ji,jj))*qns(ji,jj)
[825]	269	zcmo(ji,jj,19) = sprecip(ji,jj)
	270	zcmo(ji,jj,20) = smt_i(ji,jj)
	271	zcmo(ji,jj,21) = ot_i(ji,jj)
	272	zcmo(ji,jj,25) = et_i(ji,jj)
	273	zcmo(ji,jj,26) = et_s(ji,jj)
	274	zcmo(ji,jj,28) = fsbri(ji,jj)
	275	zcmo(ji,jj,29) = fseqv(ji,jj)
[921]	276
[825]	277	zcmo(ji,jj,30) = bv_i(ji,jj)
	278	zcmo(ji,jj,31) = hicol(ji,jj)
	279	zcmo(ji,jj,32) = strength(ji,jj)
	280	zcmo(ji,jj,33) = SQRT( zcmo(ji,jj,7)*zcmo(ji,jj,7) + &
[921]	281	zcmo(ji,jj,8)*zcmo(ji,jj,8) )
[825]	282	zcmo(ji,jj,34) = diag_sur_me(ji,jj) * &
[921]	283	86400.0 * zinda ! Surface melt
[825]	284	zcmo(ji,jj,35) = diag_bot_me(ji,jj) * &
[921]	285	86400.0 * zinda ! Bottom melt
[825]	286	zcmo(ji,jj,36) = divu_i(ji,jj)
	287	zcmo(ji,jj,37) = shear_i(ji,jj)
	288	END DO
	289	END DO
	290
	291	!
	292	! ecriture d'un fichier netcdf
	293	!
	294	niter = niter + 1
	295	DO jf = 1 , noumef
	296	DO jj = 1 , jpj
	297	DO ji = 1 , jpi
	298	zfield(ji,jj) = zcmo(ji,jj,jf) * cmulti(jf) + cadd(jf)
	299	END DO
	300	END DO
[921]	301
[888]	302	IF ( jf == 7 .OR. jf == 8 .OR. jf == 15 .OR. jf == 16 ) THEN
[825]	303	CALL lbc_lnk( zfield, 'T', -1. )
	304	ELSE
	305	CALL lbc_lnk( zfield, 'T', 1. )
	306	ENDIF
	307
[1055]	308	IF( ln_nicep ) THEN
	309	WRITE(numout,*)
	310	WRITE(numout,*) 'nc(jf), nice, nam(jf), niter, ndim'
	311	WRITE(numout,*) nc(jf), nice, nam(jf), niter, ndim
	312	ENDIF
[825]	313	IF ( nc(jf) == 1 ) CALL histwrite( nice, nam(jf), niter, zfield, ndim, ndex51 )
[921]	314
[825]	315	END DO
	316
[1334]	317	IF ( ( nn_fsbc * niter ) >= nitend .OR. kindic < 0 ) THEN
[1112]	318	IF( lwp) WRITE(numout,*) ' Closing the icemod file '
[921]	319	CALL histclo( nice )
[825]	320	ENDIF
	321
	322	!-----------------------------
	323	! Thickness distribution file
	324	!-----------------------------
	325	IF ( add_diag_swi .EQ. 1 ) THEN
	326
[921]	327	DO jl = 1, jpl
	328	CALL lbc_lnk( a_i(:,:,jl) , 'T' , 1. )
	329	CALL lbc_lnk( sm_i(:,:,jl) , 'T' , 1. )
	330	CALL lbc_lnk( oa_i(:,:,jl) , 'T' , 1. )
	331	CALL lbc_lnk( ht_i(:,:,jl) , 'T' , 1. )
	332	CALL lbc_lnk( ht_s(:,:,jl) , 'T' , 1. )
[825]	333	END DO
	334
[921]	335	! Compute ice age
	336	DO jl = 1, jpl
[825]	337	DO jj = 1, jpj
	338	DO ji = 1, jpi
	339	zinda = MAX( zzero , SIGN( zone , a_i(ji,jj,jl) - 1.0e-6 ) )
[921]	340	zoi(ji,jj,jl) = oa_i(ji,jj,jl) / MAX( a_i(ji,jj,jl) , 1.0e-6 ) * &
	341	zinda
[825]	342	END DO
	343	END DO
	344	END DO
	345
[921]	346	! Compute brine volume
	347	zei(:,:,:) = 0.0
	348	DO jl = 1, jpl
	349	DO jk = 1, nlay_i
	350	DO jj = 1, jpj
	351	DO ji = 1, jpi
	352	zinda = MAX( zzero , SIGN( zone , a_i(ji,jj,jl) - 1.0e-6 ) )
	353	zei(ji,jj,jl) = zei(ji,jj,jl) + 100.0* &
	354	( - tmut * s_i(ji,jj,jk,jl) / MIN( ( t_i(ji,jj,jk,jl) - rtt ), -1.0e-6 ) ) * &
	355	zinda / nlay_i
	356	END DO
	357	END DO
	358	END DO
	359	END DO
[825]	360
[921]	361	DO jl = 1, jpl
	362	CALL lbc_lnk( zei(:,:,jl) , 'T' , 1. )
	363	END DO
[825]	364
[921]	365	CALL histwrite( nicea, "iice_itd", niter, a_i , ndimitd , ndexitd ) ! area
	366	CALL histwrite( nicea, "iice_hid", niter, ht_i , ndimitd , ndexitd ) ! thickness
	367	CALL histwrite( nicea, "iice_hsd", niter, ht_s , ndimitd , ndexitd ) ! snow depth
	368	CALL histwrite( nicea, "iice_std", niter, sm_i , ndimitd , ndexitd ) ! salinity
	369	CALL histwrite( nicea, "iice_otd", niter, zoi , ndimitd , ndexitd ) ! age
	370	CALL histwrite( nicea, "iice_etd", niter, zei , ndimitd , ndexitd ) ! brine volume
	371
	372	! ! Create an output files (output.lim.abort.nc) if S < 0 or u > 20 m/s
	373	! IF( kindic < 0 ) CALL lim_wri_state( 'output.abort' )
	374	! not yet implemented
	375
[1334]	376	IF ( ( nn_fsbc * niter ) >= nitend .OR. kindic < 0 ) THEN
[1112]	377	IF(lwp) WRITE(numout,*) ' Closing the icemod file '
[921]	378	CALL histclo( nicea )
	379	ENDIF
	380
[825]	381	ENDIF
	382
	383	END SUBROUTINE lim_wri
	384	#endif
	385
	386	SUBROUTINE lim_wri_init
	387	!!-------------------------------------------------------------------
	388	!! * ROUTINE lim_wri_init *
	389	!!
	390	!! ** Purpose : ???
	391	!!
	392	!! ** Method : Read the namicewri namelist and check the parameter
	393	!! values called at the first timestep (nit000)
	394	!!
	395	!! ** input : Namelist namicewri
	396	!!
	397	!! history :
	398	!! 8.5 ! 03-08 (C. Ethe) original code
	399	!!-------------------------------------------------------------------
	400	!! * Local declarations
	401	INTEGER :: nf ! ???
	402
	403	TYPE FIELD
	404	CHARACTER(len = 35) :: ztitle
	405	CHARACTER(len = 8 ) :: zname
	406	CHARACTER(len = 8 ) :: zunit
	407	INTEGER :: znc
	408	REAL :: zcmulti
	409	REAL :: zcadd
	410	END TYPE FIELD
	411
	412	TYPE(FIELD) :: &
	413	field_1 , field_2 , field_3 , field_4 , field_5 , field_6 , &
	414	field_7 , field_8 , field_9 , field_10, field_11, field_12, &
	415	field_13, field_14, field_15, field_16, field_17, field_18, &
	416	field_19, field_20, field_21, field_22, field_23, field_24, &
	417	field_25, field_26, field_27, field_28, field_29, field_30, &
	418	field_31, field_32, field_33, field_34, field_35, field_36, &
	419	field_37
	420
	421	TYPE(FIELD) , DIMENSION(jpnoumax) :: zfield
	422
	423	NAMELIST/namiceout/ noumef, &
	424	field_1 , field_2 , field_3 , field_4 , field_5 , field_6 , &
	425	field_7 , field_8 , field_9 , field_10, field_11, field_12, &
	426	field_13, field_14, field_15, field_16, field_17, field_18, &
	427	field_19, field_20, field_21, field_22, field_23, field_24, &
	428	field_25, field_26, field_27, field_28, field_29, field_30, &
	429	field_31, field_32, field_33, field_34, field_35, field_36, &
	430	field_37, add_diag_swi
	431	!!-------------------------------------------------------------------
	432
	433	! Read Namelist namicewri
	434	REWIND ( numnam_ice )
	435	READ ( numnam_ice , namiceout )
	436
	437	zfield(1) = field_1
	438	zfield(2) = field_2
	439	zfield(3) = field_3
	440	zfield(4) = field_4
	441	zfield(5) = field_5
	442	zfield(6) = field_6
	443	zfield(7) = field_7
	444	zfield(8) = field_8
	445	zfield(9) = field_9
	446	zfield(10) = field_10
	447	zfield(11) = field_11
	448	zfield(12) = field_12
	449	zfield(13) = field_13
	450	zfield(14) = field_14
	451	zfield(15) = field_15
	452	zfield(16) = field_16
	453	zfield(17) = field_17
	454	zfield(18) = field_18
	455	zfield(19) = field_19
	456	zfield(20) = field_20
	457	zfield(21) = field_21
	458	zfield(22) = field_22
	459	zfield(23) = field_23
	460	zfield(24) = field_24
	461	zfield(25) = field_25
	462	zfield(26) = field_26
	463	zfield(27) = field_27
	464	zfield(28) = field_28
	465	zfield(29) = field_29
	466	zfield(30) = field_30
	467	zfield(31) = field_31
	468	zfield(32) = field_32
	469	zfield(33) = field_33
	470	zfield(34) = field_34
	471	zfield(35) = field_35
	472	zfield(36) = field_36
	473	zfield(37) = field_37
[921]	474
[825]	475	DO nf = 1, noumef
	476	titn (nf) = zfield(nf)%ztitle
	477	nam (nf) = zfield(nf)%zname
	478	uni (nf) = zfield(nf)%zunit
	479	nc (nf) = zfield(nf)%znc
	480	cmulti(nf) = zfield(nf)%zcmulti
	481	cadd (nf) = zfield(nf)%zcadd
	482	END DO
	483
	484	IF(lwp) THEN
	485	WRITE(numout,*)
	486	WRITE(numout,*) 'lim_wri_init : Ice parameters for outputs'
	487	WRITE(numout,*) '~~~~~~~~~~~~'
	488	WRITE(numout,*) ' number of fields to be stored noumef = ', noumef
	489	WRITE(numout,*) ' title name unit Saving (1/0) ', &
	490	& ' multiplicative constant additive constant '
	491	DO nf = 1 , noumef
	492	WRITE(numout,*) ' ', titn(nf), ' ', nam(nf),' ', uni(nf),' ', nc(nf),' ', cmulti(nf), &
	493	' ', cadd(nf)
	494	END DO
	495	WRITE(numout,*) ' add_diag_swi ', add_diag_swi
	496	ENDIF
[921]	497
[825]	498	END SUBROUTINE lim_wri_init
	499
	500	#else
	501	!!----------------------------------------------------------------------
	502	!! Default option : Empty module NO LIM sea-ice model
	503	!!----------------------------------------------------------------------
	504	CONTAINS
	505	SUBROUTINE lim_wri ! Empty routine
	506	END SUBROUTINE lim_wri
	507	#endif
	508
	509	!!======================================================================
	510	END MODULE limwri

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

Download in other formats: