subroutine sfc_sice::sfc_sice_run	(	integer, intent(in)	im,
		integer, intent(in)	km,
		real (kind=kind_phys), dimension(im), intent(in)	ps,
		real (kind=kind_phys), dimension(im), intent(in)	u1,
		real (kind=kind_phys), dimension(im), intent(in)	v1,
		real (kind=kind_phys), dimension(im), intent(in)	t1,
		real (kind=kind_phys), dimension(im), intent(in)	q1,
		real (kind=kind_phys), intent(in)	delt,
		real (kind=kind_phys), dimension(im), intent(in)	sfcemis,
		real (kind=kind_phys), dimension(im), intent(in)	dlwflx,
		real (kind=kind_phys), dimension(im), intent(in)	sfcnsw,
		real (kind=kind_phys), dimension(im), intent(in)	sfcdsw,
		real (kind=kind_phys), dimension(im), intent(in)	srflag,
		real (kind=kind_phys), dimension(im), intent(in)	cm,
		real (kind=kind_phys), dimension(im), intent(in)	ch,
		real (kind=kind_phys), dimension(im), intent(in)	prsl1,
		real (kind=kind_phys), dimension(im), intent(in)	prslki,
		integer, dimension(im), intent(in)	islimsk,
		real (kind=kind_phys), dimension(im), intent(in)	ddvel,
		logical, dimension(im), intent(in)	flag_iter,
		logical, intent(in)	mom4ice,
		integer, intent(in)	lsm,
		logical, intent(in)	lprnt,
		integer, intent(in)	ipr,
		real (kind=kind_phys), dimension(im), intent(inout)	hice,
		real (kind=kind_phys), dimension(im), intent(inout)	fice,
		real (kind=kind_phys), dimension(im), intent(inout)	tice,
		real (kind=kind_phys), dimension(im), intent(inout)	weasd,
		real (kind=kind_phys), dimension(im), intent(inout)	tskin,
		real (kind=kind_phys), dimension(im), intent(inout)	tprcp,
		real (kind=kind_phys), dimension(im,km), intent(inout)	stc,
		real (kind=kind_phys), dimension(im), intent(inout)	ep,
		real (kind=kind_phys), dimension(im), intent(inout)	snwdph,
		real (kind=kind_phys), dimension(im), intent(inout)	qsurf,
		real (kind=kind_phys), dimension(im), intent(inout)	snowmt,
		real (kind=kind_phys), dimension(im), intent(inout)	gflux,
		real (kind=kind_phys), dimension(im), intent(inout)	cmm,
		real (kind=kind_phys), dimension(im), intent(inout)	chh,
		real (kind=kind_phys), dimension(im), intent(inout)	evap,
		real (kind=kind_phys), dimension(im), intent(inout)	hflx,
		character(len=*), intent(out)	errmsg,
		integer, intent(out)	errflg
	)

• Set flag for sea-ice.
• Update/read sea ice temperature from soil temperature and initialize variables.
• Convert snow depth in water equivalent from mm to m unit.
• Calculate sensible and latent heat flux over open water & sea ice.
• Calculate net solar incoming at top sneti.
• Calculate net non-solar and upir heat flux @ ice surface hfi.
• Calculate heat flux derivative at surface hfd.
• Assigin heat flux from ocean focn and snowfall rate as constants, which should be from ocean model and other physics.
• Initialize snow depth snowd.
• Call the three-layer thermodynamics sea ice model ice3lay().

References physcons::con_cp, physcons::con_eps, physcons::con_epsm1, physcons::con_fvirt, physcons::con_g, physcons::con_hvap, physcons::con_rd, physcons::con_sbc, physcons::con_t0c, physcons::con_tice, funcphys::fpvs(), and ice3lay().

Here is the call graph for this function: