add module for icb_thermo and icb_coupling

src/fesom_module.F90

@@ -36,6 +36,7 @@ module fesom_main_storage_module
   use age_tracer_init_interface
   use iceberg_params
   use iceberg_step
+  use iceberg_ocean_coupling
   ! Define icepack module
 
 #if defined (__icepack)

src/icb_coupling.F90

@@ -1,3 +1,17 @@
+module iceberg_ocean_coupling
+ USE MOD_MESH
+ use MOD_PARTIT
+ use MOD_ICE
+ USE MOD_DYN
+ use iceberg_params
+
+ public :: reset_ib_fluxes
+ public :: prepare_icb2fesom
+ public :: icb2fesom
+
+ contains
+
+
 subroutine reset_ib_fluxes()
     use o_param
 
@@ -190,3 +204,4 @@ subroutine icb2fesom(mesh, partit, ice)
     end if
 !---wiso-code-end
 end subroutine icb2fesom
+end module

src/icb_elem.F90

@@ -3,6 +3,8 @@ module iceberg_element
  USE MOD_MESH
  USE MOD_DYN
  use iceberg_params
+ use iceberg_thermodynamics
+ use iceberg_ocean_coupling
 ! use iceberg_dynamics
 ! use iceberg_step
 

src/icb_step.F90

@@ -6,6 +6,8 @@ module iceberg_step
  use iceberg_params
  use iceberg_dynamics
  use iceberg_element
+ use iceberg_thermodynamics
+ use iceberg_ocean_coupling
 
 implicit none
 

src/icb_thermo.F90

@@ -1,3 +1,23 @@
+module iceberg_thermodynamics
+ USE MOD_MESH
+ use MOD_PARTIT
+ use MOD_ICE
+ USE MOD_DYN
+ use iceberg_params
+ !use iceberg_element
+ !use iceberg_step
+ use iceberg_ocean_coupling
+
+ implicit none
+
+ public ::  iceberg_meltrates
+ public ::  iceberg_newdimensions
+ public ::  iceberg_heat_water_fluxes_3eq
+ public ::  potit_ib
+ public ::  fcn_density
+
+ contains 
+
 !!=============================================================================
 ! calculates the empirical melt rates of the iceberg as in 
 ! Martin: 'Parameterizing the fresh-water flux from land ice to ocean
@@ -596,3 +616,4 @@ subroutine fcn_density(t,s,z,rho)
       + 4.8314e-4*s**2)
  rho = rhopot / (1.0 + 0.1*z/bulk)
 end subroutine fcn_density
+end module

src/ice_oce_coupling.F90

@@ -265,6 +265,7 @@ subroutine oce_fluxes(ice, dynamics, tracers, partit, mesh)
                             init_flux_atm_ocn
 #endif
     use iceberg_params
+    use iceberg_ocean_coupling
     use cavity_interfaces
     !---fwf-code
     use g_clock