Index: LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.f90
===================================================================
--- LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.f90	(revision 5508)
+++ LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.f90	(revision 5509)
@@ -1589,5 +1589,5 @@
 !=======================================================================
   SUBROUTINE guide_read(timestep)
-    USE netcdf, ONLY: nf90_put_var
+    USE netcdf, ONLY: nf90_get_var
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h
@@ -1803,9 +1803,15 @@
          endif
 ! Coefs ap, bp pour calcul de la pression aux differents niveaux
+!function nf90_get_var(ncid, varid, values, start, count, stride, map)
+!  integer,                         intent( in) :: ncid, varid
+!  any valid type, scalar or array of any rank, &
+!                                   intent(out) :: values
+!  integer, dimension(:), optional, intent( in) :: start, count, stride, map
+!  integer                                      :: nf90_get_var
          IF (guide_plevs.EQ.1) THEN
-             status = nf90_put_var(ncidpl, varidap, apnc, [1], [nlevnc])
-             status = nf90_put_var(ncidpl, varidbp, bpnc, [1], [nlevnc])
+             status = nf90_get_var(ncidpl, varidap, apnc, [1], [nlevnc])
+             status = nf90_get_var(ncidpl, varidbp, bpnc, [1], [nlevnc])
          ELSEIF (guide_plevs.EQ.0) THEN
-             status = nf90_put_var(ncidpl, varidpl, apnc, [1], [nlevnc])
+             status = nf90_get_var(ncidpl, varidpl, apnc, [1], [nlevnc])
 !FC Pour les corrections la pression est deja en Pascals on commente la ligne ci-dessous
              IF(convert_Pa) apnc=apnc*100.! conversion en Pascals
@@ -1833,5 +1839,5 @@
 ! Pression
      if (guide_plevs.EQ.2) then
-         status = nf90_put_var(ncidp, varidp, pnat2, start, count)
+         status = nf90_get_var(ncidp, varidp, pnat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1843,5 +1849,5 @@
 !  Vent zonal
      if (guide_u) then
-         status = nf90_put_var(ncidu, varidu, unat2, start, count)
+         status = nf90_get_var(ncidu, varidu, unat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1855,5 +1861,5 @@
 !  Temperature
      if (guide_T) then
-         status = nf90_put_var(ncidt, varidt, tnat2, start, count)
+         status = nf90_get_var(ncidt, varidt, tnat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1865,5 +1871,5 @@
 !  Humidite
      if (guide_Q) then
-         status = nf90_put_var(ncidQ, varidQ, qnat2, start, count)
+         status = nf90_get_var(ncidQ, varidQ, qnat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1879,5 +1885,5 @@
          count(2)=jjnb_v
          IF (invert_y) start(2)=jjm-jje_v+1
-         status = nf90_put_var(ncidv, varidv, vnat2, start, count)
+         status = nf90_get_var(ncidv, varidv, vnat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1896,5 +1902,5 @@
          count(4)=0
          IF (invert_y) start(2)=jjp1-jje_u+1
-         status = nf90_put_var(ncidps, varidps, psnat2, start, count)
+         status = nf90_get_var(ncidps, varidps, psnat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1908,5 +1914,5 @@
 !=======================================================================
   SUBROUTINE guide_read2D(timestep)
-    USE netcdf, ONLY: nf90_put_var
+    USE netcdf, ONLY: nf90_get_var
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h
@@ -2057,8 +2063,8 @@
 ! Coefs ap, bp pour calcul de la pression aux differents niveaux
          if (guide_plevs.EQ.1) then
-             status = nf90_put_var(ncidpl, varidap, apnc, [1], [nlevnc])
-             status = nf90_put_var(ncidpl, varidbp, bpnc, [1], [nlevnc])
+             status = nf90_get_var(ncidpl, varidap, apnc, [1], [nlevnc])
+             status = nf90_get_var(ncidpl, varidbp, bpnc, [1], [nlevnc])
          elseif (guide_plevs.EQ.0) THEN
-             status = nf90_put_var(ncidpl, varidpl, apnc, [1], [nlevnc])
+             status = nf90_get_var(ncidpl, varidpl, apnc, [1], [nlevnc])
              apnc=apnc*100.! conversion en Pascals
              bpnc(:)=0.
@@ -2085,5 +2091,5 @@
 !  Pression
      if (guide_plevs.EQ.2) then
-         status = nf90_put_var(ncidp, varidp, zu, start, count)
+         status = nf90_get_var(ncidp, varidp, zu, start, count)
          DO i=1,iip1
              pnat2(i,:,:)=zu(:,:)
@@ -2098,5 +2104,5 @@
 !  Vent zonal
      if (guide_u) then
-         status = nf90_put_var(ncidu, varidu, zu, start, count)
+         status = nf90_get_var(ncidu, varidu, zu, start, count)
          DO i=1,iip1
              unat2(i,:,:)=zu(:,:)
@@ -2113,5 +2119,5 @@
 !  Temperature
      if (guide_T) then
-         status = nf90_put_var(ncidt, varidt, zu, start, count)
+         status = nf90_get_var(ncidt, varidt, zu, start, count)
          DO i=1,iip1
              tnat2(i,:,:)=zu(:,:)
@@ -2127,5 +2133,5 @@
 !  Humidite
      if (guide_Q) then
-         status = nf90_put_var(ncidQ, varidQ, zu, start, count)
+         status = nf90_get_var(ncidQ, varidQ, zu, start, count)
          DO i=1,iip1
              qnat2(i,:,:)=zu(:,:)
@@ -2144,5 +2150,5 @@
          count(2)=jjnb_v
          IF (invert_y) start(2)=jjm-jje_v+1
-         status = nf90_put_var(ncidv, varidv, zv, start, count)
+         status = nf90_get_var(ncidv, varidv, zv, start, count)
          DO i=1,iip1
              vnat2(i,:,:)=zv(:,:)
@@ -2166,5 +2172,5 @@
          count(4)=0
          IF (invert_y) start(2)=jjp1-jje_u+1
-         status = nf90_put_var(ncidps, varidps, zu(:, 1), start, count)
+         status = nf90_get_var(ncidps, varidps, zu(:, 1), start, count)
          DO i=1,iip1
              psnat2(i,:)=zu(:,1)