/[MITgcm]/MITgcm/pkg/seaice/SEAICE_FFIELDS.h
ViewVC logotype

Diff of /MITgcm/pkg/seaice/SEAICE_FFIELDS.h

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

revision 1.4 by dimitri, Sat Dec 28 10:11:11 2002 UTC revision 1.8 by edhill, Thu Oct 9 04:19:20 2003 UTC
# Line 1  Line 1 
1  C  C $Header$
2    C $Name$
3    
4  #ifdef ALLOW_SEAICE  #ifdef ALLOW_SEAICE
5    
# Line 8  C     | o Sea ice model forcing fields Line 9  C     | o Sea ice model forcing fields
9  C     |==========================================================|  C     |==========================================================|
10  C     \==========================================================/  C     \==========================================================/
11  C  C
12  C     uwind  - Surface (10-m) zonal wind velocity in m/s  C
13  C              at North-East B-grid U point  C     uwind     :: Surface (10-m) zonal wind velocity in m/s
14  C              >0 from West to East  C                  > 0 for increase in uVel, which is west to
15  C     vwind  - Surface (10-m) meridional wind velocity in m/s  C                      east for cartesian and spherical polar grids
16  C              at North-East B-grid V point  C                  Typical range: -10 < uwind < 10
17  C              >0 from South to North  C
18  C     atemp  - Surface (2-m) air temperature in deg K  C     vwind     :: Surface (10-m) meridional wind velocity in m/s
19  C              at North-East B-grid tracer point  C                  > 0 for increase in vVel, which is south to
20  C     aqh    - Surface (2m) specific humidity in kg/kg  C                      north for cartesian and spherical polar grids
21  C              at North-East B-grid tracer point  C                  Typical range: -10 < vwind < 10
22  C     lwflux - Downward longwave radiation in W/m^2  C
23  C              at North-East B-grid tracer point  C     evap      :: Evaporation in m/s
24  C              >0 for ocean warming  C                  > 0 for increase in salt (ocean salinity)
25  C     swflux - Downward shortwave radiation in W/m^2  C                  Typical range: 0 < evap < 2.5e-7
26  C              at North-East B-grid tracer point  C
27  C              >0 for ocean warming  C     precip    :: Precipitation in m/s
28  C     precip - Precipitation in m/s  C                  > 0 for decrease in salt (ocean salinity)
29  C              at North-East B-grid tracer point  C                  Typical range: 0 < precip < 5e-7
30  C              >0 decreases salinity  C
31  C     evap   - Evaporation in m/s  C     runoff    :: River and glacier runoff in m/s
32  C              at North-East B-grid tracer point  C                  > 0 for decrease in salt (ocean salinity)
33  C              >0 increases salinity  C                  Typical range: 0 < runoff < ????
34  C     runoff - River and glacier runoff in m/s  C
35  C              at North-East B-grid tracer point  C     atemp     :: Surface (2-m) air temperature in deg K
36  C              >0 decreases salinity  C                  Typical range: 200 < atemp < 300
37    C
38    C     aqh       :: Surface (2m) specific humidity in kg/kg
39    C                  Typical range: 0 < aqh < 0.02
40    C
41    C     swdown    :: Downward shortwave radiation in W/m^2
42    C                  > 0 for increase in theta (ocean warming)
43    C                  Typical range: 0 < swdown < 450
44    C
45    C     lwdown    :: Downward longwave radiation in W/m^2
46    C                  > 0 for increase in theta (ocean warming)
47    C                  Typical range: 50 < lwdown < 450
48    C
49    C
50    C     NOTES:
51    C     ======
52    C
53    C     All surface forcing fields are defined at the center of
54    C     each grid (the rVel location in model/inc/GRID.h).
55    C
56    C     #ifdef SEAICE_EXTERNAL_FORCING, sea-ice forcing fields
57    C     are defined in exf_fields.h
58    C
59    C     #ifdef SEAICE_EXTERNAL_FLUXES, additional forcing fields, fu,
60    C     fv, Qnet, Qsw, and EmPmR, as defined in FFIELDS.h, are required.
61    C
62    C     Downward and net radiative fluxes are not the same !!!
63    C     Downward radiative fluxes, swdown and lwdown, are used by
64    C     radiation bulk formulae to compute net radiative fluxes.
65  C  C
66    
67  #ifdef SEAICE_EXTERNAL_FORCING  #ifdef SEAICE_EXTERNAL_FORCING
# Line 40  C Line 69  C
69  C--   Define forcing fields outside pkg/seaice.  C--   Define forcing fields outside pkg/seaice.
70  #include "exf_fields.h"  #include "exf_fields.h"
71    
72  #else SEAICE_EXTERNAL_FORCING  #else /* SEAICE_EXTERNAL_FORCING undefined */
73    
74  C--   Define forcing fields internally.  C--   Define forcing fields internally.
75        COMMON /SEAICE_FFIELDS/ uwind, vwind, atemp,        COMMON /SEAICE_FFIELDS/ uwind, vwind, atemp,
76       &     aqh, lwflux, swflux, precip, evap, runoff       &     aqh, lwdown, swdown, precip, evap, runoff
77        _RL  uwind    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _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)        _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)        _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)        _RL  aqh      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
81        _RL  lwflux   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  lwdown   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
82        _RL  swflux   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _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)        _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)        _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)        _RL  runoff   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
86    
87        COMMON /SEAICE_TDFIELDS/ uwind0, uwind1, vwind0, vwind1,        COMMON /SEAICE_TDFIELDS/ uwind0, uwind1, vwind0, vwind1,
88       &     atemp0, atemp1, aqh0, aqh1, lwflux0, lwflux1,       &     atemp0, atemp1, aqh0, aqh1, lwdown0, lwdown1,
89       &     swflux0, swflux1, precip0, precip1, evap0, evap1,       &     swdown0, swdown1, precip0, precip1, evap0, evap1,
90       &     runoff0, runoff1, SSSsi0, SSSsi1, SSTsi0, SSTsi1       &     runoff0, runoff1, SSSsi0, SSSsi1, SSTsi0, SSTsi1
91        _RS  uwind0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _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)        _RS  uwind1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 67  C--   Define forcing fields internally. Line 96  C--   Define forcing fields internally.
96        _RS  atemp1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _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)        _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)        _RS  aqh1     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
99        _RS  lwflux0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  lwdown0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
100        _RS  lwflux1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  lwdown1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
101        _RS  swflux0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  swdown0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
102        _RS  swflux1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _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)        _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)        _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)        _RS  evap0    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 82  C--   Define forcing fields internally. Line 111  C--   Define forcing fields internally.
111        _RS  SSTsi0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _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)        _RS  SSTsi1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
113    
114  #endif SEAICE_EXTERNAL_FORCING  #endif /* SEAICE_EXTERNAL_FORCING */
115    
116  #endif ALLOW_SEAICE  #endif /* ALLOW_SEAICE */
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22