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	jscott	1.2	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	jscott	1.1
4			#include "ctrparam.h"
5			subroutine kvcarbon(iyr)
6	jscott	1.2	#include "BD2G04.COM"
7	jscott	1.1	#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