!
! $Id$
!
      SUBROUTINE dynetat0_loc(fichnom,vcov,ucov,
     .                    teta,q,masse,ps,phis,time)
      USE infotrac
      use control_mod, only : planet_type
      USE parallel_lmdz
      IMPLICIT NONE

c=======================================================================
c
c   Auteur:  P. Le Van / L.Fairhead
c   -------
c
c   objet:
c   ------
c
c   Lecture de l'etat initial
c
c=======================================================================
c-----------------------------------------------------------------------
c   Declarations:
c   -------------

#include "dimensions.h"
#include "paramet.h"
#include "temps.h"
#include "comconst.h"
#include "comvert.h"
#include "comgeom.h"
#include "ener.h"
#include "netcdf.inc"
#include "description.h"
#include "serre.h"
#include "logic.h"
#include "iniprint.h"

c   Arguments:
c   ----------

      CHARACTER*(*) fichnom
      REAL vcov(ijb_v:ije_v,llm),ucov(ijb_u:ije_u,llm)
      REAL teta(ijb_u:ije_u,llm)
      REAL q(ijb_u:ije_u,llm,nqtot),masse(ijb_u:ije_u,llm)
      REAL ps(ijb_u:ije_u),phis(ijb_u:ije_u)

      REAL time

c   Variables 
c
      INTEGER length,iq
      PARAMETER (length = 100)
      REAL tab_cntrl(length) ! tableau des parametres du run
      INTEGER ierr, nid, nvarid
      REAL,ALLOCATABLE :: vcov_glo(:,:),ucov_glo(:,:),teta_glo(:,:)
      REAL,ALLOCATABLE :: q_glo(:,:),masse_glo(:,:),ps_glo(:)
      REAL,ALLOCATABLE :: phis_glo(:)

      INTEGER idecal

c-----------------------------------------------------------------------
c  Ouverture NetCDF du fichier etat initial

      ierr = NF_OPEN (fichnom, NF_NOWRITE,nid)
      IF (ierr.NE.NF_NOERR) THEN
        write(lunout,*)
     &  'dynetat0_loc: Pb d''ouverture du fichier start.nc'
        write(lunout,*)' ierr = ', ierr
        CALL ABORT_GCM("DYNETAT0", "", 1)
      ENDIF

c
      ierr = NF_INQ_VARID (nid, "controle", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <controle> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tab_cntrl)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echoue pour <controle>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF

      !!! AS: idecal is a hack to be able to read planeto starts...
      !!!     .... while keeping everything OK for LMDZ EARTH
      if (planet_type.eq."generic") then
          print*,'NOTE NOTE NOTE : Planeto-like start files'
          idecal = 4
          annee_ref  = 2000
      else
          print*,'NOTE NOTE NOTE : Earth-like start files'
          idecal = 5
          annee_ref  = tab_cntrl(5)
      endif


      im         = tab_cntrl(1)
      jm         = tab_cntrl(2)
      lllm       = tab_cntrl(3)
      day_ref    = tab_cntrl(4)
      rad        = tab_cntrl(idecal+1)
      omeg       = tab_cntrl(idecal+2)
      g          = tab_cntrl(idecal+3)
      cpp        = tab_cntrl(idecal+4)
      kappa      = tab_cntrl(idecal+5)
      daysec     = tab_cntrl(idecal+6)
      dtvr       = tab_cntrl(idecal+7)
      etot0      = tab_cntrl(idecal+8)
      ptot0      = tab_cntrl(idecal+9)
      ztot0      = tab_cntrl(idecal+10)
      stot0      = tab_cntrl(idecal+11)
      ang0       = tab_cntrl(idecal+12)
      pa         = tab_cntrl(idecal+13)
      preff      = tab_cntrl(idecal+14)
c
      clon       = tab_cntrl(idecal+15)
      clat       = tab_cntrl(idecal+16)
      grossismx  = tab_cntrl(idecal+17)
      grossismy  = tab_cntrl(idecal+18)
c
      IF ( tab_cntrl(idecal+19).EQ.1. )  THEN
        fxyhypb  = . TRUE .
c        dzoomx   = tab_cntrl(25)
c        dzoomy   = tab_cntrl(26)
c        taux     = tab_cntrl(28)
c        tauy     = tab_cntrl(29)
      ELSE
        fxyhypb = . FALSE .
        ysinus  = . FALSE .
        IF( tab_cntrl(idecal+22).EQ.1. ) ysinus = . TRUE. 
      ENDIF

      day_ini = tab_cntrl(30)
      itau_dyn = tab_cntrl(31)
c   .................................................................
c
c
      write(lunout,*)'dynetat0_loc: rad,omeg,g,cpp,kappa',
     &               rad,omeg,g,cpp,kappa

      IF(   im.ne.iim           )  THEN
          PRINT 1,im,iim
          STOP
      ELSE  IF( jm.ne.jjm       )  THEN
          PRINT 2,jm,jjm
          STOP
      ELSE  IF( lllm.ne.llm     )  THEN
          PRINT 3,lllm,llm
          STOP
      ENDIF

      ierr = NF_INQ_VARID (nid, "rlonu", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <rlonu> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlonu)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, rlonu)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <rlonu>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF

      ierr = NF_INQ_VARID (nid, "rlatu", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <rlatu> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlatu)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, rlatu)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <rlatu>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF

      ierr = NF_INQ_VARID (nid, "rlonv", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <rlonv> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlonv)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, rlonv)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <rlonv>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF

      ierr = NF_INQ_VARID (nid, "rlatv", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <rlatv> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlatv)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, rlatv)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour rlatv"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF

      ierr = NF_INQ_VARID (nid, "cu", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <cu> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, cu)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, cu)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <cu>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF

      ierr = NF_INQ_VARID (nid, "cv", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <cv> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, cv)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, cv)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <cv>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF

      ierr = NF_INQ_VARID (nid, "aire", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <aire> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, aire)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, aire)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <aire>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
      
      ALLOCATE(phis_glo(ip1jmp1))
      
      ierr = NF_INQ_VARID (nid, "phisinit", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <phisinit> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, phis_glo)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, phis_glo)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <phisinit>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
      phis(ijb_u:ije_u)=phis_glo(ijb_u:ije_u)
      DEALLOCATE(phis_glo)

      ierr = NF_INQ_VARID (nid, "temps", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0: Le champ <temps> est absent"
         write(lunout,*)"dynetat0: J essaie <Time>"
         ierr = NF_INQ_VARID (nid, "Time", nvarid)
         IF (ierr .NE. NF_NOERR) THEN
            write(lunout,*)"dynetat0: Le champ <Time> est absent"
            CALL abort_gcm("dynetat0", "", 1)
         ENDIF
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, time)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, time)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee <temps>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF

      ierr = NF_INQ_VARID (nid, "ucov", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <ucov> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
      
      ALLOCATE(ucov_glo(ip1jmp1,llm))
      
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ucov_glo)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, ucov_glo)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <ucov>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF

      ucov(ijb_u:ije_u,:)=ucov_glo(ijb_u:ije_u,:)
      DEALLOCATE(ucov_glo)
      ALLOCATE(vcov_glo(ip1jm,llm))
      
      ierr = NF_INQ_VARID (nid, "vcov", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <vcov> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, vcov_glo)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, vcov_glo)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <vcov>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
      vcov(ijb_v:ije_v,:)=vcov_glo(ijb_v:ije_v,:)
      DEALLOCATE(vcov_glo)
      ALLOCATE(teta_glo(ip1jmp1,llm))

      ierr = NF_INQ_VARID (nid, "teta", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <teta> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, teta_glo)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, teta_glo)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <teta>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF

      teta(ijb_u:ije_u,:)=teta_glo(ijb_u:ije_u,:)
      DEALLOCATE(teta_glo)
      ALLOCATE(q_glo(ip1jmp1,llm))


      DO iq=1,nqtot
        ierr =  NF_INQ_VARID (nid, tname(iq), nvarid)
        IF (ierr .NE. NF_NOERR) THEN
           write(lunout,*)"dynetat0_loc: Le traceur <"                  &
     &     //trim(tname(iq))//"> est absent"
           write(lunout,*)"Il est donc initialise a zero"
           q(:,:,iq)=0.
        ELSE
#ifdef NC_DOUBLE
          ierr = NF_GET_VAR_DOUBLE(nid, nvarid, q_glo)
#else
          ierr = NF_GET_VAR_REAL(nid, nvarid, q_glo)
#endif
          IF (ierr .NE. NF_NOERR) THEN
            write(lunout,*)
     &      "dynetat0_loc: Lecture echouee pour "//tname(iq)
            CALL abort_gcm("dynetat0", "", 1)
          ENDIF
        q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)

        ENDIF
      ENDDO

      DEALLOCATE(q_glo)
      ALLOCATE(masse_glo(ip1jmp1,llm))

      ierr = NF_INQ_VARID (nid, "masse", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <masse> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, masse_glo)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, masse_glo)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <masse>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
      masse(ijb_u:ije_u,:)=masse_glo(ijb_u:ije_u,:)
      DEALLOCATE(masse_glo)
      ALLOCATE(ps_glo(ip1jmp1))

      ierr = NF_INQ_VARID (nid, "ps", nvarid)
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Le champ <ps> est absent"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF
#ifdef NC_DOUBLE
      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ps_glo)
#else
      ierr = NF_GET_VAR_REAL(nid, nvarid, ps_glo)
#endif
      IF (ierr .NE. NF_NOERR) THEN
         write(lunout,*)"dynetat0_loc: Lecture echouee pour <ps>"
         CALL abort_gcm("dynetat0", "", 1)
      ENDIF

      ps(ijb_u:ije_u)=ps_glo(ijb_u:ije_u)
      DEALLOCATE(ps_glo)

      ierr = NF_CLOSE(nid)

       day_ini=day_ini+INT(time)
       time=time-INT(time)

  1   FORMAT(//10x,'la valeur de im =',i4,2x,'lue sur le fichier de dem
     *arrage est differente de la valeur parametree iim =',i4//)
   2  FORMAT(//10x,'la valeur de jm =',i4,2x,'lue sur le fichier de dem
     *arrage est differente de la valeur parametree jjm =',i4//)
   3  FORMAT(//10x,'la valeur de lmax =',i4,2x,'lue sur le fichier dema
     *rrage est differente de la valeur parametree llm =',i4//)
   4  FORMAT(//10x,'la valeur de dtrv =',i4,2x,'lue sur le fichier dema
     *rrage est differente de la valeur  dtinteg =',i4//)

      RETURN
      END