Context Navigation

← Previous Revision
Next Revision →
Normal
Revision Log

icestp.F90

Last change on this file was 8892, checked in by frrh, 6 years ago
Commit updates with debugging write statements.
File size: 26.0 KB

Rev	Line
[8321]	1	MODULE icestp
	2	!!======================================================================
	3	!! * MODULE icestp *
	4	!! Surface module : update the ocean surface boundary condition over ice
[8534]	5	!! covered area using ESIM sea-ice model
[8321]	6	!!=====================================================================
	7	!! History : 2.0 ! 2006-12 (M. Vancoppenolle) Original code
	8	!! 3.0 ! 2008-02 (C. Talandier) Surface module from icestp.F90
[8411]	9	!! - ! 2008-04 (G. Madec) sltyle and ice_ctl routine
[8321]	10	!! 3.3 ! 2010-11 (G. Madec) ice-ocean stress always computed at each ocean time-step
	11	!! 3.4 ! 2011-01 (A Porter) dynamical allocation
[8426]	12	!! - ! 2012-10 (C. Rousset) add ice_dia
[8321]	13	!! 3.6 ! 2014-07 (M. Vancoppenolle, G. Madec, O. Marti) revise coupled interface
	14	!! 4.0 ! 2016-06 (L. Brodeau) new unified bulk routine (based on AeroBulk)
	15	!!----------------------------------------------------------------------
	16	#if defined key_lim3
	17	!!----------------------------------------------------------------------
[8534]	18	!! 'key_lim3' ESIM sea-ice model
[8321]	19	!!----------------------------------------------------------------------
[8534]	20	!! ice_stp : sea-ice model time-stepping and update ocean SBC over ice-covered area
	21	!! ice_init : initialize sea-ice
[8321]	22	!!----------------------------------------------------------------------
[8486]	23	USE oce ! ocean dynamics and tracers
	24	USE dom_oce ! ocean space and time domain
	25	USE c1d ! 1D vertical configuration
	26	USE ice ! sea-ice: variables
	27	USE ice1D ! sea-ice: thermodynamical 1D variables
[8321]	28	!
[8534]	29	USE phycst ! Define parameters for the routines
	30	USE eosbn2 ! equation of state
[8486]	31	USE sbc_oce ! Surface boundary condition: ocean fields
	32	USE sbc_ice ! Surface boundary condition: ice fields
[8534]	33	!
[8486]	34	USE iceforcing ! sea-ice: Surface boundary condition !!gm why not icesbc module name
[8516]	35	USE icedyn ! sea-ice: dynamics
[8486]	36	USE icethd ! sea-ice: thermodynamics
[8534]	37	USE icecor ! sea-ice: corrections
[8486]	38	USE iceupdate ! sea-ice: sea surface boundary condition update
	39	USE icedia ! sea-ice: budget diagnostics
	40	USE icewri ! sea-ice: outputs
	41	USE icerst ! sea-ice: restarts
	42	USE icevar ! sea-ice: operations
	43	USE icectl ! sea-ice: control
[8514]	44	USE iceistate ! sea-ice: initial state
[8505]	45	USE iceitd ! sea-ice: remapping thickness distribution
	46	USE icealb ! sea-ice: albedo
[8321]	47	!
[8486]	48	USE bdy_oce , ONLY : ln_bdy ! flag for bdy
	49	USE bdyice ! unstructured open boundary data for sea-ice
[8321]	50	# if defined key_agrif
	51	USE agrif_ice
	52	USE agrif_lim3_update
	53	USE agrif_lim3_interp
	54	# endif
[8486]	55	!
	56	USE in_out_manager ! I/O manager
	57	USE iom ! I/O manager library
	58	USE lib_mpp ! MPP library
[8534]	59	USE lib_fortran ! fortran utilities (glob_sum + no signed zero)
[8486]	60	USE timing ! Timing
[8534]	61	USE prtctl ! Print control
[8321]	62
	63	IMPLICIT NONE
	64	PRIVATE
	65
[8486]	66	PUBLIC ice_stp ! called by sbcmod.F90
	67	PUBLIC ice_init ! called by sbcmod.F90
[8321]	68
	69	!! * Substitutions
	70	# include "vectopt_loop_substitute.h90"
	71	!!----------------------------------------------------------------------
[8486]	72	!! NEMO/ICE 4.0 , NEMO Consortium (2017)
[8321]	73	!! $Id: icestp.F90 8319 2017-07-11 15:00:44Z clem $
	74	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
	75	!!----------------------------------------------------------------------
	76	CONTAINS
	77
	78	SUBROUTINE ice_stp( kt, ksbc )
	79	!!---------------------------------------------------------------------
	80	!! * ROUTINE ice_stp *
	81	!!
[8486]	82	!! ** Purpose : sea-ice model time-stepping and update ocean surface
	83	!! boundary condition over ice-covered area
[8321]	84	!!
	85	!! ** Method : ice model time stepping
	86	!! - call the ice dynamics routine
	87	!! - call the ice advection/diffusion routine
	88	!! - call the ice thermodynamics routine
	89	!! - call the routine that computes mass and
	90	!! heat fluxes at the ice/ocean interface
	91	!! - save the outputs
	92	!! - save the outputs for restart when necessary
	93	!!
	94	!! ** Action : - time evolution of the LIM sea-ice model
	95	!! - update all sbc variables below sea-ice:
	96	!! utau, vtau, taum, wndm, qns , qsr, emp , sfx
	97	!!---------------------------------------------------------------------
	98	INTEGER, INTENT(in) :: kt ! ocean time step
[8486]	99	INTEGER, INTENT(in) :: ksbc ! flux formulation (user defined, bulk, or Pure Coupled)
	100	!
	101	INTEGER :: jl ! dummy loop index
[8321]	102	!!----------------------------------------------------------------------
[8486]	103	!
[8321]	104	IF( nn_timing == 1 ) CALL timing_start('ice_stp')
[8486]	105	!
	106	! !-----------------------!
	107	IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN ! --- Ice time step --- !
	108	! !-----------------------!
[8892]	109	!meaning
[8498]	110	kt_ice = kt ! -- Ice model time step
[8486]	111	!
[8321]	112	# if defined key_agrif
	113	IF( .NOT. Agrif_Root() ) lim_nbstep = MOD( lim_nbstep, Agrif_irhot() * Agrif_Parent(nn_fsbc) / nn_fsbc ) + 1
	114	# endif
[8498]	115	u_oce(:,:) = ssu_m(:,:) ! -- mean surface ocean current
	116	v_oce(:,:) = ssv_m(:,:)
	117	!
[8892]	118	write(numout,*) "RSRH ice_stp 00 before eos_stp ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
	119	write(numout,*) "RSRH ice_stp 00 qns,sfx,emp before eos_stp ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
[8498]	120	CALL eos_fzp( sss_m(:,:) , t_bo(:,:) ) ! -- freezing temperature [Kelvin] (set to rt0 over land)
[8321]	121	t_bo(:,:) = ( t_bo(:,:) + rt0 ) * tmask(:,:,1) + rt0 * ( 1._wp - tmask(:,:,1) )
	122	!
[8498]	123	CALL store_fields ! -- Store now ice values
[8486]	124	!
[8321]	125	!------------------------------------------------!
	126	! --- Dynamical coupling with the atmosphere --- !
	127	!------------------------------------------------!
[8486]	128	! It provides the following fields used in sea ice model:
[8404]	129	! utau_ice, vtau_ice = surface ice stress [N/m2]
[8486]	130	!------------------------------------------------!
[8562]	131	CALL ice_forcing_tau( kt, ksbc, utau_ice, vtau_ice )
[8486]	132	!-------------------------------------!
	133	! --- ice dynamics and advection --- !
	134	!-------------------------------------!
[8531]	135	CALL diag_set0 ! set diag of mass, heat and salt fluxes to 0
[8486]	136	CALL ice_rst_opn( kt ) ! Open Ice restart file (if necessary)
[8321]	137	!
[8892]	138	write(numout,*) "RSRH ice_stp 00 before ice_dyn ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
	139	write(numout,*) "RSRH ice_stp 00 qns,sfx,emp before ice_dyn ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
[8562]	140	IF( ln_icedyn .AND. .NOT.lk_c1d ) &
	141	& CALL ice_dyn( kt ) ! -- Ice dynamics
[8892]	142	write(numout,*) "RSRH ice_stp 00 after ice_dyn ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
	143	write(numout,*) "RSRH ice_stp 00 qns,sfx,emp aft ice_dyn ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
[8486]	144	!
	145	! !== lateral boundary conditions ==!
[8321]	146	#if defined key_agrif
[8562]	147	IF( .NOT. Agrif_Root() ) CALL agrif_interp_lim3('T') ! -- AGRIF
[8321]	148	#endif
[8562]	149	IF( ln_icethd .AND. ln_bdy ) CALL bdy_ice( kt ) ! -- bdy ice thermo
[8486]	150	!
	151	!
	152	! !== previous lead fraction and ice volume for flux calculations
[8892]	153	write(numout,*) "RSRH ice_stp 00 be4 glo2eqv ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
	154	write(numout,*) "RSRH ice_stp 00 qns,sfx,emp b4 glo2eqv ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
[8563]	155	CALL ice_var_glo2eqv ! h_i and h_s for ice albedo calculation
[8892]	156	write(numout,*) "RSRH ice_stp 00 be4 ice_var_agg ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8486]	157	CALL ice_var_agg(1) ! at_i for coupling
[8892]	158	write(numout,*) "RSRH ice_stp 00 qns,sfx,emp b4 store_filds ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	159	write(numout,*) "RSRH ice_stp 00 be4 store_filds ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8498]	160	CALL store_fields ! Store now ice values
[8321]	161
	162	!------------------------------------------------------!
	163	! --- Thermodynamical coupling with the atmosphere --- !
	164	!------------------------------------------------------!
[8404]	165	! It provides the following fields used in sea ice model:
[8486]	166	! emp_oce , emp_ice = E-P over ocean and sea ice [Kg/m2/s]
	167	! sprecip = solid precipitation [Kg/m2/s]
	168	! evap_ice = sublimation [Kg/m2/s]
	169	! qsr_tot , qns_tot = solar & non solar heat flux (total) [W/m2]
	170	! qsr_ice , qns_ice = solar & non solar heat flux over ice [W/m2]
	171	! dqns_ice = non solar heat sensistivity [W/m2]
	172	! qemp_oce, qemp_ice, = sensible heat (associated with evap & precip) [W/m2]
	173	! qprec_ice, qevap_ice
	174	!------------------------------------------------------!
[8892]	175	write(numout,*) "RSRH ice_stp 00 qns,sfx,emp before forc ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	176	write(numout,*) "RSRH ice_stp 00 be4 forc ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8562]	177	CALL ice_forcing_flx( kt, ksbc )
[8321]	178	!----------------------------!
	179	! --- ice thermodynamics --- !
	180	!----------------------------!
[8892]	181	write(numout,*) "RSRH ice_stp 00 qns,sfx,emp before thd ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	182	write(numout,*) "RSRH ice_stp 00 be4 thd", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8562]	183	IF( ln_icethd ) CALL ice_thd( kt ) ! -- Ice thermodynamics
	184	!
[8892]	185	write(numout,*) "RSRH ice_stp 00 qns,sfx,emp before icecor ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	186	write(numout,*) "RSRH ice_stp 00 before icecor ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8562]	187	!
	188	IF( ln_icethd ) CALL ice_cor( kt , 2 ) ! -- Corrections
[8892]	189	write(numout,*) "RSRH ice_stp AA qns,sfx,emp ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	190	write(numout,*) "RSRH ice_stp AA ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8321]	191	# if defined key_agrif
[8562]	192	IF( .NOT. Agrif_Root() ) CALL agrif_update_lim3( kt )
[8321]	193	# endif
[8562]	194	CALL ice_var_glo2eqv ! necessary calls (at least for coupling)
[8892]	195	write(numout,*) "RSRH ice_stp BB qns,sfx,emp ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	196	write(numout,*) "RSRH ice_stp BB ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8562]	197	CALL ice_var_agg( 2 ) ! necessary calls (at least for coupling)
[8892]	198	write(numout,*) "RSRH ice_stp CC qns,sfx,emp ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	199	write(numout,*) "RSRH ice_stp CC ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8562]	200	!
[8321]	201	!! clem: one should switch the calculation of the fluxes onto the parent grid but the following calls do not work
	202	!! moreover it should only be called at the update frequency (as in agrif_lim3_update.F90)
	203	!!# if defined key_agrif
[8562]	204	!! IF( .NOT. Agrif_Root() ) CALL Agrif_ChildGrid_To_ParentGrid()
[8321]	205	!!# endif
[8562]	206	CALL ice_update_flx( kt ) ! -- Update ocean surface mass, heat and salt fluxes
[8892]	207	write(numout,*) "RSRH ice_stp DD qns,sfx,emp ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	208	write(numout,*) "RSRH ice_stp DD ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8321]	209	!!# if defined key_agrif
[8562]	210	!! IF( .NOT. Agrif_Root() ) CALL Agrif_ParentGrid_To_ChildGrid()
[8321]	211	!!# endif
[8562]	212	IF( ln_icediahsb ) CALL ice_dia( kt ) ! -- Diagnostics outputs
[8892]	213	write(numout,*) "RSRH ice_stp EE qns,sfx,emp ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	214	write(numout,*) "RSRH ice_stp EE ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8321]	215	!
[8562]	216	CALL ice_wri( kt ) ! -- Ice outputs
[8892]	217	write(numout,*) "RSRH ice_stp FF qns,sfx,emp ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	218	write(numout,*) "RSRH ice_stp FF ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8321]	219	!
[8562]	220	IF( lrst_ice ) CALL ice_rst_write( kt ) ! -- Ice restart file
[8892]	221	write(numout,*) "RSRH ice_stp HH qns,sfx,emp ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	222	write(numout,*) "RSRH ice_stp HH ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8321]	223	!
[8562]	224	IF( ln_icectl ) CALL ice_ctl( kt ) ! -- alerts in case of model crash
[8892]	225	write(numout,*) "RSRH ice_stp II qns,sfx,emp ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	226	write(numout,*) "RSRH ice_stp II ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8321]	227	!
	228	ENDIF ! End sea-ice time step only
	229
	230	!-------------------------!
	231	! --- Ocean time step --- !
	232	!-------------------------!
[8562]	233	IF( ln_icedyn ) CALL ice_update_tau( kt, ub(:,:,1), vb(:,:,1) ) ! -- update surface ocean stresses
[8892]	234	write(numout,*) "RSRH ice_stp JJ qns,sfx,emp ", kt , glob_sum(qns(:,:)), glob_sum(sfx(:,:)),glob_sum(emp(:,:)) ; flush(numout)
	235	write(numout,*) "RSRH ice_stp JJ ", kt , glob_sum(ssv_m(:,:)), glob_sum(v_ice(:,:)),glob_sum(ssu_m(:,:)), glob_sum(u_ice(:,:)) ; flush(numout)
[8321]	236	!!gm remark, the ocean-ice stress is not saved in ice diag call above ..... find a solution!!!
	237	!
[8562]	238	IF( nn_timing == 1 ) CALL timing_stop('ice_stp')
[8321]	239	!
	240	END SUBROUTINE ice_stp
	241
	242
	243	SUBROUTINE ice_init
	244	!!----------------------------------------------------------------------
	245	!! * ROUTINE ice_init *
	246	!!
[8534]	247	!! ** purpose : Initialize sea-ice parameters
[8321]	248	!!----------------------------------------------------------------------
	249	INTEGER :: ji, jj, ierr
	250	!!----------------------------------------------------------------------
	251	IF(lwp) WRITE(numout,*)
[8514]	252	IF(lwp) WRITE(numout,*) 'ice_init: Arrays allocation & Initialization off all routines & init state'
	253	IF(lwp) WRITE(numout,*) '~~~~~~~~'
[8321]	254	!
	255	! ! Open the reference and configuration namelist files and namelist output file
	256	CALL ctl_opn( numnam_ice_ref, 'namelist_ice_ref', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
	257	CALL ctl_opn( numnam_ice_cfg, 'namelist_ice_cfg', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
	258	IF(lwm) CALL ctl_opn( numoni, 'output.namelist.ice', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, 1 )
	259	!
[8531]	260	CALL par_init ! set some ice run parameters
[8321]	261	!
[8324]	262	! ! Allocate the ice arrays (sbc_ice already allocated in sbc_init)
[8321]	263	ierr = ice_alloc () ! ice variables
[8331]	264	ierr = ierr + sbc_ice_alloc () ! surface forcing
[8420]	265	ierr = ierr + ice1D_alloc () ! thermodynamics
[8321]	266	!
	267	IF( lk_mpp ) CALL mpp_sum( ierr )
	268	IF( ierr /= 0 ) CALL ctl_stop('STOP', 'ice_init : unable to allocate ice arrays')
	269	!
	270	CALL ice_itd_init ! ice thickness distribution initialization
	271	!
[8752]	272	IF( ln_icethd ) THEN
	273	CALL ice_thd_init ! set ice thermodynics parameters (clem: important to call it first for melt ponds)
	274	ENDIF
[8321]	275	! ! Initial sea-ice state
	276	IF( .NOT. ln_rstart ) THEN ! start from rest: sea-ice deduced from sst
[8515]	277	CALL ice_istate_init
[8514]	278	CALL ice_istate
[8321]	279	ELSE ! start from a restart file
[8413]	280	CALL ice_rst_read
[8321]	281	ENDIF
[8534]	282	CALL ice_var_glo2eqv
[8424]	283	CALL ice_var_agg(2)
[8321]	284	!
[8531]	285	CALL ice_forcing_init ! set ice-ocean and ice-atm. coupling parameters
	286	!
[8518]	287	IF( ln_icedyn ) THEN
	288	CALL ice_dyn_init ! set ice dynamics parameters
	289	ENDIF
[8531]	290	!
[8517]	291	CALL ice_update_init ! ice surface boundary condition
[8321]	292	!
[8517]	293	CALL ice_alb_init ! ice surface albedo
[8321]	294	!
[8517]	295	CALL ice_dia_init ! initialization for diags
[8505]	296	!
[8771]	297	fr_i (:,:) = at_i(:,:) ! initialisation of sea-ice fraction
	298	tn_ice (:,:,:) = t_su(:,:,:) ! initialisation of surface temp for coupled simu
[8321]	299	!
[8517]	300	! ! set max concentration in both hemispheres
[8515]	301	WHERE( gphit(:,:) > 0._wp ) ; rn_amax_2d(:,:) = rn_amax_n ! NH
	302	ELSEWHERE ; rn_amax_2d(:,:) = rn_amax_s ! SH
	303	END WHERE
[8518]	304
	305	IF( ln_rstart ) CALL iom_close( numrir ) ! close input ice restart file
[8321]	306	!
	307	END SUBROUTINE ice_init
	308
	309
[8531]	310	SUBROUTINE par_init
[8321]	311	!!-------------------------------------------------------------------
[8531]	312	!! * ROUTINE par_init *
[8321]	313	!!
[8534]	314	!! ** Purpose : Definition generic parameters for ice model
[8321]	315	!!
[8534]	316	!! ** Method : Read namelist and check the parameter
[8517]	317	!! values called at the first timestep (nit000)
[8321]	318	!!
[8531]	319	!! ** input : Namelist nampar
[8321]	320	!!-------------------------------------------------------------------
	321	INTEGER :: ios ! Local integer output status for namelist read
[8531]	322	NAMELIST/nampar/ jpl, nlay_i, nlay_s, nn_monocat, ln_icedyn, ln_icethd, rn_amax_n, rn_amax_s, &
	323	& cn_icerst_in, cn_icerst_indir, cn_icerst_out, cn_icerst_outdir
[8321]	324	!!-------------------------------------------------------------------
	325	!
[8531]	326	REWIND( numnam_ice_ref ) ! Namelist nampar in reference namelist : Parameters for ice
	327	READ ( numnam_ice_ref, nampar, IOSTAT = ios, ERR = 901)
	328	901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampar in reference namelist', lwp )
[8321]	329
[8531]	330	REWIND( numnam_ice_cfg ) ! Namelist nampar in configuration namelist : Parameters for ice
	331	READ ( numnam_ice_cfg, nampar, IOSTAT = ios, ERR = 902 )
	332	902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampar in configuration namelist', lwp )
	333	IF(lwm) WRITE ( numoni, nampar )
[8321]	334	!
[8486]	335	IF(lwp) THEN ! control print
[8321]	336	WRITE(numout,*)
[8531]	337	WRITE(numout,*) 'par_init: ice parameters shared among all the routines'
	338	WRITE(numout,*) ' ~~~~~~~'
	339	WRITE(numout,*) ' Namelist nampar: '
[8486]	340	WRITE(numout,*) ' number of ice categories jpl = ', jpl
	341	WRITE(numout,*) ' number of ice layers nlay_i = ', nlay_i
	342	WRITE(numout,*) ' number of snow layers nlay_s = ', nlay_s
	343	WRITE(numout,*) ' virtual ITD mono-category param (1-4) or not (0) nn_monocat = ', nn_monocat
[8531]	344	WRITE(numout,*) ' Ice dynamics (T) or not (F) ln_icedyn = ', ln_icedyn
	345	WRITE(numout,*) ' Ice thermodynamics (T) or not (F) ln_icethd = ', ln_icethd
[8486]	346	WRITE(numout,*) ' maximum ice concentration for NH = ', rn_amax_n
	347	WRITE(numout,*) ' maximum ice concentration for SH = ', rn_amax_s
[8321]	348	ENDIF
	349	!
[8517]	350	! !--- check consistency
[8486]	351	IF ( jpl > 1 .AND. nn_monocat == 1 ) THEN
[8321]	352	nn_monocat = 0
	353	IF(lwp) WRITE(numout,*)
	354	IF(lwp) WRITE(numout,*) ' nn_monocat forced to 0 as jpl>1, i.e. multi-category case is chosen'
	355	ENDIF
[8752]	356	! IF ( jpl == 1 .AND. nn_monocat == 0 ) THEN
	357	! CALL ctl_stop( 'STOP', 'par_init : if jpl=1 then nn_monocat should be between 1 and 4' )
	358	! ENDIF
[8321]	359	!
[8531]	360	IF( ln_bdy .AND. ln_icediachk ) CALL ctl_warn('par_init: online conservation check does not work with BDY')
[8321]	361	!
[8517]	362	rdt_ice = REAL(nn_fsbc) * rdt !--- sea-ice timestep and inverse
[8486]	363	r1_rdtice = 1._wp / rdt_ice
[8517]	364	IF( lwp ) WRITE(numout,*) ' ice timestep rdt_ice = ', rdt_ice
[8321]	365	!
[8517]	366	r1_nlay_i = 1._wp / REAL( nlay_i, wp ) !--- inverse of nlay_i and nlay_s
[8486]	367	r1_nlay_s = 1._wp / REAL( nlay_s, wp )
	368	!
[8531]	369	END SUBROUTINE par_init
[8321]	370
	371
[8498]	372	SUBROUTINE store_fields
[8321]	373	!!----------------------------------------------------------------------
[8498]	374	!! * ROUTINE store_fields *
[8321]	375	!!
	376	!! ** purpose : store ice variables at "before" time step
	377	!!----------------------------------------------------------------------
[8360]	378	INTEGER :: ji, jj, jl ! dummy loop index
	379	!!----------------------------------------------------------------------
[8321]	380	!
[8564]	381	a_i_b (:,:,:) = a_i (:,:,:) ! ice area
	382	v_i_b (:,:,:) = v_i (:,:,:) ! ice volume
	383	v_s_b (:,:,:) = v_s (:,:,:) ! snow volume
	384	sv_i_b(:,:,:) = sv_i(:,:,:) ! salt content
	385	oa_i_b(:,:,:) = oa_i(:,:,:) ! areal age content
	386	e_s_b (:,:,:,:) = e_s (:,:,:,:) ! snow thermal energy
	387	e_i_b (:,:,:,:) = e_i (:,:,:,:) ! ice thermal energy
[8515]	388	WHERE( a_i_b(:,:,:) >= epsi20 )
[8563]	389	h_i_b(:,:,:) = v_i_b (:,:,:) / a_i_b(:,:,:) ! ice thickness
	390	h_s_b(:,:,:) = v_s_b (:,:,:) / a_i_b(:,:,:) ! snw thickness
[8515]	391	ELSEWHERE
[8563]	392	h_i_b(:,:,:) = 0._wp
	393	h_s_b(:,:,:) = 0._wp
[8515]	394	END WHERE
[8321]	395
[8360]	396	! ice velocities & total concentration
[8892]	397
	398
	399	write(numout,*) "RSRH store_fields a_i_b sum", glob_sum(a_i_b(:,:,:));flush(numout)
	400
[8515]	401	at_i_b(:,:) = SUM( a_i_b(:,:,:), dim=3 )
	402	u_ice_b(:,:) = u_ice(:,:)
	403	v_ice_b(:,:) = v_ice(:,:)
[8486]	404	!
[8498]	405	END SUBROUTINE store_fields
[8321]	406
	407
[8531]	408	SUBROUTINE diag_set0
[8321]	409	!!----------------------------------------------------------------------
[8531]	410	!! * ROUTINE diag_set0 *
[8321]	411	!!
	412	!! ** purpose : set ice-ocean and ice-atm. fluxes to zeros at the beggining
	413	!! of the time step
	414	!!----------------------------------------------------------------------
[8360]	415	INTEGER :: ji, jj ! dummy loop index
	416	!!----------------------------------------------------------------------
[8515]	417	sfx (:,:) = 0._wp ;
	418	sfx_bri(:,:) = 0._wp ; sfx_lam(:,:) = 0._wp
	419	sfx_sni(:,:) = 0._wp ; sfx_opw(:,:) = 0._wp
	420	sfx_bog(:,:) = 0._wp ; sfx_dyn(:,:) = 0._wp
	421	sfx_bom(:,:) = 0._wp ; sfx_sum(:,:) = 0._wp
	422	sfx_res(:,:) = 0._wp ; sfx_sub(:,:) = 0._wp
	423	!
	424	wfx_snw(:,:) = 0._wp ; wfx_ice(:,:) = 0._wp
	425	wfx_sni(:,:) = 0._wp ; wfx_opw(:,:) = 0._wp
	426	wfx_bog(:,:) = 0._wp ; wfx_dyn(:,:) = 0._wp
	427	wfx_bom(:,:) = 0._wp ; wfx_sum(:,:) = 0._wp
	428	wfx_res(:,:) = 0._wp ; wfx_sub(:,:) = 0._wp
	429	wfx_spr(:,:) = 0._wp ; wfx_lam(:,:) = 0._wp
	430	wfx_snw_dyn(:,:) = 0._wp ; wfx_snw_sum(:,:) = 0._wp
	431	wfx_snw_sub(:,:) = 0._wp ; wfx_ice_sub(:,:) = 0._wp
	432	wfx_snw_sni(:,:) = 0._wp
	433	wfx_pnd(:,:) = 0._wp
	434
	435	hfx_thd(:,:) = 0._wp ;
	436	hfx_snw(:,:) = 0._wp ; hfx_opw(:,:) = 0._wp
	437	hfx_bog(:,:) = 0._wp ; hfx_dyn(:,:) = 0._wp
	438	hfx_bom(:,:) = 0._wp ; hfx_sum(:,:) = 0._wp
	439	hfx_res(:,:) = 0._wp ; hfx_sub(:,:) = 0._wp
	440	hfx_spr(:,:) = 0._wp ; hfx_dif(:,:) = 0._wp
[8518]	441	hfx_err_rem(:,:) = 0._wp
[8515]	442	hfx_err_dif(:,:) = 0._wp
	443	wfx_err_sub(:,:) = 0._wp
	444	!
	445	afx_tot(:,:) = 0._wp ;
	446	!
[8564]	447	diag_heat(:,:) = 0._wp ; diag_sice(:,:) = 0._wp
[8515]	448	diag_vice(:,:) = 0._wp ; diag_vsnw(:,:) = 0._wp
	449
	450	! SIMIP diagnostics
	451	diag_fc_bo(:,:) = 0._wp ; diag_fc_su(:,:) = 0._wp
	452
	453	tau_icebfr(:,:) = 0._wp; ! landfast ice param only (clem: important to keep the init here)
[8782]	454	cnd_ice(:,:,:) = 0._wp ! initialisation of the effective conductivity at the top of ice/snow (Jules coupling)
[8321]	455
[8531]	456	END SUBROUTINE diag_set0
[8321]	457
	458	#else
	459	!!----------------------------------------------------------------------
[8534]	460	!! Default option Dummy module NO ESIM sea-ice model
[8321]	461	!!----------------------------------------------------------------------
	462	CONTAINS
	463	SUBROUTINE ice_stp ( kt, ksbc ) ! Dummy routine
	464	INTEGER, INTENT(in) :: kt, ksbc
	465	WRITE(,) 'ice_stp: You should not have seen this print! error?', kt
	466	END SUBROUTINE ice_stp
	467	SUBROUTINE ice_init ! Dummy routine
	468	END SUBROUTINE ice_init
	469	#endif
	470
	471	!!======================================================================
	472	END MODULE icestp

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/dev_r8126_LIM3_couple/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90

Download in other formats: