/[MITgcm]/MITgcm_contrib/jscott/igsm/src/kvcarbon.F
ViewVC logotype

Contents of /MITgcm_contrib/jscott/igsm/src/kvcarbon.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (show annotations) (download)
Tue Sep 1 22:08:08 2009 UTC (15 years, 10 months ago) by jscott
Branch: MAIN
new routine for carbon uptake

1 C $Header$
2 C $Name$
3
4 #include "ctrparam.h"
5 subroutine kvcarbon(iyr)
6 #include "ATMSIZE.h"
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.5*pi-atan(yy))/(0.5*pi))
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

  ViewVC Help
Powered by ViewVC 1.1.22