Index: /LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90	(revision 5071)
+++ /LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90	(revision 5072)
@@ -9,5 +9,5 @@
   USE Write_Field_loc
   use lmdz_netcdf, only: nf90_nowrite, nf90_open, nf90_inq_varid, nf90_close, &
-          nf90_inq_dimid, nf90_inquire_dimension, nf_get_vara_rd, nf_inq_dimid, &
+          nf90_inq_dimid, nf90_inquire_dimension, nf_inq_dimid, &
           nf_inq_dimlen, nf_enddef, nf_def_dim, nf90_put_var, nf_noerr, nf_close, nf_inq_varid, &
           nf_redef, nf_write, nf_unlimited, nf_float, nf_clobber, nf_64bit_offset, nf90_float, &
@@ -1786,8 +1786,8 @@
 ! Coefs ap, bp pour calcul de la pression aux differents niveaux
          IF (guide_plevs==1) THEN
-             status=nf_get_vara_rd(ncidpl,varidap,[1],[nlevnc],apnc)
-             status=nf_get_vara_rd(ncidpl,varidbp,[1],[nlevnc],bpnc)
+             status=nf90_put_var(ncidpl,varidap,apnc,[1],[nlevnc])
+             status=nf90_put_var(ncidpl,varidbp,bpnc,[1],[nlevnc])
          ELSEIF (guide_plevs==0) THEN
-             status=nf_get_vara_rd(ncidpl,varidpl,[1],[nlevnc],apnc)
+             status=nf90_put_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
@@ -1815,5 +1815,5 @@
 ! Pression 
      if (guide_plevs==2) then
-         status=nf_get_vara_rd(ncidp,varidp,start,count,pnat2)
+         status=nf90_put_var(ncidp,varidp,pnat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1825,5 +1825,5 @@
 !  Vent zonal
      if (guide_u) then
-         status=nf_get_vara_rd(ncidu,varidu,start,count,unat2)
+         status=nf90_put_var(ncidu,varidu,unat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1837,5 +1837,5 @@
 !  Temperature
      if (guide_T) then
-         status=nf_get_vara_rd(ncidt,varidt,start,count,tnat2)
+         status=nf90_put_var(ncidt,varidt,tnat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1847,5 +1847,5 @@
 !  Humidite
      if (guide_Q) then
-         status=nf_get_vara_rd(ncidQ,varidQ,start,count,qnat2)
+         status=nf90_put_var(ncidQ,varidQ,qnat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1862,5 +1862,5 @@
          IF (invert_y) start(2)=jjm-jje_v+1
 
-         status=nf_get_vara_rd(ncidv,varidv,start,count,vnat2)
+         status=nf90_put_var(ncidv,varidv,vnat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1879,5 +1879,5 @@
          count(4)=0
          IF (invert_y) start(2)=jjp1-jje_u+1
-         status=nf_get_vara_rd(ncidps,varidps,start,count,psnat2)
+         status=nf90_put_var(ncidps,varidps,psnat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -2039,8 +2039,8 @@
 ! Coefs ap, bp pour calcul de la pression aux differents niveaux
          if (guide_plevs==1) then
-             status=nf_get_vara_rd(ncidpl,varidap,[1],[nlevnc],apnc)
-             status=nf_get_vara_rd(ncidpl,varidbp,[1],[nlevnc],bpnc)
+             status=nf90_put_var(ncidpl,varidap,apnc,[1],[nlevnc])
+             status=nf90_put_var(ncidpl,varidbp,bpnc,[1],[nlevnc])
          elseif (guide_plevs==0) THEN
-             status=nf_get_vara_rd(ncidpl,varidpl,[1],[nlevnc],apnc)
+             status=nf90_put_var(ncidpl,varidpl,apnc,[1],[nlevnc])
              apnc=apnc*100.! conversion en Pascals
              bpnc(:)=0.
@@ -2067,5 +2067,5 @@
 !  Pression
      if (guide_plevs==2) then
-         status=nf_get_vara_rd(ncidp,varidp,start,count,zu)
+         status=nf90_put_var(ncidp,varidp,zu,start,count)
          DO i=1,iip1
              pnat2(i,:,:)=zu(:,:)
@@ -2080,5 +2080,5 @@
 !  Vent zonal
      if (guide_u) then
-         status=nf_get_vara_rd(ncidu,varidu,start,count,zu)
+         status=nf90_put_var(ncidu,varidu,zu,start,count)
          DO i=1,iip1
              unat2(i,:,:)=zu(:,:)
@@ -2095,5 +2095,5 @@
 !  Temperature
      if (guide_T) then
-         status=nf_get_vara_rd(ncidt,varidt,start,count,zu)
+         status=nf90_put_var(ncidt,varidt,zu,start,count)
          DO i=1,iip1
              tnat2(i,:,:)=zu(:,:)
@@ -2109,5 +2109,5 @@
 !  Humidite
      if (guide_Q) then
-         status=nf_get_vara_rd(ncidQ,varidQ,start,count,zu)
+         status=nf90_put_var(ncidQ,varidQ,zu,start,count)
          DO i=1,iip1
              qnat2(i,:,:)=zu(:,:)
@@ -2126,5 +2126,5 @@
          count(2)=jjnb_v
          IF (invert_y) start(2)=jjm-jje_v+1
-         status=nf_get_vara_rd(ncidv,varidv,start,count,zv)
+         status=nf90_put_var(ncidv,varidv,zv,start,count)
          DO i=1,iip1
              vnat2(i,:,:)=zv(:,:)
@@ -2148,5 +2148,5 @@
          count(4)=0
          IF (invert_y) start(2)=jjp1-jje_u+1
-         status=nf_get_vara_rd(ncidps,varidps,start,count,zu(:,1))
+         status=nf90_put_var(ncidps,varidps,zu(:,1),start,count)
          DO i=1,iip1
              psnat2(i,:)=zu(:,1)