!
! $Id: $
!
MODULE phyaqua_mod

  IMPLICIT NONE

CONTAINS

  SUBROUTINE iniaqua(nlon, latfi, lonfi, iflag_phys)

  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  !  Create an initial state (startphy.nc) for the physics
  !  Usefull for idealised cases (e.g. aquaplanets or testcases)
  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  USE phys_state_var_mod, ONLY: rlat, rlon, phys_state_var_init
  USE mod_phys_lmdz_para, ONLY: klon_omp
  USE comgeomphy, ONLY: rlond, rlatd
  IMPLICIT NONE
      
  INTEGER,INTENT(IN) :: nlon,iflag_phys
  REAL,INTENT(IN) :: lonfi(nlon),latfi(nlon)

  ! local variables
  REAL :: pi

  ! initializations:
  pi=2.*ASIN(1.)

  CALL phys_state_var_init()

  rlat(1:klon_omp)=rlatd(1:klon_omp)*180./pi
  rlon(1:klon_omp)=rlond(1:klon_omp)*180./pi


  ! Here you could create an initial condition for the physics
  ! ...
  ! ... fill in the fields...
  ! ...
  ! ... and create a "startphy.nc" file
      CALL phyredem ("startphy.nc")

  END SUBROUTINE iniaqua

END MODULE phyaqua_mod