-
Notifications
You must be signed in to change notification settings - Fork 0
/
cregin.h
65 lines (65 loc) · 3.21 KB
/
cregin.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
c====================== include file "cregin.h" ========================
c
c variables used for computing regional tracer averages (see
c "reg1st.F" & sub "diagt1" in "tracer.F") and for computing term
c balances for tracer and momentum equations (see "clinic.F",
c "tracer.F" and "diag2.F")
c
c
c mskhr = mask field defining regions in the horizontal
c (eg: mskhr(i,j) = n indicates point (i,j) is in the
c "nth" horizontal region where n=1..nhreg)
c The "mskhr" masks are used in "diagt1" when
c computing volume weighted tracer averages and in
c "clinic.F", "tracer.F" and "diag2.F" when computing
c term balances for tracers and momentum.
c mskvr = mask field defining regions in the vertical
c (eg: mskvr(k) = m indicates all points under a horizontal
c mask at level "k" are in the "mth" vertical region
c where m=1..nvreg)
c The "mskvr" masks are used in "diag.F", but not
c in "diagt1", where tracer averages are calculated
c for each k-level.
c
c hregnm = horizontal region name
c vregnm = vertical region name
c volbt = total volume under a given horizontal region
c volbk = volume contained in a horizontal region at level "k"
c volgt = total ocean volume
c volgk = total ocean volume at level "k"
c areab = total ocean surface area for a given horizontal region
c areag = total ocean surface area
c
c volt = ocean volume within a particular horizontal & vertical
c region (on the "t" grid) for tracer term balances
c volt(0) represents the sum of all regions
c rvolt = 1/volt ( 0.0 if volt = 0.0)
c areat = horizontal ocean surface area corresponding to "volt"
c areat(0) represents the sum of all regions
c rareat = 1/areat ( 0.0 if areat = 0.0)
c volu = ocean volume within a particular horizontal & vertical
c region (on the "u" grid) for momentum term balances
c volu(0) represents the sum of all regions
c rvolu = 1/volu ( 0.0 if volu = 0.0)
c areau = horizontal ocean area corresponding to "volu"
c areau(0) represents the sum of all regions
c rareau = 1/areau ( 0.0 if areau = 0.0)
c llvreg = level limits for defining vertical regions in term
c balance calculations (not used in computing volume
c weighted tracer averages)
c (eg: llvreg(3,1) = 4... means that starting level for
c the third region in the vertical is 4. similarly,
c llvreg(3,2) = 6 means the ending level is 6 for that
c region. note regions should not overlap.)
c
character*40 hregnm
character*20 vregnm
c
common /cregni/ mskhr(imt,jmt), mskvr(km), llvreg(numreg,2)
common /cregnr/ volbk(nhreg,km), volbt(nhreg), volgk(km)
common /cregnr/ volgt, areab(nhreg), areag, volt(0:numreg)
common /cregnr/ volu(0:numreg), areat(0:numreg), areau(0:numreg)
common /cregnr/ rvolt(0:numreg), rvolu(0:numreg), rareat(0:numreg)
common /cregnr/ rareau(0:numreg)
common /cregnc/ hregnm(nhreg), vregnm(nvreg)
c