SUBROUTINE surfini(ngrid,nq,qsurf,psolaralb)

      USE surfdat_h, only: albedodat, albedice
      USE tracer_h, only: igcm_co2_ice
      use comgeomfi_h, only: lati
      use planetwide_mod, only: planetwide_maxval, planetwide_minval

      IMPLICIT NONE
c=======================================================================
c
c   creation des calottes pour l'etat initial
c
c=======================================================================
c-----------------------------------------------------------------------
c   Declarations:
c   -------------
!#include "dimensions.h"
!#include "dimphys.h"
c
      INTEGER,INTENT(IN) :: ngrid
      INTEGER,INTENT(IN) :: nq
      REAL,INTENT(OUT) :: psolaralb(ngrid)
      REAL,INTENT(IN) :: qsurf(ngrid,nq) !tracer on surface (kg/m2)

      INTEGER :: ig,icap
      REAL :: min_albedo,max_albedo
c
c=======================================================================


      DO ig=1,ngrid
         psolaralb(ig)=albedodat(ig)
      ENDDO

      call planetwide_minval(albedodat,min_albedo)
      call planetwide_maxval(albedodat,max_albedo)
      write(*,*) 'surfini: minimum corrected albedo',min_albedo
      write(*,*) 'surfini: maximum corrected albedo',max_albedo

      if (igcm_co2_ice.ne.0) then
! Change Albedo if there is CO2 ice on the surface
        DO ig=1,ngrid
          IF (qsurf(ig,igcm_co2_ice) .GT. 0.) THEN
             IF(lati(ig).LT.0.) THEN
                icap=2 ! Southern hemisphere
             ELSE
                icap=1 ! Northern hemisphere
             ENDIF
             psolaralb(ig) = albedice(icap)
          END IF
        ENDDO ! of DO ig=1,ngrid     
      else
        write(*,*) "surfini: No CO2 ice tracer on surface  ..."
        write(*,*) "         and therefore no albedo change."
      endif

      call planetwide_minval(psolaralb,min_albedo)
      call planetwide_maxval(psolaralb,max_albedo)
      write(*,*) 'surfini: minimum corrected albedo',min_albedo
      write(*,*) 'surfini: maximum corrected albedo',max_albedo

      END