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.

limthd_sal.F90 in branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/LIM_SRC_3 – NEMO

Context Navigation

source: branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/LIM_SRC_3/limthd_sal.F90 @ 2612

Last change on this file since 2612 was 2612, checked in by gm, 13 years ago
dynamic mem: #785 ; LIM-3 case: changes required for compilation (continuation)
Property svn:keywords set to `Id`
File size: 12.7 KB

Rev	Line
[825]	1	MODULE limthd_sal
	2	!!======================================================================
	3	!! * MODULE limthd_sal *
[2528]	4	!! LIM-3 sea-ice : computation of salinity variations in the ice
[825]	5	!!======================================================================
[2528]	6	!! History : - ! 2003-05 (M. Vancoppenolle) UCL-ASTR first coding for LIM3-1D
	7	!! 3.0 ! 2005-12 (M. Vancoppenolle) adapted to the 3-D version
[2612]	8	!! 4.0 ! 2011-02 (G. Madec) dynamical allocation
[2528]	9	!!---------------------------------------------------------------------
[888]	10	#if defined key_lim3
[825]	11	!!----------------------------------------------------------------------
[2528]	12	!! 'key_lim3' LIM-3 sea-ice model
	13	!!----------------------------------------------------------------------
[825]	14	!! lim_thd_sal : salinity variations in the ice
[2528]	15	!!----------------------------------------------------------------------
[825]	16	USE par_oce ! ocean parameters
	17	USE phycst ! physical constants (ocean directory)
[888]	18	USE sbc_oce ! Surface boundary condition: ocean fields
[2612]	19	USE ice ! LIM variables
	20	USE par_ice ! LIM parameters
	21	USE thd_ice ! LIM thermodynamics
	22	USE limvar ! LIM variables
	23	USE wrk_nemo ! workspace manager
[2528]	24	USE in_out_manager ! I/O manager
[825]	25
	26	IMPLICIT NONE
	27	PRIVATE
	28
[2528]	29	PUBLIC lim_thd_sal ! called by limthd module
	30	PUBLIC lim_thd_sal_init ! called by iceini module
[825]	31
[1156]	32	!!----------------------------------------------------------------------
[2612]	33	!! NEMO/LIM3 4.0 , UCL - NEMO Consortium (2011)
[1156]	34	!! $Id$
[2528]	35	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
[1156]	36	!!----------------------------------------------------------------------
[921]	37	CONTAINS
[825]	38
[2528]	39	SUBROUTINE lim_thd_sal( kideb, kiut )
[825]	40	!!-------------------------------------------------------------------
[2528]	41	!! * ROUTINE lim_thd_sal *
	42	!!
	43	!! ** Purpose : computes new salinities in the ice
[825]	44	!!
	45	!! ** Method : 4 possibilities
	46	!! -> num_sal = 1 -> constant salinity for z,t
	47	!! -> num_sal = 2 -> S = S(z,t) [simple Vancoppenolle et al 2005]
	48	!! -> num_sal = 3 -> S = S(z) [multiyear ice]
	49	!! -> num_sal = 4 -> S = S(h) [Cox and Weeks 74]
	50	!!---------------------------------------------------------------------
[2528]	51	INTEGER, INTENT(in) :: kideb, kiut ! thickness category index
	52	!
	53	INTEGER :: ji, jk ! dummy loop indices
	54	INTEGER :: zji, zjj ! local integers
[2612]	55	REAL(wp) :: zsold, iflush, iaccrbo, igravdr, isnowic, i_ice_switch, ztmelts ! local scalars
[2528]	56	REAL(wp) :: zaaa, zbbb, zccc, zdiscrim ! local scalars
[2612]	57	!
	58	REAL(wp), POINTER, DIMENSION(:) :: ze_init, zhiold, zsiold
[2528]	59	!!---------------------------------------------------------------------
[825]	60
[2612]	61	IF( .NOT. wrk_use(1, 1,2,3) ) THEN
	62	CALL ctl_stop('lim_thd_dh : requestead workspace arrays unavailable.') ; RETURN
	63	END IF
	64	! Set-up pointers to sub-arrays of workspace arrays
	65	ze_init => wrk_1d_1 (1:jpij)
	66	zhiold => wrk_1d_2 (1:jpij)
	67	zsiold => wrk_1d_3 (1:jpij)
[825]	68
[921]	69	!------------------------------------------------------------------------------\|
	70	! 1) Constant salinity, constant in time \|
	71	!------------------------------------------------------------------------------\|
[2612]	72	!!gm comment: if num_sal = 1 s_i_b and sm_i_b can be set to bulk_sal one for all in the initialisation phase !!
[2528]	73	IF( num_sal == 1 ) THEN
[2612]	74	!
[825]	75	DO jk = 1, nlay_i
	76	DO ji = kideb, kiut
	77	s_i_b(ji,jk) = bulk_sal
	78	END DO ! ji
	79	END DO ! jk
[2528]	80	!
[825]	81	DO ji = kideb, kiut
	82	sm_i_b(ji) = bulk_sal
	83	END DO ! ji
[2528]	84	!
	85	ENDIF
[825]	86
[921]	87	!------------------------------------------------------------------------------\|
	88	! Module 2 : Constant salinity varying in time \|
	89	!------------------------------------------------------------------------------\|
[825]	90
[2612]	91	IF( num_sal == 2 .OR. num_sal == 4 ) THEN
[825]	92
	93	!---------------------------------
	94	! Thickness at previous time step
	95	!---------------------------------
	96	DO ji = kideb, kiut
[2612]	97	zhiold(ji) = ht_i_b(ji) - dh_i_bott(ji) - dh_snowice(ji) - dh_i_surf(ji)
	98	END DO
[825]	99
	100	!---------------------
	101	! Global heat content
	102	!---------------------
[2612]	103	ze_init(:) = 0._wp
[825]	104	DO jk = 1, nlay_i
	105	DO ji = kideb, kiut
	106	ze_init(ji) = ze_init(ji) + q_i_b(ji,jk) * ht_i_b(ji) / nlay_i
[2612]	107	END DO
	108	END DO
[825]	109
	110	DO ji = kideb, kiut
[2612]	111	!
[825]	112	! Switches
	113	!----------
[2612]	114	iflush = MAX( 0._wp , SIGN( 1.0 , t_su_b(ji) - rtt ) ) ! =1 if summer
	115	igravdr = MAX( 0._wp , SIGN( 1.0 , t_bo_b(ji) - t_su_b(ji) ) ) ! =1 if t_su < t_bo
	116	iaccrbo = MAX( 0._wp , SIGN( 1.0 , dh_i_bott(ji) ) ) ! =1 if bottom accretion
	117	i_ice_switch = 1._wp - MAX ( 0._wp , SIGN( 1._wp , - ht_i_b(ji) + 1.e-2 ) )
	118	isnowic = 1._wp - MAX ( 0._wp , SIGN( 1._wp , - dh_snowice(ji) ) ) * i_ice_switch ! =1 if snow ice formation
[825]	119
	120	!---------------------
	121	! Salinity tendencies
	122	!---------------------
[2612]	123	! ! drainage by gravity drainage
[2528]	124	dsm_i_gd_1d(ji) = - igravdr * MAX( sm_i_b(ji) - sal_G , 0._wp ) / time_G * rdt_ice
[2612]	125	! ! drainage by flushing
	126	dsm_i_fl_1d(ji) = - iflush * MAX( sm_i_b(ji) - sal_F , 0._wp ) / time_F * rdt_ice
[825]	127
	128	!-----------------
	129	! Update salinity
	130	!-----------------
	131	! only drainage terms ( gravity drainage and flushing )
[2612]	132	! snow ice / bottom sources are added in lim_thd_ent to conserve energy
[825]	133	zsiold(ji) = sm_i_b(ji)
[2528]	134	sm_i_b(ji) = sm_i_b(ji) + dsm_i_fl_1d(ji) + dsm_i_gd_1d(ji)
[825]	135
[2612]	136	! if no ice, salinity = 0.1
[2528]	137	i_ice_switch = 1._wp - MAX ( 0._wp, SIGN( 1._wp , - ht_i_b(ji) ) )
[2612]	138	sm_i_b(ji) = i_ice_switch * sm_i_b(ji) + s_i_min * ( 1._wp - i_ice_switch )
[825]	139	END DO ! ji
	140
	141	! Salinity profile
[2528]	142	CALL lim_var_salprof1d( kideb, kiut )
[825]	143
	144	!----------------------------
	145	! Heat flux - brine drainage
	146	!----------------------------
	147
	148	DO ji = kideb, kiut
[2612]	149	!!gm useless
[825]	150	! iflush : 1 if summer
[2528]	151	iflush = MAX( 0._wp , SIGN ( 1._wp , t_su_b(ji) - rtt ) )
[825]	152	! igravdr : 1 if t_su lt t_bo
[2528]	153	igravdr = MAX( 0._wp , SIGN ( 1._wp , t_bo_b(ji) - t_su_b(ji) ) )
[825]	154	! iaccrbo : 1 if bottom accretion
[2528]	155	iaccrbo = MAX( 0._wp , SIGN ( 1._wp , dh_i_bott(ji) ) )
[2612]	156	!!gm end useless
[2528]	157	!
	158	fhbri_1d(ji) = 0._wp
[825]	159	END DO ! ji
	160
	161	!----------------------------
	162	! Salt flux - brine drainage
	163	!----------------------------
	164	DO ji = kideb, kiut
[2528]	165	i_ice_switch = 1._wp - MAX ( 0._wp, SIGN( 1._wp , - ht_i_b(ji) ) )
	166	fsbri_1d(ji) = fsbri_1d(ji) - i_ice_switch * rhoic * a_i_b(ji) * ht_i_b(ji) &
	167	& * ( MAX(dsm_i_gd_1d(ji) + dsm_i_fl_1d(ji), sm_i_b(ji) - zsiold(ji) ) ) / rdt_ice
	168	IF( num_sal == 4 ) fsbri_1d(ji) = 0._wp
[825]	169	END DO ! ji
	170
	171	! Only necessary for conservation check since salinity is modified
	172	!--------------------
	173	! Temperature update
	174	!--------------------
	175	DO jk = 1, nlay_i
	176	DO ji = kideb, kiut
	177	ztmelts = -tmut*s_i_b(ji,jk) + rtt
	178	!Conversion q(S,T) -> T (second order equation)
	179	zaaa = cpic
[2528]	180	zbbb = ( rcp - cpic ) * ( ztmelts - rtt ) + q_i_b(ji,jk) / rhoic - lfus
[825]	181	zccc = lfus * ( ztmelts - rtt )
[2612]	182	zdiscrim = SQRT( MAX( zbbbzbbb - 4.0zaaa*zccc, 0._wp ) )
[2528]	183	t_i_b(ji,jk) = rtt - ( zbbb + zdiscrim ) / ( 2.0 *zaaa )
[2612]	184	END DO
	185	END DO
[2528]	186	!
[825]	187	ENDIF ! num_sal .EQ. 2
	188
[921]	189	!------------------------------------------------------------------------------\|
	190	! Module 3 : Profile of salinity, constant in time \|
	191	!------------------------------------------------------------------------------\|
[825]	192
[2612]	193	IF( num_sal == 3 ) CALL lim_var_salprof1d( kideb, kiut )
[825]	194
[921]	195	!------------------------------------------------------------------------------\|
	196	! Module 4 : Constant salinity varying in time \|
	197	!------------------------------------------------------------------------------\|
[825]	198
[2612]	199	IF( num_sal == 5 ) THEN ! Cox and Weeks, 1974
	200	!
[825]	201	DO ji = kideb, kiut
	202	zsold = sm_i_b(ji)
[2612]	203	IF( ht_i_b(ji) < 0.4 ) THEN
	204	sm_i_b(ji) = 14.24 - 19.39 * ht_i_b(ji)
[825]	205	ELSE
[2612]	206	sm_i_b(ji) = 7.88 - 1.59 * ht_i_b(ji)
	207	sm_i_b(ji) = MIN( sm_i_b(ji) , zsold )
[825]	208	ENDIF
[2612]	209	IF( ht_i_b(ji) > 3.06918239 ) THEN
	210	sm_i_b(ji) = 3._wp
[825]	211	ENDIF
	212	DO jk = 1, nlay_i
	213	s_i_b(ji,jk) = sm_i_b(ji)
	214	END DO
[2612]	215	END DO
	216	!
[825]	217	ENDIF ! num_sal
	218
[921]	219	!------------------------------------------------------------------------------\|
	220	! 5) Computation of salt flux due to Bottom growth
	221	!------------------------------------------------------------------------------\|
[825]	222
[2612]	223	IF ( num_sal == 4 ) THEN
[825]	224	DO ji = kideb, kiut
[2612]	225	zji = MOD( npb(ji) - 1 , jpi ) + 1
	226	zjj = ( npb(ji) - 1 ) / jpi + 1
[2528]	227	fseqv_1d(ji) = fseqv_1d(ji) + ( sss_m(zji,zjj) - bulk_sal ) &
	228	& * rhoic * a_i_b(ji) * MAX( dh_i_bott(ji) , 0.0 ) / rdt_ice
[825]	229	END DO
	230	ELSE
	231	DO ji = kideb, kiut
[2612]	232	zji = MOD( npb(ji) - 1 , jpi ) + 1
	233	zjj = ( npb(ji) - 1 ) / jpi + 1
[2528]	234	fseqv_1d(ji) = fseqv_1d(ji) + ( sss_m(zji,zjj) - s_i_new(ji) ) &
	235	& * rhoic * a_i_b(ji) * MAX( dh_i_bott(ji) , 0.0 ) / rdt_ice
[2612]	236	END DO
[825]	237	ENDIF
[2528]	238	!
[2612]	239	IF( .NOT. wrk_release(1, 1,2,3) ) CALL ctl_stop( 'lim_thd_lac : failed to release workspace arrays.' )
	240	!
[825]	241	END SUBROUTINE lim_thd_sal
	242
[2528]	243
[825]	244	SUBROUTINE lim_thd_sal_init
	245	!!-------------------------------------------------------------------
	246	!! * ROUTINE lim_thd_sal_init *
	247	!!
	248	!! ** Purpose : initialization of ice salinity parameters
	249	!!
[2528]	250	!! ** Method : Read the namicesal namelist and check the parameter
	251	!! values called at the first timestep (nit000)
[825]	252	!!
	253	!! ** input : Namelist namicesal
	254	!!-------------------------------------------------------------------
[2528]	255	NAMELIST/namicesal/ num_sal, bulk_sal, sal_G, time_G, sal_F, time_F, &
	256	& s_i_max, s_i_min, s_i_0, s_i_1
[825]	257	!!-------------------------------------------------------------------
[2528]	258	!
	259	REWIND( numnam_ice ) ! Read Namelist namicesal
	260	READ ( numnam_ice , namicesal )
	261	!
	262	IF(lwp) THEN ! control print
[825]	263	WRITE(numout,*)
	264	WRITE(numout,*) 'lim_thd_sal_init : Ice parameters for salinity '
	265	WRITE(numout,*) '~~~~~~~~~~~~~~~~'
	266	WRITE(numout,*) ' switch for salinity num_sal : ', num_sal
	267	WRITE(numout,*) ' bulk salinity value if num_sal = 1 : ', bulk_sal
	268	WRITE(numout,*) ' restoring salinity for GD : ', sal_G
	269	WRITE(numout,*) ' restoring time for GD : ', time_G
	270	WRITE(numout,*) ' restoring salinity for flushing : ', sal_F
	271	WRITE(numout,*) ' restoring time for flushing : ', time_F
	272	WRITE(numout,*) ' Maximum tolerated ice salinity : ', s_i_max
	273	WRITE(numout,*) ' Minimum tolerated ice salinity : ', s_i_min
	274	WRITE(numout,*) ' 1st salinity for salinity profile : ', s_i_0
	275	WRITE(numout,*) ' 2nd salinity for salinity profile : ', s_i_1
	276	ENDIF
[2528]	277	!
[825]	278	END SUBROUTINE lim_thd_sal_init
	279
	280	#else
	281	!!----------------------------------------------------------------------
[2528]	282	!! Default option Dummy Module No LIM-3 sea-ice model
[825]	283	!!----------------------------------------------------------------------
	284	#endif
	285	!!======================================================================
	286	END MODULE limthd_sal

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

Download in other formats: