MODULE mod_1D_cases_read

  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  !Declarations specifiques au cas standard
  CHARACTER*80 :: fich_cas
  ! Discr?tisation
  INTEGER nlev_cas, nt_cas


  !       integer year_ini_cas, day_ini_cas, mth_ini_cas
  !       real heure_ini_cas
  !       real day_ju_ini_cas   ! Julian day of case first day
  !       parameter (year_ini_cas=2011)
  !       parameter (year_ini_cas=1969)
  !       parameter (mth_ini_cas=10)
  !       parameter (mth_ini_cas=6)
  !       parameter (day_ini_cas=1)  ! 10 = 10Juil2006
  !       parameter (day_ini_cas=24)  ! 24 = 24 juin 1969
  !       parameter (heure_ini_cas=0.) !0h en secondes
  !       real pdt_cas
  !       parameter (pdt_cas=3.*3600)

  !CR ATTENTION TEST AMMA
  !        parameter (year_ini_cas=2006)
  !        parameter (mth_ini_cas=7)
  !        parameter (day_ini_cas=10)  ! 10 = 10Juil2006
  !        parameter (heure_ini_cas=0.) !0h en secondes
  !        parameter (pdt_cas=1800.)

  !profils environnementaux
  REAL, ALLOCATABLE :: plev_cas(:, :)

  REAL, ALLOCATABLE :: z_cas(:, :)
  REAL, ALLOCATABLE :: t_cas(:, :), q_cas(:, :), rh_cas(:, :)
  REAL, ALLOCATABLE :: th_cas(:, :), rv_cas(:, :)
  REAL, ALLOCATABLE :: u_cas(:, :)
  REAL, ALLOCATABLE :: v_cas(:, :)

  !forcing
  REAL, ALLOCATABLE :: ht_cas(:, :), vt_cas(:, :), dt_cas(:, :), dtrad_cas(:, :)
  REAL, ALLOCATABLE :: hth_cas(:, :), vth_cas(:, :), dth_cas(:, :)
  REAL, ALLOCATABLE :: hq_cas(:, :), vq_cas(:, :), dq_cas(:, :)
  REAL, ALLOCATABLE :: hr_cas(:, :), vr_cas(:, :), dr_cas(:, :)
  REAL, ALLOCATABLE :: hu_cas(:, :), vu_cas(:, :), du_cas(:, :)
  REAL, ALLOCATABLE :: hv_cas(:, :), vv_cas(:, :), dv_cas(:, :)
  REAL, ALLOCATABLE :: vitw_cas(:, :)
  REAL, ALLOCATABLE :: ug_cas(:, :), vg_cas(:, :)
  REAL, ALLOCATABLE :: lat_cas(:), sens_cas(:), ts_cas(:), ustar_cas(:)
  REAL, ALLOCATABLE :: uw_cas(:, :), vw_cas(:, :), q1_cas(:, :), q2_cas(:, :)

  !champs interpoles
  REAL, ALLOCATABLE :: plev_prof_cas(:)
  REAL, ALLOCATABLE :: t_prof_cas(:)
  REAL, ALLOCATABLE :: q_prof_cas(:)
  REAL, ALLOCATABLE :: u_prof_cas(:)
  REAL, ALLOCATABLE :: v_prof_cas(:)

  REAL, ALLOCATABLE :: vitw_prof_cas(:)
  REAL, ALLOCATABLE :: ug_prof_cas(:)
  REAL, ALLOCATABLE :: vg_prof_cas(:)
  REAL, ALLOCATABLE :: ht_prof_cas(:)
  REAL, ALLOCATABLE :: hq_prof_cas(:)
  REAL, ALLOCATABLE :: vt_prof_cas(:)
  REAL, ALLOCATABLE :: vq_prof_cas(:)
  REAL, ALLOCATABLE :: dt_prof_cas(:)
  REAL, ALLOCATABLE :: dtrad_prof_cas(:)
  REAL, ALLOCATABLE :: dq_prof_cas(:)
  REAL, ALLOCATABLE :: hu_prof_cas(:)
  REAL, ALLOCATABLE :: hv_prof_cas(:)
  REAL, ALLOCATABLE :: vu_prof_cas(:)
  REAL, ALLOCATABLE :: vv_prof_cas(:)
  REAL, ALLOCATABLE :: du_prof_cas(:)
  REAL, ALLOCATABLE :: dv_prof_cas(:)
  REAL, ALLOCATABLE :: uw_prof_cas(:)
  REAL, ALLOCATABLE :: vw_prof_cas(:)
  REAL, ALLOCATABLE :: q1_prof_cas(:)
  REAL, ALLOCATABLE :: q2_prof_cas(:)

  REAL lat_prof_cas, sens_prof_cas, ts_prof_cas, ustar_prof_cas


CONTAINS

  SUBROUTINE read_1D_cas

    INTEGER nid, rid, ierr
    INTEGER ii, jj

    fich_cas = 'setup/cas.nc'
    PRINT*, 'fich_cas ', fich_cas
    ierr = nf90_open(fich_cas, nf90_nowrite, nid)
    PRINT*, 'fich_cas,nf90_nowrite,nid ', fich_cas, nf90_nowrite, nid
    IF (ierr/=nf90_noerr) THEN
      WRITE(*, *) 'ERROR: GROS Pb opening forcings nc file '
      WRITE(*, *) nf90_strerror(ierr)
      stop ""
    endif
    !.......................................................................
    ierr = nf90_inq_dimid(nid, 'lat', rid)
    IF (ierr/=nf90_noerr) THEN
      PRINT*, 'Oh probleme lecture dimension lat'
    ENDIF
    ierr = nf90_inquire_dimension(nid, rid, len = ii)
    PRINT*, 'OK1 nid,rid,lat', nid, rid, ii
    !.......................................................................
    ierr = nf90_inq_dimid(nid, 'lon', rid)
    IF (ierr/=nf90_noerr) THEN
      PRINT*, 'Oh probleme lecture dimension lon'
    ENDIF
    ierr = nf90_inquire_dimension(nid, rid, len = jj)
    PRINT*, 'OK2 nid,rid,lat', nid, rid, jj
    !.......................................................................
    ierr = nf90_inq_dimid(nid, 'lev', rid)
    IF (ierr/=nf90_noerr) THEN
      PRINT*, 'Oh probleme lecture dimension zz'
    ENDIF
    ierr = nf90_inquire_dimension(nid, rid, len = nlev_cas)
    PRINT*, 'OK3 nid,rid,nlev_cas', nid, rid, nlev_cas
    !.......................................................................
    ierr = nf90_inq_dimid(nid, 'time', rid)
    PRINT*, 'nid,rid', nid, rid
    nt_cas = 0
    IF (ierr/=nf90_noerr) THEN
      stop 'probleme lecture dimension sens'
    ENDIF
    ierr = nf90_inquire_dimension(nid, rid, len = nt_cas)
    PRINT*, 'OK4 nid,rid,nt_cas', nid, rid, nt_cas

    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !profils moyens:
    allocate(plev_cas(nlev_cas, nt_cas))
    allocate(z_cas(nlev_cas, nt_cas))
    allocate(t_cas(nlev_cas, nt_cas), q_cas(nlev_cas, nt_cas), rh_cas(nlev_cas, nt_cas))
    allocate(th_cas(nlev_cas, nt_cas), rv_cas(nlev_cas, nt_cas))
    allocate(u_cas(nlev_cas, nt_cas))
    allocate(v_cas(nlev_cas, nt_cas))

    !forcing
    allocate(ht_cas(nlev_cas, nt_cas), vt_cas(nlev_cas, nt_cas), dt_cas(nlev_cas, nt_cas), dtrad_cas(nlev_cas, nt_cas))
    allocate(hq_cas(nlev_cas, nt_cas), vq_cas(nlev_cas, nt_cas), dq_cas(nlev_cas, nt_cas))
    allocate(hth_cas(nlev_cas, nt_cas), vth_cas(nlev_cas, nt_cas), dth_cas(nlev_cas, nt_cas))
    allocate(hr_cas(nlev_cas, nt_cas), vr_cas(nlev_cas, nt_cas), dr_cas(nlev_cas, nt_cas))
    allocate(hu_cas(nlev_cas, nt_cas), vu_cas(nlev_cas, nt_cas), du_cas(nlev_cas, nt_cas))
    allocate(hv_cas(nlev_cas, nt_cas), vv_cas(nlev_cas, nt_cas), dv_cas(nlev_cas, nt_cas))
    allocate(vitw_cas(nlev_cas, nt_cas))
    allocate(ug_cas(nlev_cas, nt_cas))
    allocate(vg_cas(nlev_cas, nt_cas))
    allocate(lat_cas(nt_cas), sens_cas(nt_cas), ts_cas(nt_cas), ustar_cas(nt_cas))
    allocate(uw_cas(nlev_cas, nt_cas), vw_cas(nlev_cas, nt_cas), q1_cas(nlev_cas, nt_cas), q2_cas(nlev_cas, nt_cas))


    !champs interpoles
    allocate(plev_prof_cas(nlev_cas))
    allocate(t_prof_cas(nlev_cas))
    allocate(q_prof_cas(nlev_cas))
    allocate(u_prof_cas(nlev_cas))
    allocate(v_prof_cas(nlev_cas))

    allocate(vitw_prof_cas(nlev_cas))
    allocate(ug_prof_cas(nlev_cas))
    allocate(vg_prof_cas(nlev_cas))
    allocate(ht_prof_cas(nlev_cas))
    allocate(hq_prof_cas(nlev_cas))
    allocate(hu_prof_cas(nlev_cas))
    allocate(hv_prof_cas(nlev_cas))
    allocate(vt_prof_cas(nlev_cas))
    allocate(vq_prof_cas(nlev_cas))
    allocate(vu_prof_cas(nlev_cas))
    allocate(vv_prof_cas(nlev_cas))
    allocate(dt_prof_cas(nlev_cas))
    allocate(dtrad_prof_cas(nlev_cas))
    allocate(dq_prof_cas(nlev_cas))
    allocate(du_prof_cas(nlev_cas))
    allocate(dv_prof_cas(nlev_cas))
    allocate(uw_prof_cas(nlev_cas))
    allocate(vw_prof_cas(nlev_cas))
    allocate(q1_prof_cas(nlev_cas))
    allocate(q2_prof_cas(nlev_cas))

    PRINT*, 'Allocations OK'
    CALL read_cas(nid, nlev_cas, nt_cas                                       &
            , z_cas, plev_cas, t_cas, q_cas, rh_cas, th_cas, rv_cas, u_cas, v_cas         &
            , ug_cas, vg_cas, vitw_cas, du_cas, hu_cas, vu_cas, dv_cas, hv_cas, vv_cas    &
            , dt_cas, dtrad_cas, ht_cas, vt_cas, dq_cas, hq_cas, vq_cas                 &
            , dth_cas, hth_cas, vth_cas, dr_cas, hr_cas, vr_cas, sens_cas, lat_cas, ts_cas&
            , ustar_cas, uw_cas, vw_cas, q1_cas, q2_cas)
    PRINT*, 'Read cas OK'

  END SUBROUTINE read_1D_cas


  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  SUBROUTINE deallocate_1D_cases
    !profils environnementaux:
    deallocate(plev_cas)

    deallocate(z_cas)
    deallocate(t_cas, q_cas, rh_cas)
    deallocate(th_cas, rv_cas)
    deallocate(u_cas)
    deallocate(v_cas)

    !forcing
    deallocate(ht_cas, vt_cas, dt_cas, dtrad_cas)
    deallocate(hq_cas, vq_cas, dq_cas)
    deallocate(hth_cas, vth_cas, dth_cas)
    deallocate(hr_cas, vr_cas, dr_cas)
    deallocate(hu_cas, vu_cas, du_cas)
    deallocate(hv_cas, vv_cas, dv_cas)
    deallocate(vitw_cas)
    deallocate(ug_cas)
    deallocate(vg_cas)
    deallocate(lat_cas, sens_cas, ts_cas, ustar_cas, uw_cas, vw_cas, q1_cas, q2_cas)

    !champs interpoles
    deallocate(plev_prof_cas)
    deallocate(t_prof_cas)
    deallocate(q_prof_cas)
    deallocate(u_prof_cas)
    deallocate(v_prof_cas)

    deallocate(vitw_prof_cas)
    deallocate(ug_prof_cas)
    deallocate(vg_prof_cas)
    deallocate(ht_prof_cas)
    deallocate(hq_prof_cas)
    deallocate(hu_prof_cas)
    deallocate(hv_prof_cas)
    deallocate(vt_prof_cas)
    deallocate(vq_prof_cas)
    deallocate(vu_prof_cas)
    deallocate(vv_prof_cas)
    deallocate(dt_prof_cas)
    deallocate(dtrad_prof_cas)
    deallocate(dq_prof_cas)
    deallocate(du_prof_cas)
    deallocate(dv_prof_cas)
    deallocate(t_prof_cas)
    deallocate(q_prof_cas)
    deallocate(u_prof_cas)
    deallocate(v_prof_cas)
    deallocate(uw_prof_cas)
    deallocate(vw_prof_cas)
    deallocate(q1_prof_cas)
    deallocate(q2_prof_cas)

  END SUBROUTINE deallocate_1D_cases

  !=====================================================================
  SUBROUTINE read_cas(nid, nlevel, ntime                          &
          , zz, pp, temp, qv, rh, theta, rv, u, v, ug, vg, w, &
          du, hu, vu, dv, hv, vv, dt, dtrad, ht, vt, dq, hq, vq, &
          dth, hth, vth, dr, hr, vr, sens, flat, ts, ustar, uw, vw, q1, q2)

    !program reading forcing of the case study
    USE netcdf, ONLY: nf90_get_var

    INTEGER ntime, nlevel

    REAL zz(nlevel, ntime)
    REAL pp(nlevel, ntime)
    REAL temp(nlevel, ntime), qv(nlevel, ntime), rh(nlevel, ntime)
    REAL theta(nlevel, ntime), rv(nlevel, ntime)
    REAL u(nlevel, ntime)
    REAL v(nlevel, ntime)
    REAL ug(nlevel, ntime)
    REAL vg(nlevel, ntime)
    REAL w(nlevel, ntime)
    REAL du(nlevel, ntime), hu(nlevel, ntime), vu(nlevel, ntime)
    REAL dv(nlevel, ntime), hv(nlevel, ntime), vv(nlevel, ntime)
    REAL dt(nlevel, ntime), ht(nlevel, ntime), vt(nlevel, ntime)
    REAL dtrad(nlevel, ntime)
    REAL dq(nlevel, ntime), hq(nlevel, ntime), vq(nlevel, ntime)
    REAL dth(nlevel, ntime), hth(nlevel, ntime), vth(nlevel, ntime)
    REAL dr(nlevel, ntime), hr(nlevel, ntime), vr(nlevel, ntime)
    REAL flat(ntime), sens(ntime), ts(ntime), ustar(ntime)
    REAL uw(nlevel, ntime), vw(nlevel, ntime), q1(nlevel, ntime), q2(nlevel, ntime)

    INTEGER nid, ierr, rid
    INTEGER nbvar3d
    parameter(nbvar3d = 39)
    INTEGER var3didin(nbvar3d)

    ierr = nf90_inq_varid(nid, "zz", var3didin(1))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'lev'
    endif

    ierr = nf90_inq_varid(nid, "pp", var3didin(2))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'plev'
    endif

    ierr = nf90_inq_varid(nid, "temp", var3didin(3))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'temp'
    endif

    ierr = nf90_inq_varid(nid, "qv", var3didin(4))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'qv'
    endif

    ierr = nf90_inq_varid(nid, "rh", var3didin(5))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'rh'
    endif

    ierr = nf90_inq_varid(nid, "theta", var3didin(6))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'theta'
    endif

    ierr = nf90_inq_varid(nid, "rv", var3didin(7))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'rv'
    endif

    ierr = nf90_inq_varid(nid, "u", var3didin(8))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'u'
    endif

    ierr = nf90_inq_varid(nid, "v", var3didin(9))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'v'
    endif

    ierr = nf90_inq_varid(nid, "ug", var3didin(10))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'ug'
    endif

    ierr = nf90_inq_varid(nid, "vg", var3didin(11))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'vg'
    endif

    ierr = nf90_inq_varid(nid, "w", var3didin(12))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'w'
    endif

    ierr = nf90_inq_varid(nid, "advu", var3didin(13))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'advu'
    endif

    ierr = nf90_inq_varid(nid, "hu", var3didin(14))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'hu'
    endif

    ierr = nf90_inq_varid(nid, "vu", var3didin(15))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'vu'
    endif

    ierr = nf90_inq_varid(nid, "advv", var3didin(16))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'advv'
    endif

    ierr = nf90_inq_varid(nid, "hv", var3didin(17))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'hv'
    endif

    ierr = nf90_inq_varid(nid, "vv", var3didin(18))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'vv'
    endif

    ierr = nf90_inq_varid(nid, "advT", var3didin(19))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'advT'
    endif

    ierr = nf90_inq_varid(nid, "hT", var3didin(20))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'hT'
    endif

    ierr = nf90_inq_varid(nid, "vT", var3didin(21))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'vT'
    endif

    ierr = nf90_inq_varid(nid, "advq", var3didin(22))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'advq'
    endif

    ierr = nf90_inq_varid(nid, "hq", var3didin(23))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'hq'
    endif

    ierr = nf90_inq_varid(nid, "vq", var3didin(24))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'vq'
    endif

    ierr = nf90_inq_varid(nid, "advth", var3didin(25))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'advth'
    endif

    ierr = nf90_inq_varid(nid, "hth", var3didin(26))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'hth'
    endif

    ierr = nf90_inq_varid(nid, "vth", var3didin(27))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'vth'
    endif

    ierr = nf90_inq_varid(nid, "advr", var3didin(28))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'advr'
    endif

    ierr = nf90_inq_varid(nid, "hr", var3didin(29))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'hr'
    endif

    ierr = nf90_inq_varid(nid, "vr", var3didin(30))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'vr'
    endif

    ierr = nf90_inq_varid(nid, "radT", var3didin(31))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'radT'
    endif

    ierr = nf90_inq_varid(nid, "sens", var3didin(32))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'sens'
    endif

    ierr = nf90_inq_varid(nid, "flat", var3didin(33))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'flat'
    endif

    ierr = nf90_inq_varid(nid, "ts", var3didin(34))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'ts'
    endif

    ierr = nf90_inq_varid(nid, "ustar", var3didin(35))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'ustar'
    endif

    ierr = nf90_inq_varid(nid, "uw", var3didin(36))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'uw'
    endif

    ierr = nf90_inq_varid(nid, "vw", var3didin(37))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'vw'
    endif

    ierr = nf90_inq_varid(nid, "q1", var3didin(38))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'q1'
    endif

    ierr = nf90_inq_varid(nid, "q2", var3didin(39))
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop 'q2'
    endif

    ierr = nf90_get_var(nid, var3didin(1), zz)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture z ok',zz

    ierr = nf90_get_var(nid, var3didin(2), pp)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture pp ok',pp

    ierr = nf90_get_var(nid, var3didin(3), temp)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture T ok',temp

    ierr = nf90_get_var(nid, var3didin(4), qv)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture qv ok',qv

    ierr = nf90_get_var(nid, var3didin(5), rh)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture rh ok',rh

    ierr = nf90_get_var(nid, var3didin(6), theta)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture theta ok',theta

    ierr = nf90_get_var(nid, var3didin(7), rv)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture rv ok',rv

    ierr = nf90_get_var(nid, var3didin(8), u)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture u ok',u

    ierr = nf90_get_var(nid, var3didin(9), v)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture v ok',v

    ierr = nf90_get_var(nid, var3didin(10), ug)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture ug ok',ug

    ierr = nf90_get_var(nid, var3didin(11), vg)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture vg ok',vg

    ierr = nf90_get_var(nid, var3didin(12), w)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture w ok',w

    ierr = nf90_get_var(nid, var3didin(13), du)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture du ok',du

    ierr = nf90_get_var(nid, var3didin(14), hu)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture hu ok',hu

    ierr = nf90_get_var(nid, var3didin(15), vu)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture vu ok',vu

    ierr = nf90_get_var(nid, var3didin(16), dv)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture dv ok',dv

    ierr = nf90_get_var(nid, var3didin(17), hv)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture hv ok',hv

    ierr = nf90_get_var(nid, var3didin(18), vv)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture vv ok',vv

    ierr = nf90_get_var(nid, var3didin(19), dt)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture dt ok',dt

    ierr = nf90_get_var(nid, var3didin(20), ht)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture ht ok',ht

    ierr = nf90_get_var(nid, var3didin(21), vt)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture vt ok',vt

    ierr = nf90_get_var(nid, var3didin(22), dq)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture dq ok',dq

    ierr = nf90_get_var(nid, var3didin(23), hq)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture hq ok',hq

    ierr = nf90_get_var(nid, var3didin(24), vq)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture vq ok',vq

    ierr = nf90_get_var(nid, var3didin(25), dth)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture dth ok',dth

    ierr = nf90_get_var(nid, var3didin(26), hth)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture hth ok',hth

    ierr = nf90_get_var(nid, var3didin(27), vth)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture vth ok',vth

    ierr = nf90_get_var(nid, var3didin(28), dr)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture dr ok',dr

    ierr = nf90_get_var(nid, var3didin(29), hr)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture hr ok',hr

    ierr = nf90_get_var(nid, var3didin(30), vr)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture vr ok',vr

    ierr = nf90_get_var(nid, var3didin(31), dtrad)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture dtrad ok',dtrad

    ierr = nf90_get_var(nid, var3didin(32), sens)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture sens ok',sens

    ierr = nf90_get_var(nid, var3didin(33), flat)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture flat ok',flat

    ierr = nf90_get_var(nid, var3didin(34), ts)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !          WRITE(*,*)'lecture ts ok',ts

    ierr = nf90_get_var(nid, var3didin(35), ustar)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !         WRITE(*,*)'lecture ustar ok',ustar

    ierr = nf90_get_var(nid, var3didin(36), uw)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !         WRITE(*,*)'lecture uw ok',uw

    ierr = nf90_get_var(nid, var3didin(37), vw)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !         WRITE(*,*)'lecture vw ok',vw

    ierr = nf90_get_var(nid, var3didin(38), q1)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !         WRITE(*,*)'lecture q1 ok',q1

    ierr = nf90_get_var(nid, var3didin(39), q2)
    IF(ierr/=nf90_noerr) THEN
      WRITE(*, *) nf90_strerror(ierr)
      stop "getvarup"
    endif
    !         WRITE(*,*)'lecture q2 ok',q2

  END SUBROUTINE  read_cas
  !======================================================================
  SUBROUTINE interp_case_time(day, day1, annee_ref                &
          !    &         ,year_cas,day_cas,nt_cas,pdt_forc,nlev_cas      &
          , nt_cas, nlev_cas                                       &
          , ts_cas, plev_cas, t_cas, q_cas, u_cas, v_cas               &
          , ug_cas, vg_cas, vitw_cas, du_cas, hu_cas, vu_cas           &
          , dv_cas, hv_cas, vv_cas, dt_cas, ht_cas, vt_cas, dtrad_cas   &
          , dq_cas, hq_cas, vq_cas, lat_cas, sens_cas, ustar_cas       &
          , uw_cas, vw_cas, q1_cas, q2_cas                           &
          , ts_prof_cas, plev_prof_cas, t_prof_cas, q_prof_cas       &
          , u_prof_cas, v_prof_cas, ug_prof_cas, vg_prof_cas         &
          , vitw_prof_cas, du_prof_cas, hu_prof_cas, vu_prof_cas     &
          , dv_prof_cas, hv_prof_cas, vv_prof_cas, dt_prof_cas       &
          , ht_prof_cas, vt_prof_cas, dtrad_prof_cas, dq_prof_cas    &
          , hq_prof_cas, vq_prof_cas, lat_prof_cas, sens_prof_cas    &
          , ustar_prof_cas, uw_prof_cas, vw_prof_cas, q1_prof_cas, q2_prof_cas)

    USE lmdz_compar1d
    USE lmdz_date_cas, ONLY: year_ini_cas, mth_ini_cas, day_deb, heure_ini_cas, pdt_cas, day_ju_ini_cas

    IMPLICIT NONE

    !---------------------------------------------------------------------------------------
    ! Time interpolation of a 2D field to the timestep corresponding to day

    ! day: current julian day (e.g. 717538.2)
    ! day1: first day of the simulation
    ! nt_cas: total nb of data in the forcing
    ! pdt_cas: total time interval (in sec) between 2 forcing data
    !---------------------------------------------------------------------------------------

    ! inputs:
    INTEGER annee_ref
    INTEGER nt_cas, nlev_cas
    REAL day, day1, day_cas
    REAL ts_cas(nt_cas)
    REAL plev_cas(nlev_cas, nt_cas)
    REAL t_cas(nlev_cas, nt_cas), q_cas(nlev_cas, nt_cas)
    REAL u_cas(nlev_cas, nt_cas), v_cas(nlev_cas, nt_cas)
    REAL ug_cas(nlev_cas, nt_cas), vg_cas(nlev_cas, nt_cas)
    REAL vitw_cas(nlev_cas, nt_cas)
    REAL du_cas(nlev_cas, nt_cas), hu_cas(nlev_cas, nt_cas), vu_cas(nlev_cas, nt_cas)
    REAL dv_cas(nlev_cas, nt_cas), hv_cas(nlev_cas, nt_cas), vv_cas(nlev_cas, nt_cas)
    REAL dt_cas(nlev_cas, nt_cas), ht_cas(nlev_cas, nt_cas), vt_cas(nlev_cas, nt_cas)
    REAL dtrad_cas(nlev_cas, nt_cas)
    REAL dq_cas(nlev_cas, nt_cas), hq_cas(nlev_cas, nt_cas), vq_cas(nlev_cas, nt_cas)
    REAL lat_cas(nt_cas)
    REAL sens_cas(nt_cas)
    REAL ustar_cas(nt_cas), uw_cas(nlev_cas, nt_cas), vw_cas(nlev_cas, nt_cas)
    REAL q1_cas(nlev_cas, nt_cas), q2_cas(nlev_cas, nt_cas)

    ! outputs:
    REAL plev_prof_cas(nlev_cas)
    REAL t_prof_cas(nlev_cas), q_prof_cas(nlev_cas)
    REAL u_prof_cas(nlev_cas), v_prof_cas(nlev_cas)
    REAL ug_prof_cas(nlev_cas), vg_prof_cas(nlev_cas)
    REAL vitw_prof_cas(nlev_cas)
    REAL du_prof_cas(nlev_cas), hu_prof_cas(nlev_cas), vu_prof_cas(nlev_cas)
    REAL dv_prof_cas(nlev_cas), hv_prof_cas(nlev_cas), vv_prof_cas(nlev_cas)
    REAL dt_prof_cas(nlev_cas), ht_prof_cas(nlev_cas), vt_prof_cas(nlev_cas)
    REAL dtrad_prof_cas(nlev_cas)
    REAL dq_prof_cas(nlev_cas), hq_prof_cas(nlev_cas), vq_prof_cas(nlev_cas)
    REAL lat_prof_cas, sens_prof_cas, ts_prof_cas, ustar_prof_cas
    REAL uw_prof_cas(nlev_cas), vw_prof_cas(nlev_cas), q1_prof_cas(nlev_cas), q2_prof_cas(nlev_cas)
    ! local:
    INTEGER it_cas1, it_cas2, k
    REAL timeit, time_cas1, time_cas2, frac

    PRINT*, 'Check time', day1, day_ju_ini_cas, day_deb + 1, pdt_cas

    ! On teste si la date du cas AMMA est correcte.
    ! C est pour memoire car en fait les fichiers .def
    ! sont censes etre corrects.
    ! A supprimer a terme (MPL 20150623)
    !     if ((forcing_type.EQ.10).AND.(1.EQ.0)) THEN
    ! Check that initial day of the simulation consistent with AMMA case:
    !      if (annee_ref.NE.2006) THEN
    !       PRINT*,'Pour AMMA, annee_ref doit etre 2006'
    !       PRINT*,'Changer annee_ref dans run.def'
    !       stop
    !      endif
    !      if (annee_ref.EQ.2006 .AND. day1.lt.day_cas) THEN
    !       PRINT*,'AMMA a debute le 10 juillet 2006',day1,day_cas
    !       PRINT*,'Changer dayref dans run.def'
    !       stop
    !      endif
    !      if (annee_ref.EQ.2006 .AND. day1.gt.day_cas+1) THEN
    !       PRINT*,'AMMA a fini le 11 juillet'
    !       PRINT*,'Changer dayref ou nday dans run.def'
    !       stop
    !      endif
    !      endif

    ! Determine timestep relative to the 1st day:
    !       timeit=(day-day1)*86400.
    !       if (annee_ref.EQ.1992) THEN
    !        timeit=(day-day_cas)*86400.
    !       else
    !        timeit=(day+61.-1.)*86400. ! 61 days between Nov01 and Dec31 1992
    !       endif
    timeit = (day - day_ju_ini_cas) * 86400
    PRINT *, 'day=', day
    PRINT *, 'day_ju_ini_cas=', day_ju_ini_cas
    PRINT *, 'pdt_cas=', pdt_cas
    PRINT *, 'timeit=', timeit
    PRINT *, 'nt_cas=', nt_cas

    ! Determine the closest observation times:
    !       it_cas1=INT(timeit/pdt_cas)+1
    !       it_cas2=it_cas1 + 1
    !       time_cas1=(it_cas1-1)*pdt_cas
    !       time_cas2=(it_cas2-1)*pdt_cas

    it_cas1 = INT(timeit / pdt_cas) + 1
    IF (it_cas1 == nt_cas) THEN
      it_cas2 = it_cas1
    ELSE
      it_cas2 = it_cas1 + 1
    ENDIF
    time_cas1 = (it_cas1 - 1) * pdt_cas
    time_cas2 = (it_cas2 - 1) * pdt_cas
    PRINT *, 'it_cas1=', it_cas1
    PRINT *, 'it_cas2=', it_cas2
    PRINT *, 'time_cas1=', time_cas1
    PRINT *, 'time_cas2=', time_cas2

    IF (it_cas1 > nt_cas) THEN
      WRITE(*, *) 'PB-stop: day, day_ju_ini_cas,it_cas1, it_cas2, timeit: '            &
              , day, day_ju_ini_cas, it_cas1, it_cas2, timeit
      stop
    endif

    ! time interpolation:
    IF (it_cas1 == it_cas2) THEN
      frac = 0.
    ELSE
      frac = (time_cas2 - timeit) / (time_cas2 - time_cas1)
      frac = max(frac, 0.0)
    ENDIF

    lat_prof_cas = lat_cas(it_cas2)                                       &
            - frac * (lat_cas(it_cas2) - lat_cas(it_cas1))
    sens_prof_cas = sens_cas(it_cas2)                                     &
            - frac * (sens_cas(it_cas2) - sens_cas(it_cas1))
    ts_prof_cas = ts_cas(it_cas2)                                         &
            - frac * (ts_cas(it_cas2) - ts_cas(it_cas1))
    ustar_prof_cas = ustar_cas(it_cas2)                                   &
            - frac * (ustar_cas(it_cas2) - ustar_cas(it_cas1))

    DO k = 1, nlev_cas
      plev_prof_cas(k) = plev_cas(k, it_cas2)                               &
              - frac * (plev_cas(k, it_cas2) - plev_cas(k, it_cas1))
      t_prof_cas(k) = t_cas(k, it_cas2)                               &
              - frac * (t_cas(k, it_cas2) - t_cas(k, it_cas1))
      q_prof_cas(k) = q_cas(k, it_cas2)                               &
              - frac * (q_cas(k, it_cas2) - q_cas(k, it_cas1))
      u_prof_cas(k) = u_cas(k, it_cas2)                               &
              - frac * (u_cas(k, it_cas2) - u_cas(k, it_cas1))
      v_prof_cas(k) = v_cas(k, it_cas2)                               &
              - frac * (v_cas(k, it_cas2) - v_cas(k, it_cas1))
      ug_prof_cas(k) = ug_cas(k, it_cas2)                               &
              - frac * (ug_cas(k, it_cas2) - ug_cas(k, it_cas1))
      vg_prof_cas(k) = vg_cas(k, it_cas2)                               &
              - frac * (vg_cas(k, it_cas2) - vg_cas(k, it_cas1))
      vitw_prof_cas(k) = vitw_cas(k, it_cas2)                               &
              - frac * (vitw_cas(k, it_cas2) - vitw_cas(k, it_cas1))
      du_prof_cas(k) = du_cas(k, it_cas2)                                   &
              - frac * (du_cas(k, it_cas2) - du_cas(k, it_cas1))
      hu_prof_cas(k) = hu_cas(k, it_cas2)                                   &
              - frac * (hu_cas(k, it_cas2) - hu_cas(k, it_cas1))
      vu_prof_cas(k) = vu_cas(k, it_cas2)                                   &
              - frac * (vu_cas(k, it_cas2) - vu_cas(k, it_cas1))
      dv_prof_cas(k) = dv_cas(k, it_cas2)                                   &
              - frac * (dv_cas(k, it_cas2) - dv_cas(k, it_cas1))
      hv_prof_cas(k) = hv_cas(k, it_cas2)                                   &
              - frac * (hv_cas(k, it_cas2) - hv_cas(k, it_cas1))
      vv_prof_cas(k) = vv_cas(k, it_cas2)                                   &
              - frac * (vv_cas(k, it_cas2) - vv_cas(k, it_cas1))
      dt_prof_cas(k) = dt_cas(k, it_cas2)                                   &
              - frac * (dt_cas(k, it_cas2) - dt_cas(k, it_cas1))
      ht_prof_cas(k) = ht_cas(k, it_cas2)                                   &
              - frac * (ht_cas(k, it_cas2) - ht_cas(k, it_cas1))
      vt_prof_cas(k) = vt_cas(k, it_cas2)                                   &
              - frac * (vt_cas(k, it_cas2) - vt_cas(k, it_cas1))
      dtrad_prof_cas(k) = dtrad_cas(k, it_cas2)                                   &
              - frac * (dtrad_cas(k, it_cas2) - dtrad_cas(k, it_cas1))
      dq_prof_cas(k) = dq_cas(k, it_cas2)                                   &
              - frac * (dq_cas(k, it_cas2) - dq_cas(k, it_cas1))
      hq_prof_cas(k) = hq_cas(k, it_cas2)                                   &
              - frac * (hq_cas(k, it_cas2) - hq_cas(k, it_cas1))
      vq_prof_cas(k) = vq_cas(k, it_cas2)                                   &
              - frac * (vq_cas(k, it_cas2) - vq_cas(k, it_cas1))
      uw_prof_cas(k) = uw_cas(k, it_cas2)                                   &
              - frac * (uw_cas(k, it_cas2) - uw_cas(k, it_cas1))
      vw_prof_cas(k) = vw_cas(k, it_cas2)                                   &
              - frac * (vw_cas(k, it_cas2) - vw_cas(k, it_cas1))
      q1_prof_cas(k) = q1_cas(k, it_cas2)                                   &
              - frac * (q1_cas(k, it_cas2) - q1_cas(k, it_cas1))
      q2_prof_cas(k) = q2_cas(k, it_cas2)                                   &
              - frac * (q2_cas(k, it_cas2) - q2_cas(k, it_cas1))
    enddo

    RETURN
  END

  !**********************************************************************************************
END MODULE mod_1D_cases_read