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.

sbc_ice.F90 in NEMO/trunk/src/OCE/SBC – NEMO

Context Navigation

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

Last change on this file since 9656 was 9656, checked in by clem, 6 years ago
remove the remaining references to LIM
Property svn:keywords set to `Id`
File size: 13.1 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
[3294]	22	USE ice_domain_size, only: ncat
	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
[3294]	34	LOGICAL , PUBLIC, PARAMETER :: lk_cice = .FALSE. !: no CICE
[2528]	35	CHARACTER(len=1), PUBLIC, PARAMETER :: cp_ice_msh = 'C' !: 'C'-grid ice-velocity
[1465]	36	# endif
[3294]	37	# if defined key_cice
[9656]	38	LOGICAL , PUBLIC, PARAMETER :: lk_si3 = .FALSE. !: no SI3
[3294]	39	LOGICAL , PUBLIC, PARAMETER :: lk_cice = .TRUE. !: CICE ice model
	40	CHARACTER(len=1), PUBLIC :: cp_ice_msh = 'F' !: 'F'-grid ice-velocity
	41	# endif
[1465]	42
[4990]	43	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qns_ice !: non solar heat flux over ice [W/m2]
	44	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qsr_ice !: solar heat flux over ice [W/m2]
	45	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qla_ice !: latent flux over ice [W/m2]
	46	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: dqla_ice !: latent sensibility over ice [W/m2/K]
	47	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: dqns_ice !: non solar heat flux over ice (LW+SEN+LA) [W/m2/K]
	48	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: tn_ice !: ice surface temperature [K]
	49	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: alb_ice !: ice albedo [-]
[888]	50
[9019]	51	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qml_ice !: heat available for snow / ice surface melting [W/m2]
	52	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qcn_ice !: heat conduction flux in the layer below surface [W/m2]
	53	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qsr_ice_tr !: solar flux transmitted below the ice surface [W/m2]
	54
[4990]	55	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: utau_ice !: atmos-ice u-stress. VP: I-pt ; EVP: U,V-pts [N/m2]
	56	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: vtau_ice !: atmos-ice v-stress. VP: I-pt ; EVP: U,V-pts [N/m2]
[5407]	57	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: emp_ice !: sublimation - precip over sea ice [kg/m2/s]
[888]	58
[4990]	59	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: topmelt !: category topmelt
	60	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: botmelt !: category botmelt
[1465]	61
[9570]	62	#if defined key_si3
[5407]	63	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: evap_ice !: sublimation [kg/m2/s]
	64	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: devap_ice !: sublimation sensitivity [kg/m2/s/K]
	65	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qns_oce !: non solar heat flux over ocean [W/m2]
	66	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qsr_oce !: non solar heat flux over ocean [W/m2]
	67	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qemp_oce !: heat flux of precip and evap over ocean [W/m2]
	68	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qemp_ice !: heat flux of precip and evap over ice [W/m2]
[6416]	69	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qevap_ice !: heat flux of evap over ice [W/m2]
	70	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qprec_ice !: enthalpy of precip over ice [J/m3]
[5407]	71	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: emp_oce !: evap - precip over ocean [kg/m2/s]
	72	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: wndm_ice !: wind speed module at T-point [m/s]
[9019]	73	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: sstfrz !: wind speed module at T-point [m/s]
	74	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: tsfc_ice !: sea ice surface skin temperature (on categories)
[5407]	75	#endif
	76
[4990]	77	#if defined key_cice
[3294]	78	!
[9656]	79	! for consistency with SI3, these are declared with three dimensions
[3294]	80	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qlw_ice !: incoming long-wave
	81	!
	82	! other forcing arrays are two dimensional
	83	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ss_iou !: x ice-ocean surface stress at NEMO U point
	84	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ss_iov !: y ice-ocean surface stress at NEMO V point
	85	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qatm_ice !: specific humidity
	86	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: wndi_ice !: i wind at T point
	87	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: wndj_ice !: j wind at T point
	88	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: nfrzmlt !: NEMO frzmlt
	89	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fr_iu !: ice fraction at NEMO U point
	90	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fr_iv !: ice fraction at NEMO V point
[4990]	91
	92	! variables used in the coupled interface
	93	INTEGER , PUBLIC, PARAMETER :: jpl = ncat
	94	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: u_ice, v_ice ! jpi, jpj
	95
[9656]	96	! already defined in ice.F90 for SI3
[4990]	97	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: a_i
[9019]	98	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: h_i, h_s
[3294]	99
[4990]	100	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tatm_ice !: air temperature [K]
	101	#endif
	102
	103	REAL(wp), PUBLIC, SAVE :: cldf_ice = 0.81 !: cloud fraction over sea ice, summer CLIO value [-]
	104
[9019]	105	!! arrays relating to embedding ice in the ocean
	106	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_mass !: mass of snow and ice at current ice time step [Kg/m2]
	107	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_mass_b !: mass of snow and ice at previous ice time step [Kg/m2]
	108	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_fmass !: time evolution of mass of snow+ice [Kg/m2/s]
	109
[2715]	110	!!----------------------------------------------------------------------
[9598]	111	!! NEMO/OCE 4.0 , NEMO Consortium (2018)
[2715]	112	!! $Id$
[9598]	113	!! Software governed by the CeCILL licence (./LICENSE)
[2715]	114	!!----------------------------------------------------------------------
	115	CONTAINS
	116
	117	INTEGER FUNCTION sbc_ice_alloc()
	118	!!----------------------------------------------------------------------
	119	!! * FUNCTION sbc_ice_alloc *
	120	!!----------------------------------------------------------------------
[9019]	121	INTEGER :: ierr(4)
[4306]	122	!!----------------------------------------------------------------------
	123	ierr(:) = 0
	124
[9019]	125	ALLOCATE( snwice_mass(jpi,jpj) , snwice_mass_b(jpi,jpj), snwice_fmass(jpi,jpj) , STAT=ierr(1) )
	126
[9570]	127	#if defined key_si3
[9019]	128	ALLOCATE( qns_ice (jpi,jpj,jpl) , qsr_ice (jpi,jpj,jpl) , &
	129	& qla_ice (jpi,jpj,jpl) , dqla_ice (jpi,jpj,jpl) , &
	130	& dqns_ice(jpi,jpj,jpl) , tn_ice (jpi,jpj,jpl) , alb_ice (jpi,jpj,jpl) , &
	131	& qml_ice (jpi,jpj,jpl) , qcn_ice (jpi,jpj,jpl) , qsr_ice_tr(jpi,jpj,jpl) , &
	132	& utau_ice(jpi,jpj) , vtau_ice (jpi,jpj) , wndm_ice (jpi,jpj) , &
	133	& evap_ice(jpi,jpj,jpl) , devap_ice(jpi,jpj,jpl) , qprec_ice (jpi,jpj) , &
	134	& qemp_ice(jpi,jpj) , qevap_ice(jpi,jpj,jpl) , qemp_oce (jpi,jpj) , &
	135	& qns_oce (jpi,jpj) , qsr_oce (jpi,jpj) , emp_oce (jpi,jpj) , &
	136	& emp_ice (jpi,jpj) , tsfc_ice (jpi,jpj,jpl) , sstfrz (jpi,jpj) , STAT= ierr(2) )
[5407]	137	#endif
[4990]	138
	139	#if defined key_cice
[3294]	140	ALLOCATE( qla_ice(jpi,jpj,1) , qlw_ice(jpi,jpj,1) , qsr_ice(jpi,jpj,1) , &
	141	wndi_ice(jpi,jpj) , tatm_ice(jpi,jpj) , qatm_ice(jpi,jpj) , &
	142	wndj_ice(jpi,jpj) , nfrzmlt(jpi,jpj) , ss_iou(jpi,jpj) , &
	143	ss_iov(jpi,jpj) , fr_iu(jpi,jpj) , fr_iv(jpi,jpj) , &
[4990]	144	a_i(jpi,jpj,ncat) , topmelt(jpi,jpj,ncat) , botmelt(jpi,jpj,ncat) , &
[9019]	145	STAT= ierr(2) )
	146	IF( ln_cpl ) ALLOCATE( u_ice(jpi,jpj) , tn_ice (jpi,jpj,1) , &
	147	& v_ice(jpi,jpj) , alb_ice(jpi,jpj,1) , &
[4990]	148	& emp_ice(jpi,jpj) , qns_ice(jpi,jpj,1) , dqns_ice(jpi,jpj,1) , &
[9019]	149	& STAT= ierr(3) )
	150	IF( ln_cpl ) ALLOCATE( h_i(jpi,jpj,jpl) , h_s(jpi,jpj,jpl) , STAT=ierr(4) )
[3294]	151	#endif
[4990]	152
[4306]	153	sbc_ice_alloc = MAXVAL( ierr )
[2715]	154	IF( lk_mpp ) CALL mpp_sum ( sbc_ice_alloc )
	155	IF( sbc_ice_alloc > 0 ) CALL ctl_warn('sbc_ice_alloc: allocation of arrays failed')
	156	END FUNCTION sbc_ice_alloc
	157
[888]	158	#else
	159	!!----------------------------------------------------------------------
[9654]	160	!! Default option NO SI3 or CICE sea-ice model
[888]	161	!!----------------------------------------------------------------------
[9019]	162	USE lib_mpp ! MPP library
[4990]	163	USE in_out_manager ! I/O manager
[9019]	164
	165	IMPLICIT NONE
	166	PRIVATE
	167
	168	PUBLIC sbc_ice_alloc !
	169
[9656]	170	LOGICAL , PUBLIC, PARAMETER :: lk_si3 = .FALSE. !: no SI3 ice model
	171	LOGICAL , PUBLIC, PARAMETER :: lk_cice = .FALSE. !: no CICE ice model
[2528]	172	CHARACTER(len=1), PUBLIC, PARAMETER :: cp_ice_msh = '-' !: no grid ice-velocity
[9019]	173	REAL(wp) , PUBLIC, PARAMETER :: cldf_ice = 0.81 !: cloud fraction over sea ice, summer CLIO value [-]
[4990]	174	INTEGER , PUBLIC, PARAMETER :: jpl = 1
[9019]	175	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: u_ice, v_ice ! jpi, jpj
[4990]	176	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: tn_ice, alb_ice, qns_ice, dqns_ice ! (jpi,jpj,jpl)
	177	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: a_i
	178	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: emp_ice
	179	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qsr_ice
[9019]	180	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: h_i, h_s
[4990]	181	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: topmelt, botmelt
[9019]	182	!
	183	!! arrays related to embedding ice in the ocean.
	184	!! These arrays need to be declared even if no ice model is required.
	185	!! In the no ice model or traditional levitating ice cases they contain only zeros
	186	!! ---------------------
	187	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_mass !: mass of snow and ice at current ice time step [Kg/m2]
	188	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_mass_b !: mass of snow and ice at previous ice time step [Kg/m2]
	189	REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: snwice_fmass !: time evolution of mass of snow+ice [Kg/m2/s]
	190	!!----------------------------------------------------------------------
	191	CONTAINS
	192
	193	INTEGER FUNCTION sbc_ice_alloc()
	194	!!----------------------------------------------------------------------
	195	!! * FUNCTION sbc_ice_alloc *
	196	!!----------------------------------------------------------------------
	197	INTEGER :: ierr(1)
	198	!!----------------------------------------------------------------------
	199	ierr(:) = 0
	200	ALLOCATE( snwice_mass(jpi,jpj) , snwice_mass_b(jpi,jpj), snwice_fmass(jpi,jpj) , STAT=ierr(1) )
	201	sbc_ice_alloc = MAXVAL( ierr )
	202	IF( lk_mpp ) CALL mpp_sum ( sbc_ice_alloc )
	203	IF( sbc_ice_alloc > 0 ) CALL ctl_warn('sbc_ice_alloc: allocation of arrays failed')
	204	END FUNCTION sbc_ice_alloc
[888]	205	#endif
	206
[1469]	207	!!======================================================================
[888]	208	END MODULE sbc_ice

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

Download in other formats: