Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

sbcssm.F90 @ 12286

Last change on this file since 12286 was 12286, checked in by jcastill, 4 years ago
Remove svn keywords
File size: 13.2 KB

Rev	Line
[888]	1	MODULE sbcssm
	2	!!======================================================================
	3	!! * MODULE sbcssm *
	4	!! Surface module : provide time-mean ocean surface variables
	5	!!======================================================================
[2715]	6	!! History : 9.0 ! 2006-07 (G. Madec) Original code
[2528]	7	!! 3.3 ! 2010-10 (C. Bricaud, G. Madec) add the Patm forcing for sea-ice
[6140]	8	!! 3.7 ! 2015-11 (G. Madec) non linear free surface by default: e3t_m always computed
[888]	9	!!----------------------------------------------------------------------
	10
	11	!!----------------------------------------------------------------------
[6140]	12	!! sbc_ssm : calculate sea surface mean currents, temperature,
	13	!! and salinity over nn_fsbc time-step
[888]	14	!!----------------------------------------------------------------------
[6140]	15	USE oce ! ocean dynamics and tracers
	16	USE dom_oce ! ocean space and time domain
	17	USE sbc_oce ! surface boundary condition: ocean fields
	18	USE sbcapr ! surface boundary condition: atmospheric pressure
	19	USE eosbn2 ! equation of state and related derivatives
[4990]	20	!
[6140]	21	USE in_out_manager ! I/O manager
	22	USE prtctl ! Print control
	23	USE iom ! IOM library
[888]	24
	25	IMPLICIT NONE
	26	PRIVATE
	27
[6140]	28	PUBLIC sbc_ssm ! routine called by step.F90
	29	PUBLIC sbc_ssm_init ! routine called by sbcmod.F90
[4292]	30
[6140]	31	LOGICAL, SAVE :: l_ssm_mean = .FALSE. ! keep track of whether means have been read from restart file
[4292]	32
[888]	33	!!----------------------------------------------------------------------
[2528]	34	!! NEMO/OPA 3.3 , NEMO Consortium (2010)
[1156]	35	!! $Id$
[2715]	36	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
[888]	37	!!----------------------------------------------------------------------
	38	CONTAINS
	39
	40	SUBROUTINE sbc_ssm( kt )
	41	!!---------------------------------------------------------------------
	42	!! * ROUTINE sbc_oce *
	43	!!
	44	!! ** Purpose : provide ocean surface variable to sea-surface boundary
	45	!! condition computation
	46	!!
	47	!! ** Method : compute mean surface velocity (2 components at U and
[7646]	48	!! V-points) [m/s], temperature [Celsius] and salinity [psu] over
[888]	49	!! the periode (kt - nn_fsbc) to kt
[2528]	50	!! Note that the inverse barometer ssh (i.e. ssh associated with Patm)
	51	!! is add to ssh_m when ln_apr_dyn = T. Required for sea-ice dynamics.
[888]	52	!!---------------------------------------------------------------------
[2715]	53	INTEGER, INTENT(in) :: kt ! ocean time step
[888]	54	!
[4990]	55	INTEGER :: ji, jj ! loop index
[2715]	56	REAL(wp) :: zcoef, zf_sbc ! local scalar
[4990]	57	REAL(wp), DIMENSION(jpi,jpj,jpts) :: zts
[888]	58	!!---------------------------------------------------------------------
[6140]	59	!
[5120]	60	! !* surface T-, U-, V- ocean level variables (T, S, depth, velocity)
[4990]	61	DO jj = 1, jpj
	62	DO ji = 1, jpi
	63	zts(ji,jj,jp_tem) = tsn(ji,jj,mikt(ji,jj),jp_tem)
	64	zts(ji,jj,jp_sal) = tsn(ji,jj,mikt(ji,jj),jp_sal)
	65	END DO
	66	END DO
	67	!
[3607]	68	IF( nn_fsbc == 1 ) THEN ! Instantaneous surface fields !
[888]	69	! ! ---------------------------------------- !
[7753]	70	ssu_m(:,:) = ub(:,:,1)
	71	ssv_m(:,:) = vb(:,:,1)
	72	IF( l_useCT ) THEN ; sst_m(:,:) = eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) )
	73	ELSE ; sst_m(:,:) = zts(:,:,jp_tem)
[4990]	74	ENDIF
[7753]	75	sss_m(:,:) = zts(:,:,jp_sal)
[2528]	76	! ! removed inverse barometer ssh when Patm forcing is used (for sea-ice dynamics)
[7753]	77	IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) )
	78	ELSE ; ssh_m(:,:) = sshn(:,:)
[2528]	79	ENDIF
[888]	80	!
[7753]	81	e3t_m(:,:) = e3t_n(:,:,1)
[4292]	82	!
[7753]	83	frq_m(:,:) = fraqsr_1lev(:,:)
[5407]	84	!
[888]	85	ELSE
[3607]	86	! ! ----------------------------------------------- !
[3614]	87	IF( kt == nit000 .AND. .NOT. l_ssm_mean ) THEN ! Initialisation: 1st time-step, no input means !
[3607]	88	! ! ----------------------------------------------- !
[888]	89	IF(lwp) WRITE(numout,*)
[7646]	90	IF(lwp) WRITE(numout,*) 'sbc_ssm : mean fields initialised to instantaneous values'
	91	IF(lwp) WRITE(numout,*) '~~~~~~~ '
[3607]	92	zcoef = REAL( nn_fsbc - 1, wp )
[7753]	93	ssu_m(:,:) = zcoef * ub(:,:,1)
	94	ssv_m(:,:) = zcoef * vb(:,:,1)
	95	IF( l_useCT ) THEN ; sst_m(:,:) = zcoef * eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) )
	96	ELSE ; sst_m(:,:) = zcoef * zts(:,:,jp_tem)
[4990]	97	ENDIF
[7753]	98	sss_m(:,:) = zcoef * zts(:,:,jp_sal)
[4990]	99	! ! removed inverse barometer ssh when Patm forcing is used (for sea-ice dynamics)
[7753]	100	IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = zcoef * ( sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) )
	101	ELSE ; ssh_m(:,:) = zcoef * sshn(:,:)
[888]	102	ENDIF
[4990]	103	!
[7753]	104	e3t_m(:,:) = zcoef * e3t_n(:,:,1)
	105	!
	106	frq_m(:,:) = zcoef * fraqsr_1lev(:,:)
[888]	107	! ! ---------------------------------------- !
	108	ELSEIF( MOD( kt - 2 , nn_fsbc ) == 0 ) THEN ! Initialisation: New mean computation !
	109	! ! ---------------------------------------- !
[7753]	110	ssu_m(:,:) = 0._wp ! reset to zero ocean mean sbc fields
	111	ssv_m(:,:) = 0._wp
	112	sst_m(:,:) = 0._wp
	113	sss_m(:,:) = 0._wp
	114	ssh_m(:,:) = 0._wp
	115	e3t_m(:,:) = 0._wp
	116	frq_m(:,:) = 0._wp
[888]	117	ENDIF
	118	! ! ---------------------------------------- !
	119	! ! Cumulate at each time step !
	120	! ! ---------------------------------------- !
[7753]	121	ssu_m(:,:) = ssu_m(:,:) + ub(:,:,1)
	122	ssv_m(:,:) = ssv_m(:,:) + vb(:,:,1)
	123	IF( l_useCT ) THEN ; sst_m(:,:) = sst_m(:,:) + eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) )
	124	ELSE ; sst_m(:,:) = sst_m(:,:) + zts(:,:,jp_tem)
[4990]	125	ENDIF
[7753]	126	sss_m(:,:) = sss_m(:,:) + zts(:,:,jp_sal)
[2528]	127	! ! removed inverse barometer ssh when Patm forcing is used (for sea-ice dynamics)
[7753]	128	IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = ssh_m(:,:) + sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) )
	129	ELSE ; ssh_m(:,:) = ssh_m(:,:) + sshn(:,:)
[2528]	130	ENDIF
[4990]	131	!
[7753]	132	e3t_m(:,:) = e3t_m(:,:) + e3t_n(:,:,1)
	133	!
	134	frq_m(:,:) = frq_m(:,:) + fraqsr_1lev(:,:)
[888]	135
	136	! ! ---------------------------------------- !
	137	IF( MOD( kt - 1 , nn_fsbc ) == 0 ) THEN ! Mean value at each nn_fsbc time-step !
	138	! ! ---------------------------------------- !
	139	zcoef = 1. / REAL( nn_fsbc, wp )
[7753]	140	sst_m(:,:) = sst_m(:,:) * zcoef ! mean SST [Celsius]
	141	sss_m(:,:) = sss_m(:,:) * zcoef ! mean SSS [psu]
	142	ssu_m(:,:) = ssu_m(:,:) * zcoef ! mean suface current [m/s]
	143	ssv_m(:,:) = ssv_m(:,:) * zcoef !
	144	ssh_m(:,:) = ssh_m(:,:) * zcoef ! mean SSH [m]
	145	e3t_m(:,:) = e3t_m(:,:) * zcoef ! mean vertical scale factor [m]
	146	frq_m(:,:) = frq_m(:,:) * zcoef ! mean fraction of solar net radiation absorbed in the 1st T level [-]
[888]	147	!
	148	ENDIF
	149	! ! ---------------------------------------- !
	150	IF( lrst_oce ) THEN ! Write in the ocean restart file !
	151	! ! ---------------------------------------- !
	152	IF(lwp) WRITE(numout,*)
	153	IF(lwp) WRITE(numout,*) 'sbc_ssm : sea surface mean fields written in ocean restart file ', &
	154	& 'at it= ', kt,' date= ', ndastp
	155	IF(lwp) WRITE(numout,*) '~~~~~~~'
	156	zf_sbc = REAL( nn_fsbc, wp )
	157	CALL iom_rstput( kt, nitrst, numrow, 'nn_fsbc', zf_sbc ) ! sbc frequency
	158	CALL iom_rstput( kt, nitrst, numrow, 'ssu_m' , ssu_m ) ! sea surface mean fields
	159	CALL iom_rstput( kt, nitrst, numrow, 'ssv_m' , ssv_m )
	160	CALL iom_rstput( kt, nitrst, numrow, 'sst_m' , sst_m )
	161	CALL iom_rstput( kt, nitrst, numrow, 'sss_m' , sss_m )
[1244]	162	CALL iom_rstput( kt, nitrst, numrow, 'ssh_m' , ssh_m )
[6140]	163	CALL iom_rstput( kt, nitrst, numrow, 'e3t_m' , e3t_m )
[5407]	164	CALL iom_rstput( kt, nitrst, numrow, 'frq_m' , frq_m )
[888]	165	!
	166	ENDIF
	167	!
	168	ENDIF
	169	!
[5407]	170	IF( MOD( kt - 1 , nn_fsbc ) == 0 ) THEN ! Mean value at each nn_fsbc time-step !
	171	CALL iom_put( 'ssu_m', ssu_m )
	172	CALL iom_put( 'ssv_m', ssv_m )
	173	CALL iom_put( 'sst_m', sst_m )
	174	CALL iom_put( 'sss_m', sss_m )
	175	CALL iom_put( 'ssh_m', ssh_m )
[6140]	176	CALL iom_put( 'e3t_m', e3t_m )
[5407]	177	CALL iom_put( 'frq_m', frq_m )
	178	ENDIF
	179	!
[888]	180	END SUBROUTINE sbc_ssm
	181
[6140]	182
[3607]	183	SUBROUTINE sbc_ssm_init
	184	!!----------------------------------------------------------------------
	185	!! * ROUTINE sbc_ssm_init *
	186	!!
	187	!! ** Purpose : Initialisation of the sbc data
	188	!!
	189	!! ** Action : - read parameters
	190	!!----------------------------------------------------------------------
[6140]	191	REAL(wp) :: zcoef, zf_sbc ! local scalar
[3607]	192	!!----------------------------------------------------------------------
[6140]	193	!
[3607]	194	IF( nn_fsbc == 1 ) THEN
	195	!
	196	IF(lwp) WRITE(numout,*)
[7646]	197	IF(lwp) WRITE(numout,*) 'sbc_ssm_init : sea surface mean fields, nn_fsbc=1 : instantaneous values'
	198	IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ '
[3607]	199	!
	200	ELSE
	201	!
	202	IF(lwp) WRITE(numout,*)
[7646]	203	IF(lwp) WRITE(numout,*) 'sbc_ssm_init : sea surface mean fields'
	204	IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~ '
[3607]	205	!
	206	IF( ln_rstart .AND. iom_varid( numror, 'nn_fsbc', ldstop = .FALSE. ) > 0 ) THEN
	207	l_ssm_mean = .TRUE.
[6140]	208	CALL iom_get( numror , 'nn_fsbc', zf_sbc ) ! sbc frequency of previous run
	209	CALL iom_get( numror, jpdom_autoglo, 'ssu_m' , ssu_m ) ! sea surface mean velocity (U-point)
	210	CALL iom_get( numror, jpdom_autoglo, 'ssv_m' , ssv_m ) ! " " velocity (V-point)
	211	CALL iom_get( numror, jpdom_autoglo, 'sst_m' , sst_m ) ! " " temperature (T-point)
	212	CALL iom_get( numror, jpdom_autoglo, 'sss_m' , sss_m ) ! " " salinity (T-point)
	213	CALL iom_get( numror, jpdom_autoglo, 'ssh_m' , ssh_m ) ! " " height (T-point)
	214	CALL iom_get( numror, jpdom_autoglo, 'e3t_m' , e3t_m ) ! 1st level thickness (T-point)
[5407]	215	! fraction of solar net radiation absorbed in 1st T level
	216	IF( iom_varid( numror, 'frq_m', ldstop = .FALSE. ) > 0 ) THEN
	217	CALL iom_get( numror, jpdom_autoglo, 'frq_m' , frq_m )
	218	ELSE
[7753]	219	frq_m(:,:) = 1._wp ! default definition
[5407]	220	ENDIF
[3607]	221	!
	222	IF( zf_sbc /= REAL( nn_fsbc, wp ) ) THEN ! nn_fsbc has changed between 2 runs
[7646]	223	IF(lwp) WRITE(numout,*) ' restart with a change in the frequency of mean from ', zf_sbc, ' to ', nn_fsbc
[3607]	224	zcoef = REAL( nn_fsbc - 1, wp ) / zf_sbc
[7753]	225	ssu_m(:,:) = zcoef * ssu_m(:,:)
	226	ssv_m(:,:) = zcoef * ssv_m(:,:)
	227	sst_m(:,:) = zcoef * sst_m(:,:)
	228	sss_m(:,:) = zcoef * sss_m(:,:)
	229	ssh_m(:,:) = zcoef * ssh_m(:,:)
	230	e3t_m(:,:) = zcoef * e3t_m(:,:)
	231	frq_m(:,:) = zcoef * frq_m(:,:)
[3607]	232	ELSE
[7646]	233	IF(lwp) WRITE(numout,*) ' mean fields read in the ocean restart file'
[3607]	234	ENDIF
	235	ENDIF
	236	ENDIF
	237	!
[5407]	238	IF( .NOT. l_ssm_mean ) THEN ! default initialisation. needed by lim_istate
	239	!
[7646]	240	IF(lwp) WRITE(numout,*) ' default initialisation of ss._m arrays'
[7753]	241	ssu_m(:,:) = ub(:,:,1)
	242	ssv_m(:,:) = vb(:,:,1)
[6489]	243	IF( l_useCT ) THEN ; sst_m(:,:) = eos_pt_from_ct( tsn(:,:,1,jp_tem), tsn(:,:,1,jp_sal) )
[7646]	244	ELSE ; sst_m(:,:) = tsn(:,:,1,jp_tem)
[5407]	245	ENDIF
[7753]	246	sss_m(:,:) = tsn (:,:,1,jp_sal)
	247	ssh_m(:,:) = sshn (:,:)
	248	e3t_m(:,:) = e3t_n(:,:,1)
	249	frq_m(:,:) = 1._wp
[5407]	250	!
	251	ENDIF
	252	!
[3607]	253	END SUBROUTINE sbc_ssm_init
	254
[888]	255	!!======================================================================
	256	END MODULE sbcssm

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

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

Context Navigation

source: branches/UKMO/r8395_cpl_tauwav/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90 @ 12286

Download in other formats: