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