Context Navigation

sbc_ice.F90 @ 15720

Last change on this file since 15720 was 14072, checked in by laurent, 4 years ago
Merging branch "2020/dev_r13648_ASINTER-04_laurent_bulk_ice", ticket #2369
Property svn:keywords set to `Id`
File size: 12.6 KB

Rev	Line
[888]	1	MODULE sbc_ice
	2	!!======================================================================
	3	!! * MODULE sbc_ice *
[9656]	4	!! Surface module - SI3 & CICE: parameters & variables defined in memory
[888]	5	!!======================================================================
[9656]	6	!! History : 3.0 ! 2006-08 (G. Madec) Surface module
	7	!! 3.2 ! 2009-06 (S. Masson) merge with ice_oce
	8	!! 3.3.1 ! 2011-01 (A. R. Porter, STFC Daresbury) dynamical allocation
	9	!! 3.4 ! 2011-11 (C. Harris) CICE added as an option
	10	!! 4.0 ! 2018 (many people) SI3 compatibility
[888]	11	!!----------------------------------------------------------------------
[9570]	12	#if defined key_si3 \|\| defined key_cice
[888]	13	!!----------------------------------------------------------------------
[9654]	14	!! 'key_si3' or 'key_cice' : SI3 or CICE sea-ice model
[888]	15	!!----------------------------------------------------------------------
	16	USE par_oce ! ocean parameters
[4990]	17	USE sbc_oce ! surface boundary condition: ocean
[9570]	18	# if defined key_si3
[9656]	19	USE ice ! SI3 parameters
[1465]	20	# endif
[4990]	21	# if defined key_cice
[14072]	22	USE ice_domain_size, only: ncat
[3294]	23	#endif
[2715]	24	USE lib_mpp ! MPP library
	25	USE in_out_manager ! I/O manager
[888]	26
	27	IMPLICIT NONE
	28	PRIVATE
	29
[9019]	30	PUBLIC sbc_ice_alloc ! called in sbcmod.F90 or sbcice_cice.F90
[2715]	31
[9570]	32	# if defined key_si3
[9656]	33	LOGICAL , PUBLIC, PARAMETER :: lk_si3 = .TRUE. !: SI3 ice model
[14072]	34	LOGICAL , PUBLIC, PARAMETER :: lk_cice = .FALSE. !: no CICE
[1465]	35	# endif
[3294]	36	# if defined key_cice
[9656]	37	LOGICAL , PUBLIC, PARAMETER :: lk_si3 = .FALSE. !: no SI3
[3294]	38	LOGICAL , PUBLIC, PARAMETER :: lk_cice = .TRUE. !: CICE ice model
	39	# endif
[1465]	40
[4990]	41	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qns_ice !: non solar heat flux over ice [W/m2]
	42	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qsr_ice !: solar heat flux over ice [W/m2]
	43	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qla_ice !: latent flux over ice [W/m2]
	44	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: dqla_ice !: latent sensibility over ice [W/m2/K]
	45	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: dqns_ice !: non solar heat flux over ice (LW+SEN+LA) [W/m2/K]
	46	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: tn_ice !: ice surface temperature [K]
	47	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: alb_ice !: ice albedo [-]
[888]	48
[14072]	49	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qml_ice !: heat available for snow / ice surface melting [W/m2]
	50	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qcn_ice !: heat conduction flux in the layer below surface [W/m2]
[9910]	51	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qtr_ice_top !: solar flux transmitted below the ice surface [W/m2]
[9019]	52
[4990]	53	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: utau_ice !: atmos-ice u-stress. VP: I-pt ; EVP: U,V-pts [N/m2]
	54	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: vtau_ice !: atmos-ice v-stress. VP: I-pt ; EVP: U,V-pts [N/m2]
[5407]	55	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: emp_ice !: sublimation - precip over sea ice [kg/m2/s]
[888]	56
[4990]	57	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: topmelt !: category topmelt
	58	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: botmelt !: category botmelt
[1465]	59
[9570]	60	#if defined key_si3
[5407]	61	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: evap_ice !: sublimation [kg/m2/s]
	62	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: devap_ice !: sublimation sensitivity [kg/m2/s/K]
	63	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qns_oce !: non solar heat flux over ocean [W/m2]
	64	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qsr_oce !: non solar heat flux over ocean [W/m2]
	65	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qemp_oce !: heat flux of precip and evap over ocean [W/m2]
	66	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qemp_ice !: heat flux of precip and evap over ice [W/m2]
[6416]	67	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qevap_ice !: heat flux of evap over ice [W/m2]
	68	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qprec_ice !: enthalpy of precip over ice [J/m3]
[5407]	69	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: emp_oce !: evap - precip over ocean [kg/m2/s]
	70	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: wndm_ice !: wind speed module at T-point [m/s]
[13472]	71	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: sstfrz !: sea surface freezing temperature [degC]
	72	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: rCdU_ice !: ice-ocean drag at T-point (<0) [m/s]
[5407]	73	#endif
	74
[4990]	75	#if defined key_cice
[3294]	76	!
[9656]	77	! for consistency with SI3, these are declared with three dimensions
[3294]	78	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qlw_ice !: incoming long-wave
	79	!
	80	! other forcing arrays are two dimensional
	81	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ss_iou !: x ice-ocean surface stress at NEMO U point
	82	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ss_iov !: y ice-ocean surface stress at NEMO V point
	83	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qatm_ice !: specific humidity
	84	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: wndi_ice !: i wind at T point
	85	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: wndj_ice !: j wind at T point
	86	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: nfrzmlt !: NEMO frzmlt
	87	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fr_iu !: ice fraction at NEMO U point
	88	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fr_iv !: ice fraction at NEMO V point
[14072]	89
[4990]	90	! variables used in the coupled interface
	91	INTEGER , PUBLIC, PARAMETER :: jpl = ncat
[14072]	92	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: u_ice, v_ice
	93
[9656]	94	! already defined in ice.F90 for SI3
[4990]	95	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: a_i
[9019]	96	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: h_i, h_s
[3294]	97
[4990]	98	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tatm_ice !: air temperature [K]
	99	#endif
	100
[9019]	101	!! arrays relating to embedding ice in the ocean
	102	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_mass !: mass of snow and ice at current ice time step [Kg/m2]
	103	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_mass_b !: mass of snow and ice at previous ice time step [Kg/m2]
	104	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_fmass !: time evolution of mass of snow+ice [Kg/m2/s]
	105
[2715]	106	!!----------------------------------------------------------------------
[9598]	107	!! NEMO/OCE 4.0 , NEMO Consortium (2018)
[14072]	108	!! $Id$
[10068]	109	!! Software governed by the CeCILL license (see ./LICENSE)
[2715]	110	!!----------------------------------------------------------------------
	111	CONTAINS
	112
	113	INTEGER FUNCTION sbc_ice_alloc()
	114	!!----------------------------------------------------------------------
	115	!! * FUNCTION sbc_ice_alloc *
	116	!!----------------------------------------------------------------------
[9019]	117	INTEGER :: ierr(4)
[4306]	118	!!----------------------------------------------------------------------
	119	ierr(:) = 0
	120
[9019]	121	ALLOCATE( snwice_mass(jpi,jpj) , snwice_mass_b(jpi,jpj), snwice_fmass(jpi,jpj) , STAT=ierr(1) )
	122
[9570]	123	#if defined key_si3
[9019]	124	ALLOCATE( qns_ice (jpi,jpj,jpl) , qsr_ice (jpi,jpj,jpl) , &
	125	& qla_ice (jpi,jpj,jpl) , dqla_ice (jpi,jpj,jpl) , &
[9910]	126	& dqns_ice(jpi,jpj,jpl) , tn_ice (jpi,jpj,jpl) , alb_ice (jpi,jpj,jpl) , &
	127	& qml_ice (jpi,jpj,jpl) , qcn_ice (jpi,jpj,jpl) , qtr_ice_top(jpi,jpj,jpl) , &
	128	& utau_ice(jpi,jpj) , vtau_ice (jpi,jpj) , wndm_ice (jpi,jpj) , &
	129	& evap_ice(jpi,jpj,jpl) , devap_ice(jpi,jpj,jpl) , qprec_ice (jpi,jpj) , &
	130	& qemp_ice(jpi,jpj) , qevap_ice(jpi,jpj,jpl) , qemp_oce (jpi,jpj) , &
	131	& qns_oce (jpi,jpj) , qsr_oce (jpi,jpj) , emp_oce (jpi,jpj) , &
[13472]	132	& emp_ice (jpi,jpj) , sstfrz (jpi,jpj) , rCdU_ice (jpi,jpj) , STAT= ierr(2) )
[5407]	133	#endif
[4990]	134
	135	#if defined key_cice
[3294]	136	ALLOCATE( qla_ice(jpi,jpj,1) , qlw_ice(jpi,jpj,1) , qsr_ice(jpi,jpj,1) , &
	137	wndi_ice(jpi,jpj) , tatm_ice(jpi,jpj) , qatm_ice(jpi,jpj) , &
	138	wndj_ice(jpi,jpj) , nfrzmlt(jpi,jpj) , ss_iou(jpi,jpj) , &
	139	ss_iov(jpi,jpj) , fr_iu(jpi,jpj) , fr_iv(jpi,jpj) , &
[4990]	140	a_i(jpi,jpj,ncat) , topmelt(jpi,jpj,ncat) , botmelt(jpi,jpj,ncat) , &
[9019]	141	STAT= ierr(2) )
	142	IF( ln_cpl ) ALLOCATE( u_ice(jpi,jpj) , tn_ice (jpi,jpj,1) , &
	143	& v_ice(jpi,jpj) , alb_ice(jpi,jpj,1) , &
[4990]	144	& emp_ice(jpi,jpj) , qns_ice(jpi,jpj,1) , dqns_ice(jpi,jpj,1) , &
[14072]	145	& STAT= ierr(3) )
[9019]	146	IF( ln_cpl ) ALLOCATE( h_i(jpi,jpj,jpl) , h_s(jpi,jpj,jpl) , STAT=ierr(4) )
[3294]	147	#endif
[4990]	148
[4306]	149	sbc_ice_alloc = MAXVAL( ierr )
[10425]	150	CALL mpp_sum ( 'sbc_ice', sbc_ice_alloc )
[2715]	151	IF( sbc_ice_alloc > 0 ) CALL ctl_warn('sbc_ice_alloc: allocation of arrays failed')
	152	END FUNCTION sbc_ice_alloc
	153
[888]	154	#else
	155	!!----------------------------------------------------------------------
[9654]	156	!! Default option NO SI3 or CICE sea-ice model
[888]	157	!!----------------------------------------------------------------------
[9019]	158	USE lib_mpp ! MPP library
[4990]	159	USE in_out_manager ! I/O manager
[9019]	160
	161	IMPLICIT NONE
	162	PRIVATE
	163
	164	PUBLIC sbc_ice_alloc !
	165
[9656]	166	LOGICAL , PUBLIC, PARAMETER :: lk_si3 = .FALSE. !: no SI3 ice model
	167	LOGICAL , PUBLIC, PARAMETER :: lk_cice = .FALSE. !: no CICE ice model
[14072]	168
	169	INTEGER , PUBLIC, PARAMETER :: jpl = 1
[9019]	170	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: u_ice, v_ice ! jpi, jpj
[4990]	171	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: tn_ice, alb_ice, qns_ice, dqns_ice ! (jpi,jpj,jpl)
	172	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: a_i
	173	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: emp_ice
	174	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qsr_ice
[9019]	175	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: h_i, h_s
[4990]	176	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: topmelt, botmelt
[9019]	177	!
[14072]	178	!! arrays related to embedding ice in the ocean.
	179	!! These arrays need to be declared even if no ice model is required.
[9019]	180	!! In the no ice model or traditional levitating ice cases they contain only zeros
	181	!! ---------------------
	182	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_mass !: mass of snow and ice at current ice time step [Kg/m2]
	183	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_mass_b !: mass of snow and ice at previous ice time step [Kg/m2]
	184	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_fmass !: time evolution of mass of snow+ice [Kg/m2/s]
	185	!!----------------------------------------------------------------------
	186	CONTAINS
	187
	188	INTEGER FUNCTION sbc_ice_alloc()
	189	!!----------------------------------------------------------------------
	190	!! * FUNCTION sbc_ice_alloc *
	191	!!----------------------------------------------------------------------
	192	INTEGER :: ierr(1)
	193	!!----------------------------------------------------------------------
	194	ierr(:) = 0
	195	ALLOCATE( snwice_mass(jpi,jpj) , snwice_mass_b(jpi,jpj), snwice_fmass(jpi,jpj) , STAT=ierr(1) )
	196	sbc_ice_alloc = MAXVAL( ierr )
[10425]	197	CALL mpp_sum ( 'sbc_ice', sbc_ice_alloc )
[9019]	198	IF( sbc_ice_alloc > 0 ) CALL ctl_warn('sbc_ice_alloc: allocation of arrays failed')
	199	END FUNCTION sbc_ice_alloc
[888]	200	#endif
	201
[1469]	202	!!======================================================================
[888]	203	END MODULE sbc_ice

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

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

Context Navigation

source: NEMO/trunk/src/OCE/SBC/sbc_ice.F90 @ 15720

Download in other formats: