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.

sbcabl.F90 in NEMO/trunk/src/ABL – NEMO

Context Navigation

source: NEMO/trunk/src/ABL/sbcabl.F90 @ 12549

Last change on this file since 12549 was 12549, checked in by smasson, 4 years ago
trunk: fix compilation bug for ABL/sbcabl.F90, see #2411
File size: 20.1 KB

Rev	Line
[11305]	1	MODULE sbcabl
	2	!!======================================================================
	3	!! * MODULE sbcabl *
[12199]	4	!! Ocean forcing: momentum, heat and freshwater flux formulation
	5	!! derived from an ABL model
[11305]	6	!!=====================================================================
	7	!! History : 4.0 ! 2019-03 (F. Lemarié & G. Samson) Original code
	8	!!----------------------------------------------------------------------
	9
	10	!!----------------------------------------------------------------------
[12199]	11	!! sbc_abl_init : Initialization of ABL model based on namelist options
	12	!! sbc_abl : driver for the computation of momentum, heat and freshwater
[11305]	13	!! fluxes over ocean via the ABL model
	14	!!----------------------------------------------------------------------
	15	USE abl ! ABL
	16	USE par_abl ! abl parameters
	17	USE ablmod
[11858]	18	USE ablrst
[11305]	19
	20	USE phycst ! physical constants
	21	USE fldread ! read input fields
	22	USE sbc_oce ! Surface boundary condition: ocean fields
	23	USE sbcblk ! Surface boundary condition: bulk formulae
[12021]	24	USE sbcblk_phy ! Surface boundary condition: bulk formulae
[11305]	25	USE dom_oce, ONLY : tmask
	26	!
	27	USE iom ! I/O manager library
	28	USE in_out_manager ! I/O manager
	29	USE lib_mpp ! distribued memory computing library
	30	USE lib_fortran ! to use key_nosignedzero
	31	USE timing ! Timing
	32	USE lbclnk ! ocean lateral boundary conditions (or mpp link)
	33	USE prtctl ! Print control
[11334]	34	#if defined key_si3
	35	USE ice , ONLY : u_ice, v_ice, tm_su, ato_i ! ato_i = total open water fractional area
	36	USE sbc_ice, ONLY : wndm_ice, utau_ice, vtau_ice
[12199]	37	#endif
[11334]	38	#if ! defined key_iomput
	39	USE diawri , ONLY : dia_wri_alloc_abl
	40	#endif
[11305]	41	IMPLICIT NONE
	42	PRIVATE
	43
	44	PUBLIC sbc_abl_init ! routine called in sbcmod module
	45	PUBLIC sbc_abl ! routine called in sbcmod module
	46
	47	!!----------------------------------------------------------------------
	48	!! NEMO/OPA 3.7 , NEMO-consortium (2014)
	49	!! $Id: sbcabl.F90 6416 2016-04-01 12:22:17Z clem $
	50	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
	51	!!----------------------------------------------------------------------
	52	CONTAINS
	53
	54	SUBROUTINE sbc_abl_init
	55	!!---------------------------------------------------------------------
	56	!! * ROUTINE sbc_abl_init *
	57	!!
	58	!! ** Purposes : - read namelist section namsbc_abl
[12199]	59	!! - initialize and check parameter values
	60	!! - initialize variables of ABL model
[11305]	61	!!
	62	!!----------------------------------------------------------------------
	63	INTEGER :: ji, jj, jk, jbak, jbak_dta ! dummy loop indices
	64	INTEGER :: ios, ierror, ioptio ! Local integer
	65	INTEGER :: inum, indims, idimsz(4), id
	66	CHARACTER(len=100) :: cn_dir, cn_dom ! Atmospheric grid directory
	67	REAL(wp) :: zcff,zcff1
	68	LOGICAL :: lluldl
[11858]	69	NAMELIST/namsbc_abl/ cn_dir, cn_dom, cn_ablrst_in, cn_ablrst_out, &
	70	& cn_ablrst_indir, cn_ablrst_outdir, &
	71	& ln_hpgls_frc, ln_geos_winds, nn_dyn_restore, &
[11305]	72	& rn_ldyn_min , rn_ldyn_max, rn_ltra_min, rn_ltra_max, &
[11322]	73	& nn_amxl, rn_cm, rn_ct, rn_ce, rn_ceps, rn_Rod, rn_Ric, &
[12199]	74	& ln_smth_pblh
[11305]	75	!!---------------------------------------------------------------------
	76
[12549]	77	! Namelist namsbc_abl in reference namelist : ABL parameters
[11305]	78	READ ( numnam_ref, namsbc_abl, IOSTAT = ios, ERR = 901 )
[11348]	79	901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_abl in reference namelist' )
[12549]	80	! Namelist namsbc_abl in configuration namelist : ABL parameters
[11305]	81	READ ( numnam_cfg, namsbc_abl, IOSTAT = ios, ERR = 902 )
[11348]	82	902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_abl in configuration namelist' )
[11305]	83	!
	84	IF(lwm) WRITE( numond, namsbc_abl )
	85	!
[12199]	86	! Check ABL mixing length option
[11305]	87	IF( nn_amxl < 0 .OR. nn_amxl > 2 ) &
[12199]	88	& CALL ctl_stop( 'abl_init : bad flag, nn_amxl must be 0, 1 or 2 ' )
[11305]	89	!
[12199]	90	! Check ABL dyn restore option
[11305]	91	IF( nn_dyn_restore < 0 .OR. nn_dyn_restore > 2 ) &
[12199]	92	& CALL ctl_stop( 'abl_init : bad flag, nn_dyn_restore must be 0, 1 or 2 ' )
[12063]	93
[11305]	94	!!---------------------------------------------------------------------
	95	!! Control prints
[12199]	96	!!---------------------------------------------------------------------
[11305]	97	IF(lwp) THEN ! Control print (other namelist variable)
	98	WRITE(numout,*)
	99	WRITE(numout,*) ' ABL -- cn_dir = ', cn_dir
	100	WRITE(numout,*) ' ABL -- cn_dom = ', cn_dom
	101	IF( ln_hpgls_frc ) THEN
	102	WRITE(numout,*) ' ABL -- winds forced by large-scale pressure gradient'
	103	IF(ln_geos_winds) THEN
	104	ln_geos_winds = .FALSE.
[12199]	105	WRITE(numout,*) ' ABL -- geostrophic guide disabled (not compatible with ln_hpgls_frc = .T.)'
[11305]	106	END IF
	107	ELSE IF( ln_geos_winds ) THEN
[12199]	108	WRITE(numout,*) ' ABL -- winds forced by geostrophic winds'
[11305]	109	ELSE
[12199]	110	WRITE(numout,*) ' ABL -- Geostrophic winds and large-scale pressure gradient are ignored'
[11305]	111	END IF
	112	!
	113	SELECT CASE ( nn_dyn_restore )
[12199]	114	CASE ( 0 )
	115	WRITE(numout,*) ' ABL -- No restoring for ABL winds'
[11305]	116	CASE ( 1 )
[12199]	117	WRITE(numout,*) ' ABL -- Restoring of ABL winds only in the equatorial region '
[11305]	118	CASE ( 2 )
[12199]	119	WRITE(numout,*) ' ABL -- Restoring of ABL winds activated everywhere '
[11305]	120	END SELECT
	121	!
[12199]	122	IF( ln_smth_pblh ) WRITE(numout,*) ' ABL -- Smoothing of PBL height is activated'
	123	!
[11305]	124	ENDIF
	125
	126	!!---------------------------------------------------------------------
	127	!! Convert nudging coefficient from hours to 1/sec
[12199]	128	!!---------------------------------------------------------------------
[11305]	129	zcff = 1._wp / 3600._wp
	130	rn_ldyn_min = zcff / rn_ldyn_min
[12199]	131	rn_ldyn_max = zcff / rn_ldyn_max
[11305]	132	rn_ltra_min = zcff / rn_ltra_min
[12199]	133	rn_ltra_max = zcff / rn_ltra_max
[11305]	134
	135	!!---------------------------------------------------------------------
[12199]	136	!! ABL grid initialization
	137	!!---------------------------------------------------------------------
	138	CALL iom_open( TRIM(cn_dir)//TRIM(cn_dom), inum )
[11305]	139	id = iom_varid( inum, 'e3t_abl', kdimsz=idimsz, kndims=indims, lduld=lluldl )
[12199]	140	jpka = idimsz(indims - COUNT( (/lluldl/) ) )
[11305]	141	jpkam1 = jpka - 1
	142
	143	IF( abl_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'abl_init : unable to allocate arrays' )
	144	CALL iom_get( inum, jpdom_unknown, 'e3t_abl', e3t_abl(:) )
	145	CALL iom_get( inum, jpdom_unknown, 'e3w_abl', e3w_abl(:) )
	146	CALL iom_get( inum, jpdom_unknown, 'ght_abl', ght_abl(:) )
	147	CALL iom_get( inum, jpdom_unknown, 'ghw_abl', ghw_abl(:) )
	148	CALL iom_close( inum )
[11348]	149
[11334]	150	#if ! defined key_iomput
[12199]	151	IF( dia_wri_alloc_abl() /= 0 ) CALL ctl_stop( 'STOP', 'abl_init : unable to allocate arrays' )
[11334]	152	#endif
[11305]	153
	154	IF(lwp) THEN
	155	WRITE(numout,*)
	156	WRITE(numout,*) ' sbc_abl_init : ABL Reference vertical grid'
[12199]	157	WRITE(numout,*) ' ~~~~~~~'
[11305]	158	WRITE(numout, "(9x,' level ght_abl ghw_abl e3t_abl e3w_abl ')" )
	159	WRITE(numout, "(10x, i4, 4f9.2)" ) ( jk, ght_abl(jk), ghw_abl(jk), e3t_abl(jk), e3w_abl(jk), jk = 1, jpka )
	160	END IF
	161
	162	!!---------------------------------------------------------------------
[12199]	163	!! Check TKE closure parameters
	164	!!---------------------------------------------------------------------
[11305]	165	rn_Sch = rn_ce / rn_cm
	166	mxl_min = (avm_bak / rn_cm) / sqrt( tke_min )
	167
	168	IF(lwp) THEN
	169	WRITE(numout,*)
	170	WRITE(numout,*) ' abl_zdf_tke : ABL TKE turbulent closure'
[12199]	171	WRITE(numout,*) ' ~~~~~~~~~~~'
[11305]	172	IF(nn_amxl==0) WRITE(numout,*) 'Deardorff 80 length-scale '
	173	IF(nn_amxl==1) WRITE(numout,*) 'length-scale based on the distance to the PBL height '
	174	WRITE(numout,*) ' Minimum value of atmospheric TKE = ',tke_min,' m^2 s^-2'
	175	WRITE(numout,*) ' Minimum value of atmospheric mixing length = ',mxl_min,' m'
[12199]	176	WRITE(numout,*) ' Constant for turbulent viscosity = ',rn_Cm
	177	WRITE(numout,*) ' Constant for turbulent diffusivity = ',rn_Ct
	178	WRITE(numout,*) ' Constant for Schmidt number = ',rn_Sch
	179	WRITE(numout,*) ' Constant for TKE dissipation = ',rn_Ceps
[11305]	180	END IF
	181
	182	!!-------------------------------------------------------------------------------------------
	183	!! Compute parameters to build the vertical profile for the nudging term (used in abl_stp())
[12199]	184	!!-------------------------------------------------------------------------------------------
[11305]	185	zcff1 = 1._wp / ( jp_bmax - jp_bmin )**3
	186	! for active tracers
	187	jp_alp3_tra = -2._wp * zcff1 * ( rn_ltra_max - rn_ltra_min )
	188	jp_alp2_tra = 3._wp * zcff1 * (jp_bmax + jp_bmin) * ( rn_ltra_max - rn_ltra_min )
	189	jp_alp1_tra = -6._wp * zcff1 * jp_bmax * jp_bmin * ( rn_ltra_max - rn_ltra_min )
	190	jp_alp0_tra = zcff1 * ( rn_ltra_max * jp_bminjp_bmin (3._wp*jp_bmax - jp_bmin) &
[12199]	191	& - rn_ltra_min * jp_bmaxjp_bmax (3._wp*jp_bmin - jp_bmax) )
[11305]	192	! for dynamics
	193	jp_alp3_dyn = -2._wp * zcff1 * ( rn_ldyn_max - rn_ldyn_min )
	194	jp_alp2_dyn = 3._wp * zcff1 * (jp_bmax + jp_bmin) * ( rn_ldyn_max - rn_ldyn_min )
	195	jp_alp1_dyn = -6._wp * zcff1 * jp_bmax * jp_bmin * ( rn_ldyn_max - rn_ldyn_min )
	196	jp_alp0_dyn = zcff1 * ( rn_ldyn_max * jp_bminjp_bmin (3._wp*jp_bmax - jp_bmin) &
[12199]	197	& - rn_ldyn_min * jp_bmaxjp_bmax (3._wp*jp_bmin - jp_bmax) )
[11305]	198
	199	jp_pblh_min = ghw_abl( 4) / jp_bmin !<-- at least 3 grid points at the bottom have value rn_ltra_min
	200	jp_pblh_max = ghw_abl(jpka-3) / jp_bmax !<-- at least 3 grid points at the top have value rn_ltra_max
	201
[11363]	202	! ABL timestep
[12489]	203	rDt_abl = nn_fsbc * rn_Dt
[11363]	204
[11305]	205	! Check parameters for dynamics
[11363]	206	zcff = ( jp_alp3_dyn * jp_bmin*3 + jp_alp2_dyn jp_bmin**2 &
[12489]	207	& + jp_alp1_dyn * jp_bmin + jp_alp0_dyn ) * rDt_abl
[11363]	208	zcff1 = ( jp_alp3_dyn * jp_bmax*3 + jp_alp2_dyn jp_bmax**2 &
[12489]	209	& + jp_alp1_dyn * jp_bmax + jp_alp0_dyn ) * rDt_abl
[11305]	210	IF(lwp) THEN
	211	IF(nn_dyn_restore > 0) THEN
[11363]	212	WRITE(numout,*) ' ABL Minimum value for dynamics restoring = ',zcff
	213	WRITE(numout,*) ' ABL Maximum value for dynamics restoring = ',zcff1
[11305]	214	! Check that restoring coefficients are between 0 and 1
[11586]	215	IF( zcff1 - nn_fsbc > 0.001_wp .OR. zcff1 < 0._wp ) &
[11305]	216	& CALL ctl_stop( 'abl_init : wrong value for rn_ldyn_max' )
[11586]	217	IF( zcff - nn_fsbc > 0.001_wp .OR. zcff < 0._wp ) &
[11305]	218	& CALL ctl_stop( 'abl_init : wrong value for rn_ldyn_min' )
[11363]	219	IF( zcff > zcff1 ) &
[11586]	220	& CALL ctl_stop( 'abl_init : rn_ldyn_max must be smaller than rn_ldyn_min' )
[11305]	221	END IF
	222	END IF
	223
	224	! Check parameters for active tracers
[11363]	225	zcff = ( jp_alp3_tra * jp_bmin*3 + jp_alp2_tra jp_bmin**2 &
[12489]	226	& + jp_alp1_tra * jp_bmin + jp_alp0_tra ) * rDt_abl
[11363]	227	zcff1 = ( jp_alp3_tra * jp_bmax*3 + jp_alp2_tra jp_bmax**2 &
[12489]	228	& + jp_alp1_tra * jp_bmax + jp_alp0_tra ) * rDt_abl
[11305]	229	IF(lwp) THEN
[11363]	230	WRITE(numout,*) ' ABL Minimum value for tracers restoring = ',zcff
	231	WRITE(numout,*) ' ABL Maximum value for tracers restoring = ',zcff1
[11305]	232	! Check that restoring coefficients are between 0 and 1
[11586]	233	IF( zcff1 - nn_fsbc > 0.001_wp .OR. zcff1 < 0._wp ) &
[11305]	234	& CALL ctl_stop( 'abl_init : wrong value for rn_ltra_max' )
[11586]	235	IF( zcff - nn_fsbc > 0.001_wp .OR. zcff < 0._wp ) &
[11305]	236	& CALL ctl_stop( 'abl_init : wrong value for rn_ltra_min' )
[11363]	237	IF( zcff > zcff1 ) &
[11586]	238	& CALL ctl_stop( 'abl_init : rn_ltra_max must be smaller than rn_ltra_min' )
[11305]	239	END IF
	240
	241	!!-------------------------------------------------------------------------------------------
	242	!! Initialize Coriolis frequency, equatorial restoring and land/sea mask
[12199]	243	!!-------------------------------------------------------------------------------------------
[11305]	244	fft_abl(:,:) = 2._wp * omega * SIN( rad * gphit(:,:) )
	245
	246	! Equatorial restoring
	247	IF( nn_dyn_restore == 1 ) THEN
[12199]	248	zcff = 2._wp * omega * SIN( rad * 90._wp ) !++ fmax
[11363]	249	rest_eq(:,:) = SIN( 0.5_wprpi( (fft_abl(:,:) - zcff) / zcff ) )**8
	250	!!GS: alternative shape
	251	!rest_eq(:,:) = SIN( 0.5_wprpi(zcff - ABS(ff_t(:,:))) / (zcff - 3.e-5) )**8
	252	!WHERE(ABS(ff_t(:,:)).LE.3.e-5) rest_eq(:,:) = 1._wp
[11305]	253	ELSE
	254	rest_eq(:,:) = 1._wp
	255	END IF
	256	! T-mask
	257	msk_abl(:,:) = tmask(:,:,1)
	258
[12199]	259	!!-------------------------------------------------------------------------------------------
[11305]	260
	261	! initialize 2D bulk fields AND 3D abl data
	262	CALL sbc_blk_init
	263
	264	! Initialize the time index for now time (nt_n) and after time (nt_a)
	265	nt_n = 1 + MOD( nit000 , 2)
	266	nt_a = 1 + MOD( nit000+1, 2)
	267
[11858]	268	! initialize ABL from data or restart
	269	IF( ln_rstart ) THEN
	270	CALL abl_rst_read
	271	ELSE
	272	CALL fld_read( nit000, nn_fsbc, sf ) ! input fields provided at the first time-step
[11305]	273
[12199]	274	u_abl(:,:,:,nt_n ) = sf(jp_wndi)%fnow(:,:,:)
	275	v_abl(:,:,:,nt_n ) = sf(jp_wndj)%fnow(:,:,:)
[11858]	276	tq_abl(:,:,:,nt_n,jp_ta) = sf(jp_tair)%fnow(:,:,:)
	277	tq_abl(:,:,:,nt_n,jp_qa) = sf(jp_humi)%fnow(:,:,:)
[12199]	278
[11858]	279	tke_abl(:,:,:,nt_n ) = tke_min
	280	avm_abl(:,:,: ) = avm_bak
	281	avt_abl(:,:,: ) = avt_bak
[12199]	282	mxl_abl(:,:,: ) = mxl_min
	283	pblh (:,: ) = ghw_abl( 3 ) !<-- assume that the pbl contains 3 grid points
[11858]	284	u_abl (:,:,:,nt_a ) = 0._wp
	285	v_abl (:,:,:,nt_a ) = 0._wp
	286	tq_abl (:,:,:,nt_a,: ) = 0._wp
	287	tke_abl(:,:,:,nt_a ) = 0._wp
	288	ENDIF
[12021]	289
	290	rhoa(:,:) = rho_air( tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), sf(jp_slp)%fnow(:,:,1) ) !!GS: rhoa must be (re)computed here here to avoid division by zero in blk_ice_1 (TBI)
[12199]	291
[11305]	292	END SUBROUTINE sbc_abl_init
	293
	294
	295	SUBROUTINE sbc_abl( kt )
	296	!!---------------------------------------------------------------------
	297	!! * ROUTINE sbc_abl *
	298	!!
	299	!! ** Purpose : provide the momentum, heat and freshwater fluxes at
[12199]	300	!! the ocean surface from an ABL calculation at each oceanic time step
[11305]	301	!!
[12199]	302	!! ** Method :
[11305]	303	!! - Pre-compute part of turbulent fluxes in blk_oce_1
[12199]	304	!! - Perform 1 time-step of the ABL model
[11305]	305	!! - Finalize flux computation in blk_oce_2
	306	!!
	307	!! ** Outputs : - utau : i-component of the stress at U-point (N/m2)
	308	!! - vtau : j-component of the stress at V-point (N/m2)
	309	!! - taum : Wind stress module at T-point (N/m2)
	310	!! - wndm : Wind speed module at T-point (m/s)
	311	!! - qsr : Solar heat flux over the ocean (W/m2)
	312	!! - qns : Non Solar heat flux over the ocean (W/m2)
	313	!! - emp : evaporation minus precipitation (kg/m2/s)
	314	!!
	315	!!---------------------------------------------------------------------
	316	INTEGER , INTENT(in) :: kt ! ocean time step
	317	!!
	318	REAL(wp), DIMENSION(jpi,jpj) :: zssq, zcd_du, zsen, zevp
[11334]	319	#if defined key_si3
	320	REAL(wp), DIMENSION(jpi,jpj) :: zssqi, zcd_dui, zseni, zevpi
[12199]	321	#endif
[11305]	322	INTEGER :: jbak, jbak_dta, ji, jj
	323	!!---------------------------------------------------------------------
	324	!
	325	!!-------------------------------------------------------------------------------------------
	326	!! 1 - Read Atmospheric 3D data for large-scale forcing
	327	!!-------------------------------------------------------------------------------------------
[12199]	328
[11305]	329	CALL fld_read( kt, nn_fsbc, sf ) ! input fields provided at the current time-step
[12199]	330
[11305]	331	!!-------------------------------------------------------------------------------------------
	332	!! 2 - Compute Cd x \|\|U\|\|, Ch x \|\|U\|\|, Ce x \|\|U\|\|, and SSQ using now fields
[12199]	333	!!-------------------------------------------------------------------------------------------
[11334]	334
[11348]	335	CALL blk_oce_1( kt, u_abl(:,:,2,nt_n ), v_abl(:,:,2,nt_n ), & ! <<= in
	336	& tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), & ! <<= in
	337	& sf(jp_slp )%fnow(:,:,1) , sst_m, ssu_m, ssv_m , & ! <<= in
[12015]	338	& sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1) , & ! <<= in
[12199]	339	& tsk_m, zssq, zcd_du, zsen, zevp ) ! =>> out
	340
[11334]	341	#if defined key_si3
[11360]	342	CALL blk_ice_1( u_abl(:,:,2,nt_n ), v_abl(:,:,2,nt_n ), & ! <<= in
	343	& tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), & ! <<= in
	344	& sf(jp_slp)%fnow(:,:,1) , u_ice, v_ice, tm_su , & ! <<= in
	345	& pseni=zseni, pevpi=zevpi, pssqi=zssqi, pcd_dui=zcd_dui ) ! <<= out
[12199]	346	#endif
[11334]	347
[11305]	348	!!-------------------------------------------------------------------------------------------
	349	!! 3 - Advance ABL variables from now (n) to after (n+1)
[12199]	350	!!-------------------------------------------------------------------------------------------
	351
	352	CALL abl_stp( kt, tsk_m, ssu_m, ssv_m, zssq, & ! <<= in
[11348]	353	& sf(jp_wndi)%fnow(:,:,:), sf(jp_wndj)%fnow(:,:,:), & ! <<= in
	354	& sf(jp_tair)%fnow(:,:,:), sf(jp_humi)%fnow(:,:,:), & ! <<= in
	355	& sf(jp_slp )%fnow(:,:,1), & ! <<= in
	356	& sf(jp_hpgi)%fnow(:,:,:), sf(jp_hpgj)%fnow(:,:,:), & ! <<= in
	357	& zcd_du, zsen, zevp, & ! <=> in/out
	358	& wndm, utau, vtau, taum & ! =>> out
[12199]	359	#if defined key_si3
[11348]	360	& , tm_su, u_ice, v_ice, zssqi, zcd_dui & ! <<= in
	361	& , zseni, zevpi, wndm_ice, ato_i & ! <<= in
	362	& , utau_ice, vtau_ice & ! =>> out
[12199]	363	#endif
[11348]	364	& )
[11305]	365	!!-------------------------------------------------------------------------------------------
[12199]	366	!! 4 - Finalize flux computation using ABL variables at (n+1), nt_n corresponds to (n+1) since
[11305]	367	!! time swap is done in abl_stp
[12199]	368	!!-------------------------------------------------------------------------------------------
[11305]	369
[11360]	370	CALL blk_oce_2( tq_abl(:,:,2,nt_n,jp_ta), &
	371	& sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1), &
	372	& sf(jp_prec)%fnow(:,:,1) , sf(jp_snow)%fnow(:,:,1), &
[12199]	373	& tsk_m, zsen, zevp )
[11305]	374
[12199]	375	CALL abl_rst_opn( kt ) ! Open abl restart file (if necessary)
	376	IF( lrst_abl ) CALL abl_rst_write( kt ) ! -- abl restart file
[11858]	377
[11334]	378	#if defined key_si3
[11348]	379	! Avoid a USE abl in icesbc module
[11334]	380	sf(jp_tair)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_ta); sf(jp_humi)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_qa)
[12199]	381	#endif
[11334]	382
[11305]	383	END SUBROUTINE sbc_abl
	384
	385	!!======================================================================
	386	END MODULE sbcabl

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

Download in other formats: