/[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.8 - (show annotations) (download)
Thu Oct 9 04:19:20 2003 UTC (20 years, 7 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint53f_post, checkpoint54a_pre, checkpoint55c_post, checkpoint53b_pre, checkpoint51l_post, checkpoint52l_pre, checkpoint52e_pre, hrcube4, checkpoint52n_post, checkpoint52j_post, checkpoint51o_pre, checkpoint52e_post, checkpoint51n_pre, checkpoint57b_post, checkpoint52d_pre, checkpoint53c_post, checkpoint53d_post, checkpoint55d_pre, checkpoint57a_post, checkpoint55j_post, checkpoint56b_post, checkpoint52j_pre, checkpoint54a_post, branch-netcdf, checkpoint52l_post, checkpoint55h_post, checkpoint51r_post, checkpoint52k_post, checkpoint52b_pre, checkpoint54b_post, checkpoint51i_post, checkpoint54d_post, checkpoint56c_post, checkpoint54e_post, checkpoint55b_post, checkpoint51l_pre, checkpoint52m_post, checkpoint55, checkpoint53a_post, checkpoint55a_post, hrcube5, checkpoint53b_post, checkpoint55g_post, checkpoint51o_post, checkpoint55f_post, checkpoint51q_post, checkpoint57a_pre, checkpoint54, checkpoint57, checkpoint56, checkpoint53, checkpoint52, checkpoint52d_post, checkpoint52a_post, checkpoint53g_post, checkpoint52f_post, checkpoint52c_post, ecco_c52_e35, checkpoint54f_post, checkpoint52a_pre, checkpoint52b_post, checkpoint51m_post, checkpoint51t_post, checkpoint53d_pre, checkpoint55e_post, checkpoint54c_post, checkpoint52i_post, checkpoint51p_post, checkpoint51n_post, checkpoint55i_post, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint56a_post, checkpoint51s_post, checkpoint55d_post
Branch point for: branch-nonh, netcdf-sm0, tg2-branch, checkpoint51n_branch
Changes since 1.7: +2 -1 lines
File MIME type: text/plain 
 o first check-in for the "branch-genmake2" merge
 o verification suite as run on shelley (gcc 3.2.2):

Wed Oct  8 23:42:29 EDT 2003
                T           S           U           V
G D M    c        m  s        m  s        m  s        m  s
E p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
N n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
2 d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .

OPTFILE=NONE

Y Y Y Y 13 16 16 16  0 16 16 16 16 16 16 16 16 13 12  0  0 pass  adjustment.128x64x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16  0  0 16 16  0  0 pass  adjustment.cs-32x32x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16 22  0 16 16 22  0 pass  adjust_nlfs.cs-32x32x1
Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O   advect_cs
Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O   advect_xy
Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O   advect_xz
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  aim.5l_cs
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass  aim.5l_Equatorial_Channel
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass  aim.5l_LatLon
Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass  exp0
Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass  exp1
Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass  exp2
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  exp4
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass  exp5
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  front_relax
Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass  global_ocean.90x40x15
Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL  global_ocean.cs32x15
Y Y Y Y  6 11 12 13 13 12 13 16 13  9  9  9  9 10  9  9 11 FAIL  global_ocean_pressure
Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass  global_with_exf
Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16  9 16 pass  hs94.128x64x5
Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass  hs94.1x64x5
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass  hs94.cs-32x32x5
Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL  ideal_2D_oce
Y Y Y Y  8 16 16 16 16 16 16 16 16 13 13  8 16 16 16 16 16 FAIL  internal_wave
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass  inverted_barometer
Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL  lab_sea
Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL  natl_box
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  plume_on_slope
Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass  solid-body.cs-32x32x1
1 C $Header: /u/u3/gcmpack/MITgcm/pkg/seaice/SEAICE_FFIELDS.h,v 1.7.2.1 2003/10/02 18:18:33 adcroft Exp $
2 C $Name: $
3
4 #ifdef ALLOW_SEAICE
5
6 C /==========================================================\
7 C | SEAICE_FFIELDS.h |
8 C | o Sea ice model forcing fields |
9 C |==========================================================|
10 C \==========================================================/
11 C
12 C
13 C uwind :: Surface (10-m) zonal wind velocity in m/s
14 C > 0 for increase in uVel, which is west to
15 C east for cartesian and spherical polar grids
16 C Typical range: -10 < uwind < 10
17 C
18 C vwind :: Surface (10-m) meridional wind velocity in m/s
19 C > 0 for increase in vVel, which is south to
20 C north for cartesian and spherical polar grids
21 C Typical range: -10 < vwind < 10
22 C
23 C evap :: Evaporation in m/s
24 C > 0 for increase in salt (ocean salinity)
25 C Typical range: 0 < evap < 2.5e-7
26 C
27 C precip :: Precipitation in m/s
28 C > 0 for decrease in salt (ocean salinity)
29 C Typical range: 0 < precip < 5e-7
30 C
31 C runoff :: River and glacier runoff in m/s
32 C > 0 for decrease in salt (ocean salinity)
33 C Typical range: 0 < runoff < ????
34 C
35 C atemp :: Surface (2-m) air temperature in deg K
36 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
66
67 #ifdef SEAICE_EXTERNAL_FORCING
68
69 C-- Define forcing fields outside pkg/seaice.
70 #include "exf_fields.h"
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 #endif /* ALLOW_SEAICE */
  ViewVC Help
Powered by ViewVC 1.1.22