pkg/seaice/SEAICE_FFIELDS.h

C $Header: /u/gcmpack/MITgcm/pkg/seaice/SEAICE_FFIELDS.h,v 1.12 2007/05/14 14:17:25 jmc Exp $
C $Name:  $

C     /==========================================================\
C     | SEAICE_FFIELDS.h                                         |
C     | o Sea ice model forcing fields                           |
C     |==========================================================|
C     \==========================================================/
C
C
C     uwind     :: Surface (10-m) zonal wind velocity in m/s
C                  > 0 for increase in uVel, which is west to
C                      east for cartesian and spherical polar grids
C                  Typical range: -10 < uwind < 10
C
C     vwind     :: Surface (10-m) meridional wind velocity in m/s
C                  > 0 for increase in vVel, which is south to
C                      north for cartesian and spherical polar grids
C                  Typical range: -10 < vwind < 10
C
C     evap      :: Evaporation in m/s
C                  > 0 for increase in salt (ocean salinity)
C                  Typical range: 0 < evap < 2.5e-7
C
C     precip    :: Precipitation in m/s
C                  > 0 for decrease in salt (ocean salinity)
C                  Typical range: 0 < precip < 5e-7
C
C     runoff    :: River and glacier runoff in m/s
C                  > 0 for decrease in salt (ocean salinity)
C                  Typical range: 0 < runoff < ????
C
C     atemp     :: Surface (2-m) air temperature in deg K
C                  Typical range: 200 < atemp < 300
C
C     aqh       :: Surface (2m) specific humidity in kg/kg
C                  Typical range: 0 < aqh < 0.02
C
C     swdown    :: Downward shortwave radiation in W/m^2
C                  > 0 for increase in theta (ocean warming)
C                  Typical range: 0 < swdown < 450
C
C     lwdown    :: Downward longwave radiation in W/m^2
C                  > 0 for increase in theta (ocean warming)
C                  Typical range: 50 < lwdown < 450
C
C
C     NOTES:
C     ======
C
C     All surface forcing fields are defined at the center of
C     each grid (the rVel location in model/inc/GRID.h).
C
C     #ifdef SEAICE_EXTERNAL_FORCING, sea-ice forcing fields
C     are defined in EXF_FIELDS.h
C
C     #ifdef SEAICE_EXTERNAL_FLUXES, additional forcing fields, fu,
C     fv, Qnet, Qsw, and EmPmR, as defined in FFIELDS.h, are required.
C
C     Downward and net radiative fluxes are not the same !!!
C     Downward radiative fluxes, swdown and lwdown, are used by
C     radiation bulk formulae to compute net radiative fluxes.
C

#ifdef SEAICE_EXTERNAL_FORCING

C--   Define forcing fields outside pkg/seaice.
#ifdef ALLOW_EXF
#include "EXF_FIELDS.h"
#endif

#else /* SEAICE_EXTERNAL_FORCING undefined */

C--   Define forcing fields internally.
      COMMON /SEAICE_FFIELDS/ uwind, vwind, atemp,
     &     aqh, lwdown, swdown, precip, evap, runoff
      _RL  uwind    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  vwind    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  atemp    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  aqh      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  lwdown   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  swdown   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  precip   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  evap     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  runoff   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)

      COMMON /SEAICE_TDFIELDS/ uwind0, uwind1, vwind0, vwind1,
     &     atemp0, atemp1, aqh0, aqh1, lwdown0, lwdown1,
     &     swdown0, swdown1, precip0, precip1, evap0, evap1,
     &     runoff0, runoff1, SSSsi0, SSSsi1, SSTsi0, SSTsi1
      _RS  uwind0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  uwind1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  vwind0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  vwind1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  atemp0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  atemp1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  aqh0     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  aqh1     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  lwdown0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  lwdown1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  swdown0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  swdown1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  precip0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  precip1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  evap0    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  evap1    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  runoff0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  runoff1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  SSSsi0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  SSSsi1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  SSTsi0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS  SSTsi1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)

#endif /* SEAICE_EXTERNAL_FORCING */

CEH3 ;;; Local Variables: ***
CEH3 ;;; mode:fortran ***
CEH3 ;;; End: ***
1	C $Header: /u/gcmpack/MITgcm/pkg/seaice/SEAICE_FFIELDS.h,v 1.12 2007/05/14 14:17:25 jmc Exp $
2	C $Name: $
3
4	C /==========================================================\
5	C \| SEAICE_FFIELDS.h \|
6	C \| o Sea ice model forcing fields \|
7	C \|==========================================================\|
8	C \==========================================================/
9	C
10	C
11	C uwind :: Surface (10-m) zonal wind velocity in m/s
12	C > 0 for increase in uVel, which is west to
13	C east for cartesian and spherical polar grids
14	C Typical range: -10 < uwind < 10
15	C
16	C vwind :: Surface (10-m) meridional wind velocity in m/s
17	C > 0 for increase in vVel, which is south to
18	C north for cartesian and spherical polar grids
19	C Typical range: -10 < vwind < 10
20	C
21	C evap :: Evaporation in m/s
22	C > 0 for increase in salt (ocean salinity)
23	C Typical range: 0 < evap < 2.5e-7
24	C
25	C precip :: Precipitation in m/s
26	C > 0 for decrease in salt (ocean salinity)
27	C Typical range: 0 < precip < 5e-7
28	C
29	C runoff :: River and glacier runoff in m/s
30	C > 0 for decrease in salt (ocean salinity)
31	C Typical range: 0 < runoff < ????
32	C
33	C atemp :: Surface (2-m) air temperature in deg K
34	C Typical range: 200 < atemp < 300
35	C
36	C aqh :: Surface (2m) specific humidity in kg/kg
37	C Typical range: 0 < aqh < 0.02
38	C
39	C swdown :: Downward shortwave radiation in W/m^2
40	C > 0 for increase in theta (ocean warming)
41	C Typical range: 0 < swdown < 450
42	C
43	C lwdown :: Downward longwave radiation in W/m^2
44	C > 0 for increase in theta (ocean warming)
45	C Typical range: 50 < lwdown < 450
46	C
47	C
48	C NOTES:
49	C ======
50	C
51	C All surface forcing fields are defined at the center of
52	C each grid (the rVel location in model/inc/GRID.h).
53	C
54	C #ifdef SEAICE_EXTERNAL_FORCING, sea-ice forcing fields
55	C are defined in EXF_FIELDS.h
56	C
57	C #ifdef SEAICE_EXTERNAL_FLUXES, additional forcing fields, fu,
58	C fv, Qnet, Qsw, and EmPmR, as defined in FFIELDS.h, are required.
59	C
60	C Downward and net radiative fluxes are not the same !!!
61	C Downward radiative fluxes, swdown and lwdown, are used by
62	C radiation bulk formulae to compute net radiative fluxes.
63	C
64
65	#ifdef SEAICE_EXTERNAL_FORCING
66
67	C-- Define forcing fields outside pkg/seaice.
68	#ifdef ALLOW_EXF
69	#include "EXF_FIELDS.h"
70	#endif
71
72	#else /* SEAICE_EXTERNAL_FORCING undefined */
73
74	C-- Define forcing fields internally.
75	COMMON /SEAICE_FFIELDS/ uwind, vwind, atemp,
76	& aqh, lwdown, swdown, precip, evap, runoff
77	_RL uwind (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
78	_RL vwind (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
79	_RL atemp (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
80	_RL aqh (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
81	_RL lwdown (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
82	_RL swdown (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
83	_RL precip (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
84	_RL evap (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
85	_RL runoff (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
86
87	COMMON /SEAICE_TDFIELDS/ uwind0, uwind1, vwind0, vwind1,
88	& atemp0, atemp1, aqh0, aqh1, lwdown0, lwdown1,
89	& swdown0, swdown1, precip0, precip1, evap0, evap1,
90	& runoff0, runoff1, SSSsi0, SSSsi1, SSTsi0, SSTsi1
91	_RS uwind0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
92	_RS uwind1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
93	_RS vwind0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
94	_RS vwind1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
95	_RS atemp0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
96	_RS atemp1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
97	_RS aqh0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
98	_RS aqh1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
99	_RS lwdown0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
100	_RS lwdown1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
101	_RS swdown0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
102	_RS swdown1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
103	_RS precip0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
104	_RS precip1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
105	_RS evap0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
106	_RS evap1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
107	_RS runoff0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
108	_RS runoff1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
109	_RS SSSsi0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
110	_RS SSSsi1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
111	_RS SSTsi0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
112	_RS SSTsi1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
113
114	#endif /* SEAICE_EXTERNAL_FORCING */
115
116	CEH3 ;;; Local Variables: ***
117	CEH3 ;;; mode:fortran ***
118	CEH3 ;;; End: ***