Index: LMDZ4/trunk/libf/dyn3d/wrgrads.F
===================================================================
--- LMDZ4/trunk/libf/dyn3d/wrgrads.F	(revision 1024)
+++ LMDZ4/trunk/libf/dyn3d/wrgrads.F	(revision 1025)
@@ -17,4 +17,8 @@
       integer if,nl
       real field(imx*jmx*lmx)
+
+      integer, parameter:: wp = selected_real_kind(p=6, r=36)
+      real(wp) field4(imx*jmx*lmx)
+
       character*10 name,file
       character*10 titlevar
@@ -83,4 +87,5 @@
 c     print*,'ivar(if),nvar(if),var(ivar(if),if),writectl'
 c     print*,ivar(if),nvar(if),var(ivar(if),if),writectl
+      field4(1:imd(if)*jmd(if)*nl)=field(1:imd(if)*jmd(if)*nl)
       do l=1,nl
          irec(if)=irec(if)+1
@@ -89,5 +94,5 @@
 c    s ,(l-1)*imd(if)*jmd(if)+(ijf-1)*imd(if)+iif
          write(unit(if)+1,rec=irec(if))
-     s   ((field((l-1)*imd(if)*jmd(if)+(j-1)*imd(if)+i)
+     s   ((field4((l-1)*imd(if)*jmd(if)+(j-1)*imd(if)+i)
      s   ,i=iii,iif),j=iji,ijf)
       enddo