/[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.5 - (show annotations) (download)
Tue Feb 18 05:33:55 2003 UTC (21 years, 2 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint48f_post, checkpoint48i_post, checkpoint48h_post, checkpoint50b_pre, checkpoint49, checkpoint48g_post, checkpoint50, checkpoint50b_post, checkpoint50a_post
Changes since 1.4: +74 -37 lines
File MIME type: text/plain 
Merging from release1_p12:
o Modifications for using pkg/exf with pkg/seaice
  - improved description of the various forcing configurations
  - added basic radiation bulk formulae to pkg/exf
  - units/sign fix for evap computation in exf_getffields.F
  - updated verification/global_with_exf/results/output.txt
o Added pkg/sbo for computing IERS Special Bureau for the Oceans
  (SBO) core products, including oceanic mass, center-of-mass,
  angular, and bottom pressure (see pkg/sbo/README.sbo).
o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
  to avoid negative values in shallow regions.
  - updated verification/natl_box/results/output.txt
  - updated verification/lab_sea/results/output.txt
o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
o Added useSingleCpuIO option (see PARAMS.h).
o Updated useSingleCpuIO option in mdsio_writefield.F to
  work with multi-field files, e.g., for single-file pickup.
o pkg/seaice:
  - bug fix in growth.F: QNET for no shortwave case
  - added HeffFile for specifying initial sea-ice thickness
  - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
  kpp_transport_t.F, forward_step.F, and the_main_loop.F
o pkg/seaice:
  - adjoint-friendly modifications
  - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
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 Northeast B-grid U point
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 Northeast B-grid V point
23 C
24 C evap :: Evaporation in m/s
25 C > 0 for increase in salt (ocean salinity)
26 C Typical range: 0 < evap < 2.5e-7
27 C Northeast B-grid tracer point
28 C
29 C precip :: Precipitation in m/s
30 C > 0 for decrease in salt (ocean salinity)
31 C Typical range: 0 < precip < 5e-7
32 C Northeast B-grid tracer point
33 C
34 C runoff :: River and glacier runoff in m/s
35 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
74
75 #ifdef SEAICE_EXTERNAL_FORCING
76
77 C-- Define forcing fields outside pkg/seaice.
78 #include "exf_fields.h"
79
80 #else /* SEAICE_EXTERNAL_FORCING undefined */
81
82 C-- Define forcing fields internally.
83 COMMON /SEAICE_FFIELDS/ uwind, vwind, atemp,
84 & aqh, lwdown, swdown, precip, evap, runoff
85 _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)
87 _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)
89 _RL lwdown (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
90 _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)
92 _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)
94
95 COMMON /SEAICE_TDFIELDS/ uwind0, uwind1, vwind0, vwind1,
96 & atemp0, atemp1, aqh0, aqh1, lwdown0, lwdown1,
97 & swdown0, swdown1, precip0, precip1, evap0, evap1,
98 & runoff0, runoff1, SSSsi0, SSSsi1, SSTsi0, SSTsi1
99 _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)
101 _RS vwind0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
102 _RS vwind1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
103 _RS atemp0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
104 _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)
106 _RS aqh1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
107 _RS lwdown0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
108 _RS lwdown1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
109 _RS swdown0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
110 _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)
112 _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)
114 _RS evap1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
115 _RS runoff0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
116 _RS runoff1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
117 _RS SSSsi0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
118 _RS SSSsi1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
119 _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)
121
122 #endif /* SEAICE_EXTERNAL_FORCING */
123
124 #endif /* ALLOW_SEAICE */
  ViewVC Help
Powered by ViewVC 1.1.22