/[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.5 by dimitri, Tue Feb 18 05:33:55 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                  Northeast B-grid U point
17  C     atemp  - Surface (2-m) air temperature in deg K  C
18  C              at North-East B-grid tracer point  C     vwind     :: Surface (10-m) meridional wind velocity in m/s
19  C     aqh    - Surface (2m) specific humidity in kg/kg  C                  > 0 for increase in vVel, which is south to
20  C              at North-East B-grid tracer point  C                      north for cartesian and spherical polar grids
21  C     lwflux - Downward longwave radiation in W/m^2  C                  Typical range: -10 < vwind < 10
22  C              at North-East B-grid tracer point  C                  Northeast B-grid V point
23  C              >0 for ocean warming  C
24  C     swflux - Downward shortwave radiation in W/m^2  C     evap      :: Evaporation in m/s
25  C              at North-East B-grid tracer point  C                  > 0 for increase in salt (ocean salinity)
26  C              >0 for ocean warming  C                  Typical range: 0 < evap < 2.5e-7
27  C     precip - Precipitation in m/s  C                  Northeast B-grid tracer point
28  C              at North-East B-grid tracer point  C
29  C              >0 decreases salinity  C     precip    :: Precipitation in m/s
30  C     evap   - Evaporation in m/s  C                  > 0 for decrease in salt (ocean salinity)
31  C              at North-East B-grid tracer point  C                  Typical range: 0 < precip < 5e-7
32  C              >0 increases salinity  C                  Northeast B-grid tracer point
33  C     runoff - River and glacier runoff in m/s  C
34  C              at North-East B-grid tracer point  C     runoff    :: River and glacier runoff in m/s
35  C              >0 decreases salinity  C                  > 0 for decrease in salt (ocean salinity)
36    C                  Typical range: 0 < runoff < ????
37    C                  Northeast B-grid tracer point
38    C
39    C     atemp     :: Surface (2-m) air temperature in deg K
40    C                  Typical range: 200 < atemp < 300
41    C                  Northeast B-grid tracer point
42    C
43    C     aqh       :: Surface (2m) specific humidity in kg/kg
44    C                  Typical range: 0 < aqh < 0.02
45    C                  Northeast B-grid tracer point
46    C
47    C     swdown    :: Downward shortwave radiation in W/m^2
48    C                  > 0 for increase in theta (ocean warming)
49    C                  Typical range: 0 < swdown < 450
50    C                  Northeast B-grid tracer point
51    C
52    C     lwdown    :: Downward longwave radiation in W/m^2
53    C                  > 0 for increase in theta (ocean warming)
54    C                  Typical range: 50 < lwdown < 450
55    C                  Northeast B-grid tracer point
56    C
57    C
58    C     NOTES:
59    C     ======
60    C
61    C     #ifdef SEAICE_EXTERNAL_FORCING, sea-ice forcing fields
62    C     are defined in exf_fields.h
63    C
64    C     #ifdef SEAICE_EXTERNAL_FLUXES, additional forcing fields, fu,
65    C     fv, Qnet, Qsw, and EmPmR, as defined in FFIELDS.h, are required.
66    C
67    C     If (SEAICEwindOnCgrid .EQ. .TRUE.), uwind and vwind are
68    C     defined on southwest C-grid U and V points, respectively.
69    C
70    C     Downward and net radiative fluxes are not the same !!!
71    C     Downward radiative fluxes, swdown and lwdown, are used by
72    C     radiation bulk formulae to compute net radiative fluxes.
73  C  C
74    
75  #ifdef SEAICE_EXTERNAL_FORCING  #ifdef SEAICE_EXTERNAL_FORCING
# Line 40  C Line 77  C
77  C--   Define forcing fields outside pkg/seaice.  C--   Define forcing fields outside pkg/seaice.
78  #include "exf_fields.h"  #include "exf_fields.h"
79    
80  #else SEAICE_EXTERNAL_FORCING  #else /* SEAICE_EXTERNAL_FORCING undefined */
81    
82  C--   Define forcing fields internally.  C--   Define forcing fields internally.
83        COMMON /SEAICE_FFIELDS/ uwind, vwind, atemp,        COMMON /SEAICE_FFIELDS/ uwind, vwind, atemp,
84       &     aqh, lwflux, swflux, precip, evap, runoff       &     aqh, lwdown, swdown, precip, evap, runoff
85        _RL  uwind    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  uwind    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
86        _RL  vwind    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  vwind    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
87        _RL  atemp    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  atemp    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
88        _RL  aqh      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  aqh      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
89        _RL  lwflux   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  lwdown   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
90        _RL  swflux   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  swdown   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
91        _RL  precip   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  precip   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
92        _RL  evap     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  evap     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
93        _RL  runoff   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  runoff   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
94    
95        COMMON /SEAICE_TDFIELDS/ uwind0, uwind1, vwind0, vwind1,        COMMON /SEAICE_TDFIELDS/ uwind0, uwind1, vwind0, vwind1,
96       &     atemp0, atemp1, aqh0, aqh1, lwflux0, lwflux1,       &     atemp0, atemp1, aqh0, aqh1, lwdown0, lwdown1,
97       &     swflux0, swflux1, precip0, precip1, evap0, evap1,       &     swdown0, swdown1, precip0, precip1, evap0, evap1,
98       &     runoff0, runoff1, SSSsi0, SSSsi1, SSTsi0, SSTsi1       &     runoff0, runoff1, SSSsi0, SSSsi1, SSTsi0, SSTsi1
99        _RS  uwind0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  uwind0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
100        _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 104  C--   Define forcing fields internally.
104        _RS  atemp1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  atemp1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
105        _RS  aqh0     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  aqh0     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
106        _RS  aqh1     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  aqh1     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
107        _RS  lwflux0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  lwdown0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
108        _RS  lwflux1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  lwdown1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
109        _RS  swflux0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  swdown0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
110        _RS  swflux1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  swdown1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
111        _RS  precip0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  precip0  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
112        _RS  precip1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  precip1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
113        _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 119  C--   Define forcing fields internally.
119        _RS  SSTsi0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  SSTsi0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
120        _RS  SSTsi1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  SSTsi1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
121    
122  #endif SEAICE_EXTERNAL_FORCING  #endif /* SEAICE_EXTERNAL_FORCING */
123    
124  #endif ALLOW_SEAICE  #endif /* ALLOW_SEAICE */

Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5

  ViewVC Help
Powered by ViewVC 1.1.22