!       if (sorties) then
      print*,'Debut des wrgradsfi'

!      print*,'16 OK convect8'
         call wrgradsfi(1,nlay,fraca,'fraca     ','fraca     ')
         call wrgradsfi(1,nlay,pt,'pt        ','pt        ')
         call wrgradsfi(1,nlay,zh,'zh        ','zh        ')
         call wrgradsfi(1,nlay,zha,'zha        ','zha        ')
         call wrgradsfi(1,nlay,zua,'zua        ','zua        ')
         call wrgradsfi(1,nlay,zva,'zva        ','zva        ')
         call wrgradsfi(1,nlay,zu,'zu        ','zu        ')
         call wrgradsfi(1,nlay,zv,'zv        ','zv        ')
         call wrgradsfi(1,nlay,zo,'zo        ','zo        ')
         call wrgradsfi(1,1,zmax,'zmax      ','zmax      ')
         call wrgradsfi(1,nlay,zdhadj,'zdhadj    ','zdhadj    ')
         call wrgradsfi(1,nlay,pduadj,'pduadj    ','pduadj    ')
         call wrgradsfi(1,nlay,pdvadj,'pdvadj    ','pdvadj    ')
         call wrgradsfi(1,nlay,pdoadj,'pdoadj    ','pdoadj    ')
         call wrgradsfi(1,nlay,entr,'entr      ','entr      ')
         call wrgradsfi(1,nlay,detr,'detr      ','detr      ')
         call wrgradsfi(1,nlay,fm,'fm        ','fm        ')
         call wrgradsfi(1,nlay,zw2,'zw2       ','zw2       ')
         call wrgradsfi(1,nlay,w_est,'w_est      ','w_est      ')
!on sort les moments
         call wrgradsfi(1,nlay,thetath2,'zh2       ','zh2       ')
         call wrgradsfi(1,nlay,wth2,'w2       ','w2       ')
         call wrgradsfi(1,nlay,wth3,'w3       ','w3       ')
         call wrgradsfi(1,nlay,q2,'q2       ','q2       ')
!
         call wrgradsfi(1,nlay,ztva,'ztva      ','ztva      ')
         call wrgradsfi(1,nlay,ztv,'ztv       ','ztv       ')

         call wrgradsfi(1,nlay,zo,'zo        ','zo        ')
         call wrgradsfi(1,nlay,zoa,'zoa        ','zoa        ')

!nouveaux diagnostiques
         call wrgradsfi(1,nlay,zthl,'zthl        ','zthl        ')
         call wrgradsfi(1,nlay,zta,'zta        ','zta        ')
         call wrgradsfi(1,nlay,zl,'zl        ','zl        ')
         call wrgradsfi(1,nlay,zdthladj,'zdthladj    ',  &
     &        'zdthladj    ')
         call wrgradsfi(1,nlay,ztla,'ztla      ','ztla      ')
         call wrgradsfi(1,nlay,zqta,'zqta      ','zqta      ')
         call wrgradsfi(1,nlay,zqla,'zqla      ','zqla      ')
         call wrgradsfi(1,nlay,deltaz,'deltaz      ','deltaz      ')
!nouveaux diagnostiques
      call wrgradsfi(1,nlay,entr_star  ,'entr_star   ','entr_star   ')
      call wrgradsfi(1,nlay,detr_star  ,'detr_star   ','detr_star   ')     
      call wrgradsfi(1,nlay,f_star    ,'f_star   ','f_star   ')
      call wrgradsfi(1,nlay,zqsat    ,'zqsat   ','zqsat   ')
      call wrgradsfi(1,nlay,zqsatth    ,'qsath   ','qsath   ')
      call wrgradsfi(1,nlay,alim_star    ,'alim_star   ','alim_star   ')
      call wrgradsfi(1,nlay,alim    ,'alim   ','alim   ')
      call wrgradsfi(1,1,f,'f      ','f      ')
      call wrgradsfi(1,1,alim_star_tot,'a_s_t      ','a_s_t      ')
      call wrgradsfi(1,1,alim_star2,'a_2      ','a_2      ')
      call wrgradsfi(1,1,zmax,'zmax      ','zmax      ')
      call wrgradsfi(1,1,zmax_sec,'z_sec      ','z_sec      ')
      call wrgradsfi(1,1,zmix,'zmix      ','zmix      ') 
      call wrgradsfi(1,1,nivcon,'nivcon      ','nivcon      ')
      call wrgradsfi(1,1,zcon,'zcon      ','zcon      ')
      call wrgradsfi(1,1,zcon2,'zcon2      ','zcon2      ')
      zsortie1d(:)=lmax(:)
      call wrgradsfi(1,1,zsortie1d,'lmax      ','lmax      ')
      call wrgradsfi(1,1,wmax,'wmax      ','wmax      ')
      call wrgradsfi(1,1,wmax_sec,'w_sec      ','w_sec      ')
      zsortie1d(:)=lmix(:)
      call wrgradsfi(1,1,zsortie1d,'lmix      ','lmix      ')
      zsortie1d(:)=lentr(:)
      call wrgradsfi(1,1,zsortie1d,'lentr      ','lentr     ')

      print*,'Fin des wrgradsfi'