/[MITgcm]/MITgcm_contrib/MPMice/beaufort/code/EXF_OPTIONS.h

Diff of /MITgcm_contrib/MPMice/beaufort/code/EXF_OPTIONS.h

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

-revision 1.2 by dimitri,
Wed Dec 21 23:06:07 2011 UTC
+revision 1.3 by dimitri,
Sat Oct  4 03:24:19 2014 UTC
 Line 1
  C $Header$
  C $Name$
+ CBOP
+ C !ROUTINE: EXF_OPTIONS.h
+ C !INTERFACE:
+ C #include "EXF_OPTIONS.h"
+ C !DESCRIPTION:
+ C *==================================================================*
+ C | CPP options file for EXternal Forcing (EXF) package:
+ C | Control which optional features to compile in this package code.
+ C *==================================================================*
+ CEOP
  #ifndef EXF_OPTIONS_H
  #define EXF_OPTIONS_H
  #include "PACKAGES_CONFIG.h"
- #include "AD_CONFIG.h"
+ #include "CPP_OPTIONS.h"
  #ifdef ALLOW_EXF
+ #ifdef ECCO_CPPOPTIONS_H
- #include "CPP_OPTIONS.h"
+ C-- When multi-package option-file ECCO_CPPOPTIONS.h is used (directly included
+ C    in CPP_OPTIONS.h), this option file is left empty since all options that
+ C   are specific to this package are assumed to be set in ECCO_CPPOPTIONS.h
+ #else /* ndef ECCO_CPPOPTIONS_H */
- #ifdef ALLOW_AUTODIFF_TAMC
+ C-- Package-specific Options & Macros go here
- C When compile for AD mode,
+ C   pkg/exf CPP options:
- C this package is currently configured by ECCO_CPPOPTIONS.h
+ C   --------------------
- C which is directly included in CPP_OPTIONS.h
+ C
+ C   > ( EXF_VERBOSE ) < replaced with run-time, logical parameter "exf_verbose".
- #else /* ndef ALLOW_AUTODIFF_TAMC */
+ C
+ C   >>> ALLOW_ATM_WIND <<<
- C CPP flags controlling which code is included in the files that
+ C       If defined, 10-m wind fields can be read-in from files.
- C will be compiled.
+ C
- C
+ C   >>> ALLOW_ATM_TEMP <<<
+ C       If defined, atmospheric temperature and specific
- c   pkg/exf CPP options:
+ C       humidity fields can be read-in from files.
- c   --------------------
+ C
- c
+ C   >>> ALLOW_DOWNWARD_RADIATION <<<
- c   > ( EXF_VERBOSE ) < replaced with run-time, logical parameter "exf_verbose".
+ C       If defined, downward long-wave and short-wave radiation
- c
+ C       can be read-in form files or computed from lwflux and swflux.
- c   >>> ALLOW_ATM_WIND <<<
+ C
- c       If defined, 10-m wind fields can be read-in from files.
+ C   >>> ALLOW_ZENITHANGLE <<<
- c
+ C       If defined, ocean albedo varies with the zenith angle, and
- c   >>> ALLOW_ATM_TEMP <<<
+ C       incoming fluxes at the top of the atmosphere are computed
- c       If defined, atmospheric temperature and specific
+ C
- c       humidity fields can be read-in from files.
+ C   >>> ALLOW_BULKFORMULAE <<<
- c
+ C       Allows the use of bulk formulae in order to estimate
- c   >>> ALLOW_DOWNWARD_RADIATION <<<
+ C       turbulent and radiative fluxes at the ocean surface.
- c       If defined, downward long-wave and short-wave radiation
+ C
- c       can be read-in form files or computed from lwflux and swflux.
+ C   >>> EXF_READ_EVAP <<<
- c
+ C       If defined, evaporation fields are read-in, rather than
- c   >>> ALLOW_ZENITHANGLE <<<
+ C       computed from atmospheric state.
- c       If defined, ocean albedo varies with the zenith angle, and
+ C
- c       incoming fluxes at the top of the atmosphere are computed
+ C   >>> ALLOW_RUNOFF <<<
- c
+ C       If defined, river and glacier runoff can be read-in from files.
- c   >>> ALLOW_BULKFORMULAE <<<
+ C
- c       Allows the use of bulk formulae in order to estimate
+ C   >>> ALLOW_RUNOFTEMP <<<
- c       turbulent and radiative fluxes at the ocean surface.
+ C       If defined, river and glacier runoff temperature
- c
+ C       can be read-in from files.
- c   >>> EXF_READ_EVAP <<<
+ C
- c       If defined, evaporation fields are read-in, rather than
+ C   >>> ATMOSPHERIC_LOADING <<<
- c       computed from atmospheric state.
+ C       If defined, atmospheric pressure can be read-in from files.
- c
+ C   WARNING: this flag is set (define/undef) in CPP_OPTIONS.h
- c   >>> ALLOW_RUNOFF <<<
+ C            and cannot be changed here (in EXF_OPTIONS)
- c       If defined, river and glacier runoff can be read-in from files.
+ C
- c
+ C   >>> EXF_SEAICE_FRACTION <<<
- c   >>> ATMOSPHERIC_LOADING <<<
+ C       If defined, seaice fraction can be read-in from files (areaMaskFile)
- c       If defined, atmospheric pressure can be read-in from files.
+ C
- c   WARNING: this flag is set (define/undef) in CPP_OPTIONS.h
+ C   >>> ALLOW_CLIMSST_RELAXATION <<<
- c            and cannot be changed here (in EXF_OPTIONS)
+ C       Allow the relaxation to a monthly climatology of sea surface
- c
+ C       temperature, e.g. the Reynolds climatology.
- c   >>> ICE_AREAMASK <<<
+ C
- c       If defined, fractional ice-covered area MASK can be read-in from files.
+ C   >>> ALLOW_CLIMSSS_RELAXATION <<<
- c
+ C       Allow the relaxation to a monthly climatology of sea surface
- c   >>> ALLOW_CLIMSST_RELAXATION <<<
+ C       salinity, e.g. the Levitus climatology.
- c       Allow the relaxation to a monthly climatology of sea surface
+ C
- c       temperature, e.g. the Reynolds climatology.
+ C   >>> USE_EXF_INTERPOLATION <<<
- c
+ C       Allows specification of arbitrary Cartesian input grids.
- c   >>> ALLOW_CLIMSSS_RELAXATION <<<
+ C
- c       Allow the relaxation to a monthly climatology of sea surface
+ C   >>> EXF_CALC_ATMRHO
- c       salinity, e.g. the Levitus climatology.
+ C       Calculate the local atm density as function of temp, humidity
- c
+ C       and pressure
- c   >>> USE_EXF_INTERPOLATION <<<
+ C
- c       Allows specification of arbitrary Cartesian input grids.
+ C   ====================================================================
- c
+ C
- c   ====================================================================
+ C       The following CPP options:
- c
+ C
- c       The following CPP options:
+ C          ALLOW_ATM_WIND              (WIND)
- c
+ C          ALLOW_ATM_TEMP              (TEMP)
- c          ALLOW_ATM_WIND              (WIND)
+ C          ALLOW_DOWNWARD_RADIATION    (DOWN)
- c          ALLOW_ATM_TEMP              (TEMP)
+ C          ALLOW_BULKFORMULAE          (BULK)
- c          ALLOW_DOWNWARD_RADIATION    (DOWN)
+ C          EXF_READ_EVAP               (EVAP)
- c          ALLOW_BULKFORMULAE          (BULK)
+ C
- c          EXF_READ_EVAP               (EVAP)
+ C       permit the ocean-model forcing configurations listed in the
- c
+ C       table below.  The first configuration is the default,
- c       permit the ocean-model forcing configurations listed in the
+ C       flux-forced, ocean model.  The next four are stand-alone
- c       table below.  The first configuration is the default,
+ C       configurations that use pkg/exf, open-water bulk formulae to
- c       flux-forced, ocean model.  The next four are stand-alone
+ C       compute the missing surface fluxes from atmospheric variables.
- c       configurations that use pkg/exf, open-water bulk formulae to
+ C       The last four configurations can be used in conjunction with
- c       compute the missing surface fluxes from atmospheric variables.
+ C       pkg/seaice to model ice-covered regions.  The forcing fields
- c       The last four configurations can be used in conjunction with
+ C       in the rightmost column are defined in exf_fields.
- c       pkg/seaice to model ice-covered regions.  The forcing fields
+ C
- c       in the rightmost column are defined in exf_fields.
+ C
- c
+ C    WIND |TEMP |DOWN |BULK |EVAP |            actions
- c
+ C    -----|-----|-----|-----|-----|-------------------------------------
- c    WIND |TEMP |DOWN |BULK |EVAP |            actions
+ C         |     |     |     |     |
- c    -----|-----|-----|-----|-----|-------------------------------------
+ C      -  |  -  |  -  |  -  |  -  | Read-in ustress, vstress, hflux,
- c         |     |     |     |     |
+ C         |     |     |     |     | swflux, and sflux.
- c      -  |  -  |  -  |  -  |  -  | Read-in ustress, vstress, hflux,
+ C         |     |     |     |     |
- c         |     |     |     |     | swflux, and sflux.
+ C     def | def | def | def |  -  | Read-in uwind, vwind, atemp, aqh,
- c         |     |     |     |     |
+ C         |     |     |     |     | swdown, lwdown, precip, and runoff.
- c     def | def | def | def |  -  | Read-in uwind, vwind, atemp, aqh,
+ C         |     |     |     |     | Compute ustress, vstress, hflux,
- c         |     |     |     |     | swdown, lwdown, precip, and runoff.
+ C         |     |     |     |     | swflux, and sflux.
- c         |     |     |     |     | Compute ustress, vstress, hflux,
+ C         |     |     |     |     |
- c         |     |     |     |     | swflux, and sflux.
+ C     def | def |  -  | def |  -  | Read-in uwind, vwind, atemp, aqh,
- c         |     |     |     |     |
+ C         |     |     |     |     | swflux, lwflux, precip, and runoff.
- c     def | def |  -  | def |  -  | Read-in uwind, vwind, atemp, aqh,
+ C         |     |     |     |     | Compute ustress, vstress, hflux,
- c         |     |     |     |     | swflux, lwflux, precip, and runoff.
+ C         |     |     |     |     | and sflux.
- c         |     |     |     |     | Compute ustress, vstress, hflux,
+ C         |     |     |     |     |
- c         |     |     |     |     | and sflux.
+ C     def |  -  |  -  | def |  -  | Read-in uwind, vwind, hflux,
- c         |     |     |     |     |
+ C         |     |     |     |     | swflux, and sflux.
- c     def |  -  |  -  | def |  -  | Read-in uwind, vwind, hflux,
+ C         |     |     |     |     | Compute ustress and vstress.
- c         |     |     |     |     | swflux, and sflux.
+ C         |     |     |     |     |
- c         |     |     |     |     | Compute ustress and vstress.
+ C      -  | def |  -  | def |  -  | Read-in ustress, vstress, atemp,
- c         |     |     |     |     |
+ C         |     |     |     |     | aqh, swflux, lwflux, precip, and
- c      -  | def |  -  | def |  -  | Read-in ustress, vstress, atemp,
+ C         |     |     |     |     | runoff.  Compute hflux and sflux.
- c         |     |     |     |     | aqh, swflux, lwflux, precip, and
+ C         |     |     |     |     |
- c         |     |     |     |     | runoff.  Compute hflux and sflux.
+ C     def | def |  -  |  -  | def | Read-in uwind, vwind, atemp, aqh,
- c         |     |     |     |     |
+ C         |     |     |     |     | swflux, lwflux, precip, runoff,
- c     def | def |  -  |  -  | def | Read-in uwind, vwind, atemp, aqh,
+ C         |     |     |     |     | and evap.
- c         |     |     |     |     | swflux, lwflux, precip, runoff,
+ C         |     |     |     |     |
- c         |     |     |     |     | and evap.
+ C     def | def |  -  | def |  -  | Read-in uwind, vwind, atemp, aqh,
- c         |     |     |     |     |
+ C         |     |     |     |     | swflux, lwflux, precip, and runoff.
- c     def | def |  -  | def |  -  | Read-in uwind, vwind, atemp, aqh,
+ C         |     |     |     |     | Compute open-water ustress, vstress,
- c         |     |     |     |     | swflux, lwflux, precip, and runoff.
+ C         |     |     |     |     | hflux, swflux, and evap.
- c         |     |     |     |     | Compute open-water ustress, vstress,
+ C         |     |     |     |     |
- c         |     |     |     |     | hflux, swflux, and evap.
+ C     def | def | def |  -  | def | Read-in uwind, vwind, atemp, aqh,
- c         |     |     |     |     |
+ C         |     |     |     |     | swdown, lwdown, precip, runoff,
- c     def | def | def |  -  | def | Read-in uwind, vwind, atemp, aqh,
+ C         |     |     |     |     | and evap.
- c         |     |     |     |     | swdown, lwdown, precip, runoff,
+ C         |     |     |     |     |
- c         |     |     |     |     | and evap.
+ C     def | def | def | def |  -  | Read-in uwind, vwind, atemp, aqh,
- c         |     |     |     |     |
+ C         |     |     |     |     | swdown, lwdown, precip, and runoff.
- c     def | def | def | def |  -  | Read-in uwind, vwind, atemp, aqh,
+ C         |     |     |     |     | Compute open-water ustress, vstress,
- c         |     |     |     |     | swdown, lwdown, precip, and runoff.
+ C         |     |     |     |     | hflux, swflux, and evap.
- c         |     |     |     |     | Compute open-water ustress, vstress,
+ C
- c         |     |     |     |     | hflux, swflux, and evap.
+ C   ====================================================================
- c
- c   ====================================================================
  C   Bulk formulae related flags.
  #define  ALLOW_ATM_TEMP
  #define  ALLOW_ATM_WIND
  #define  ALLOW_DOWNWARD_RADIATION
  #define  ALLOW_RUNOFF
- #if (defined (ALLOW_ATM_TEMP) || \
+ #undef   ALLOW_RUNOFTEMP
-      defined (ALLOW_ATM_WIND))
+ #if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND))
  # define ALLOW_BULKFORMULAE
- # undef ALLOW_BULK_LARGEYEAGER04
+ # undef  ALLOW_BULK_LARGEYEAGER04
+ #endif
+ #if (defined (ALLOW_ATM_TEMP) && defined (ATMOSPHERIC_LOADING))
+ # undef  EXF_CALC_ATMRHO
  #endif
  C   Zenith Angle/Albedo related flags.
-Line 152 
 C   Zenith Angle/Albedo related flags.
+Line 170 
 C   Zenith Angle/Albedo related flags.
  # undef ALLOW_ZENITHANGLE
  #endif
- C   Use ocean_emissivity*lwdwon in lwFlux. This flag should always be define
+ C   Use ocean_emissivity*lwdown in lwFlux. This flag should be defined
- C   but is not because of backward compatibility
+ C   unless to reproduce old results (obtained with inconsistent old code)
  #ifdef ALLOW_DOWNWARD_RADIATION
- # undef EXF_LWDOWN_WITH_EMISSIVITY
+ # define EXF_LWDOWN_WITH_EMISSIVITY
  #endif
  C   Relaxation to monthly climatologies.
-Line 165 
 C   Relaxation to monthly climatologies.
+Line 183 
 C   Relaxation to monthly climatologies.
  C   Use spatial interpolation to interpolate
  C   forcing files from input grid to model grid.
  #define USE_EXF_INTERPOLATION
- C   runoff is a special case for which one might want to bypass
+ C   for interpolated vector fields, rotate towards model-grid axis
- C   interpolation from an input grid
+ C   using old rotation formulae (instead of grid-angles)
- #ifdef USE_EXF_INTERPOLATION
+ #undef EXF_USE_OLD_VEC_ROTATION
- # undef USE_NO_INTERP_RUNOFF
+ C   for interpolation around N & S pole, use the old formulation
- #endif
+ C   (no pole symmetry, single vector-comp interp, reset to 0 zonal-comp @ N.pole)
+ #undef EXF_USE_OLD_INTERP_POLE
  #define EXF_INTERP_USE_DYNALLOC
- #if ( defined (EXF_INTERP_USE_DYNALLOC) & defined (USING_THREADS) )
+ #if ( defined (EXF_INTERP_USE_DYNALLOC) && defined (USING_THREADS) )
  # define EXF_IREAD_USE_GLOBAL_POINTER
  #endif
- #endif /* ndef ALLOW_AUTODIFF_TAMC */
+ #endif /* ndef ECCO_CPPOPTIONS_H */
  #endif /* ALLOW_EXF */
  #endif /* EXF_OPTIONS_H */

 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.3
 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.3
-Removed from v.1.2
+Added in v.1.3

	ViewVC Help
Powered by ViewVC 1.1.22