IF (ok_isccp) THEN
c
       ndex2d = 0
       ndex3d = 0
c
       itau_w = itau_phy + itap
c
        DO k=1,kmaxm1
         zx_tmp_fi3d(1:klon, 1:lmaxm1)=fq_isccp(1:klon,k,1:lmaxm1)*100.
         CALL gr_fi_ecrit(lmaxm1,klon,iim,jjmp1,zx_tmp_fi3d,
     .                    zx_tmp_3d)
c
cIM: champ 3d : (lon,lat,pres) pour un tau fixe
c
         CALL histwrite(nid_isccp,"cldISCCP_"//taulev(k),itau_w,
     .                  zx_tmp_3d,iim*jjmp1*lmaxm1,ndex3d)
        ENDDO !k
c
        CALL gr_fi_ecrit(1, klon,iim,jjmp1,nbsunlit(1,:),zx_tmp_2d)
        CALL histwrite(nid_isccp,"nsunlit",itau_w,
     .                 zx_tmp_2d,iim*jjmp1,ndex2d)
c
       if (ok_sync) then
        call histsync(nid_isccp)
       endif

      ENDIF !ok_isccp