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 @ 9654

Last change on this file since 9654 was 9654, checked in by clem, 6 years ago
remove lk_lim3 from the code and replace by lk_si3
Property svn:keywords set to `Id`
File size: 13.0 KB

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

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

Download in other formats: