MODULE iniphyparam_mod
#include "use_logging.h"
  IMPLICIT NONE
  PRIVATE

  REAL, PARAMETER :: perfect_gas_const = 8314.46261815324 ! NB using g instead of kg for mass

  PUBLIC :: iniphyparam

CONTAINS

  SUBROUTINE iniphyparam(ngrid,nlayer, &
       &           punjours,              &
       &           pdayref,ptimestep,     &
       &           prad,pg,pr,pcpp)
    USE callkeys
    USE phys_const, ONLY : planet_rad,g,r,cpp,rcp,dtphys,unjours,mugaz
    USE planet, ONLY : coefir, coefvis
    USE astronomy
    USE turbulence, ONLY : lmixmin, emin_turb
    USE surface
    USE read_param_mod

    INTEGER, INTENT(IN) :: &
         ngrid, &       ! Size of the horizontal grid
         nlayer, &      ! Number of vertical layers.
         pdayref        ! Day of reference for the simulation
    REAL, INTENT(IN)  :: ptimestep, prad, pg, pr, pcpp, punjours

    CALL read_param('unjours', 86400., unjours,'unjours')
    CALL read_param('planet_rad',prad,planet_rad,'planet_rad')
    CALL read_param('g',9.8           ,g,'g')
    CALL read_param('cpp',1004.       ,cpp,'cpp')
    CALL read_param('mugaz',28.       ,mugaz,'mugaz')
    r=perfect_gas_const/mugaz
    rcp=r/cpp

    CALL read_param('year_day',360.   ,year_day,'year_day')
    CALL read_param('periheli',150.   ,periheli,'periheli')
    CALL read_param('aphelie',150.    ,aphelie,'aphelie')
    CALL read_param('peri_day',0.     ,peri_day,'peri_day')
    CALL read_param('obliquit',23.    ,obliquit,'obliquit')
    CALL read_param('Cd_mer',.01      ,Cd_mer,'Cd_mer')
    CALL read_param('Cd_ter',.01      ,Cd_ter,'Cd_ter')
    CALL read_param('I_mer',30000.    ,I_mer,'I_mer')
    CALL read_param('I_ter',30000.    ,I_ter,'I_ter')
    CALL read_param('alb_ter',.112    ,alb_ter,'alb_ter')
    CALL read_param('alb_mer',.112    ,alb_mer,'alb_mer')
    CALL read_param('emi_mer',1.      ,emi_mer,'emi_mer')
    CALL read_param('emi_mer',1.      ,emi_mer,'emi_mer')
    CALL read_param('emin_turb',1.e-16 ,emin_turb,'emin_turb')
    CALL read_param('lmixmin',100.    ,lmixmin,'lmixmin')
    CALL read_param('coefvis',.99     ,coefvis,'coefvis')
    CALL read_param('coefir',.08      ,coefir,'coefir')

    CALL read_param('callrad',.true.,callrad,'appel rayonnemen')
    CALL read_param('calldifv',.true.,calldifv,'appel difv')
    CALL read_param('calladj',.true.,calladj,'appel adj')
    CALL read_param('callcond',.true.,callcond,'appel cond')
    CALL read_param('callsoil',.true.,callsoil,'appel soil')
    CALL read_param('season',.true.,season,'appel soil')
    CALL read_param('diurnal',.false.,diurnal,'appel soil')
    CALL read_param('lverbose',.true.,lverbose,'appel soil')
    CALL read_param('period_sort',1.,period_sort,'period sorties en jour')

    CALL check_mismatch('unjours', punjours, unjours)
    CALL check_mismatch('rad', prad, planet_rad)
    CALL check_mismatch('g', pg, g)
    CALL check_mismatch('R', pr, r)
    CALL check_mismatch('cpp', pcpp, cpp)

    WRITELOG(*,*) 'Activation de la physique:'
    WRITELOG(*,*) ' R=',r
    WRITELOG(*,*) ' Rayonnement ',callrad
    WRITELOG(*,*) ' Diffusion verticale turbulente ', calldifv
    WRITELOG(*,*) ' Ajustement convectif ',calladj
    WRITELOG(*,*) ' Sol ',callsoil
    WRITELOG(*,*) ' Cycle diurne ',diurnal

    !   choice of the frequency of the computation of radiations
    IF(diurnal) THEN
       iradia=NINT(punjours/(20.*ptimestep))
    ELSE
       iradia=NINT(punjours/(4.*ptimestep))
    ENDIF
    iradia=1
    WRITELOG(*,*) 'unjours',punjours
    WRITELOG(*,*) 'The radiative transfer is computed each ', &
         &   iradia,' physical time-step or each ', &
         &   iradia*ptimestep,' seconds'

    dtphys=ptimestep

    LOG_INFO('iniphyparama')
  END SUBROUTINE iniphyparam

  SUBROUTINE check_mismatch(name, a,b)
    CHARACTER(*), INTENT(IN) :: name
    REAL, INTENT(IN) :: a,b
    IF(a /= b) THEN
       PRINT *, 'Phys/dyn mismatch for ', name, ' : ',a,b
    END IF
  END SUBROUTINE check_mismatch

END MODULE iniphyparam_mod