fabm_calendar_year: provide calendar_year and maximum_time_step to FABM (#64)

Author: knutaros

src/fabm/gotm_fabm.F90

@@ -115,7 +115,7 @@ module gotm_fabm
    REALTYPE,allocatable,dimension(:),target :: par,pres,swr,k_par,z,nuh_ct
 
    ! External variables
-   REALTYPE :: dt,dt_eff   ! External and internal time steps
+   REALTYPE, target :: dt,dt_eff   ! External and internal time steps
    integer  :: w_adv_ctr   ! Scheme for vertical advection (0 if not used)
    REALTYPE,pointer,dimension(:) :: nuh,h,bioshade,w,rho
    REALTYPE,pointer,dimension(:) :: SRelaxTau,sProf,salt
@@ -124,6 +124,7 @@ module gotm_fabm
    REALTYPE,pointer :: I_0,A,g1,g2
    integer,pointer  :: yearday,secondsofday
    REALTYPE, target :: decimal_yearday
+   REALTYPE, target :: decimal_year
    logical          :: fabm_ready
 
    REALTYPE,pointer :: fabm_airp
@@ -836,6 +837,13 @@ subroutine init_gotm_fabm_state(nlev)
 
 !-----------------------------------------------------------------------
 !BOC
+
+   if ( model%variable_needs_values( standard_variables%calendar_year ) ) then
+      if ( associated(fabm_calendar_date) .and. associated(fabm_julianday) ) then
+         call model%link_scalar( standard_variables%calendar_year, decimal_year )
+      end if
+   end if
+
    call fabm_check_ready(model)
    fabm_ready = .true.
 
@@ -1032,6 +1040,7 @@ subroutine set_env_gotm_fabm(latitude,longitude,dt_,w_adv_method_,w_adv_ctr_,tem
 
    ! Calculate and save internal time step.
    dt_eff = dt/split_factor
+   call model%link_scalar( standard_variables%maximum_time_step, dt_eff )
 
    I_0 => I_0_
    A => A_
@@ -1229,6 +1238,7 @@ subroutine do_gotm_fabm(nlev,itime)
    do split=1,split_factor
       if (has_date) then
          call fabm_update_time(model, itime, yyyy, mm, dd, real(secondsofday, gotmrk))
+         decimal_year = real( yyyy, gotmrk )
       else
          call fabm_update_time(model, itime)
       end if