/[MITgcm]/manual/s_phys_pkgs/text/exf.tex
ViewVC logotype

Diff of /manual/s_phys_pkgs/text/exf.tex

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

revision 1.2 by heimbach, Thu Jul 14 22:04:52 2005 UTC revision 1.3 by heimbach, Fri Jul 15 20:08:42 2005 UTC
# Line 59  Table \ref{tab:pkg:exf:cpp} summarizes t Line 59  Table \ref{tab:pkg:exf:cpp} summarizes t
59          \texttt{ALLOW\_ATM\_WIND} &          \texttt{ALLOW\_ATM\_WIND} &
60            compute wind stress from wind speed input\\            compute wind stress from wind speed input\\
61          \texttt{ALLOW\_BULKFORMULAE} &          \texttt{ALLOW\_BULKFORMULAE} &
62            is used if either ALLOW\_ATM\_TEMP or ALLOW\_ATM\_WIND is enabled \\            is used if \texttt{ALLOW\_ATM\_TEMP} or
63              \texttt{ALLOW\_ATM\_WIND} is enabled \\
64          \texttt{EXF\_READ\_EVAP} & read evaporation instead of computing it \\          \texttt{EXF\_READ\_EVAP} & read evaporation instead of computing it \\
65          \texttt{ALLOW\_RUNOFF} & read time-constant river/glacier run-off field \\          \texttt{ALLOW\_RUNOFF} & read time-constant river/glacier run-off field \\
66          \texttt{ALLOW\_DOWNWARD\_RADIATION} & compute net from downward or downward from net radiation \\          \texttt{ALLOW\_DOWNWARD\_RADIATION} & compute net from downward or downward from net radiation \\
67          \texttt{USE\_EXF\_INTERPOLATION} & enable on-the-fly bilinear or bicubic interpolation of input fields \\          \texttt{USE\_EXF\_INTERPOLATION} & enable on-the-fly bilinear or bicubic interpolation of input fields \\
68        \hline        \hline
69             \multicolumn{2}{|c|}{\textit{used in conjunction with relaxation to prescribed (climatological) fields}} \\
70             \hline
71          \texttt{ALLOW\_CLIMTEMP\_RELAXATION} &          \texttt{ALLOW\_CLIMTEMP\_RELAXATION} &
72            relaxation to 3-D potential temperature climatology \\            relaxation to 3-D temperature climatology \\
73          \texttt{ALLOW\_CLIMSALT\_RELAXATION} &          \texttt{ALLOW\_CLIMSALT\_RELAXATION} &
74            relaxation to 3-D salinity climatology \\            relaxation to 3-D salinity climatology \\
75          \texttt{ALLOW\_CLIMSST\_RELAXATION} &          \texttt{ALLOW\_CLIMSST\_RELAXATION} &
# Line 74  Table \ref{tab:pkg:exf:cpp} summarizes t Line 77  Table \ref{tab:pkg:exf:cpp} summarizes t
77          \texttt{ALLOW\_CLIMSSS\_RELAXATION} &          \texttt{ALLOW\_CLIMSSS\_RELAXATION} &
78            relaxation to 2-D SSS climatology  \\            relaxation to 2-D SSS climatology  \\
79        \hline        \hline
80          \texttt{SHORTWAVE\_HEATING} & in \texttt{CPP\_OPTIONS.h}: enable shortwave radiation \\           \multicolumn{2}{|c|}{\textit{these are set outside of EXF in} \texttt{CPP\_OPTIONS.h}} \\
81          \texttt{ATMOSPHERIC\_LOADING} &  in \texttt{CPP\_OPTIONS.h}: enable surface pressure forcing \\           \hline
82            \texttt{SHORTWAVE\_HEATING} & enable shortwave radiation \\
83            \texttt{ATMOSPHERIC\_LOADING} & enable surface pressure forcing \\
84        \hline        \hline
85      \end{tabular}      \end{tabular}
86    }    }
# Line 89  Table \ref{tab:pkg:exf:cpp} summarizes t Line 94  Table \ref{tab:pkg:exf:cpp} summarizes t
94  \label{sec:pkg:exf:runtime}}  \label{sec:pkg:exf:runtime}}
95    
96  Run-time parameters are set in files \texttt{data.pkg},  Run-time parameters are set in files \texttt{data.pkg},
97  and \texttt{data.pkg\_clim} (for relaxation/climatological fields).  and \texttt{data.pkg\_clim} (for relaxation/climatological fields)
98    which are read in \texttt{exf\_readparms.F}.
99  Run-time parameters may be broken into 2 categories:  Run-time parameters may be broken into 2 categories:
100  (i) general flags and parameters, and  (i) general flags and parameters, and
101  (ii) attributes for each forcing and climatological field.  (ii) attributes for each forcing and climatological field.
102    
103  \paragraph{General flags and parameters}  \paragraph{General flags and parameters}
104    
105    ~
106    
107  \begin{table}[h!]  \begin{table}[h!]
108    \label{tab:pkg:exf:runtime_flags}    \label{tab:pkg:exf:runtime_flags}
109    {\footnotesize    {\footnotesize
110      \begin{tabular}{|l|cl|}      \begin{tabular}{|l|c|l|}
111        \hline        \hline
112        \textbf{Flag/parameter} & \textbf{default} &  \textbf{Description}  \\        \textbf{Flag/parameter} & \textbf{default} &  \textbf{Description}  \\
113        \hline \hline        \hline \hline
114          useExfCheckRange & \texttt{.TRUE.} & ~ \\          useExfCheckRange & \texttt{.TRUE.} &
115          useExfYearlyFields & \texttt{.FALSE.} & ~ \\             check range of input fields and stop if out of range \\
116          twoDigitYear & \texttt{.FALSE.} & ~ \\          useExfYearlyFields & \texttt{.FALSE.} &
117          repeatPeriod & \texttt{0.0} & ~ \\             append current year postfix of form \texttt{\_YYYY} on filename \\
118          windstressmax & \texttt{2.0} & ~ \\          twoDigitYear & \texttt{.FALSE.} &
119          exf\_albedo & \texttt{0.1} & ~ \\             instead of appending \texttt{\_YYYY} append  \texttt{YY} \\
120          exf\_iprec  & \texttt{32} & ~ \\          repeatPeriod & \texttt{0.0} & $ > 0 $ :
121          exf\_yftype & \texttt{'RL'} & ~ \\             cycle through all input fields at the same period (in seconds) \\
122            ~            & ~            & $ = 0 $ :
123               use period assigned to each field \\
124            exf\_offset\_atemp & \texttt{0.0} & set to 273.16 to convert from deg. Kelvin (assumed input) to Celsius \\
125            windstressmax & \texttt{2.0} &
126               max. allowed wind stress $N/m^2$ \\
127            exf\_albedo & \texttt{0.1} &
128              surface albedo used to compute downward vs. net radiative fluxes \\
129            exf\_iprec  & \texttt{32} &
130              precision of input fields (32-bit or 64-bit) \\
131            exf\_yftype & \texttt{'RL'} &
132              precision of arrays ('RL' vs. 'RS') \\
133        \hline        \hline
134      \end{tabular}      \end{tabular}
135    }    }
# Line 118  Run-time parameters may be broken into 2 Line 137  Run-time parameters may be broken into 2
137  \end{table}  \end{table}
138    
139    
140  \paragraph{Field attributes}  \paragraph{Field attributes} ~ \\
141    %
142    All EXF fields are listed in Section \ref{sec:pkg:exf:fields_units}.
143    Each field has a number of attributes which can be customized.
144    They are summarized in
145    Table \ref{tab:pkg:exf:runtime_attributes}.
146    To obtain an attribute for a specific field, e.g. \texttt{uwind}
147    prepend the field name to the listed attribute, e.g. for attribute
148    \texttt{period} this yields \texttt{uwindperiod}:
149    %
150    \begin{eqnarray*}
151      \begin{array}{cccccc}
152        ~ & \texttt{field} & \& & \texttt{attribute} & \longrightarrow & \texttt{parameter} \\
153        \text{e.g.} & \text{uwind} & \& & \text{period} & \longrightarrow & \text{uwindperiod} \\
154      \end{array}
155    \end{eqnarray*}
156    %
157    
158    \begin{table}[h!]
159      \label{tab:pkg:exf:runtime_attributes}
160      {\footnotesize
161        \begin{tabular}{|l|c|l|}
162          \hline
163          \textbf{attribute} &  \textbf{Default} &  \textbf{Description}  \\
164          \hline \hline
165             \textit{field}\texttt{file} & ' ' &
166               filename; if left empty no file will be read; \texttt{const} will be used instead \\
167             \textit{field}\texttt{const} & 0. &
168               constant that will be used if no file is read  \\
169             \textit{field}\texttt{startdate1} & 0. &
170               format: \texttt{YYYYMMDD}; start year (YYYY), month (MM), day (YY) \\
171               ~&~& of field to determine record number \\
172             \textit{field}\texttt{startdate2} & 0. &
173               format: \texttt{HHMMSS}; start hour (HH), minute (MM), second(SS) \\
174               ~&~& of field to determine record number\\
175             \textit{field}\texttt{period} & 0. &
176               interval in seconds between two records \\
177             \texttt{exf\_inscal\_}\textit{field}& ~ &
178               optional rescaling of input fields to comply with EXF units \\
179             \texttt{exf\_outscal\_}\textit{field}& ~ &
180               optional rescaling of EXF fields when mapped onto MITgcm fields \\
181             \hline
182             \multicolumn{3}{|c|}{\textit{used in conjunction with}
183                                  \texttt{EXF\_USE\_INTERPOLATION}} \\
184             \hline
185             \textit{field}\texttt{\_lon0} & $thetaMin+delX/2$  &
186               starting longitude of input \\
187             \textit{field}\texttt{\_lon\_inc} & $delX$ &
188               increment in longitude of input \\
189             \textit{field}\texttt{\_lat0} &  $phiMin+delY/2$ &
190               starting latitude of input \\
191             \textit{field}\texttt{\_lat\_inc} & $delY$ &
192               increment in latitude of input \\
193             \textit{field}\texttt{\_nlon} & $Nx$ &
194               number of grid points in longitude of input \\
195             \textit{field}\texttt{\_nlat} & $Ny$ &
196               number of grid points in longitude of input \\
197          \hline
198        \end{tabular}
199       }
200       \caption{\newline
201                Note one exception for the default of
202                \texttt{atempconst} = celsius2K = 273.16}
203    \end{table}
204    
205    \paragraph{Example configuration} ~ \\
206    %
207    The following block is taken from the \texttt{data.exf} file
208    of the veification experiment \texttt{global\_with\_exf/}.
209    It defines attributes for the heat flux variable \texttt{hflux}:
210    
211    \begin{verbatim}
212     hfluxfile       = 'ncep_qnet.bin',
213     hfluxstartdate1 = 19920101,
214     hfluxstartdate2 = 000000,
215     hfluxperiod     = 2592000.0,
216     hflux_lon0      = 2
217     hflux_lon_inc   = 4
218     hflux_lat0      = -78
219     hflux_lat_inc   = 39*4
220     hflux_nlon      = 90
221     hflux_nlat      = 40
222    \end{verbatim}
223    
224    EXF will read a file of name 'ncep\_qnet.bin'.
225    Its first record represents January 1st, 1991 at 00:00 UTC.
226    Next record is 2592000 seconds (or 30 days) later.
227    Interpolation on-the-fly is used (in the present case trivially
228    on the same grid, but included nevertheless for illustration),
229    and input field grid starting coordinates and increments are
230    supplied as well.
231    
232  %----------------------------------------------------------------------  %----------------------------------------------------------------------
233    
234  \subsection{EXF fields and units  \subsection{EXF input fields and units
235  \label{sec:pkg:exf:fields_units}}  \label{sec:pkg:exf:fields_units}}
236    
237  The following list is taken from the header file \texttt{exf\_fields.h}.  The following list is taken from the header file \texttt{exf\_fields.h}.
238    It comprises all EXF input fields.
239    
240  {\footnotesize  Output fields which EXF provides to the MITgcm are fields
241  \begin{verbatim}  \textbf{fu}, \textbf{fv}, \textbf{Qnet}, \textbf{Qsw}, \textbf{EmPmR},
242    and \textbf{pload}. They are defined in \texttt{FFIELDS.h}.
243    
244    {\scriptsize
245    \begin{verbatim}
246    
247    c----------------------------------------------------------------------
248    c               |
249    c     field     :: Description
250    c               |
251    c----------------------------------------------------------------------
252  c     ustress   :: Zonal surface wind stress in N/m^2  c     ustress   :: Zonal surface wind stress in N/m^2
253  c                  > 0 for increase in uVel, which is west to  c               |  > 0 for increase in uVel, which is west to
254  c                      east for cartesian and spherical polar grids  c               |      east for cartesian and spherical polar grids
255  c                  Typical range: -0.5 < ustress < 0.5  c               |  Typical range: -0.5 < ustress < 0.5
256  c                  Southwest C-grid U point  c               |  Southwest C-grid U point
257  c                  Input field  c               |  Input field
258  c  c----------------------------------------------------------------------
259  c     vstress   :: Meridional surface wind stress in N/m^2  c     vstress   :: Meridional surface wind stress in N/m^2
260  c                  > 0 for increase in vVel, which is south to  c               |  > 0 for increase in vVel, which is south to
261  c                      north for cartesian and spherical polar grids  c               |      north for cartesian and spherical polar grids
262  c                  Typical range: -0.5 < vstress < 0.5  c               |  Typical range: -0.5 < vstress < 0.5
263  c                  Southwest C-grid V point  c               |  Southwest C-grid V point
264  c                  Input field  c               |  Input field
265  c  c----------------------------------------------------------------------
266  c     hflux     :: Net upward surface heat flux in W/m^2  c     hflux     :: Net upward surface heat flux in W/m^2
267  c                  excluding shortwave (on input)  c               |  excluding shortwave (on input)
268  c                  hflux = latent + sensible + lwflux  c               |  hflux = latent + sensible + lwflux
269  c                  > 0 for decrease in theta (ocean cooling)  c               |  > 0 for decrease in theta (ocean cooling)
270  c                  Typical range: -250 < hflux < 600  c               |  Typical range: -250 < hflux < 600
271  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
272  c                  Input field  c               |  Input field
273  c  c----------------------------------------------------------------------
274  c     sflux     :: Net upward freshwater flux in m/s  c     sflux     :: Net upward freshwater flux in m/s
275  c                  sflux = evap - precip - runoff  c               |  sflux = evap - precip - runoff
276  c                  > 0 for increase in salt (ocean salinity)  c               |  > 0 for increase in salt (ocean salinity)
277  c                  Typical range: -1e-7 < sflux < 1e-7  c               |  Typical range: -1e-7 < sflux < 1e-7
278  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
279  c                  Input field  c               |  Input field
280  c  c----------------------------------------------------------------------
281  c     swflux    :: Net upward shortwave radiation in W/m^2  c     swflux    :: Net upward shortwave radiation in W/m^2
282  c                  swflux = - ( swdown - ice and snow absorption - reflected )  c               |  swflux = - ( swdown - ice and snow absorption - reflected )
283  c                  > 0 for decrease in theta (ocean cooling)  c               |  > 0 for decrease in theta (ocean cooling)
284  c                  Typical range: -350 < swflux < 0  c               |  Typical range: -350 < swflux < 0
285  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
286  c                  Input field  c               |  Input field
287  c  c----------------------------------------------------------------------
288  c     uwind     :: Surface (10-m) zonal wind velocity in m/s  c     uwind     :: Surface (10-m) zonal wind velocity in m/s
289  c                  > 0 for increase in uVel, which is west to  c               |  > 0 for increase in uVel, which is west to
290  c                      east for cartesian and spherical polar grids  c               |      east for cartesian and spherical polar grids
291  c                  Typical range: -10 < uwind < 10  c               |  Typical range: -10 < uwind < 10
292  c                  Southwest C-grid U point  c               |  Southwest C-grid U point
293  c                  Input or input/output field  c               |  Input or input/output field
294  c  c----------------------------------------------------------------------
295  c     vwind     :: Surface (10-m) meridional wind velocity in m/s  c     vwind     :: Surface (10-m) meridional wind velocity in m/s
296  c                  > 0 for increase in vVel, which is south to  c               |  > 0 for increase in vVel, which is south to
297  c                      north for cartesian and spherical polar grids  c               |      north for cartesian and spherical polar grids
298  c                  Typical range: -10 < vwind < 10  c               |  Typical range: -10 < vwind < 10
299  c                  Southwest C-grid V point  c               |  Southwest C-grid V point
300  c                  Input or input/output field  c               |  Input or input/output field
301  c  c----------------------------------------------------------------------
302  c     atemp     :: Surface (2-m) air temperature in deg K  c     atemp     :: Surface (2-m) air temperature in deg K
303  c                  Typical range: 200 < atemp < 300  c               |  Typical range: 200 < atemp < 300
304  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
305  c                  Input or input/output field  c               |  Input or input/output field
306  c  c----------------------------------------------------------------------
307  c     aqh       :: Surface (2m) specific humidity in kg/kg  c     aqh       :: Surface (2m) specific humidity in kg/kg
308  c                  Typical range: 0 < aqh < 0.02  c               |  Typical range: 0 < aqh < 0.02
309  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
310  c                  Input or input/output field  c               |  Input or input/output field
311  c  c----------------------------------------------------------------------
312  c     lwflux    :: Net upward longwave radiation in W/m^2  c     lwflux    :: Net upward longwave radiation in W/m^2
313  c                  lwflux = - ( lwdown - ice and snow absorption - emitted )  c               |  lwflux = - ( lwdown - ice and snow absorption - emitted )
314  c                  > 0 for decrease in theta (ocean cooling)  c               |  > 0 for decrease in theta (ocean cooling)
315  c                  Typical range: -20 < lwflux < 170  c               |  Typical range: -20 < lwflux < 170
316  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
317  c                  Input field  c               |  Input field
318  c  c----------------------------------------------------------------------
319  c     evap      :: Evaporation in m/s  c     evap      :: Evaporation in m/s
320  c                  > 0 for increase in salt (ocean salinity)  c               |  > 0 for increase in salt (ocean salinity)
321  c                  Typical range: 0 < evap < 2.5e-7  c               |  Typical range: 0 < evap < 2.5e-7
322  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
323  c                  Input, input/output, or output field  c               |  Input, input/output, or output field
324  c  c----------------------------------------------------------------------
325  c     precip    :: Precipitation in m/s  c     precip    :: Precipitation in m/s
326  c                  > 0 for decrease in salt (ocean salinity)  c               |  > 0 for decrease in salt (ocean salinity)
327  c                  Typical range: 0 < precip < 5e-7  c               |  Typical range: 0 < precip < 5e-7
328  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
329  c                  Input or input/output field  c               |  Input or input/output field
330  c  c----------------------------------------------------------------------
331  c     runoff    :: River and glacier runoff in m/s  c     runoff    :: River and glacier runoff in m/s
332  c                  > 0 for decrease in salt (ocean salinity)  c               |  > 0 for decrease in salt (ocean salinity)
333  c                  Typical range: 0 < runoff < ????  c               |  Typical range: 0 < runoff < ????
334  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
335  c                  Input or input/output field  c               |  Input or input/output field
336  c                  !!! WATCH OUT: Default exf_inscal_runoff !!!  c               |  !!! WATCH OUT: Default exf_inscal_runoff !!!
337  c                  !!! in exf_readparms.F is not 1.0        !!!  c               |  !!! in exf_readparms.F is not 1.0        !!!
338  c  c----------------------------------------------------------------------
339  c     swdown    :: Downward shortwave radiation in W/m^2  c     swdown    :: Downward shortwave radiation in W/m^2
340  c                  > 0 for increase in theta (ocean warming)  c               |  > 0 for increase in theta (ocean warming)
341  c                  Typical range: 0 < swdown < 450  c               |  Typical range: 0 < swdown < 450
342  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
343  c                  Input/output field  c               |  Input/output field
344  c  c----------------------------------------------------------------------
345  c     lwdown    :: Downward longwave radiation in W/m^2  c     lwdown    :: Downward longwave radiation in W/m^2
346  c                  > 0 for increase in theta (ocean warming)  c               |  > 0 for increase in theta (ocean warming)
347  c                  Typical range: 50 < lwdown < 450  c               |  Typical range: 50 < lwdown < 450
348  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
349  c                  Input/output field  c               |  Input/output field
350  c  c----------------------------------------------------------------------
351  c     apressure :: Atmospheric pressure field in N/m^2  c     apressure :: Atmospheric pressure field in N/m^2
352  c                  > 0 for ????  c               |  > 0 for ????
353  c                  Typical range: ???? < apressure < ????  c               |  Typical range: ???? < apressure < ????
354  c                  Southwest C-grid tracer point  c               |  Southwest C-grid tracer point
355  c                  Input field  c               |  Input field
356  C  c----------------------------------------------------------------------
357  C  
 c     NOTES:  
 c     ======  
 c  
 c     Input and output units and sign conventions can be customized  
 c     using variables exf_inscal_* and exf_outscal_*, which are set  
 c     by exf_readparms.F  
 c  
 c     Output fields fu, fv, Qnet, Qsw, and EmPmR are  
 c     defined in FFIELDS.h  
 c  
 c     #ifndef SHORTWAVE_HEATING, hflux includes shortwave,  
 c     that is, hflux = latent + sensible + lwflux +swflux  
 c  
 c     If (EXFwindOnBgrid .EQ. .TRUE.), uwind and vwind are  
 c     defined on northeast B-grid U and V points, respectively.  
 c  
 c     Arrays *0 and *1 below are used for temporal interpolation.  
358  \end{verbatim}  \end{verbatim}
359  }  }
360    
# Line 265  c     Arrays *0 and *1 below are used fo Line 363  c     Arrays *0 and *1 below are used fo
363  \subsection{Key subroutines  \subsection{Key subroutines
364  \label{sec:pkg:exf:subroutines}}  \label{sec:pkg:exf:subroutines}}
365    
366    \subsubsection{Top-level routine: \texttt{exf\_getforcing.F}}
367    
368    {\scriptsize
369    \begin{verbatim}
370    
371    C     !CALLING SEQUENCE:
372    c ...
373    c  exf_getforcing (TOP LEVEL ROUTINE)
374    c  |
375    c  |-- exf_getclim (get climatological fields used e.g. for relax.)
376    c  |   |--- exf_set_climtemp (relax. to 3-D temperature field)
377    c  |   |--- exf_set_climsalt (relax. to 3-D salinity field)
378    c  |   |--- exf_set_climsst  (relax. to 2-D SST field)
379    c  |   |--- exf_set_climsss  (relax. to 2-D SSS field)
380    c  |   o
381    c  |
382    c  |-- exf_getffields <- this one does almost everything
383    c  |   |   1. reads in fields, either flux or atmos. state,
384    c  |   |      depending on CPP options (for each variable two fields
385    c  |   |      consecutive in time are read in and interpolated onto
386    c  |   |      current time step).
387    c  |   |   2. If forcing is atmos. state and control is atmos. state,
388    c  |   |      then the control variable anomalies are read here
389    c  |   |          * ctrl_getatemp
390    c  |   |          * ctrl_getaqh
391    c  |   |          * ctrl_getuwind
392    c  |   |          * ctrl_getvwind
393    c  |   |      If forcing and control are fluxes, then
394    c  |   |      controls are added later.
395    c  |   o
396    c  |
397    c  |-- exf_check_range
398    c  |   |   1. Check whether read fields are within assumed range
399    c  |   |      (may capture mismatches in units)
400    c  |   o
401    c  |
402    c  |-- exf_bulkformulae
403    c  |   |   1. Compute net or downwelling radiative fluxes via
404    c  |   |      Stefan-Boltzmann law in case only one is known.
405    c  |   |   2. Compute air-sea momentum and buoyancy fluxes from
406    c  |   |      atmospheric state following Large and Pond, JPO, 1981/82
407    c  |   o
408    c  |
409    c  |-- < add time-mean river runoff here, if available >
410    c  |
411    c  |-- < update tile edges here >
412    c  |
413    c  |-- exf_getsurfacefluxes
414    c  |   |   1. If forcing and control are fluxes, then
415    c  |   |      controls are added here.
416    c  |   o
417    c  |
418    c  |-- < treatment of hflux w.r.t. swflux >
419    c  |
420    c  |-- exf_diagnostics_fill
421    c  |   |   1. Do EXF-related diagnostics output here.
422    c  |   o
423    c  |
424    c  |-- exf_mapfields
425    c  |   |   1. Map the EXF variables onto the core MITgcm
426    c  |   |      forcing fields.
427    c  |   o
428    c  |
429    c  |-- exf_bulkformulae
430    c  |   If ALLOW_BULKFORMULAE, compute fluxes via bulkformulae
431    c  |
432    c  |-- exf_getsurfacefluxes
433    c  |   If forcing and control is flux, then the
434    c  |   control vector anomalies are read here
435    c  |      * ctrl_getheatflux
436    c  |      * ctrl_getsaltflux
437    c  |      * ctrl_getzonstress
438    c  |      * call ctrl_getmerstress
439    c  |
440    c  |-- exf_mapfields
441    c  |   Forcing fields from exf package are mapped onto
442    c  |   mitgcm forcing arrays.
443    c  |   Mapping enables a runtime rescaling of fields
444    
445    \end{verbatim}
446    }
447    
448    \subsubsection{Bulk formula routine: \texttt{exf\_bulkformulae.F}}
449    
450    \subsubsection{Generic I/O routine: \texttt{exf\_set\_gen.F}}
451    
452    \subsubsection{Interpolation routine: \texttt{exf\_interp.F}}
453    
454    \subsubsection{Header routines}
455    
456  %----------------------------------------------------------------------  %----------------------------------------------------------------------
457    
458  \subsection{EXF diagnostics  \subsection{EXF diagnostics
# Line 275  Diagnostics output is available via the Line 463  Diagnostics output is available via the
463  Available output fields are summarized in  Available output fields are summarized in
464  Table \ref{tab:pkg:exf:diagnostics}.  Table \ref{tab:pkg:exf:diagnostics}.
465    
466  \begin{table}  \begin{table}[h!]
467  \label{tab:pkg:exf:diagnostics}  \label{tab:pkg:exf:diagnostics}
 \caption{~}  
468  {\footnotesize  {\footnotesize
469  \begin{verbatim}  \begin{verbatim}
470  ------------------------------------------------------  ------------------------------------------------------
471   <-Name->|Levs|grid|<--  Units   -->|<- Tile (max=80c)   <-Name->|Levs|grid|<--  Units   -->|<- Tile (max=80c)
472  ------------------------------------------------------  ------------------------------------------------------
473   EXFlwdn |  1 | SM |W/m^2           |Downward longwave radiation, >0 increases theta   EXFlwdn |  1 | SM | W/m^2          | Downward longwave radiation, >0 increases theta
474   EXFswdn |  1 | SM |W/m^2           |Downward shortwave radiation, >0 increases theta   EXFswdn |  1 | SM | W/m^2          | Downward shortwave radiation, >0 increases theta
475   EXFqnet |  1 | SM |W/m^2           |Net upward heat flux (turb+rad), >0 decreases theta   EXFqnet |  1 | SM | W/m^2          | Net upward heat flux (turb+rad), >0 decreases theta
476   EXFtaux |  1 | SU |N/m^2           |zonal surface wind stress, >0 increases uVel   EXFtaux |  1 | SU | N/m^2          | zonal surface wind stress, >0 increases uVel
477   EXFtauy |  1 | SV |N/m^2           |meridional surface wind stress, >0 increases vVel   EXFtauy |  1 | SV | N/m^2          | meridional surface wind stress, >0 increases vVel
478   EXFuwind|  1 | SM |m/s             |zonal 10-m wind speed, >0 increases uVel   EXFuwind|  1 | SM | m/s            | zonal 10-m wind speed, >0 increases uVel
479   EXFvwind|  1 | SM |m/s             |meridional 10-m wind speed, >0 increases uVel   EXFvwind|  1 | SM | m/s            | meridional 10-m wind speed, >0 increases uVel
480   EXFatemp|  1 | SM |degK            |surface (2-m) air temperature   EXFatemp|  1 | SM | degK           | surface (2-m) air temperature
481   EXFaqh  |  1 | SM |kg/kg           |surface (2-m) specific humidity   EXFaqh  |  1 | SM | kg/kg          | surface (2-m) specific humidity
482   EXFevap |  1 | SM |m/s             |evaporation, > 0 increases salinity   EXFevap |  1 | SM | m/s            | evaporation, > 0 increases salinity
483   EXFpreci|  1 | SM |m/s             |evaporation, > 0 decreases salinity   EXFpreci|  1 | SM | m/s            | evaporation, > 0 decreases salinity
484   EXFempmr|  1 | SM |m/s             |net upward freshwater flux, > 0 increases salinity   EXFempmr|  1 | SM | m/s            | net upward freshwater flux, > 0 increases salinity
485   EXFpress|  1 | SM |N/m^2           |atmospheric pressure field   EXFpress|  1 | SM | N/m^2          | atmospheric pressure field
486  \end{verbatim}  \end{verbatim}
487  }  }
488    \caption{~}
489  \end{table}  \end{table}
490    
491  %----------------------------------------------------------------------  %----------------------------------------------------------------------
492    
493  \subsection{Reference experiments}  \subsection{Reference experiments}
494    
495    \subsubsection{global\_with\_exf}
496    
497    \subsubsection{lab\_sea}
498    
499  %----------------------------------------------------------------------  %----------------------------------------------------------------------
500    
501  \subsection{References}  \subsection{References}

Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3

  ViewVC Help
Powered by ViewVC 1.1.22