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

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

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

Revision 1.5 - (hide 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 heimbach 1.2 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 dimitri 1.5 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 heimbach 1.2 C
74 dimitri 1.4
75     #ifdef SEAICE_EXTERNAL_FORCING
76    
77     C-- Define forcing fields outside pkg/seaice.
78     #include "exf_fields.h"
79    
80 dimitri 1.5 #else /* SEAICE_EXTERNAL_FORCING undefined */
81 dimitri 1.4
82     C-- Define forcing fields internally.
83     COMMON /SEAICE_FFIELDS/ uwind, vwind, atemp,
84 dimitri 1.5 & aqh, lwdown, swdown, precip, evap, runoff
85 dimitri 1.4 _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 dimitri 1.5 _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 dimitri 1.4 _RL precip (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
92 dimitri 1.3 _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 heimbach 1.2
95 dimitri 1.4 COMMON /SEAICE_TDFIELDS/ uwind0, uwind1, vwind0, vwind1,
96 dimitri 1.5 & atemp0, atemp1, aqh0, aqh1, lwdown0, lwdown1,
97     & swdown0, swdown1, precip0, precip1, evap0, evap1,
98 heimbach 1.2 & runoff0, runoff1, SSSsi0, SSSsi1, SSTsi0, SSTsi1
99 dimitri 1.4 _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 dimitri 1.5 _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 dimitri 1.4 _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 heimbach 1.2 _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 dimitri 1.4
122 dimitri 1.5 #endif /* SEAICE_EXTERNAL_FORCING */
123 heimbach 1.2
124 dimitri 1.5 #endif /* ALLOW_SEAICE */
  ViewVC Help
Powered by ViewVC 1.1.22