!
! $Id: create_etat0_limit.F 1279 2009-12-10 09:02:56Z lguez $
!
       PROGRAM create_etat0_limit
#ifdef CPP_EARTH
! This prog. is designed to work for Earth
       USE dimphy
       USE comgeomphy
       USE infotrac
#ifdef CPP_IOIPSL
       use ioipsl, only: ioconf_calendar
#endif
       IMPLICIT NONE
c
c
c     Programme d'appel a etat0, creation des etats initiaux et limit_netcdf
c   
c
c     interbar = .T . si appel a  interpol. barycentrique inter_barxy
c
c     extrap   = .T . si on fait une extrapolation de donnees , comme pour
c       les  SST  lorsque  le fichier ne contient pas uniquement  des points 
c     oceaniques .
c
c     oldice   = .T. si l'on veut garder les anciennes glaces , obtenues
c     par  grille_m  ( grid_atob ) .
c
c     on cree le masque dans etat0 que l'on passe ensuite dans limit pour 
c     garder les coherences

      LOGICAL interbar, extrap , oldice
      PARAMETER ( interbar = .true. , extrap = .FALSE. , oldice=.false.)
#include "dimensions.h"
#include "paramet.h"
#include "indicesol.h"
#include  "control.h"
      REAL :: masque(iip1,jjp1)
!      REAL :: pctsrf(iim*(jjm-1)+2, nbsrf)

      IF (config_inca /= 'none') THEN
#ifdef INCA
         call init_const_lmdz(
     $        nbtr,anneeref,dayref,
     $        iphysiq, day_step,nday)
#endif
         print *, 'nbtr =' , nbtr 
      END IF

      CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
      PRINT *,'---> klon=',klon
      call InitComgeomphy

#ifdef CPP_IOIPSL
      call ioconf_calendar('360d')
#endif

      WRITE(6,*) '  *********************  '
      WRITE(6,*) ' interbar = ',interbar
      CALL etat0_netcdf ( interbar, masque )
c
      WRITE(6,1)
      WRITE(6,*) '  *********************  '
      WRITE(6,*) '  ***  Limit_netcdf ***  '
      WRITE(6,*) '  *********************  '
      WRITE(6,1)
      
c     
      CALL  limit_netcdf ( interbar, extrap , oldice, masque)

1     FORMAT(//)

#endif
! of #ifdef CPP_EARTH
      STOP
      END