C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/jscott/igsm/src/kvcarbon.F,v 1.1 2009/09/01 22:08:08 jscott Exp $
C $Name:  $

#include "ctrparam.h"
      subroutine kvcarbon(iyr)
#include "ATMSIZE.h"
#include "DRIVER.h"
      common/Garydiff/depthml(jm0),edzcar(jm0)
      common /Garykvct/cfkvct,edzcart(jm0)
      real depthml,edzcar,cfkvct,edzcart
      REAL pi, yy
      pi=4.*atan(1.)
      cfkvct=1.0
      if (iyr.ge.1991) then
         yy=(iyr-1990)/40.
         cfkvct=sqrt((0.5*pi-atan(yy))/(0.5*pi))
      endif
      print *,'iyr=',iyr,' cfkvct=',cfkvct
       do j=1,jm0
         edzcart(j)=cfkvct*edzcar(j)
      enddo
      return
      end