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.

nemogcm.F90 in NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/SAO – NEMO

Context Navigation

source: NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/SAO/nemogcm.F90 @ 12586

Last change on this file since 12586 was 12586, checked in by francesca, 4 years ago
Add extra-halo support (jperio 3,4) - ticket #2366
Property svn:keywords set to `Id`
File size: 22.2 KB

Rev	Line
[2496]	1	MODULE nemogcm
[2442]	2	!!======================================================================
[2496]	3	!! * MODULE nemogcm *
	4	!! Ocean system : NEMO GCM (ocean dynamics, on-line tracers, biochemistry and sea-ice)
[2442]	5	!!======================================================================
[9598]	6	!! History : 3.6 ! 2015-12 (A. Ryan) Original code (from OCE/)
[7646]	7	!! 4.0 ! 2016-11 (G. Madec, S. Flavoni) domain configuration / user defined interface
[1593]	8	!!----------------------------------------------------------------------
[3]	9
	10	!!----------------------------------------------------------------------
[7646]	11	!! nemo_gcm : solve ocean dynamics, tracer, biogeochemistry and/or sea-ice
	12	!! nemo_init : initialization of the NEMO system
	13	!! nemo_ctl : initialisation of the contol print
	14	!! nemo_closefile: close remaining open files
	15	!! nemo_alloc : dynamical allocation
[3]	16	!!----------------------------------------------------------------------
[7646]	17	USE step_oce ! module used in the ocean time stepping module (step.F90)
	18	USE domain ! domain initialization (dom_init & dom_cfg routines)
	19	USE istate ! initial state setting (istate_init routine)
	20	USE phycst ! physical constant (par_cst routine)
	21	USE step ! NEMO time-stepping (stp routine)
	22	USE cpl_oasis3 ! OASIS3 coupling
	23	USE diaobs ! Observation diagnostics (dia_obs_init routine)
[3625]	24	#if defined key_nemocice_decomp
	25	USE ice_domain_size, only: nx_global, ny_global
	26	#endif
[7646]	27	! ! Stand Alone Observation operator modules
	28	USE sao_data
	29	USE sao_intp
	30	!
	31	USE lib_mpp ! distributed memory computing
	32	USE mppini ! shared/distributed memory setting (mpp_init routine)
[12586]	33	USE lbcnfd , ONLY : isendto, nsndto ! Setup of north fold exchanges
[7646]	34	USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)
[1412]	35	#if defined key_iomput
[7646]	36	USE xios ! xIOserver
[1359]	37	#endif
[12586]	38	USE halo_mng
[268]	39
[2715]	40	IMPLICIT NONE
[3]	41	PRIVATE
	42
[4829]	43	PUBLIC nemo_gcm ! called by model.F90
[2496]	44	PUBLIC nemo_init ! needed by AGRIF
[3764]	45	PUBLIC nemo_alloc ! needed by TAM
[467]	46
[2498]	47	CHARACTER(lc) :: cform_aaa="( /, 'AAAAAAAA', / ) " ! flag for output listing
[1593]	48
[3]	49	!!----------------------------------------------------------------------
[9598]	50	!! NEMO/OCE 4.0 , NEMO Consortium (2018)
[2392]	51	!! $Id$
[10068]	52	!! Software governed by the CeCILL license (see ./LICENSE)
[3]	53	!!----------------------------------------------------------------------
	54	CONTAINS
	55
[4120]	56	SUBROUTINE nemo_gcm
	57	!!----------------------------------------------------------------------
	58	!! * SUBROUTINE offline_obs_oper *
	59	!!
	60	!! ** Purpose : To use NEMO components to interpolate model fields
	61	!! to observation space.
	62	!!
	63	!! ** Method : 1. Initialise NEMO
	64	!! 2. Initialise offline obs_oper
	65	!! 3. Cycle through match ups
	66	!! 4. Write results to file
	67	!!----------------------------------------------------------------------
[7646]	68	!
	69	CALL nemo_init ! Initialise NEMO
	70	!
	71	CALL sao_data_init ! Initialise Stand Alone Observation operator data
	72	!
	73	CALL dia_obs_init ! Initialise obs_operator
	74	!
	75	CALL sao_interp ! Interpolate to observation space
	76	!
	77	CALL dia_obs_wri ! Pipe to output files
	78	!
	79	CALL dia_obs_dealloc ! Reset the obs_oper between
	80	!
	81	IF(lk_mpp) CALL mppstop ! Safely stop MPI (end mpp communications)
	82	!
[4120]	83	END SUBROUTINE nemo_gcm
[389]	84
[4829]	85
[2496]	86	SUBROUTINE nemo_init
[389]	87	!!----------------------------------------------------------------------
[2496]	88	!! * ROUTINE nemo_init *
[389]	89	!!
[2496]	90	!! ** Purpose : initialization of the NEMO GCM
[389]	91	!!----------------------------------------------------------------------
[7646]	92	INTEGER :: ios, ilocal_comm ! local integer
[5600]	93	!
[12377]	94	NAMELIST/namctl/ sn_cfctl, nn_print, nn_ictls, nn_ictle, &
[10570]	95	& nn_isplt , nn_jsplt, nn_jctls, nn_jctle, &
[9019]	96	& ln_timing, ln_diacfl
[9436]	97	NAMELIST/namcfg/ ln_read_cfg, cn_domcfg, ln_closea, ln_write_cfg, cn_domcfg_out, ln_use_jattr
[3]	98	!!----------------------------------------------------------------------
[1593]	99	!
[5600]	100	cxios_context = 'nemo'
[12586]	101	nn_hls = 1
[2496]	102	!
[11536]	103	! !-------------------------------------------------!
	104	! ! set communicator & select the local rank !
	105	! ! must be done as soon as possible to get narea !
	106	! !-------------------------------------------------!
[1593]	107	!
[1412]	108	#if defined key_iomput
[2200]	109	IF( Agrif_Root() ) THEN
[5600]	110	IF( lk_oasis ) THEN
[11536]	111	CALL cpl_init( "oceanx", ilocal_comm ) ! nemo local communicator given by oasis
	112	CALL xios_initialize( "not used" , local_comm =ilocal_comm ) ! send nemo communicator to xios
[5600]	113	ELSE
[11536]	114	CALL xios_initialize( "for_xios_mpi_id", return_comm=ilocal_comm ) ! nemo local communicator given by xios
[5600]	115	ENDIF
[2200]	116	ENDIF
[11536]	117	CALL mpp_start( ilocal_comm )
[532]	118	#else
[5600]	119	IF( lk_oasis ) THEN
	120	IF( Agrif_Root() ) THEN
[11536]	121	CALL cpl_init( "oceanx", ilocal_comm ) ! nemo local communicator given by oasis
[5600]	122	ENDIF
[11536]	123	CALL mpp_start( ilocal_comm )
[5600]	124	ELSE
[11536]	125	CALL mpp_start( )
[2236]	126	ENDIF
[532]	127	#endif
[11536]	128	!
	129	narea = mpprank + 1 ! mpprank: the rank of proc (0 --> mppsize -1 )
	130	lwm = (narea == 1) ! control of output namelists
	131	!
	132	! !---------------------------------------------------------------!
	133	! ! Open output files, reference and configuration namelist files !
	134	! !---------------------------------------------------------------!
	135	!
	136	! open ocean.output as soon as possible to get all output prints (including errors messages)
	137	IF( lwm ) CALL ctl_opn( numout, 'ocean.output', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE. )
	138	! open reference and configuration namelist files
[12377]	139	CALL load_nml( numnam_ref, 'namelist_ref', -1, lwm )
	140	CALL load_nml( numnam_cfg, 'namelist_cfg', -1, lwm )
[11536]	141	IF( lwm ) CALL ctl_opn( numond, 'output.namelist.dyn', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE. )
	142	! open /dev/null file to be able to supress output write easily
	143	CALL ctl_opn( numnul, '/dev/null', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE. )
	144	!
	145	! !--------------------!
[12377]	146	! ! Open listing units ! -> need sn_cfctl from namctl to define lwp
[11536]	147	! !--------------------!
	148	!
	149	READ ( numnam_ref, namctl, IOSTAT = ios, ERR = 901 )
	150	901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namctl in reference namelist' )
	151	READ ( numnam_cfg, namctl, IOSTAT = ios, ERR = 902 )
	152	902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namctl in configuration namelist' )
	153	!
[12377]	154	! finalize the definition of namctl variables
	155	IF( sn_cfctl%l_allon ) THEN
	156	! Turn on all options.
	157	CALL nemo_set_cfctl( sn_cfctl, .TRUE., .TRUE. )
	158	! Ensure all processors are active
	159	sn_cfctl%procmin = 0 ; sn_cfctl%procmax = 1000000 ; sn_cfctl%procincr = 1
	160	ELSEIF( sn_cfctl%l_config ) THEN
	161	! Activate finer control of report outputs
	162	! optionally switch off output from selected areas (note this only
	163	! applies to output which does not involve global communications)
	164	IF( ( narea < sn_cfctl%procmin .OR. narea > sn_cfctl%procmax ) .OR. &
	165	& ( MOD( narea - sn_cfctl%procmin, sn_cfctl%procincr ) /= 0 ) ) &
	166	& CALL nemo_set_cfctl( sn_cfctl, .FALSE., .FALSE. )
	167	ELSE
	168	! turn off all options.
	169	CALL nemo_set_cfctl( sn_cfctl, .FALSE., .TRUE. )
	170	ENDIF
[11536]	171	!
[12377]	172	lwp = (narea == 1) .OR. sn_cfctl%l_oceout ! control of all listing output print
	173	!
[11536]	174	IF(lwp) THEN ! open listing units
[1593]	175	!
[11536]	176	IF( .NOT. lwm ) & ! alreay opened for narea == 1
	177	& CALL ctl_opn( numout, 'ocean.output', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE., narea )
[1593]	178	!
[1579]	179	WRITE(numout,*)
[11536]	180	WRITE(numout,*) ' CNRS - NERC - Met OFFICE - MERCATOR-ocean - CMCC'
[1593]	181	WRITE(numout,*) ' NEMO team'
[5042]	182	WRITE(numout,*) ' Stand Alone Observation operator'
[10570]	183	WRITE(numout,*) ' NEMO version 4.0 (2019) '
[1579]	184	WRITE(numout,*)
[10570]	185	WRITE(numout,*) " ._ ._ ._ ._ ._ "
	186	WRITE(numout,*) " _.-._)`\_.-._)`\_.-._)`\_.-._)`\_.-._)`\_ "
[1579]	187	WRITE(numout,*)
[10570]	188	WRITE(numout,*) " o _, _, "
	189	WRITE(numout,*) " o .' ( .-' / "
	190	WRITE(numout,*) " o _/..._'. .' / "
	191	WRITE(numout,*) " ( o .-'` ` '-./ _.' "
	192	WRITE(numout,*) " ) ( o) ;= <_ ( "
	193	WRITE(numout,*) " ( '-.,\\__ __.-;`\ '. ) "
	194	WRITE(numout,*) " ) ) \) \|`\ \) '. \ ( ( "
	195	WRITE(numout,*) " ( ( \_/ '-._\ ) ) "
[11536]	196	WRITE(numout,*) " ) ) jgs ` ( ( "
[10570]	197	WRITE(numout,*) " ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ "
	198	WRITE(numout,*)
[1593]	199	!
[7646]	200	WRITE(numout,cform_aaa) ! Flag AAAAAAA
	201	!
[473]	202	ENDIF
[11536]	203	!
	204	IF(lwm) WRITE( numond, namctl )
	205	!
	206	! !------------------------------------!
	207	! ! Set global domain size parameters !
	208	! !------------------------------------!
	209	!
	210	READ ( numnam_ref, namcfg, IOSTAT = ios, ERR = 903 )
	211	903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcfg in reference namelist' )
	212	READ ( numnam_cfg, namcfg, IOSTAT = ios, ERR = 904 )
	213	904 IF( ios > 0 ) CALL ctl_nam ( ios , 'namcfg in configuration namelist' )
	214	!
	215	IF( ln_read_cfg ) THEN ! Read sizes in domain configuration file
	216	CALL domain_cfg ( cn_cfg, nn_cfg, jpiglo, jpjglo, jpkglo, jperio )
	217	ELSE ! user-defined namelist
	218	CALL usr_def_nam( cn_cfg, nn_cfg, jpiglo, jpjglo, jpkglo, jperio )
	219	ENDIF
	220	!
	221	IF(lwm) WRITE( numond, namcfg )
	222	!
	223	! !-----------------------------------------!
	224	! ! mpp parameters and domain decomposition !
	225	! !-----------------------------------------!
	226	CALL mpp_init
[2715]	227
[12586]	228	CALL halo_mng_init()
[7646]	229	! Now we know the dimensions of the grid and numout has been set: we can allocate arrays
[2715]	230	CALL nemo_alloc()
	231
[2496]	232	! !-------------------------------!
	233	! ! NEMO general initialization !
	234	! !-------------------------------!
[473]	235
[9436]	236	CALL nemo_ctl ! Control prints
[2382]	237	!
[9436]	238	! ! General initialization
	239	IF( ln_timing ) CALL timing_init ! timing
	240	IF( ln_timing ) CALL timing_start( 'nemo_init')
[3294]	241	!
[9019]	242	CALL phy_cst ! Physical constants
	243	CALL eos_init ! Equation of state
[12377]	244	CALL dom_init( Nbb, Nnn, Naa, 'SAO') ! Domain
[413]	245
[3294]	246
[12377]	247	IF( sn_cfctl%l_prtctl ) &
	248	& CALL prt_ctl_init ! Print control
[2027]	249
[12377]	250	CALL istate_init ! ocean initial state (Dynamics and tracers)
[2496]	251	END SUBROUTINE nemo_init
[467]	252
	253
[2496]	254	SUBROUTINE nemo_ctl
[467]	255	!!----------------------------------------------------------------------
[2496]	256	!! * ROUTINE nemo_ctl *
[467]	257	!!
[3764]	258	!! ** Purpose : control print setting
[467]	259	!!
[12377]	260	!! ** Method : - print namctl and namcfg information and check some consistencies
[467]	261	!!----------------------------------------------------------------------
[2442]	262	!
[2496]	263	IF(lwp) THEN ! control print
[531]	264	WRITE(numout,*)
[7646]	265	WRITE(numout,*) 'nemo_ctl: Control prints'
[9436]	266	WRITE(numout,*) '~~~~~~~~'
[1593]	267	WRITE(numout,*) ' Namelist namctl'
[12377]	268	WRITE(numout,*) ' sn_cfctl%l_glochk = ', sn_cfctl%l_glochk
	269	WRITE(numout,*) ' sn_cfctl%l_allon = ', sn_cfctl%l_allon
[10570]	270	WRITE(numout,*) ' finer control over o/p sn_cfctl%l_config = ', sn_cfctl%l_config
	271	WRITE(numout,*) ' sn_cfctl%l_runstat = ', sn_cfctl%l_runstat
	272	WRITE(numout,*) ' sn_cfctl%l_trcstat = ', sn_cfctl%l_trcstat
	273	WRITE(numout,*) ' sn_cfctl%l_oceout = ', sn_cfctl%l_oceout
	274	WRITE(numout,*) ' sn_cfctl%l_layout = ', sn_cfctl%l_layout
[12377]	275	WRITE(numout,*) ' sn_cfctl%l_prtctl = ', sn_cfctl%l_prtctl
	276	WRITE(numout,*) ' sn_cfctl%l_prttrc = ', sn_cfctl%l_prttrc
	277	WRITE(numout,*) ' sn_cfctl%l_oasout = ', sn_cfctl%l_oasout
[10570]	278	WRITE(numout,*) ' sn_cfctl%procmin = ', sn_cfctl%procmin
	279	WRITE(numout,*) ' sn_cfctl%procmax = ', sn_cfctl%procmax
	280	WRITE(numout,*) ' sn_cfctl%procincr = ', sn_cfctl%procincr
	281	WRITE(numout,*) ' sn_cfctl%ptimincr = ', sn_cfctl%ptimincr
[1601]	282	WRITE(numout,*) ' level of print nn_print = ', nn_print
	283	WRITE(numout,*) ' Start i indice for SUM control nn_ictls = ', nn_ictls
	284	WRITE(numout,*) ' End i indice for SUM control nn_ictle = ', nn_ictle
	285	WRITE(numout,*) ' Start j indice for SUM control nn_jctls = ', nn_jctls
	286	WRITE(numout,*) ' End j indice for SUM control nn_jctle = ', nn_jctle
	287	WRITE(numout,*) ' number of proc. following i nn_isplt = ', nn_isplt
	288	WRITE(numout,*) ' number of proc. following j nn_jsplt = ', nn_jsplt
[9019]	289	WRITE(numout,*) ' timing by routine ln_timing = ', ln_timing
	290	WRITE(numout,*) ' CFL diagnostics ln_diacfl = ', ln_diacfl
[531]	291	ENDIF
[2442]	292	!
[1601]	293	nprint = nn_print ! convert DOCTOR namelist names into OLD names
	294	nictls = nn_ictls
	295	nictle = nn_ictle
	296	njctls = nn_jctls
	297	njctle = nn_jctle
	298	isplt = nn_isplt
	299	jsplt = nn_jsplt
[4829]	300
	301	IF(lwp) THEN ! control print
	302	WRITE(numout,*)
	303	WRITE(numout,*) ' Namelist namcfg'
[7646]	304	WRITE(numout,*) ' read domain configuration file ln_read_cfg = ', ln_read_cfg
	305	WRITE(numout,*) ' filename to be read cn_domcfg = ', TRIM(cn_domcfg)
[9436]	306	WRITE(numout,*) ' keep closed seas in the domain (if exist) ln_closea = ', ln_closea
[10570]	307	WRITE(numout,*) ' create a configuration definition file ln_write_cfg = ', ln_write_cfg
[7646]	308	WRITE(numout,*) ' filename to be written cn_domcfg_out = ', TRIM(cn_domcfg_out)
	309	WRITE(numout,*) ' use file attribute if exists as i/p j-start ln_use_jattr = ', ln_use_jattr
[4829]	310	ENDIF
[9436]	311	IF( .NOT.ln_read_cfg ) ln_closea = .false. ! dealing possible only with a domcfg file
	312	!
[2442]	313	! ! Parameter control
[1593]	314	!
[12377]	315	IF( sn_cfctl%l_prtctl .OR. sn_cfctl%l_prttrc ) THEN ! sub-domain area indices for the control prints
[3294]	316	IF( lk_mpp .AND. jpnij > 1 ) THEN
[2496]	317	isplt = jpni ; jsplt = jpnj ; ijsplt = jpni*jpnj ! the domain is forced to the real split domain
[531]	318	ELSE
	319	IF( isplt == 1 .AND. jsplt == 1 ) THEN
[1593]	320	CALL ctl_warn( ' - isplt & jsplt are equal to 1', &
	321	& ' - the print control will be done over the whole domain' )
[531]	322	ENDIF
[1593]	323	ijsplt = isplt * jsplt ! total number of processors ijsplt
[531]	324	ENDIF
	325	IF(lwp) WRITE(numout,*)' - The total number of processors over which the'
	326	IF(lwp) WRITE(numout,*)' print control will be done is ijsplt : ', ijsplt
[1593]	327	!
	328	! ! indices used for the SUM control
	329	IF( nictls+nictle+njctls+njctle == 0 ) THEN ! print control done over the default area
[3764]	330	lsp_area = .FALSE.
[1593]	331	ELSE ! print control done over a specific area
[531]	332	lsp_area = .TRUE.
	333	IF( nictls < 1 .OR. nictls > jpiglo ) THEN
	334	CALL ctl_warn( ' - nictls must be 1<=nictls>=jpiglo, it is forced to 1' )
	335	nictls = 1
	336	ENDIF
	337	IF( nictle < 1 .OR. nictle > jpiglo ) THEN
	338	CALL ctl_warn( ' - nictle must be 1<=nictle>=jpiglo, it is forced to jpiglo' )
	339	nictle = jpiglo
	340	ENDIF
	341	IF( njctls < 1 .OR. njctls > jpjglo ) THEN
	342	CALL ctl_warn( ' - njctls must be 1<=njctls>=jpjglo, it is forced to 1' )
	343	njctls = 1
	344	ENDIF
	345	IF( njctle < 1 .OR. njctle > jpjglo ) THEN
	346	CALL ctl_warn( ' - njctle must be 1<=njctle>=jpjglo, it is forced to jpjglo' )
	347	njctle = jpjglo
	348	ENDIF
[1593]	349	ENDIF
	350	ENDIF
[2442]	351	!
[10570]	352	IF( 1._wp /= SIGN(1._wp,-0._wp) ) CALL ctl_stop( 'nemo_ctl: The intrinsec SIGN function follows f2003 standard.', &
	353	& 'Compile with key_nosignedzero enabled:', &
	354	& '--> add -Dkey_nosignedzero to the definition of %CPP in your arch file' )
[3764]	355	!
[2496]	356	END SUBROUTINE nemo_ctl
[467]	357
	358
[2496]	359	SUBROUTINE nemo_closefile
[467]	360	!!----------------------------------------------------------------------
[2496]	361	!! * ROUTINE nemo_closefile *
[467]	362	!!
	363	!! ** Purpose : Close the files
	364	!!----------------------------------------------------------------------
[1593]	365	!
	366	IF( lk_mpp ) CALL mppsync
	367	!
[1685]	368	CALL iom_close ! close all input/output files managed by iom_*
[1593]	369	!
[4829]	370	IF( numstp /= -1 ) CLOSE( numstp ) ! time-step file
[9019]	371	IF( numrun /= -1 ) CLOSE( numrun ) ! run statistics file
[4829]	372	IF( lwm.AND.numond /= -1 ) CLOSE( numond ) ! oce output namelist
	373	IF( lwm.AND.numoni /= -1 ) CLOSE( numoni ) ! ice output namelist
	374	IF( numevo_ice /= -1 ) CLOSE( numevo_ice ) ! ice variables (temp. evolution)
	375	IF( numout /= 6 ) CLOSE( numout ) ! standard model output file
	376	IF( numdct_vol /= -1 ) CLOSE( numdct_vol ) ! volume transports
	377	IF( numdct_heat /= -1 ) CLOSE( numdct_heat ) ! heat transports
	378	IF( numdct_salt /= -1 ) CLOSE( numdct_salt ) ! salt transports
[1593]	379	!
[2442]	380	numout = 6 ! redefine numout in case it is used after this point...
	381	!
[2496]	382	END SUBROUTINE nemo_closefile
[467]	383
[2715]	384
	385	SUBROUTINE nemo_alloc
	386	!!----------------------------------------------------------------------
	387	!! * ROUTINE nemo_alloc *
	388	!!
	389	!! ** Purpose : Allocate all the dynamic arrays of the OPA modules
	390	!!
	391	!! ** Method :
	392	!!----------------------------------------------------------------------
	393	USE diawri , ONLY: dia_wri_alloc
	394	USE dom_oce , ONLY: dom_oce_alloc
	395	!
	396	INTEGER :: ierr
	397	!!----------------------------------------------------------------------
	398	!
[3764]	399	ierr = oce_alloc () ! ocean
[2715]	400	ierr = ierr + dia_wri_alloc ()
	401	ierr = ierr + dom_oce_alloc () ! ocean domain
	402	!
[10425]	403	CALL mpp_sum( 'nemogcm', ierr )
[2715]	404	IF( ierr /= 0 ) CALL ctl_stop( 'STOP', 'nemo_alloc : unable to allocate standard ocean arrays' )
	405	!
	406	END SUBROUTINE nemo_alloc
	407
[10570]	408	SUBROUTINE nemo_set_cfctl(sn_cfctl, setto, for_all )
	409	!!----------------------------------------------------------------------
	410	!! * ROUTINE nemo_set_cfctl *
	411	!!
	412	!! ** Purpose : Set elements of the output control structure to setto.
	413	!! for_all should be .false. unless all areas are to be
	414	!! treated identically.
	415	!!
	416	!! ** Method : Note this routine can be used to switch on/off some
	417	!! types of output for selected areas but any output types
	418	!! that involve global communications (e.g. mpp_max, glob_sum)
	419	!! should be protected from selective switching by the
	420	!! for_all argument
	421	!!----------------------------------------------------------------------
	422	LOGICAL :: setto, for_all
[10601]	423	TYPE(sn_ctl) :: sn_cfctl
[10570]	424	!!----------------------------------------------------------------------
	425	IF( for_all ) THEN
	426	sn_cfctl%l_runstat = setto
	427	sn_cfctl%l_trcstat = setto
	428	ENDIF
	429	sn_cfctl%l_oceout = setto
	430	sn_cfctl%l_layout = setto
[12377]	431	sn_cfctl%l_prtctl = setto
	432	sn_cfctl%l_prttrc = setto
	433	sn_cfctl%l_oasout = setto
[10570]	434	END SUBROUTINE nemo_set_cfctl
	435
[3]	436	!!======================================================================
[2496]	437	END MODULE nemogcm
[4829]	438

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

Download in other formats: