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

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

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

Revision 1.7 - (show annotations) (download)
Thu Aug 7 02:31:29 2003 UTC (20 years, 9 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint51j_post, branchpoint-genmake2, checkpoint51f_pre, checkpoint51h_pre, checkpoint51g_post, checkpoint51f_post, checkpoint51i_pre
Branch point for: branch-genmake2
Changes since 1.6: +3 -12 lines
File MIME type: text/plain 
o Added on-the-fly spatial interpolation capability
    "USE_EXF_INTERPOLATION" to pkg/exf.
  - This is a temporary Cartesian-grid hack until
    the super-duper ESMF coupler becomes available.
  - See verification/global_with_exf/README for usage example.
  - Removed obsolete EXFwindOnBgrid and SEAICEwindOnCgrid
    flags and modified pkg/seaice accordingly.
o Bug fix to pkg/ptracers, pkg/generic_advdiff/gad_calc_rhs.F,
    and pkg/kpp/kpp_transport_ptr.F for dealing with tracer
    non-local transport term.
1 C $Header:
2
3 #ifdef ALLOW_SEAICE
4
5 C /==========================================================\
6 C | SEAICE_FFIELDS.h |
7 C | o Sea ice model forcing fields |
8 C |==========================================================|
9 C \==========================================================/
10 C
11 C
12 C uwind :: Surface (10-m) zonal wind velocity in m/s
13 C > 0 for increase in uVel, which is west to
14 C east for cartesian and spherical polar grids
15 C Typical range: -10 < uwind < 10
16 C
17 C vwind :: Surface (10-m) meridional wind velocity in m/s
18 C > 0 for increase in vVel, which is south to
19 C north for cartesian and spherical polar grids
20 C Typical range: -10 < vwind < 10
21 C
22 C evap :: Evaporation in m/s
23 C > 0 for increase in salt (ocean salinity)
24 C Typical range: 0 < evap < 2.5e-7
25 C
26 C precip :: Precipitation in m/s
27 C > 0 for decrease in salt (ocean salinity)
28 C Typical range: 0 < precip < 5e-7
29 C
30 C runoff :: River and glacier runoff in m/s
31 C > 0 for decrease in salt (ocean salinity)
32 C Typical range: 0 < runoff < ????
33 C
34 C atemp :: Surface (2-m) air temperature in deg K
35 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
65
66 #ifdef SEAICE_EXTERNAL_FORCING
67
68 C-- Define forcing fields outside pkg/seaice.
69 #include "exf_fields.h"
70
71 #else /* SEAICE_EXTERNAL_FORCING undefined */
72
73 C-- Define forcing fields internally.
74 COMMON /SEAICE_FFIELDS/ uwind, vwind, atemp,
75 & aqh, lwdown, swdown, precip, evap, runoff
76 _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)
78 _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)
80 _RL lwdown (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
81 _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)
83 _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)
85
86 COMMON /SEAICE_TDFIELDS/ uwind0, uwind1, vwind0, vwind1,
87 & atemp0, atemp1, aqh0, aqh1, lwdown0, lwdown1,
88 & swdown0, swdown1, precip0, precip1, evap0, evap1,
89 & runoff0, runoff1, SSSsi0, SSSsi1, SSTsi0, SSTsi1
90 _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)
92 _RS vwind0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
93 _RS vwind1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
94 _RS atemp0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
95 _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)
97 _RS aqh1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
98 _RS lwdown0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
99 _RS lwdown1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
100 _RS swdown0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
101 _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)
103 _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)
105 _RS evap1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
106 _RS runoff0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
107 _RS runoff1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
108 _RS SSSsi0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
109 _RS SSSsi1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
110 _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)
112
113 #endif /* SEAICE_EXTERNAL_FORCING */
114
115 #endif /* ALLOW_SEAICE */
  ViewVC Help
Powered by ViewVC 1.1.22