Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

limdia_2.F90 @ 6736

Last change on this file since 6736 was 6736, checked in by jamesharle, 8 years ago
FASTNEt code modifications
Property svn:keywords set to `Id`
File size: 12.7 KB

Rev	Line
[821]	1	MODULE limdia_2
[3]	2	!!======================================================================
[821]	3	!! * MODULE limdia_2 *
[3]	4	!! diagnostics of ice model
	5	!!======================================================================
[508]	6	!! History : 8.0 ! 97-06 (Louvain-La-Neuve) Original code
	7	!! 8.5 ! 02-09 (C. Ethe , G. Madec ) F90: Free form and module
	8	!! 9.0 ! 06-08 (S. Masson) change frequency output control
	9	!!-------------------------------------------------------------------
[821]	10	#if defined key_lim2
[3]	11	!!----------------------------------------------------------------------
[821]	12	!! 'key_lim2' : LIM 2.0 sea-ice model
[3]	13	!!----------------------------------------------------------------------
[821]	14	!! lim_dia_2 : computation of the time evolution of keys var.
	15	!! lim_dia_init_2 : initialization and namelist read
[3]	16	!!----------------------------------------------------------------------
[1715]	17	USE dom_oce ! ocean space and time domain
[77]	18	USE phycst !
[821]	19	USE par_ice_2 ! ice parameters
[888]	20	USE sbc_oce ! surface boundary condition variables
[821]	21	USE dom_ice_2 !
	22	USE ice_2 !
	23	USE limistate_2 !
[77]	24	USE in_out_manager ! I/O manager
[2715]	25	USE lib_mpp ! MPP library
[3]	26
	27	IMPLICIT NONE
	28	PRIVATE
	29
[888]	30	PUBLIC lim_dia_2 ! called by sbc_ice_lim_2
[2715]	31
[508]	32	INTEGER, PUBLIC :: ntmoy = 1 , & !: instantaneous values of ice evolution or averaging ntmoy
	33	& ninfo = 1 !: frequency of ouputs on file ice_evolu in case of averaging
[3]	34
[77]	35	INTEGER, PARAMETER :: & ! Parameters for outputs to files "evolu"
	36	jpinfmx = 100 , & ! maximum number of key variables
	37	jpchinf = 5 , & ! ???
	38	jpchsep = jpchinf + 2 ! ???
	39
	40	INTEGER :: &
	41	nfrinf = 4 , & ! number of variables written in one line
[3]	42	nferme , & ! last time step at which the var. are written on file
	43	nvinfo , & ! number of total variables
	44	nbvt , & ! number of time variables
	45	naveg ! number of step for accumulation before averaging
	46
[508]	47	CHARACTER(len= 8) :: fmtinf = '1PE13.5 ' ! format of the output values
	48	CHARACTER(len=30) :: fmtw , & ! formats
	49	& fmtr , & ! ???
	50	& fmtitr ! ???
	51	CHARACTER(len=jpchsep), DIMENSION(jpinfmx) :: titvar ! title of key variables
[77]	52
[508]	53	REAL(wp) :: epsi06 = 1.e-06 ! ???
	54	REAL(wp), DIMENSION(jpinfmx) :: vinfom ! temporary working space
[2715]	55	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: aire ! masked grid cell area
[3]	56
	57	!! * Substitutions
	58	# include "vectopt_loop_substitute.h90"
	59	!!----------------------------------------------------------------------
[2528]	60	!! NEMO/LIM2 3.3 , UCL - NEMO Consortium (2010)
[1156]	61	!! $Id$
[2528]	62	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
[3]	63	!!----------------------------------------------------------------------
	64	CONTAINS
	65
[821]	66	SUBROUTINE lim_dia_2( kt )
[3]	67	!!--------------------------------------------------------------------
[821]	68	!! * ROUTINE lim_dia_2 *
[3]	69	!!
	70	!! ** Purpose : Computation and outputs on file ice.evolu
	71	!! the temporal evolution of some key variables
[508]	72	!!-------------------------------------------------------------------
	73	INTEGER, INTENT(in) :: kt ! number of iteration
[3]	74	!!
[508]	75	INTEGER :: jv,ji, jj ! dummy loop indices
	76	INTEGER :: nv ! indice of variable
	77	REAL(wp) :: zarea , zldarea , & ! sea-ice and leads area
	78	& zextent15, zextent85, & ! sea-ice extent (15% and 85%)
	79	& zicevol , zsnwvol , & ! sea-ice and snow volume volume
	80	& zicespd ! sea-ice velocity
	81	REAL(wp), DIMENSION(jpinfmx) :: vinfor ! temporary working space
[3]	82	!!-------------------------------------------------------------------
	83
[821]	84	IF( kt == nit000 ) CALL lim_dia_init_2 ! initialisation of ice_evolu file
[3]	85
[508]	86	! computation of key variables at each time step
[3]	87
[508]	88	nv = 1
[888]	89	vinfor(nv) = REAL( kt + nn_fsbc - 1 )
[508]	90	nv = nv + 1
	91	vinfor(nv) = nyear
[3]	92
[508]	93	DO jv = nbvt + 1, nvinfo
	94	vinfor(jv) = 0.e0
	95	END DO
[3]	96
[508]	97	zextent15 = 0.e0
	98	zextent85 = 0.e0
	99	! variables in northern Hemis
	100	DO jj = njeq, jpjm1
	101	DO ji = fs_2, fs_jpim1 ! vector opt.
	102	IF( tms(ji,jj) == 1 ) THEN
	103	zarea = ( 1.0 - frld(ji,jj) ) * aire(ji,jj)
	104	IF (frld(ji,jj) <= 0.15 ) zextent15 = aire(ji,jj)
	105	IF (frld(ji,jj) <= 0.85 ) zextent85 = aire(ji,jj)
	106	zldarea = zarea / MAX( ( 1 - frld(ji,jj) ) , epsi06 )
	107	zicevol = zarea * hicif(ji,jj)
	108	zsnwvol = zarea * hsnif(ji,jj)
[1470]	109	zicespd = zicevol * ( u_ice(ji,jj) * u_ice(ji,jj) + v_ice(ji,jj) * v_ice(ji,jj) )
[508]	110	vinfor(nv+ 1) = vinfor(nv+ 1) + zarea
	111	vinfor(nv+ 3) = vinfor(nv+ 3) + zextent15
	112	vinfor(nv+ 5) = vinfor(nv+ 5) + zextent85
	113	vinfor(nv+ 7) = vinfor(nv+ 7) + zldarea
	114	vinfor(nv+ 9) = vinfor(nv+ 9) + zicevol
	115	vinfor(nv+11) = vinfor(nv+11) + zsnwvol
	116	vinfor(nv+13) = vinfor(nv+13) + zicespd
	117	ENDIF
	118	END DO
	119	END DO
	120	vinfor(nv+13) = SQRT( vinfor(nv+13) / MAX( vinfor(nv+9) , epsi06 ) )
[3]	121
	122
	123	! variables in southern Hemis
	124	nv = nv + 1
[77]	125	DO jj = 2, njeqm1
[12]	126	DO ji = fs_2, fs_jpim1 ! vector opt.
[3]	127	IF( tms(ji,jj) == 1 ) THEN
[12]	128	zarea = ( 1.0 - frld(ji,jj) ) * aire(ji,jj)
	129	IF (frld(ji,jj) <= 0.15 ) zextent15 = aire(ji,jj)
	130	IF (frld(ji,jj) <= 0.85 ) zextent85 = aire(ji,jj)
	131	zldarea = zarea / MAX( ( 1 - frld(ji,jj) ) , epsi06 )
	132	zicevol = zarea * hicif(ji,jj)
	133	zsnwvol = zarea * hsnif(ji,jj)
[1470]	134	zicespd = zicevol * ( u_ice(ji,jj) * u_ice(ji,jj) + v_ice(ji,jj) * v_ice(ji,jj) )
[12]	135	vinfor(nv+ 1) = vinfor(nv+ 1) + zarea
	136	vinfor(nv+ 3) = vinfor(nv+ 3) + zextent15
	137	vinfor(nv+ 5) = vinfor(nv+ 5) + zextent85
	138	vinfor(nv+ 7) = vinfor(nv+ 7) + zldarea
	139	vinfor(nv+ 9) = vinfor(nv+ 9) + zicevol
	140	vinfor(nv+11) = vinfor(nv+11) + zsnwvol
	141	vinfor(nv+13) = vinfor(nv+13) + zicespd
[3]	142	ENDIF
	143	END DO
	144	END DO
[12]	145	vinfor(nv+13) = SQRT( vinfor(nv+13) / MAX( vinfor(nv+9) , epsi06 ) )
[3]	146
	147	! Accumulation before averaging
[12]	148	DO jv = 1, nvinfo
[3]	149	vinfom(jv) = vinfom(jv) + vinfor(jv)
	150	END DO
	151	naveg = naveg + 1
	152
	153	! oututs on file ice_evolu
[888]	154	IF( MOD( kt + nn_fsbc - 1, ninfo ) == 0 ) THEN
[98]	155	WRITE(numevo_ice,fmtw) ( titvar(jv), vinfom(jv)/naveg, jv = 1, nvinfo )
[3]	156	naveg = 0
[12]	157	DO jv = 1, nvinfo
[77]	158	vinfom(jv) = 0.e0
[3]	159	END DO
	160	ENDIF
[508]	161	!
[821]	162	END SUBROUTINE lim_dia_2
[3]	163
	164
[821]	165	SUBROUTINE lim_dia_init_2
[3]	166	!!-------------------------------------------------------------------
[821]	167	!! * ROUTINE lim_dia_init_2 *
[3]	168	!!
	169	!! ** Purpose : Preparation of the file ice_evolu for the output of
	170	!! the temporal evolution of key variables
	171	!!
	172	!! ** input : Namelist namicedia
	173	!!-------------------------------------------------------------------
[508]	174	CHARACTER(len=jpchinf) :: titinf
[2715]	175	INTEGER :: jv ! dummy loop indice
	176	INTEGER :: ntot , ndeb, nv, ierr ! local integer
	177	REAL(wp) :: zxx0, zxx1 ! local scalars
[3]	178
[508]	179	NAMELIST/namicedia/fmtinf, nfrinf, ninfo, ntmoy
[3]	180	!!-------------------------------------------------------------------
	181
[2715]	182	REWIND( numnam_ice ) ! Read Namelist namicedia
	183	READ ( numnam_ice , namicedia )
[508]	184
[2715]	185	IF(lwp) THEN ! control print
[3]	186	WRITE(numout,*)
[821]	187	WRITE(numout,*) 'lim_dia_init_2 : ice parameters for ice diagnostics '
	188	WRITE(numout,*) '~~~~~~~~~~~~~~'
[3]	189	WRITE(numout,*) ' format of the output values fmtinf = ', fmtinf
	190	WRITE(numout,*) ' number of variables written in one line nfrinf = ', nfrinf
	191	WRITE(numout,*) ' Instantaneous values of ice evolution or averaging ntmoy = ', ntmoy
	192	WRITE(numout,*) ' frequency of ouputs on file ice_evolu in case of averaging ninfo = ', ninfo
	193	ENDIF
	194
[2715]	195	ALLOCATE( aire(jpi,jpj) , STAT=ierr ) ! masked grid cell area
	196	IF( lk_mpp ) CALL mpp_sum( ierr )
	197	IF( ierr /= 0 ) CALL ctl_stop( 'STOP', 'lim_dia_init_2 : unable to allocate standard arrays' )
[77]	198	aire(:,:) = area(:,:) * tms(:,:)
	199
[2715]	200	nv = 1 ! Titles of ice key variables
[3]	201	titvar(nv) = 'NoIt' ! iteration number
	202	nv = nv + 1
	203	titvar(nv) = 'T yr' ! time step in years
	204	nbvt = nv - 1
[508]	205	nv = nv + 1 ; titvar(nv) = 'AEFN' ! sea ice area in the northern Hemisp.(10^12 km2)
	206	nv = nv + 1 ; titvar(nv) = 'AEFS' ! sea ice area in the southern Hemisp.(10^12 km2)
	207	nv = nv + 1 ; titvar(nv) = 'A15N' ! sea ice extent (15%) in the northern Hemisp.(10^12 km2)
	208	nv = nv + 1 ; titvar(nv) = 'A15S' ! sea ice extent (15%) in the southern Hemisp.(10^12 km2)
	209	nv = nv + 1 ; titvar(nv) = 'A85N' ! sea ice extent (85%) in the northern Hemisp.(10^12 km2)
	210	nv = nv + 1 ; titvar(nv) = 'A85S' ! sea ice extent (85%) in the southern Hemisp.(10^12 km2)
	211	nv = nv + 1 ; titvar(nv) = 'ALEN' ! leads area in the northern Hemisp.(10^12 km2)
	212	nv = nv + 1 ; titvar(nv) = 'ALES' ! leads area in the southern Hemisp.(10^12 km2)
	213	nv = nv + 1 ; titvar(nv) = 'VOLN' ! sea ice volume in the northern Hemisp.(10^3 km3)
	214	nv = nv + 1 ; titvar(nv) = 'VOLS' ! sea ice volume in the southern Hemisp.(10^3 km3)
	215	nv = nv + 1 ; titvar(nv) = 'VONN' ! snow volume over sea ice in the northern Hemisp.(10^3 km3)
	216	nv = nv + 1 ; titvar(nv) = 'VONS' ! snow volume over sea ice in the southern Hemisp.(10^3 km3)
	217	nv = nv + 1 ; titvar(nv) = 'ECGN' ! mean sea ice velocity in the northern Hemisp.(m/s)
	218	nv = nv + 1 ; titvar(nv) = 'ECGS' ! mean sea ice velocity in the southern Hemisp.(m/s)
[3]	219
	220	nvinfo = nv
	221
	222	! Definition et Ecriture de l'entete : nombre d'enregistrements
[888]	223	ndeb = ( nit000 - 1 + nn_fsbc - 1 ) / ninfo
	224	IF( nit000 - 1 + nn_fsbc == 1 ) ndeb = -1
[3]	225
[888]	226	nferme = ( nitend + nn_fsbc - 1 ) / ninfo ! nit000 - 1 + nn_fsbc - 1 + nitend - nit000 + 1
[3]	227	ntot = nferme - ndeb
	228	ndeb = ninfo * ( 1 + ndeb )
	229	nferme = ninfo * nferme
	230
	231	! definition of formats
[77]	232	WRITE( fmtw , '(A,I3,A2,I1,A)' ) '(', nfrinf, '(A', jpchsep, ','//fmtinf//'))'
	233	WRITE( fmtr , '(A,I3,A,I1,A)' ) '(', nfrinf, '(', jpchsep, 'X,'//fmtinf//'))'
	234	WRITE( fmtitr, '(A,I3,A,I1,A)' ) '(', nvinfo, 'A', jpchinf, ')'
[3]	235
	236	! opening "ice_evolu" file
[1581]	237	CALL ctl_opn( numevo_ice, 'ice_evolu', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
[3]	238
	239	!- ecriture de 2 lignes d''entete :
	240	WRITE(numevo_ice,1000) fmtr, fmtw, fmtitr, nvinfo, ntot, 0, nfrinf
[77]	241	zxx0 = 0.001 * REAL( ninfo )
	242	zxx1 = 0.001 * REAL( ndeb )
	243	WRITE(numevo_ice,1111) REAL(jpchinf), 0., zxx1, zxx0, 0., 0., 0
[3]	244
	245	!- ecriture de 2 lignes de titre :
[508]	246	WRITE(numevo_ice,'(A,I8,A,I8,A,I5)') &
[3]	247	'Evolution chronologique - Experience '//cexper &
	248	//' de', ndeb, ' a', nferme, ' pas', ninfo
	249	WRITE(numevo_ice,fmtitr) ( titvar(jv), jv = 1, nvinfo )
	250
	251
	252	!--preparation de "titvar" pour l''ecriture parmi les valeurs numeriques :
	253	DO jv = 2 , nvinfo
[77]	254	titinf = titvar(jv)(:jpchinf)
[3]	255	titvar(jv) = ' '//titinf
	256	END DO
	257
	258	!--Initialisation of the arrays for the accumulation
	259	DO jv = 1, nvinfo
	260	vinfom(jv) = 0.
	261	END DO
	262	naveg = 0
	263
	264	1000 FORMAT( 3(A20),4(1x,I6) )
	265	1111 FORMAT( 3(F7.1,1X,F7.3,1X),I3,A )
[508]	266	!
[821]	267	END SUBROUTINE lim_dia_init_2
[3]	268
	269	#else
	270	!!----------------------------------------------------------------------
[821]	271	!! Default option : NO LIM 2.0 sea-ice model
[3]	272	!!----------------------------------------------------------------------
	273	CONTAINS
[821]	274	SUBROUTINE lim_dia_2 ! Empty routine
	275	END SUBROUTINE lim_dia_2
[3]	276	#endif
	277
	278	!!======================================================================
[821]	279	END MODULE limdia_2

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

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

Context Navigation

source: branches/NERC/dev_r3874_FASTNEt/NEMOGCM/NEMO/LIM_SRC_2/limdia_2.F90 @ 6736

Download in other formats: