/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/DARWIN_IO.h
ViewVC logotype

Annotation of /MITgcm_contrib/darwin2/pkg/darwin/DARWIN_IO.h

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


Revision 1.10 - (hide annotations) (download)
Tue May 19 14:30:15 2015 UTC (10 years, 2 months ago) by benw
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt65w_20160512, ctrb_darwin2_ckpt66g_20170424, ctrb_darwin2_ckpt66k_20171025, ctrb_darwin2_ckpt66n_20180118, ctrb_darwin2_ckpt65v_20160409, ctrb_darwin2_ckpt65s_20160114, ctrb_darwin2_ckpt66d_20170214, ctrb_darwin2_ckpt65m_20150615, ctrb_darwin2_ckpt65q_20151118, ctrb_darwin2_ckpt65o_20150914, ctrb_darwin2_ckpt65p_20151023, ctrb_darwin2_ckpt65z_20160929, ctrb_darwin2_ckpt65n_20150729, ctrb_darwin2_ckpt66h_20170602, ctrb_darwin2_ckpt65x_20160612, ctrb_darwin2_ckpt66f_20170407, ctrb_darwin2_ckpt66a_20161020, ctrb_darwin2_ckpt66b_20161219, ctrb_darwin2_ckpt66j_20170815, ctrb_darwin2_ckpt65y_20160801, ctrb_darwin2_ckpt66c_20170121, ctrb_darwin2_ckpt65t_20160221, ctrb_darwin2_ckpt66o_20180209, ctrb_darwin2_ckpt66e_20170314, ctrb_darwin2_ckpt65u_20160315, ctrb_darwin2_ckpt65r_20151221, ctrb_darwin2_ckpt66i_20170718, ctrb_darwin2_ckpt66l_20171025, ctrb_darwin2_ckpt66m_20171213, HEAD
Changes since 1.9: +34 -7 lines
File MIME type: text/plain
Ben Ward - some superficial structural changes allowing runs with no pfts
         - more significant structural and parameter changes to follow later

1 benw 1.10 C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/DARWIN_IO.h,v 1.9 2013/12/27 17:29:00 jahn Exp $
2 stephd 1.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 jahn 1.3 C darwin_PARunits :: units of data in PAR file; one of Ein/m2/d, uEin/m2/s, W/m2
24     C darwin_PARFileConv :: conversion factor from PAR file units to uEin/m2/s
25     C darwin_W_to_uEins :: conversion factor from W to uEin/s
26 jahn 1.4 C darwin_ironUnits :: units of data in iron file; one of mol Fe/m2/s, kg dust/m2/s
27     C darwin_dustFrac :: fraction of dust by weight that is iron
28 jahn 1.1 C darwin_forcingPeriod :: perioDARWIN forcing parameter specific for DARWIN (seconds)
29     C darwin_forcingCycle :: perioDARWIN forcing parameter specific for DARWIN (seconds)
30     C darwin_seed :: seed for the random number generator
31    
32     COMMON /DARWIN_FILENAMES/
33     & darwin_iceFile,
34     & darwin_ironFile,
35     & darwin_PARFile,
36     & darwin_nutWVelFile,
37     & darwin_PO4_relaxFile, darwin_NO3_relaxFile,
38     & darwin_FeT_relaxFile, darwin_Si_relaxFile,
39     & darwin_relaxscale,
40     & darwin_PO4_fluxFile, darwin_NO3_FluxFile,
41     & darwin_FeT_fluxFile, darwin_Si_fluxFile,
42     & darwin_oasim_edFile, darwin_oasim_esFile,
43     & darwin_waterabsorbFile, darwin_phytoabsorbFile,
44     & darwin_surfacespecFile, darwin_acdomFile,
45     & darwin_particleabsorbFile,
46 jahn 1.4 & darwin_PARunits, darwin_ironUnits,
47     & darwin_PARFileConv, darwin_ironFileConv,
48     & darwin_dustFrac,
49 jahn 1.3 & darwin_W_to_uEins,
50 jahn 1.1 & darwin_forcingPeriod, darwin_forcingCycle,
51     & darwin_seed
52    
53     CHARACTER*(MAX_LEN_FNAM) darwin_iceFile
54     CHARACTER*(MAX_LEN_FNAM) darwin_ironFile
55     CHARACTER*(MAX_LEN_FNAM) darwin_PARFile
56     CHARACTER*(MAX_LEN_FNAM) darwin_NutWVelFile
57     CHARACTER*(MAX_LEN_FNAM) darwin_PO4_relaxFile
58     CHARACTER*(MAX_LEN_FNAM) darwin_NO3_relaxFile
59     CHARACTER*(MAX_LEN_FNAM) darwin_FeT_relaxFile
60     CHARACTER*(MAX_LEN_FNAM) darwin_Si_relaxFile
61     CHARACTER*(MAX_LEN_FNAM) darwin_PO4_fluxFile
62     CHARACTER*(MAX_LEN_FNAM) darwin_NO3_fluxFile
63     CHARACTER*(MAX_LEN_FNAM) darwin_FeT_fluxFile
64     CHARACTER*(MAX_LEN_FNAM) darwin_Si_fluxFile
65     CHARACTER*(MAX_LEN_FNAM) darwin_oasim_edFile
66     CHARACTER*(MAX_LEN_FNAM) darwin_oasim_esFile
67     CHARACTER*(MAX_LEN_FNAM) darwin_waterabsorbFile
68     CHARACTER*(MAX_LEN_FNAM) darwin_phytoabsorbFile
69     CHARACTER*(MAX_LEN_FNAM) darwin_particleabsorbFile
70     CHARACTER*(MAX_LEN_FNAM) darwin_surfacespecFile
71     CHARACTER*(MAX_LEN_FNAM) darwin_acdomFile
72 jahn 1.3 CHARACTER*(MAX_LEN_FNAM) darwin_PARunits
73 jahn 1.4 CHARACTER*(MAX_LEN_FNAM) darwin_ironUnits
74 jahn 1.1
75     _RL darwin_relaxscale
76     _RL darwin_forcingPeriod
77     _RL darwin_forcingCycle
78 jahn 1.3 _RL darwin_PARFileConv
79     _RL darwin_W_to_uEins
80 jahn 1.4 _RL darwin_ironFileConv
81     _RL darwin_dustFrac
82 jahn 1.1 INTEGER darwin_seed
83    
84 jahn 1.8 COMMON/DARWIN_IO_L/
85     & darwin_useiceFile,
86     & darwin_useareamask,
87     & darwin_useiceMask,
88     & darwin_useAREA,
89     & darwin_haveIce
90     LOGICAL darwin_useiceFile
91     LOGICAL darwin_useareamask
92     LOGICAL darwin_useiceMask
93     LOGICAL darwin_useAREA
94     LOGICAL darwin_haveIce
95    
96 jahn 1.1 c INPUT:
97     c fice - ice fraction
98     c inputFe - aeolian input of iron
99     c sur_par - surface PAR
100     COMMON /DARWIN_INPUT/
101     & fice, inputFe, sur_par
102     #ifdef NUT_SUPPLY
103     & ,nut_wvel
104     #endif
105     #ifdef RELAX_NUTS
106     & ,po4_obs, no3_obs, fet_obs, si_obs
107     #endif
108     #ifdef FLUX_NUTS
109     & ,po4_flx, no3_flx, fet_flx, si_flx
110     #endif
111     c
112     _RL fice(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy)
113     _RL inputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy)
114     _RL sur_par(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy)
115     #ifdef NUT_SUPPLY
116     _RL nut_wvel(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
117     #endif
118     #ifdef RELAX_NUTS
119     _RL po4_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
120     _RL no3_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
121     _RL fet_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
122     _RL si_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
123     #endif
124     #ifdef FLUX_NUTS
125     _RL po4_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
126     _RL no3_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
127     _RL fet_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
128     _RL si_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
129     #endif
130     c
131     c OUPUT DIAGNOSTICS
132     c PPave - average primary production
133     c Nfixave - average N fixation
134     c Zoograzave - average zooplankton consumption
135     c Parave - average PAR
136     c Chlave - average chlorophyll
137 benw 1.5 c npzd_timeave - time for averaging
138     c Nlim,Flim,etc- Quota model average limitation factors
139 jahn 1.1
140     COMMON /DARWIN_OUPUT/
141     & PPave, Nfixave, Zoograzave,
142     & PARave, Chlave, Denitave,
143 benw 1.5 #ifdef QUOTA_DIAG_LIMIT
144 benw 1.10 & Tlimave,
145     & Ilimave,AP_C_ave,HP_C_ave,
146     & Nlimave,AP_N_ave,HP_N_ave,
147     #ifdef PQUOTA
148     & Plimave,AP_P_ave,HP_P_ave,
149     #endif
150     #ifdef FQUOTA
151     & Flimave,AP_F_ave,HP_F_ave,
152     #endif
153 benw 1.5 #endif
154 jahn 1.1 c ANNA_TAVE
155     #ifdef WAVES_DIAG_PCHL
156     & Pchlave,
157     #endif
158     #ifdef DAR_DIAG_ACDOM
159     & aCDOMave,
160     #endif
161 stephd 1.6 #ifdef DAR_DIAG_PARW
162     & chl2cave,
163     #endif
164     #ifdef DAR_DIAG_EK
165     & Ekave, EkoverEave, acclimave,
166     #endif
167 jahn 1.1 c ANNA end TAVE
168     #ifdef DAR_DIAG_RSTAR
169     & Rstarave, RNstarave,
170     #endif
171     #ifdef DAR_DIAG_DIVER
172     & Diver1ave, Diver2ave, Diver3ave, Diver4ave,
173     #endif
174     #ifdef DAR_DIAG_GROW
175     & Growave, Growsqave,
176     #endif
177     #ifdef ALLOW_DIAZ
178     #ifdef DAR_DIAG_NFIXP
179     & NfixPave,
180     #endif
181     #endif
182     & dar_timeave
183     c
184     _RL PPave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
185     _RL Nfixave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
186     _RL Zoograzave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
187     _RL PARave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
188     _RL Chlave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
189     _RL Denitave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
190 benw 1.10 #ifdef QUOTA_DIAG_LIMIT
191     c Carbon
192     _RL AP_C_ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
193     _RL HP_C_ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
194     c Nitrogen
195     _RL AP_N_ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
196     _RL HP_N_ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
197     _RL Nlimave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
198     c Phosphorus
199     #ifdef PQUOTA
200     _RL AP_P_ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
201     _RL HP_P_ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
202     _RL Plimave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
203     #endif
204     c Iron
205     #ifdef FQUOTA
206     _RL AP_F_ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
207     _RL HP_F_ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
208     _RL Flimave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
209     #endif
210     c Light
211     _RL Ilimave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,npmax)
212     c Temperature
213     _RL Tlimave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
214 benw 1.5 #endif
215 jahn 1.1 c ANNA_TAVE
216     #ifdef WAVES_DIAG_PCHL
217     _RL Pchlave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
218     #endif
219     #ifdef DAR_DIAG_ACDOM
220     _RL aCDOMave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
221     #endif
222 stephd 1.6 #ifdef DAR_DIAG_PARW
223     _RL chl2cave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
224     #endif
225     #ifdef DAR_DIAG_EK
226     _RL Ekave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
227     _RL EkoverEave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
228     _RL acclimave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
229     #endif
230 jahn 1.1 c ANNA end TAVE
231     #ifdef DAR_DIAG_RSTAR
232     _RL Rstarave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
233     _RL RNstarave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
234     #endif
235     #ifdef DAR_DIAG_DIVER
236     _RL Diver1ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
237     _RL Diver2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
238     _RL Diver3ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
239     _RL Diver4ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
240     #endif
241     #ifdef DAR_DIAG_GROW
242     _RL Growave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
243     _RL Growsqave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
244     #endif
245     #ifdef ALLOW_DIAZ
246     #ifdef DAR_DIAG_NFIXP
247     _RL NfixPave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
248     #endif
249     #endif
250 jahn 1.9 _RL DAR_timeave(nSx,nSy)
251 jahn 1.1
252     #ifdef CHECK_CONS
253     COMMON /DARWIN_CHECK_CONS/
254 jahn 1.7 & DAR_cons_unitC, DAR_cons_unitN,
255     & DAR_cons_unitP, DAR_cons_unitF,
256     & DAR_cons_unitS, DAR_cons_unitA,
257     & DAR_cons_unitO
258     INTEGER DAR_cons_unitC
259     INTEGER DAR_cons_unitN
260     INTEGER DAR_cons_unitP
261     INTEGER DAR_cons_unitF
262     INTEGER DAR_cons_unitS
263     INTEGER DAR_cons_unitA
264     INTEGER DAR_cons_unitO
265 jahn 1.1 #endif
266    
267 stephd 1.2 #ifdef CALC_RATE_TOTALS
268     COMMON /DARWIN_RATE_TOT/
269     & DAR_ratetot_unit, DAR_fluxtot_unit
270     INTEGER DAR_ratetot_unit
271     INTEGER DAR_fluxtot_unit
272     #endif
273    
274 jahn 1.1 #ifdef ALLOW_PAR_DAY
275     C PARday :: array for accumulating/storing daily-averaged PAR
276     COMMON /DARWIN_PAR_DAY/ PARday,
277     & darwin_PARavPeriod, darwin_PARnav
278     _RL PARday(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,2)
279     _RL darwin_PARavPeriod
280     INTEGER darwin_PARnav
281     #endif
282    
283     #ifdef ALLOW_CARBON
284     C-- COMMON /DIC_FILENAMES/
285     C DIC_windFile :: file name of wind speeds
286     C DIC_atmospFile :: file name of atmospheric pressure
287     C dic_pCO2 :: Atmospheric pCO2 to be rad in data.dic
288     C dic_int* :: place holder to read in a integer number, set at run time
289    
290     COMMON /DIC_FILENAMES/
291     & DIC_windFile, DIC_atmospFile,
292     & dic_pCO2, dic_int1, dic_int2, dic_int3, dic_int4
293     CHARACTER*(MAX_LEN_FNAM) DIC_windFile
294     CHARACTER*(MAX_LEN_FNAM) DIC_atmospFile
295     _RL dic_pCO2
296     INTEGER dic_int1
297     INTEGER dic_int2
298     INTEGER dic_int3
299     INTEGER dic_int4
300     #endif
301    

  ViewVC Help
Powered by ViewVC 1.1.22