!
! $Header$
!
! Initialisations diverses au tout debut
      IF(itap.EQ.1) THEN
         itapm1=0
!        surface terre
         DO i=1, klon
            IF(pctsrf(i,is_ter).GT.0.) THEN
               paire_ter(i)=airephy(i)*pctsrf(i,is_ter)
            ENDIF 
         ENDDO
      ENDIF

! Initialisation debut de mois
      IF(itap.EQ.itapm1+1) THEN
        nday_rain(:)=0.
!       print*,'initialisation mois suivants day_rain itap',itap
      ENDIF

! Calcul fin de journee : total_rain, nday_rain
      IF(MOD(itap,NINT(un_jour/dtime)).EQ.0) THEN
!        print*,'calcul nday_rain itap ',itap
         DO i = 1, klon
            total_rain(i)=rain_fall(i)+snow_fall(i)  
            IF(total_rain(i).GT.0.) nday_rain(i)=nday_rain(i)+1.
         ENDDO
      ENDIF

! Initialisation fin de mois
      IF(MOD(itap-itapm1,NINT(mth_len*un_jour/dtime)).EQ.0) THEN
        itapm1=itapm1+NINT(mth_len*un_jour/dtime)
!       print*,'initialisation itapm1 ',itapm1
      ENDIF