/[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.1 - (hide annotations) (download)
Wed Apr 13 18:56:23 2011 UTC (14 years, 3 months ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt62v_20110413, ctrb_darwin2_ckpt62w_20110426, ctrb_darwin2_baseline
File MIME type: text/plain
darwin2 initial checkin

1 jahn 1.1 C $Header$
2     C $Name$
3    
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     #ifdef ALLOW_PAR_DAY
190     C PARday :: array for accumulating/storing daily-averaged PAR
191     COMMON /DARWIN_PAR_DAY/ PARday,
192     & darwin_PARavPeriod, darwin_PARnav
193     _RL PARday(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,2)
194     _RL darwin_PARavPeriod
195     INTEGER darwin_PARnav
196     #endif
197    
198     #ifdef ALLOW_CARBON
199     C-- COMMON /DIC_FILENAMES/
200     C DIC_windFile :: file name of wind speeds
201     C DIC_atmospFile :: file name of atmospheric pressure
202     C dic_pCO2 :: Atmospheric pCO2 to be rad in data.dic
203     C dic_int* :: place holder to read in a integer number, set at run time
204    
205     COMMON /DIC_FILENAMES/
206     & DIC_windFile, DIC_atmospFile,
207     & dic_pCO2, dic_int1, dic_int2, dic_int3, dic_int4
208     CHARACTER*(MAX_LEN_FNAM) DIC_windFile
209     CHARACTER*(MAX_LEN_FNAM) DIC_atmospFile
210     _RL dic_pCO2
211     INTEGER dic_int1
212     INTEGER dic_int2
213     INTEGER dic_int3
214     INTEGER dic_int4
215     #endif
216    

  ViewVC Help
Powered by ViewVC 1.1.22