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

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

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

-revision 1.9 by molod,
Wed Jun 28 15:35:07 2006 UTC
+revision 1.16 by jmc,
Fri Dec  9 21:41:30 2016 UTC
 Line 22 
 forcing fields an arbitrary grid onto th
  CPP options enable or disable different aspects of the package
  (Section \ref{sec:pkg:exf:config}).
  Runtime options, flags, filenames and field-related dates/times are
- set in \texttt{data.exf} and \texttt{data.exf\_clim}
+ set in \texttt{data.exf}
  (Section \ref{sec:pkg:exf:runtime}).
  A description of key subroutines is given in Section
  \ref{sec:pkg:exf:subroutines}.
  Input fields, units and sign conventions are summarized in
  Section \ref{sec:pkg:exf:fields_units}, and available diagnostics
- output is listed in Section \ref{sec:pkg:exf:fields_diagnostics}.
+ output is listed in Section \ref{sec:pkg:exf:diagnostics}.
  %----------------------------------------------------------------------
 Line 55 
 EXF requires the calendar package \textt
  no additional CPP options are required.
  %
  \end{itemize}
- (see Section \ref{sect:buildingCode}).
+ (see Section \ref{sec:buildingCode}).
  Parts of the EXF code can be enabled or disabled at compile time
  via CPP preprocessor flags. These options are set in either
 Line 63 
 via CPP preprocessor flags. These option
  Table \ref{tab:pkg:exf:cpp} summarizes these options.
  \begin{table}[b!]
- \centering
+   \centering
-   \label{tab:pkg:exf:cpp}
    {\footnotesize
      \begin{tabular}{|l|l|}
        \hline
        \textbf{CPP option}  &  \textbf{Description}  \\
        \hline \hline
          \texttt{EXF\_VERBOSE} &
            verbose mode (recommended only for testing) \\
          \texttt{ALLOW\_ATM\_TEMP} &
            compute heat/freshwater fluxes from atmos. state input \\
          \texttt{ALLOW\_ATM\_WIND} &
            compute wind stress from wind speed input\\
          \texttt{ALLOW\_BULKFORMULAE} &
            is used if \texttt{ALLOW\_ATM\_TEMP} or
            \texttt{ALLOW\_ATM\_WIND} is enabled \\
          \texttt{EXF\_READ\_EVAP} & read evaporation instead of computing it \\
          \texttt{ALLOW\_RUNOFF} & read time-constant river/glacier run-off field \\
-Line 99 
 Table \ref{tab:pkg:exf:cpp} summarizes t
+Line 98 
 Table \ref{tab:pkg:exf:cpp} summarizes t
      \end{tabular}
    }
    \caption{~}
+   \label{tab:pkg:exf:cpp}
  \end{table}
  %----------------------------------------------------------------------
  \subsubsection{Run-time parameters
  \label{sec:pkg:exf:runtime}}
  Run-time parameters are set in files
- \texttt{data.pkg}, \texttt{data.exf}, and
+ \texttt{data.pkg} and \texttt{data.exf}
- \texttt{data.exf\_clim} (for relaxation/climatological fields)
+ %\texttt{data.exf\_clim} (for relaxation/climatological fields)
- which are read in \texttt{exf\_readparms.F}.
+ which is read in \texttt{exf\_readparms.F}.
  Run-time parameters may be broken into 3 categories:
  (i) switching on/off the package at runtime,
  (ii) general flags and parameters, and
-Line 119 
 Run-time parameters may be broken into 3
+Line 118 
 Run-time parameters may be broken into 3
  \paragraph{Enabling the package}
  ~ \\
  %
- A package is usually switched on/off at runtime by setting
+ A package is switched on/off at runtime by setting
  (e.g. for EXF) \texttt{useEXF = .TRUE.} in \texttt{data.pkg}.
- For EXF this flag is omitted, i.e. EXF is always ON if it is compiled.
  \paragraph{General flags and parameters}
  ~ \\
  %
- \begin{table}[h!]
+ \begin{table}[!ht]
- \centering
+   \centering
-   \label{tab:pkg:exf:runtime_flags}
    {\footnotesize
      \begin{tabular}{|l|c|l|}
        \hline
        \textbf{Flag/parameter} & \textbf{default} &  \textbf{Description}  \\
        \hline \hline
          useExfCheckRange & \texttt{.TRUE.} &
             check range of input fields and stop if out of range \\
          useExfYearlyFields & \texttt{.FALSE.} &
             append current year postfix of form \texttt{\_YYYY} on filename \\
          twoDigitYear & \texttt{.FALSE.} &
             instead of appending \texttt{\_YYYY} append  \texttt{YY} \\
          repeatPeriod & \texttt{0.0} & $ > 0 $ :
             cycle through all input fields at the same period (in seconds) \\
          ~            & ~            & $ = 0 $ :
             use period assigned to each field \\
          exf\_offset\_atemp & \texttt{0.0} & set to 273.16 to convert from deg. Kelvin (assumed input) to Celsius \\
          windstressmax & \texttt{2.0} &
             max. allowed wind stress $N/m^2$ \\
          exf\_albedo & \texttt{0.1} &
            surface albedo used to compute downward vs. net radiative fluxes \\
-         exf\_iprec  & \texttt{32} &
+         climtempfreeze & \texttt{-1.9} &
+           ??? \\
+         ocean\_emissivity & \texttt{} &
+           longwave ocean-surface emissivity \\
+         ice\_emissivity & \texttt{} &
+           longwave seaice emissivity \\
+         snow\_emissivity & \texttt{} &
+           longwave  snow  emissivity \\
+         exf\_iceCd & \texttt{1.63E-3} &
+           drag coefficient over sea-ice \\
+         exf\_iceCe & \texttt{1.63E-3} &
+           evaporation transfer coeff. over sea-ice \\
+         exf\_iceCh & \texttt{1.63E-3} &
+           sensible heat transfer coeff. over sea-ice \\
+         exf\_scal\_BulkCdn & \texttt{1.} &
+           overall scaling of neutral drag coeff. \\
+         useStabilityFct\_overIce  & \texttt{.FALSE.} &
+           compute turbulent transfer coeff. over sea-ice \\
+         readStressOnAgrid & \texttt{.FALSE.} &
+           read wind-streess located on model-grid, A-grid point \\
+         readStressOnCgrid & \texttt{.FALSE.} &
+           read wind-streess located on model-grid, C-grid point \\
+         useRelativeWind & \texttt{.FALSE.} &
+           subtract [U/V]VEL or [U/VICE from U/V]WIND before \\
+         ~ & ~ &  computing [U/V]STRESS \\
+         zref & \texttt{10.} &
+           reference height \\
+         hu & \texttt{10.} &
+           height of mean wind \\
+         ht & \texttt{2.} &
+           height of mean temperature and rel. humidity \\
+         umin & \texttt{0.5} &
+           minimum absolute wind speed for computing Cd \\
+         atmrho & \texttt{1.2} &
+           mean atmospheric density [kg/m\^3] \\
+         atmcp & \texttt{1005.} &
+           mean atmospheric specific heat [J/kg/K] \\
+         cdrag\_[n] & \texttt{???} &
+           n = 1,2,3; parameters for drag coeff. function \\
+         cstanton\_[n] & \texttt{???} &
+           n = 1,2; parameters for Stanton number function \\
+         cdalton & \texttt{???} &
+           parameter for Dalton number function \\
+         flamb & \texttt{2500000.} &
+           latent heat of evaporation [J/kg] \\
+         flami & \texttt{334000.} &
+           latent heat of melting of pure ice [J/kg] \\
+         zolmin & \texttt{-100.} &
+           minimum stability parameter \\
+         cvapor\_fac & \texttt{640380.} &
+           ~ \\
+         cvapor\_exp & \texttt{5107.4} &
+           ~ \\
+         cvapor\_fac\_ice & \texttt{11637800.} &
+           ~ \\
+         cvapor\_fac\_ice & \texttt{5897.8} &
+           ~ \\
+         humid\_fac & \texttt{0.606} &
+           parameter for virtual temperature calculation \\
+         gamma\_blk & \texttt{0.010} &
+           adiabatic lapse rate \\
+         saltsat & \texttt{0.980} &
+           reduction of saturation vapor pressure over salt-water \\
+         psim\_fac & \texttt{5.} &
+           ~ \\
+         exf\_monFreq & \texttt{monitorFreq} &
+           output frequency [s] \\
+         exf\_iprec  & \texttt{32} &
            precision of input fields (32-bit or 64-bit) \\
          exf\_yftype & \texttt{'RL'} &
            precision of arrays ('RL' vs. 'RS') \\
        \hline
      \end{tabular}
    }
    \caption{~}
+   \label{tab:pkg:exf:runtime_flags}
  \end{table}
+ \paragraph{Field attributes}
- \paragraph{Field attributes}
  ~ \\
  %
  All EXF fields are listed in Section \ref{sec:pkg:exf:fields_units}.
-Line 179 
 prepend the field name to the listed att
+Line 243 
 prepend the field name to the listed att
  \end{eqnarray*}
  %
- \begin{table}[h!]
+ \begin{table}[!ht]
- \centering
+   \centering
-   \label{tab:pkg:exf:runtime_attributes}
    {\footnotesize
      \begin{tabular}{|l|c|l|}
        \hline
        \textbf{attribute} &  \textbf{Default} &  \textbf{Description}  \\
        \hline \hline
           \textit{field}\texttt{file} & ' ' &
             filename; if left empty no file will be read; \texttt{const} will be used instead \\
           \textit{field}\texttt{const} & 0. &
             constant that will be used if no file is read  \\
           \textit{field}\texttt{startdate1} & 0. &
             format: \texttt{YYYYMMDD}; start year (YYYY), month (MM), day (YY) \\
             ~&~& of field to determine record number \\
           \textit{field}\texttt{startdate2} & 0. &
-Line 199 
 prepend the field name to the listed att
+Line 262 
 prepend the field name to the listed att
             ~&~& of field to determine record number\\
           \textit{field}\texttt{period} & 0. &
             interval in seconds between two records \\
           \texttt{exf\_inscal\_}\textit{field}& ~ &
             optional rescaling of input fields to comply with EXF units \\
           \texttt{exf\_outscal\_}\textit{field}& ~ &
             optional rescaling of EXF fields when mapped onto MITgcm fields \\
           \hline
           \multicolumn{3}{|c|}{\textit{used in conjunction with}
                                \texttt{EXF\_USE\_INTERPOLATION}} \\
           \hline
-          \textit{field}\texttt{\_lon0} & $thetaMin+delX/2$  &
+          \textit{field}\texttt{\_lon0} & $xgOrigin+delX/2$  &
             starting longitude of input \\
           \textit{field}\texttt{\_lon\_inc} & $delX$ &
             increment in longitude of input \\
-          \textit{field}\texttt{\_lat0} &  $phiMin+delY/2$ &
+          \textit{field}\texttt{\_lat0} &  $ygOrigin+delY/2$ &
             starting latitude of input \\
           \textit{field}\texttt{\_lat\_inc} & $delY$ &
             increment in latitude of input \\
-Line 222 
 prepend the field name to the listed att
+Line 285 
 prepend the field name to the listed att
        \hline
      \end{tabular}
     }
-    \caption{\newline
+   \caption{\newline
-             Note one exception for the default of
+            Note one exception for the default of
-             \texttt{atempconst} = celsius2K = 273.16}
+            \texttt{atempconst} = celsius2K = 273.16}
+   \label{tab:pkg:exf:runtime_attributes}
  \end{table}
  \paragraph{Example configuration} ~ \\
  %
  The following block is taken from the \texttt{data.exf} file
- of the veification experiment \texttt{global\_with\_exf/}.
+ of the verification experiment \texttt{global\_with\_exf/}.
  It defines attributes for the heat flux variable \texttt{hflux}:
  \begin{verbatim}
-Line 247 
 It defines attributes for the heat flux
+Line 311 
 It defines attributes for the heat flux
  \end{verbatim}
  EXF will read a file of name 'ncep\_qnet.bin'.
- Its first record represents January 1st, 1991 at 00:00 UTC.
+ Its first record represents January 1st, 1992 at 00:00 UTC.
  Next record is 2592000 seconds (or 30 days) later.
+ Note that the first record read and used by the EXF package corresponds to the
+ value 'startDate1' set in data.cal.  Therefore if you want to start the EXF
+ forcing from later in the 'ncep\_qnet.bin' file, it suffices to specify
+ startDate1 in data.cal as a date later than 19920101 (for example, startDate1
+ = 19940101, for starting January 1st, 1994).  For this to work,
+ 'ncep\_qnet.bin' must have at least 2 years of data because in this
+ configuration EXF will read 2 years into the file to find the 1994 starting
+ value.
  Interpolation on-the-fly is used (in the present case trivially
  on the same grid, but included nevertheless for illustration),
  and input field grid starting coordinates and increments are
  supplied as well.
  %----------------------------------------------------------------------
+ \subsubsection{EXF bulk formulae
+ \label{sec:pkg:exf:bulk_formulae}}
+ T.B.D. (cross-ref. to parameter list table)
+ %----------------------------------------------------------------------
  \subsubsection{EXF input fields and units
  \label{sec:pkg:exf:fields_units}}
- 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}.
  It comprises all EXF input fields.
  Output fields which EXF provides to the MITgcm are fields
-Line 294 
 c---------------------------------------
+Line 373 
 c---------------------------------------
  c     hl        :: latent   heat flux into ocean in W/m^2
  c               |  > 0 for increase in theta (ocean warming)
  c----------------------------------------------------------------------
  c     hflux     :: Net upward surface heat flux in W/m^2
- c               |  excluding shortwave (on input)
+ c               |  (including shortwave)
- c               |  hflux = latent + sensible + lwflux
+ c               |  hflux = latent + sensible + lwflux + swflux
  c               |  > 0 for decrease in theta (ocean cooling)
  c               |  Typical range: -250 < hflux < 600
  c               |  Southwest C-grid tracer point
-Line 492 
 Header routines
+Line 571 
 Header routines
  Diagnostics output is available via the diagnostics package
  (see Section \ref{sec:pkg:diagnostics}).
  Available output fields are summarized in
  Table \ref{tab:pkg:exf:diagnostics}.
- \begin{table}[h!]
+ \begin{table}[!ht]
  \centering
- \label{tab:pkg:exf:diagnostics}
  {\footnotesize
  \begin{verbatim}
- ------------------------------------------------------
+ ---------+----+----+----------------+-----------------
   <-Name->|Levs|grid|<--  Units   -->|<- Tile (max=80c)
- ------------------------------------------------------
+ ---------+----+----+----------------+-----------------
   EXFhs   |  1 | SM | W/m^2          | Sensible heat flux into ocean, >0 increases theta
   EXFhl   |  1 | SM | W/m^2          | Latent heat flux into ocean, >0 increases theta
   EXFlwnet|  1 | SM | W/m^2          | Net upward longwave radiation, >0 decreases theta
-Line 525 
 Table \ref{tab:pkg:exf:diagnostics}.
+Line 603 
 Table \ref{tab:pkg:exf:diagnostics}.
  \end{verbatim}
  }
  \caption{~}
+ \label{tab:pkg:exf:diagnostics}
  \end{table}
  %----------------------------------------------------------------------

 Legend:



Removed from v.1.9
 


changed lines


 
Added in v.1.16
 Legend:



Removed from v.1.9
 


changed lines


 
Added in v.1.16
-Removed from v.1.9
+Added in v.1.16

	ViewVC Help
Powered by ViewVC 1.1.22