Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

sbcssm.F90 @ 9580

Last change on this file since 9580 was 9570, checked in by nicolasmartin, 6 years ago

Global renaming for core routines (./NEMO)

Folders
- LIM_SRC_3 -> ICE_SRC
- OPA_SRC -> OCE_SRC
CPP key: key_lim3 -> key_si3
Modules, (sub)routines and variables names
- MPI: mpi_comm_opa -> mpi_comm_oce, MPI_COMM_OPA -> MPI_COMM_OCE, mpi_init_opa -> mpi_init_oce
- AGRIF: agrif_opa_* -> agrif_oce_*, agrif_lim3_* -> agrif_si3_* and few more
- TOP-PISCES: p.zlim -> p.zice, namp.zlim -> namp.zice
Comments
- NEMO/OPA -> NEMO/OCE
- ESIM|LIM3 -> SI3

Property svn:keywords set to Id

File size: 14.0 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	!!----------------------------------------------------------------------
[9570]	34	!! NEMO/OCE 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)
[9019]	123	IF( l_useCT ) THEN ; sst_m(:,:) = sst_m(:,:) + eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) )
[7753]	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 )
[9367]	157	IF( lwxios ) CALL iom_swap( cwxios_context )
	158	CALL iom_rstput( kt, nitrst, numrow, 'nn_fsbc', zf_sbc, ldxios = lwxios ) ! sbc frequency
	159	CALL iom_rstput( kt, nitrst, numrow, 'ssu_m' , ssu_m, ldxios = lwxios ) ! sea surface mean fields
	160	CALL iom_rstput( kt, nitrst, numrow, 'ssv_m' , ssv_m, ldxios = lwxios )
	161	CALL iom_rstput( kt, nitrst, numrow, 'sst_m' , sst_m, ldxios = lwxios )
	162	CALL iom_rstput( kt, nitrst, numrow, 'sss_m' , sss_m, ldxios = lwxios )
	163	CALL iom_rstput( kt, nitrst, numrow, 'ssh_m' , ssh_m, ldxios = lwxios )
	164	CALL iom_rstput( kt, nitrst, numrow, 'e3t_m' , e3t_m, ldxios = lwxios )
	165	CALL iom_rstput( kt, nitrst, numrow, 'frq_m' , frq_m, ldxios = lwxios )
[888]	166	!
[9367]	167	IF( lwxios ) CALL iom_swap( cxios_context )
[888]	168	ENDIF
	169	!
	170	ENDIF
	171	!
[5407]	172	IF( MOD( kt - 1 , nn_fsbc ) == 0 ) THEN ! Mean value at each nn_fsbc time-step !
	173	CALL iom_put( 'ssu_m', ssu_m )
	174	CALL iom_put( 'ssv_m', ssv_m )
	175	CALL iom_put( 'sst_m', sst_m )
	176	CALL iom_put( 'sss_m', sss_m )
	177	CALL iom_put( 'ssh_m', ssh_m )
[6140]	178	CALL iom_put( 'e3t_m', e3t_m )
[5407]	179	CALL iom_put( 'frq_m', frq_m )
	180	ENDIF
	181	!
[888]	182	END SUBROUTINE sbc_ssm
	183
[6140]	184
[3607]	185	SUBROUTINE sbc_ssm_init
	186	!!----------------------------------------------------------------------
	187	!! * ROUTINE sbc_ssm_init *
	188	!!
	189	!! ** Purpose : Initialisation of the sbc data
	190	!!
	191	!! ** Action : - read parameters
	192	!!----------------------------------------------------------------------
[6140]	193	REAL(wp) :: zcoef, zf_sbc ! local scalar
[3607]	194	!!----------------------------------------------------------------------
[6140]	195	!
[3607]	196	IF( nn_fsbc == 1 ) THEN
	197	!
	198	IF(lwp) WRITE(numout,*)
[7646]	199	IF(lwp) WRITE(numout,*) 'sbc_ssm_init : sea surface mean fields, nn_fsbc=1 : instantaneous values'
	200	IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ '
[3607]	201	!
	202	ELSE
	203	!
	204	IF(lwp) WRITE(numout,*)
[7646]	205	IF(lwp) WRITE(numout,*) 'sbc_ssm_init : sea surface mean fields'
	206	IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~ '
[3607]	207	!
	208	IF( ln_rstart .AND. iom_varid( numror, 'nn_fsbc', ldstop = .FALSE. ) > 0 ) THEN
	209	l_ssm_mean = .TRUE.
[9367]	210	CALL iom_get( numror , 'nn_fsbc', zf_sbc, ldxios = lrxios ) ! sbc frequency of previous run
	211	CALL iom_get( numror, jpdom_autoglo, 'ssu_m' , ssu_m, ldxios = lrxios ) ! sea surface mean velocity (U-point)
	212	CALL iom_get( numror, jpdom_autoglo, 'ssv_m' , ssv_m, ldxios = lrxios ) ! " " velocity (V-point)
	213	CALL iom_get( numror, jpdom_autoglo, 'sst_m' , sst_m, ldxios = lrxios ) ! " " temperature (T-point)
	214	CALL iom_get( numror, jpdom_autoglo, 'sss_m' , sss_m, ldxios = lrxios ) ! " " salinity (T-point)
	215	CALL iom_get( numror, jpdom_autoglo, 'ssh_m' , ssh_m, ldxios = lrxios ) ! " " height (T-point)
	216	CALL iom_get( numror, jpdom_autoglo, 'e3t_m' , e3t_m, ldxios = lrxios ) ! 1st level thickness (T-point)
[5407]	217	! fraction of solar net radiation absorbed in 1st T level
	218	IF( iom_varid( numror, 'frq_m', ldstop = .FALSE. ) > 0 ) THEN
[9367]	219	CALL iom_get( numror, jpdom_autoglo, 'frq_m' , frq_m, ldxios = lrxios )
[5407]	220	ELSE
[7753]	221	frq_m(:,:) = 1._wp ! default definition
[5407]	222	ENDIF
[3607]	223	!
	224	IF( zf_sbc /= REAL( nn_fsbc, wp ) ) THEN ! nn_fsbc has changed between 2 runs
[7646]	225	IF(lwp) WRITE(numout,*) ' restart with a change in the frequency of mean from ', zf_sbc, ' to ', nn_fsbc
[3607]	226	zcoef = REAL( nn_fsbc - 1, wp ) / zf_sbc
[7753]	227	ssu_m(:,:) = zcoef * ssu_m(:,:)
	228	ssv_m(:,:) = zcoef * ssv_m(:,:)
	229	sst_m(:,:) = zcoef * sst_m(:,:)
	230	sss_m(:,:) = zcoef * sss_m(:,:)
	231	ssh_m(:,:) = zcoef * ssh_m(:,:)
	232	e3t_m(:,:) = zcoef * e3t_m(:,:)
	233	frq_m(:,:) = zcoef * frq_m(:,:)
[3607]	234	ELSE
[7646]	235	IF(lwp) WRITE(numout,*) ' mean fields read in the ocean restart file'
[3607]	236	ENDIF
	237	ENDIF
	238	ENDIF
	239	!
[9019]	240	IF( .NOT.l_ssm_mean ) THEN ! default initialisation. needed by iceistate
[5407]	241	!
[7646]	242	IF(lwp) WRITE(numout,*) ' default initialisation of ss._m arrays'
[7753]	243	ssu_m(:,:) = ub(:,:,1)
	244	ssv_m(:,:) = vb(:,:,1)
[6489]	245	IF( l_useCT ) THEN ; sst_m(:,:) = eos_pt_from_ct( tsn(:,:,1,jp_tem), tsn(:,:,1,jp_sal) )
[7646]	246	ELSE ; sst_m(:,:) = tsn(:,:,1,jp_tem)
[5407]	247	ENDIF
[7753]	248	sss_m(:,:) = tsn (:,:,1,jp_sal)
	249	ssh_m(:,:) = sshn (:,:)
	250	e3t_m(:,:) = e3t_n(:,:,1)
	251	frq_m(:,:) = 1._wp
[5407]	252	!
	253	ENDIF
	254	!
[9367]	255	IF( lwxios.AND.nn_fsbc > 1 ) THEN
	256	CALL iom_set_rstw_var_active('nn_fsbc')
	257	CALL iom_set_rstw_var_active('ssu_m')
	258	CALL iom_set_rstw_var_active('ssv_m')
	259	CALL iom_set_rstw_var_active('sst_m')
	260	CALL iom_set_rstw_var_active('sss_m')
	261	CALL iom_set_rstw_var_active('ssh_m')
	262	CALL iom_set_rstw_var_active('e3t_m')
	263	CALL iom_set_rstw_var_active('frq_m')
	264	ENDIF
	265
[3607]	266	END SUBROUTINE sbc_ssm_init
	267
[888]	268	!!======================================================================
	269	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/2017/dev_merge_2017/NEMOGCM/NEMO/OCE_SRC/SBC/sbcssm.F90 @ 9580

Download in other formats: