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.

thd_ice.F90 in branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3 – NEMO

Context Navigation

source: branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/thd_ice.F90 @ 8373

Last change on this file since 8373 was 8371, checked in by clem, 7 years ago
minor updates
Property svn:keywords set to `Id`
File size: 13.9 KB

Rev	Line
[825]	1	MODULE thd_ice
	2	!!======================================================================
	3	!! * MODULE thd_ice *
	4	!! LIM sea-ice : Ice thermodynamics in 1D
	5	!!=====================================================================
[2715]	6	!! History : 3.0 ! 2002-11 (C. Ethe) F90: Free form and module
[825]	7	!!----------------------------------------------------------------------
[7646]	8	#if defined key_lim3
	9	!!----------------------------------------------------------------------
	10	!! 'key_lim3' LIM3 sea-ice model
	11	!!----------------------------------------------------------------------
[2715]	12	USE in_out_manager ! I/O manager
[3625]	13	USE lib_mpp ! MPP library
[8369]	14	USE ice, ONLY : nlay_i, nlay_s, jpl
[825]	15
	16	IMPLICIT NONE
	17	PRIVATE
	18
[2715]	19	PUBLIC thd_ice_alloc ! Routine called by nemogcm.F90
	20
[8369]	21	!!----------------------
	22	!! * 1D Module variables
	23	!!----------------------
[834]	24	!: In ice thermodynamics, to spare memory, the vectors are folded
[4872]	25	!: from 1D to 2D vectors. The following variables, with ending _1d
[834]	26	!: are the variables corresponding to 2d vectors
[825]	27
[8327]	28	INTEGER , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: idxice !: selected points for ice thermo
[8342]	29	INTEGER , PUBLIC :: nidx ! number of selected points
[825]	30
[4872]	31	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: qlead_1d
	32	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: ftr_ice_1d
	33	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: qsr_ice_1d
	34	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: fr1_i0_1d
	35	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: fr2_i0_1d
	36	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: qns_ice_1d
	37	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: t_bo_1d
[6416]	38	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: rn_amax_1d
[825]	39
[4688]	40	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_sum_1d
	41	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_bom_1d
	42	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_bog_1d
	43	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_dif_1d
	44	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_opw_1d
	45	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_snw_1d
	46	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_err_1d
	47	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_err_rem_1d
[5146]	48	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_err_dif_1d
[8326]	49	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_out_1d
[4688]	50
	51	! heat flux associated with ice-atmosphere mass exchange
	52	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_sub_1d
	53	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_spr_1d
	54
	55	! heat flux associated with ice-ocean mass exchange
	56	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_thd_1d
	57	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hfx_res_1d
	58
[8341]	59	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_snw_sni_1d
[8233]	60	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_snw_sum_1d
[4866]	61	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_sub_1d
[8239]	62	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_snw_sub_1d
	63	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_ice_sub_1d
[8326]	64	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_err_sub_1d
[8327]	65	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_lam_1d
[4688]	66
[4866]	67	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_bog_1d
	68	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_bom_1d
	69	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_sum_1d
	70	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_sni_1d
	71	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_opw_1d
	72	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_res_1d
	73	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: wfx_spr_1d
[4688]	74
[4866]	75	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sfx_bri_1d
	76	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sfx_bog_1d
	77	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sfx_bom_1d
	78	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sfx_sum_1d
	79	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sfx_sni_1d
	80	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sfx_opw_1d
	81	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sfx_res_1d
[6416]	82	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sfx_sub_1d
[8327]	83	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sfx_lam_1d
[6416]	84
[2715]	85	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sprecip_1d !: <==> the 2D sprecip
[4872]	86	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: at_i_1d !: <==> the 2D at_i
[4866]	87	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: fhtur_1d !: <==> the 2D fhtur
[4688]	88	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: fhld_1d !: <==> the 2D fhld
[2715]	89	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: dqns_ice_1d !: <==> the 2D dqns_ice
[5407]	90	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: evap_ice_1d !: <==> the 2D evap_ice
	91	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: qprec_ice_1d !: <==> the 2D qprec_ice
[2715]	92	! ! to reintegrate longwave flux inside the ice thermodynamics
	93	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: i0 !: fraction of radiation transmitted to the ice
[4872]	94	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: hicol_1d !: Ice collection thickness accumulated in leads
[825]	95
[4872]	96	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: t_su_1d !: <==> the 2D t_su
[8239]	97	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: t_si_1d !: <==> the 2D t_si
[4872]	98	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: a_i_1d !: <==> the 2D a_i
[8369]	99	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: a_ib_1d !: <==> the 2D a_i_b
	100	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: ht_i_1d !: <==> the 2D ht_i
	101	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: ht_ib_1d !: <==> the 2D ht_i_b
	102	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: ht_s_1d !: <==> the 2D ht_s
[4872]	103	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: fc_su !: Surface Conduction flux
	104	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: fc_bo_i !: Bottom Conduction flux
	105	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: dh_s_tot !: Snow accretion/ablation [m]
	106	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: dh_i_surf !: Ice surface accretion/ablation [m]
[6416]	107	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: dh_i_sub !: Ice surface sublimation [m]
[4872]	108	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: dh_i_bott !: Ice bottom accretion/ablation [m]
	109	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: dh_snowice !: Snow ice formation [m of ice]
	110	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sm_i_1d !: Ice bulk salinity [ppt]
	111	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: s_i_new !: Salinity of new ice at the bottom
[8371]	112	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: a_ip_1d !:
	113	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: v_ip_1d !:
[8369]	114	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: v_i_1d !:
[8371]	115	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: v_s_1d !:
	116	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: smv_i_1d !:
[825]	117
[4872]	118	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: t_s_1d !: corresponding to the 2D var t_s
	119	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: t_i_1d !: corresponding to the 2D var t_i
	120	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: s_i_1d !: profiled ice salinity
[8325]	121	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: e_i_1d !: Ice enthalpy per unit volume
	122	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: e_s_1d !: Snow enthalpy per unit volume
[825]	123
[8325]	124	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: eh_i_old !: ice heat content (q*h, J.m-2)
[4872]	125	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: h_i_old !: ice thickness layer (m)
[825]	126
[8239]	127	! Conduction flux diagnostics (SIMIP)
	128	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: diag_fc_bo_1d !: <==> the 2D diag_fc_bo
	129	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: diag_fc_su_1d !: <==> the 2D diag_fc_su
	130
[8326]	131	! surface fields from the ocean
	132	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sst_1d
	133	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sss_1d
	134
[8369]	135	!
	136	!!----------------------
	137	!! * 2D Module variables
	138	!!----------------------
	139	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: a_i_2d
	140	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: v_i_2d
	141	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: v_s_2d
	142	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: oa_i_2d
	143	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: smv_i_2d
	144	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: a_ip_2d
	145	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: v_ip_2d
	146	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: t_su_2d
	147	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ht_i_2d
	148
	149	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: a_ib_2d
	150	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ht_ib_2d
	151
[2715]	152	INTEGER , PUBLIC :: jiindex_1d ! 1D index of debugging point
[825]	153
[2715]	154	!!----------------------------------------------------------------------
[4161]	155	!! NEMO/LIM3 4.0 , UCL - NEMO Consortium (2011)
[2715]	156	!! $Id$
	157	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
	158	!!----------------------------------------------------------------------
	159	CONTAINS
	160
	161	FUNCTION thd_ice_alloc()
	162	!!---------------------------------------------------------------------!
	163	!! * ROUTINE thd_ice_alloc *
	164	!!---------------------------------------------------------------------!
	165	INTEGER :: thd_ice_alloc ! return value
[8369]	166	INTEGER :: ierr(7), ii
[2715]	167	!!---------------------------------------------------------------------!
[7646]	168	ierr(:) = 0
[2715]	169
[7646]	170	ii = 1
[8327]	171	ALLOCATE( idxice (jpij) , &
[5167]	172	& qlead_1d (jpij) , ftr_ice_1d(jpij) , qsr_ice_1d (jpij) , &
	173	& fr1_i0_1d(jpij) , fr2_i0_1d (jpij) , qns_ice_1d(jpij) , &
[5146]	174	& t_bo_1d (jpij) , &
	175	& hfx_sum_1d(jpij) , hfx_bom_1d(jpij) ,hfx_bog_1d(jpij) , &
	176	& hfx_dif_1d(jpij) , hfx_opw_1d(jpij) , &
[6403]	177	& rn_amax_1d(jpij) , &
[5146]	178	& hfx_thd_1d(jpij) , hfx_spr_1d(jpij) , &
	179	& hfx_snw_1d(jpij) , hfx_sub_1d(jpij) , hfx_err_1d(jpij) , &
[8326]	180	& hfx_res_1d(jpij) , hfx_err_rem_1d(jpij) , hfx_err_dif_1d(jpij) , hfx_out_1d(jpij), STAT=ierr(ii) )
[2715]	181	!
[7646]	182	ii = ii + 1
[8313]	183	ALLOCATE( sprecip_1d (jpij) , at_i_1d (jpij) , &
[8341]	184	& fhtur_1d (jpij) , wfx_snw_sni_1d (jpij) , wfx_spr_1d (jpij) , wfx_snw_sum_1d(jpij) , &
[5167]	185	& fhld_1d (jpij) , wfx_sub_1d (jpij) , wfx_bog_1d (jpij) , wfx_bom_1d(jpij) , &
	186	& wfx_sum_1d(jpij) , wfx_sni_1d (jpij) , wfx_opw_1d (jpij) , wfx_res_1d(jpij) , &
[8326]	187	& wfx_snw_sub_1d(jpij), wfx_ice_sub_1d(jpij), wfx_err_sub_1d(jpij) , &
[8327]	188	& wfx_lam_1d(jpij) , dqns_ice_1d(jpij) , evap_ice_1d (jpij), &
[8369]	189	& qprec_ice_1d(jpij), i0 (jpij) , &
[5167]	190	& sfx_bri_1d (jpij) , sfx_bog_1d (jpij) , sfx_bom_1d (jpij) , sfx_sum_1d (jpij), &
[6416]	191	& sfx_sni_1d (jpij) , sfx_opw_1d (jpij) , sfx_res_1d (jpij) , sfx_sub_1d (jpij), &
[8327]	192	& sfx_lam_1d (jpij) , hicol_1d (jpij) , STAT=ierr(ii) )
[2715]	193	!
[7646]	194	ii = ii + 1
[8369]	195	ALLOCATE( t_su_1d (jpij) , t_si_1d (jpij) , a_i_1d (jpij) , a_ib_1d(jpij) , &
	196	& ht_i_1d (jpij) , ht_ib_1d (jpij) , ht_s_1d (jpij) , fc_su (jpij) , fc_bo_i(jpij) , &
	197	& dh_s_tot (jpij) , dh_i_surf (jpij) , dh_i_sub(jpij) , &
	198	& dh_i_bott(jpij) , dh_snowice(jpij) , sm_i_1d (jpij) , s_i_new(jpij) , &
[8371]	199	& a_ip_1d (jpij) , v_ip_1d (jpij) , v_i_1d (jpij) , v_s_1d (jpij) , smv_i_1d(jpij) , STAT=ierr(ii) )
[2715]	200	!
[7646]	201	ii = ii + 1
	202	ALLOCATE( t_s_1d (jpij,nlay_s) , t_i_1d (jpij,nlay_i) , s_i_1d(jpij,nlay_i) , &
[8325]	203	& e_i_1d (jpij,nlay_i+1) , e_s_1d (jpij,nlay_s) , &
	204	& eh_i_old(jpij,0:nlay_i+1) , h_i_old(jpij,0:nlay_i+1) , STAT=ierr(ii) )
[7646]	205	!
[8239]	206	ii = ii + 1
	207	ALLOCATE( diag_fc_bo_1d(jpij) , diag_fc_su_1d(jpij) , STAT=ierr(ii) )
[8326]	208	!
	209	ii = ii + 1
	210	ALLOCATE( sst_1d(jpij) , sss_1d(jpij) , STAT=ierr(ii) )
[8369]	211	!
	212	ii = ii + 1
	213	ALLOCATE( a_i_2d(jpij,jpl) , a_ib_2d(jpij,jpl) , ht_i_2d(jpij,jpl) , ht_ib_2d(jpij,jpl) , &
	214	& v_i_2d(jpij,jpl) ,v_s_2d(jpij,jpl) ,oa_i_2d(jpij,jpl) ,smv_i_2d(jpij,jpl) , &
	215	& a_ip_2d(jpij,jpl) ,v_ip_2d(jpij,jpl) ,t_su_2d(jpij,jpl) , &
	216	& STAT=ierr(ii) )
[8239]	217
[7646]	218	thd_ice_alloc = MAXVAL( ierr(:) )
[2715]	219	IF( thd_ice_alloc /= 0 ) CALL ctl_warn( 'thd_ice_alloc: failed to allocate arrays.' )
	220	!
	221	END FUNCTION thd_ice_alloc
	222
[7646]	223	#else
	224	!!----------------------------------------------------------------------
	225	!! Default option : Empty module NO LIM sea-ice model
	226	!!----------------------------------------------------------------------
	227	CONTAINS
	228	SUBROUTINE thd_ice_alloc ! Empty routine
	229	END SUBROUTINE thd_ice_alloc
	230	#endif
	231
[825]	232	!!======================================================================
	233	END MODULE thd_ice

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

Download in other formats: