/[MITgcm]/MITgcm_contrib/darwin2/pkg/monod/monod_acdom.F
ViewVC logotype

Contents of /MITgcm_contrib/darwin2/pkg/monod/monod_acdom.F

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

Revision 1.3 - (show annotations) (download)
Tue Apr 16 20:21:57 2013 UTC (12 years, 3 months ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt64k_20130723, ctrb_darwin2_ckpt65w_20160512, ctrb_darwin2_ckpt65j_20150225, ctrb_darwin2_ckpt66g_20170424, ctrb_darwin2_ckpt64h_20130528, ctrb_darwin2_ckpt66k_20171025, ctrb_darwin2_ckpt66n_20180118, ctrb_darwin2_ckpt65v_20160409, ctrb_darwin2_ckpt65s_20160114, ctrb_darwin2_ckpt65_20140718, ctrb_darwin2_ckpt64m_20130820, ctrb_darwin2_ckpt66d_20170214, ctrb_darwin2_ckpt64r_20131210, ctrb_darwin2_ckpt65m_20150615, ctrb_darwin2_ckpt65q_20151118, ctrb_darwin2_ckpt65o_20150914, ctrb_darwin2_ckpt65p_20151023, ctrb_darwin2_ckpt64n_20130826, ctrb_darwin2_ckpt65e_20140929, ctrb_darwin2_ckpt64o_20131024, ctrb_darwin2_ckpt64v_20140411, ctrb_darwin2_ckpt64z_20140711, ctrb_darwin2_ckpt65l_20150504, ctrb_darwin2_ckpt65z_20160929, ctrb_darwin2_ckpt65n_20150729, ctrb_darwin2_ckpt64y_20140622, ctrb_darwin2_ckpt65d_20140915, ctrb_darwin2_ckpt64t_20140202, ctrb_darwin2_ckpt66h_20170602, ctrb_darwin2_ckpt64i_20130622, ctrb_darwin2_ckpt64s_20140105, ctrb_darwin2_ckpt64x_20140524, ctrb_darwin2_ckpt65x_20160612, ctrb_darwin2_ckpt66f_20170407, ctrb_darwin2_ckpt65g_20141120, ctrb_darwin2_ckpt65k_20150402, ctrb_darwin2_ckpt64w_20140502, ctrb_darwin2_ckpt66a_20161020, ctrb_darwin2_ckpt64g_20130503, ctrb_darwin2_ckpt64l_20130806, ctrb_darwin2_ckpt65f_20141014, ctrb_darwin2_ckpt66b_20161219, ctrb_darwin2_ckpt64u_20140308, ctrb_darwin2_ckpt64j_20130704, ctrb_darwin2_ckpt65i_20150123, ctrb_darwin2_ckpt66j_20170815, ctrb_darwin2_ckpt65y_20160801, ctrb_darwin2_ckpt66c_20170121, ctrb_darwin2_ckpt65a_20140728, ctrb_darwin2_ckpt65b_20140812, ctrb_darwin2_ckpt65t_20160221, ctrb_darwin2_ckpt64p_20131118, ctrb_darwin2_ckpt66o_20180209, ctrb_darwin2_ckpt66e_20170314, ctrb_darwin2_ckpt64q_20131118, ctrb_darwin2_ckpt64p_20131024, ctrb_darwin2_ckpt65u_20160315, ctrb_darwin2_ckpt65r_20151221, ctrb_darwin2_ckpt66i_20170718, ctrb_darwin2_ckpt65c_20140830, ctrb_darwin2_ckpt66l_20171025, ctrb_darwin2_ckpt65h_20141217, ctrb_darwin2_ckpt66m_20171213, HEAD
Changes since 1.2: +4 -4 lines 
make reference waveband for aCDOM a runtime parameter
1 C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/monod/monod_acdom.F,v 1.2 2012/05/31 21:08:25 stephd Exp $
2 C $Name: $
3
4 #include "DARWIN_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: RTWB_ACDOM
8
9 c ANNA pass extra variables if WAVEBANDS
10 C !INTERFACE: ==========================================================
11 SUBROUTINE MONOD_ACDOM(
12 #ifdef ALLOW_CDOM
13 I cdomlocal,
14 #else
15 I P_chl, aclocal, awlocal,
16 #endif
17 O acdomlocal,
18 I myThid )
19
20 C !DESCRIPTION:
21 C computes CDOM from Chlorophyll
22 C
23 c ANNA I have changed this a bit from WG code
24 c ------------- this is WG code:
25 c do k = 1,km
26 c actot450 = 0.0
27 c atot450 = 0.0
28 c do n = 1,nchl
29 c actot450 = actot450 + P(k,n)*ac(n,nl450)
30 c enddo
31 c atot450 = aw(nl450) + actot450
32 c do nl = 1,nlt
33 c acdom(k,nl) = 0.2*atot450*excdom(nl)
34 c enddo
35 c enddo
36 c do nl = 1,nlt
37 c write(6,*)'nl,lam,aw,ac,acdom = ',nl,lam(nl),
38 c *aw(nl),ac(1,nl),acdom(1,nl)
39 c enddo
40 c this implies P biomass is being scaled by the absorption of phytoplankton
41 c at 450nm. ie. the waveband on the main chl-a peak.
42 c so (i think) this is effectively getting the absorption by the phytoplankton
43 c In our version, since we have phyto chl we need to scale the
44 c chlorophyll normalised absorption spectra (the inputs)
45 c by phyt chl rather than phyto biomass
46
47 C !USES: ===============================================================
48 IMPLICIT NONE
49 C === Global variables ===
50 #include "SIZE.h"
51 #include "EEPARAMS.h"
52 #include "PARAMS.h"
53 #include "MONOD_SIZE.h"
54 #include "MONOD.h"
55 #include "SPECTRAL_SIZE.h"
56 #ifdef WAVEBANDS
57 #include "WAVEBANDS_PARAMS.h"
58 #endif
59
60 C !INPUT PARAMETERS: ===================================================
61 C P_chl :: Chlorophyll per species and level
62 c ac :: absorption spectra for each phtyo (chl normalised)
63 c awlocal :: absoprtion spectra for water (m-1)
64 C myTime :: Current time in simulation
65 C myIter :: Current iteration number in simulation
66 C myThid :: My Thread Id number
67
68 C !OUTPUT PARAMETERS: ==================================================
69 C acdom :: absortpion spectra for CDOM per level
70
71 #ifdef ALLOW_CDOM
72 _RL cdomlocal(NR)
73 #else
74 _RL P_chl(npmax,Nr)
75 _RL aclocal(npmax,tlam)
76 _RL awlocal(tlam)
77 #endif
78 _RL acdomlocal(Nr,tlam)
79 C _RL myTime
80 C INTEGER myIter
81 INTEGER myThid
82 CEOP
83
84 #ifdef WAVEBANDS
85 #ifdef DAR_CALC_ACDOM
86
87 C !LOCAL VARIABLES:
88 C == Local variables ==
89 INTEGER k, np, ilam
90 _RL actot450, atot450
91 c params in common block
92 c nl450, excdom,
93
94 #ifdef ALLOW_CDOM
95 c use cdom-like tracer
96 do k = 1,Nr
97 do ilam = 1,tlam
98 acdomlocal(k,ilam) = cdomcoeff*cdomlocal(k)*excdom(ilam)
99 enddo
100 enddo
101 #else
102 c ANNA moved cdom calculation from WG's light.f
103 c it's done for RADTRANS and WAVEBANDS_3D
104 do k = 1,Nr
105 actot450 = 0.0 _d 0
106 atot450 = 0.0 _d 0
107 do np = 1,npmax
108 actot450 = actot450 + P_chl(np,k)*aclocal(np,nlaCDOM) !nb. n,k swapped from WG
109 enddo
110 atot450 = awlocal(nlaCDOM) + actot450
111 do ilam = 1,tlam
112 acdomlocal(k,ilam) = darwin_aCDOM_fac * atot450*excdom(ilam)
113 enddo
114 enddo
115 c do nl = 1,tlam
116 c write(6,*)'nl,lam,aw,ac,acdom = ',nl,lam(nl),
117 c *aw(nl),ac(1,nl),acdom(1,nl)
118 c enddo
119 c ANNA cdom end
120 #endif
121
122 #endif /* DAR_CALC_ACDOM */
123 #endif /* WAVEBANDS */
124
125 RETURN
126 END
127
  ViewVC Help
Powered by ViewVC 1.1.22