1 |
stephd |
1.2 |
C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/DARWIN_IO.h,v 1.1 2011/04/13 18:56:23 jahn Exp $ |
2 |
|
|
C $Name: $ |
3 |
jahn |
1.1 |
|
4 |
|
|
#include "DARWIN_OPTIONS.h" |
5 |
|
|
|
6 |
|
|
c DARWIN_IO.h |
7 |
|
|
c -------------------------------------------------------------------- |
8 |
|
|
c Description: input and diagnostic output for DARWIN pkg |
9 |
|
|
c Stephanie Dutkiewicz: Spring 2006 |
10 |
|
|
c NOTE: additional diagnostics still needed |
11 |
|
|
c -------------------------------------------------------------------- |
12 |
|
|
|
13 |
|
|
C-- COMMON /DARWIN_FILENAMES/ |
14 |
|
|
C darwin_iceFile :: file name of seaice fraction |
15 |
|
|
C darwin_ironFile :: file name of aeolian iron flux |
16 |
|
|
C darwin_PARFile :: file name of Photosynthetically Active Radiation at surface |
17 |
|
|
C darwin_nutWVelFile :: file name of nutrient wvel |
18 |
|
|
C darwin_waterabsorbFile :: file name of water absorption/scattering spectra |
19 |
|
|
C darwin_phytoabsorbFile :: file name of phyto absorption/scattering spectra |
20 |
|
|
C darwin_particleabsorbFile :: file name of particle absorption/scattering spectra |
21 |
|
|
C darwin_surfacespecFile :: file name of spectrum of incident light |
22 |
|
|
C darwin_acdomFile :: file name of CDOM absorption spectrum |
23 |
|
|
C darwin_forcingPeriod :: perioDARWIN forcing parameter specific for DARWIN (seconds) |
24 |
|
|
C darwin_forcingCycle :: perioDARWIN forcing parameter specific for DARWIN (seconds) |
25 |
|
|
C darwin_seed :: seed for the random number generator |
26 |
|
|
|
27 |
|
|
COMMON /DARWIN_FILENAMES/ |
28 |
|
|
& darwin_iceFile, |
29 |
|
|
& darwin_ironFile, |
30 |
|
|
& darwin_PARFile, |
31 |
|
|
& darwin_nutWVelFile, |
32 |
|
|
& darwin_PO4_relaxFile, darwin_NO3_relaxFile, |
33 |
|
|
& darwin_FeT_relaxFile, darwin_Si_relaxFile, |
34 |
|
|
& darwin_relaxscale, |
35 |
|
|
& darwin_PO4_fluxFile, darwin_NO3_FluxFile, |
36 |
|
|
& darwin_FeT_fluxFile, darwin_Si_fluxFile, |
37 |
|
|
& darwin_oasim_edFile, darwin_oasim_esFile, |
38 |
|
|
& darwin_waterabsorbFile, darwin_phytoabsorbFile, |
39 |
|
|
& darwin_surfacespecFile, darwin_acdomFile, |
40 |
|
|
& darwin_particleabsorbFile, |
41 |
|
|
& darwin_forcingPeriod, darwin_forcingCycle, |
42 |
|
|
& darwin_seed |
43 |
|
|
|
44 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_iceFile |
45 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_ironFile |
46 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_PARFile |
47 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_NutWVelFile |
48 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_PO4_relaxFile |
49 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_NO3_relaxFile |
50 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_FeT_relaxFile |
51 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_Si_relaxFile |
52 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_PO4_fluxFile |
53 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_NO3_fluxFile |
54 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_FeT_fluxFile |
55 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_Si_fluxFile |
56 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_oasim_edFile |
57 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_oasim_esFile |
58 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_waterabsorbFile |
59 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_phytoabsorbFile |
60 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_particleabsorbFile |
61 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_surfacespecFile |
62 |
|
|
CHARACTER*(MAX_LEN_FNAM) darwin_acdomFile |
63 |
|
|
|
64 |
|
|
_RL darwin_relaxscale |
65 |
|
|
_RL darwin_forcingPeriod |
66 |
|
|
_RL darwin_forcingCycle |
67 |
|
|
INTEGER darwin_seed |
68 |
|
|
|
69 |
|
|
c INPUT: |
70 |
|
|
c fice - ice fraction |
71 |
|
|
c inputFe - aeolian input of iron |
72 |
|
|
c sur_par - surface PAR |
73 |
|
|
COMMON /DARWIN_INPUT/ |
74 |
|
|
& fice, inputFe, sur_par |
75 |
|
|
#ifdef NUT_SUPPLY |
76 |
|
|
& ,nut_wvel |
77 |
|
|
#endif |
78 |
|
|
#ifdef RELAX_NUTS |
79 |
|
|
& ,po4_obs, no3_obs, fet_obs, si_obs |
80 |
|
|
#endif |
81 |
|
|
#ifdef FLUX_NUTS |
82 |
|
|
& ,po4_flx, no3_flx, fet_flx, si_flx |
83 |
|
|
#endif |
84 |
|
|
c |
85 |
|
|
_RL fice(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy) |
86 |
|
|
_RL inputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy) |
87 |
|
|
_RL sur_par(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy) |
88 |
|
|
#ifdef NUT_SUPPLY |
89 |
|
|
_RL nut_wvel(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
90 |
|
|
#endif |
91 |
|
|
#ifdef RELAX_NUTS |
92 |
|
|
_RL po4_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
93 |
|
|
_RL no3_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
94 |
|
|
_RL fet_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
95 |
|
|
_RL si_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
96 |
|
|
#endif |
97 |
|
|
#ifdef FLUX_NUTS |
98 |
|
|
_RL po4_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
99 |
|
|
_RL no3_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
100 |
|
|
_RL fet_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
101 |
|
|
_RL si_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
102 |
|
|
#endif |
103 |
|
|
c |
104 |
|
|
c OUPUT DIAGNOSTICS |
105 |
|
|
c PPave - average primary production |
106 |
|
|
c Nfixave - average N fixation |
107 |
|
|
c Zoograzave - average zooplankton consumption |
108 |
|
|
c Parave - average PAR |
109 |
|
|
c Chlave - average chlorophyll |
110 |
|
|
c npzd_timeave - time for averaging |
111 |
|
|
|
112 |
|
|
COMMON /DARWIN_OUPUT/ |
113 |
|
|
& PPave, Nfixave, Zoograzave, |
114 |
|
|
& PARave, Chlave, Denitave, |
115 |
|
|
c ANNA_TAVE |
116 |
|
|
#ifdef WAVES_DIAG_PCHL |
117 |
|
|
& Pchlave, |
118 |
|
|
#endif |
119 |
|
|
#ifdef DAR_DIAG_ACDOM |
120 |
|
|
& aCDOMave, |
121 |
|
|
#endif |
122 |
|
|
c ANNA end TAVE |
123 |
|
|
#ifdef DAR_DIAG_RSTAR |
124 |
|
|
& Rstarave, RNstarave, |
125 |
|
|
#endif |
126 |
|
|
#ifdef DAR_DIAG_DIVER |
127 |
|
|
& Diver1ave, Diver2ave, Diver3ave, Diver4ave, |
128 |
|
|
#endif |
129 |
|
|
#ifdef DAR_DIAG_GROW |
130 |
|
|
& Growave, Growsqave, |
131 |
|
|
#endif |
132 |
|
|
#ifdef ALLOW_DIAZ |
133 |
|
|
#ifdef DAR_DIAG_NFIXP |
134 |
|
|
& NfixPave, |
135 |
|
|
#endif |
136 |
|
|
#endif |
137 |
|
|
& dar_timeave |
138 |
|
|
c |
139 |
|
|
_RL PPave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
140 |
|
|
_RL Nfixave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
141 |
|
|
_RL Zoograzave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
142 |
|
|
_RL PARave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
143 |
|
|
_RL Chlave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
144 |
|
|
_RL Denitave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
145 |
|
|
c ANNA_TAVE |
146 |
|
|
#ifdef WAVES_DIAG_PCHL |
147 |
|
|
_RL Pchlave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax) |
148 |
|
|
#endif |
149 |
|
|
#ifdef DAR_DIAG_ACDOM |
150 |
|
|
_RL aCDOMave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy) |
151 |
|
|
#endif |
152 |
|
|
c ANNA end TAVE |
153 |
|
|
#ifdef DAR_DIAG_RSTAR |
154 |
|
|
_RL Rstarave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax) |
155 |
|
|
_RL RNstarave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax) |
156 |
|
|
#endif |
157 |
|
|
#ifdef DAR_DIAG_DIVER |
158 |
|
|
_RL Diver1ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
159 |
|
|
_RL Diver2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
160 |
|
|
_RL Diver3ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
161 |
|
|
_RL Diver4ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy) |
162 |
|
|
#endif |
163 |
|
|
#ifdef DAR_DIAG_GROW |
164 |
|
|
_RL Growave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax) |
165 |
|
|
_RL Growsqave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax) |
166 |
|
|
#endif |
167 |
|
|
#ifdef ALLOW_DIAZ |
168 |
|
|
#ifdef DAR_DIAG_NFIXP |
169 |
|
|
_RL NfixPave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax) |
170 |
|
|
#endif |
171 |
|
|
#endif |
172 |
|
|
_RL DAR_timeave(nSx,nSy,nR) |
173 |
|
|
|
174 |
|
|
#ifdef CHECK_CONS |
175 |
|
|
COMMON /DARWIN_CHECK_CONS/ |
176 |
|
|
& DAR_cons_unit1, DAR_cons_unit2, |
177 |
|
|
& DAR_cons_unit3, DAR_cons_unit4, |
178 |
|
|
& DAR_cons_unit5, DAR_cons_unit6, |
179 |
|
|
& DAR_cons_unit7 |
180 |
|
|
INTEGER DAR_cons_unit1 |
181 |
|
|
INTEGER DAR_cons_unit2 |
182 |
|
|
INTEGER DAR_cons_unit3 |
183 |
|
|
INTEGER DAR_cons_unit4 |
184 |
|
|
INTEGER DAR_cons_unit5 |
185 |
|
|
INTEGER DAR_cons_unit6 |
186 |
|
|
INTEGER DAR_cons_unit7 |
187 |
|
|
#endif |
188 |
|
|
|
189 |
stephd |
1.2 |
#ifdef CALC_RATE_TOTALS |
190 |
|
|
COMMON /DARWIN_RATE_TOT/ |
191 |
|
|
& DAR_ratetot_unit, DAR_fluxtot_unit |
192 |
|
|
INTEGER DAR_ratetot_unit |
193 |
|
|
INTEGER DAR_fluxtot_unit |
194 |
|
|
#endif |
195 |
|
|
|
196 |
jahn |
1.1 |
#ifdef ALLOW_PAR_DAY |
197 |
|
|
C PARday :: array for accumulating/storing daily-averaged PAR |
198 |
|
|
COMMON /DARWIN_PAR_DAY/ PARday, |
199 |
|
|
& darwin_PARavPeriod, darwin_PARnav |
200 |
|
|
_RL PARday(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,2) |
201 |
|
|
_RL darwin_PARavPeriod |
202 |
|
|
INTEGER darwin_PARnav |
203 |
|
|
#endif |
204 |
|
|
|
205 |
|
|
#ifdef ALLOW_CARBON |
206 |
|
|
C-- COMMON /DIC_FILENAMES/ |
207 |
|
|
C DIC_windFile :: file name of wind speeds |
208 |
|
|
C DIC_atmospFile :: file name of atmospheric pressure |
209 |
|
|
C dic_pCO2 :: Atmospheric pCO2 to be rad in data.dic |
210 |
|
|
C dic_int* :: place holder to read in a integer number, set at run time |
211 |
|
|
|
212 |
|
|
COMMON /DIC_FILENAMES/ |
213 |
|
|
& DIC_windFile, DIC_atmospFile, |
214 |
|
|
& dic_pCO2, dic_int1, dic_int2, dic_int3, dic_int4 |
215 |
|
|
CHARACTER*(MAX_LEN_FNAM) DIC_windFile |
216 |
|
|
CHARACTER*(MAX_LEN_FNAM) DIC_atmospFile |
217 |
|
|
_RL dic_pCO2 |
218 |
|
|
INTEGER dic_int1 |
219 |
|
|
INTEGER dic_int2 |
220 |
|
|
INTEGER dic_int3 |
221 |
|
|
INTEGER dic_int4 |
222 |
|
|
#endif |
223 |
|
|
|