igsm/src/kvcarbon.F

C $Header: /u/gcmpack/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 "BD2G04.COM"
#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
1	C $Header: /u/gcmpack/MITgcm_contrib/jscott/igsm/src/kvcarbon.F,v 1.1 2009/09/01 22:08:08 jscott Exp $
2	C $Name: $
3
4	#include "ctrparam.h"
5	subroutine kvcarbon(iyr)
6	#include "BD2G04.COM"
7	#include "DRIVER.h"
8	common/Garydiff/depthml(jm0),edzcar(jm0)
9	common /Garykvct/cfkvct,edzcart(jm0)
10	real depthml,edzcar,cfkvct,edzcart
11	REAL pi, yy
12	pi=4.*atan(1.)
13	cfkvct=1.0
14	if (iyr.ge.1991) then
15	yy=(iyr-1990)/40.
16	cfkvct=sqrt((0.5pi-atan(yy))/(0.5pi))
17	endif
18	print *,'iyr=',iyr,' cfkvct=',cfkvct
19	do j=1,jm0
20	edzcart(j)=cfkvct*edzcar(j)
21	enddo
22	return
23	end