/[MITgcm]/MITgcm/doc/diags_changes.txt
ViewVC logotype

Contents of /MITgcm/doc/diags_changes.txt

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


Revision 1.2 - (show annotations) (download)
Mon Jun 22 03:14:37 2009 UTC (14 years, 9 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint62, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62d, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61q, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.1: +9 -5 lines
File MIME type: text/plain
Added a few more flux diagnostic equivalencies.

1 Jun 21, 2009
2 Some redundancies and relations for model diagnostics:
3
4 1) Qnet = QNETtave = - oceQnet = SIqnet
5
6 2) For open-ocean, i.e., when there is no sea ice:
7 Qnet = QNETtave = - oceQnet = SIqnet = EXFqnet = SIqneto = SIatmQnt
8
9 3) EXFqnet = EXFlwnet + EXFswnet - EXFhl - EXFhs
10
11 4) Qsw = QSWtave = - oceQsw = SIqsw
12
13 5) For open-ocean, i.e., when there is no sea ice:
14 Qsw = QSWtave = - oceQsw = SIqsw = EXFswnet
15
16 6) EmPmR = EmPmRtave = -oceFWflx = SIempmr
17
18 ================================================================
19
20 Nov 19, 2006 (after tag checkpoint58r_post)
21
22 I) Some diagnostics have been renamed (essentially, to better
23 match the content):
24
25 1) name: description (oceanic set-up):
26 PRESSURE Cell-Center Height
27 replaced by:
28 RCENTER Cell-Center Height
29
30 2) name: description (oceanic set-up):
31 TICE heat from melt/freeze of sea-ice, >0 increases theta
32 replaced by:
33 oceFreez heating from freezing of sea-water (allowFreezing=T)
34
35 3) name: description (oceanic set-up):
36 TAUX zonal surface wind stress, >0 increases uVel
37 TAUY meridional surf. wind stress, >0 increases vVel
38 replaced by:
39 oceTAUX zonal surface wind stress, >0 increases uVel
40 oceTAUY meridional surf. wind stress, >0 increases vVel
41
42 4) name: description (oceanic set-up):
43 SWFLUX net upward SW radiation, >0 increases theta
44 replaced by:
45 oceQsw net Short-Wave radiation (+=down), >0 increases theta
46
47 5) name: description (oceanic set-up):
48 DIFx_TH Zonal Diffusive Flux of Pot.Temperature
49 DIFy_TH Meridional Diffusive Flux of Pot.Temperature
50 replaced by:
51 DFxE_TH Zonal Diffusive Flux of Pot.Temperature
52 DFyE_TH Meridional Diffusive Flux of Pot.Temperature
53
54 6) name: description (oceanic set-up):
55 DIFx_SLT Zonal Diffusive Flux of Salinity
56 DIFy_SLT Meridional Diffusive Flux of Salinity
57 replaced by:
58 DFxE_SLT Zonal Diffusive Flux of Salinity
59 DFyE_SLT Meridional Diffusive Flux of Salinity
60
61 -----------------------------------------------
62 II) Change description of existing diagnostics:
63
64 name: old description (oceanic set-up):
65 TFLUX net surface heat flux, >0 increases theta
66 SFLUX net surface salt flux, >0 increases salt
67 name: new description (oceanic set-up):
68 TFLUX total heat flux (match heat-content variations), >0 increases theta
69 SFLUX total salt flux (match salt-content variations), >0 increases salt
70
71 -----------------------------------------------
72 III) New diagnostics have been added:
73
74 name: description (oceanic set-up):
75 atmPload Atmospheric pressure loading
76 sIceLoad sea-ice loading (in Mass of ice+snow / area unit)
77 oceQnet net surface heat flux into the ocean (+=down), >0 increases theta
78 oceFWflx net surface Fresh-Water flux into the ocean (+=down), >0 decreases salinity
79 oceSflux net surface Salt flux into the ocean (+=down), >0 increases salinity
80 surForcT model surface forcing for Temperature, >0 increases theta
81 surForcS model surface forcing for Salinity, >0 increases salinity
82
83 --------------------------------------------------------------------------
84 Relation between surForcT,surForcS and others surface forcing diagnostics:
85 [x] = average of model variable "x" over the diagnostic time interval
86
87 a) if useRealFreshWaterFlux=F or (nonlinFreeSurf=0 & usingZCoords=T)
88 surForcT = oceQnet + TRELAX - oceQsw
89 surForcS = oceSflux + SRELAX - [PmEpR*So]
90 (with So = local Sea-Surface Salinity (SSS) if convertFW2Salt=-1
91 and So = convertFW2Salt otherwise)
92 oceFWflx = [PmEpR]
93 TFLUX = surForcT + oceQsw + oceFreez
94 SFLUX = surForcS
95
96 b) if useRealFreshWaterFlux=T & (nonlinFreeSurf>0 or usingPCoords=T),
97 In general:
98 surForcT = oceQnet + TRELAX - oceQsw + [T_dilution_effect]*Cp
99 surForcS = oceSflux + SRELAX + [S_dilution_effect]
100 where T_dilution_effect = PmEpR*( temp_EvPrRn - SST )
101 and S_dilution_effect = PmEpR*( salt_EvPrRn - SSS )
102 oceFWflx = [PmEpR]
103 TFLUX = surForcT + oceQsw + oceFreez + [PmEpR*SST]*Cp
104 SFLUX = surForcS + [PmEpR*SSS]
105
106 And with the default value: salt_EvPrRn=0. & temp_EvPrRn=UNSET_RL
107 (=> no dilution effect on Temp.):
108 surForcT = oceQnet + TRELAX - oceQsw
109 surForcS = oceSflux + SRELAX - [PmEpR*SSS]
110
111 Notes:
112 1) here PmEpR is assumed to be the fresh-water mass flux
113 per surface area [units: kg/m^2/s] whereas the model variable
114 EmPmR is still a volume flux per surface area [units: m/s].
115 2) with Linear Free surface (nonlinFreeSurf=0), the term
116 corresponding to w_surface*SST,SSS is missing in TFLUX,SFLUX
117 [might be added later ?] to match exactly the Heat and Salt
118 budget evolution.
119 --------------------------------------------------------------------------

  ViewVC Help
Powered by ViewVC 1.1.22