[MITgcm-devel] changing some diagnostics names

[Jean-Michel Campin]

Dimitris,

I propose to send the attached text to the support list
(with some warning at the beginning if someone try
to use with the future version of the code, an old 
data.diagnostics, that contains some old diagnostics names).
Does this answer your question ?

Jean-Michel

On Tue, Nov 14, 2006 at 10:17:36AM -0800, Dimitris Menemenlis wrote:
> JM, OK with all your suggestions.
> 
> What will be really important though, is a small description
> of what these various diagnostics are and how they are related to each 
> other under differing forcing scenarios.
> 
> Right now I do not have a clear idea of what differences between oceQnet, 
> surForcT, and TFLUX will be.  When will they be identical (to within a 
> scaling factor) when will they differ?
> 
> Same question for oceSflux, SFLUX, surForcS, oceFWflx?
> 
> D.
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
-------------- next part --------------
In preparation for new diagnostics to come, and to make
it easier to identify which diagnostics to use:

I) Some diagnostics have been renamed (essentially, to better
match the content):

1) name:	description (oceanic set-up):
PRESSURE	Cell-Center Height
replaced by:
RCENTER		Cell-Center Height

2) name:	description (oceanic set-up):
TICE		heat from melt/freeze of sea-ice, >0 increases theta
replaced by:
oceFreez	heating from freezing of sea-water (allowFreezing=T)

3) name:	description (oceanic set-up):
TAUX		zonal surface wind stress, >0 increases uVel
TAUY		meridional surf. wind stress, >0 increases vVel
replaced by:
oceTAUX		zonal surface wind stress, >0 increases uVel
oceTAUY		meridional surf. wind stress, >0 increases vVel

4) name:	description (oceanic set-up):
SWFLUX		net upward SW radiation, >0 increases theta
replaced by:
oceQsw		net Short-Wave radiation (+=down), >0 increases theta

5) name:	description (oceanic set-up):
DIFx_TH		Zonal      Diffusive Flux of Pot.Temperature
DIFy_TH		Meridional Diffusive Flux of Pot.Temperature
replaced by:
DFxE_TH		Zonal      Diffusive Flux of Pot.Temperature
DFyE_TH		Meridional Diffusive Flux of Pot.Temperature

6) name:	description (oceanic set-up):
DIFx_SLT	Zonal      Diffusive Flux of Salinity
DIFy_SLT	Meridional Diffusive Flux of Salinity
replaced by:
DFxE_SLT	Zonal      Diffusive Flux of Salinity
DFyE_SLT	Meridional Diffusive Flux of Salinity

-----------------------------------------------
II) Change description of existing diagnostics:

name:		old description (oceanic set-up):
TFLUX   	net surface heat flux, >0 increases theta
SFLUX   	net surface salt flux, >0 increases salt
name:		new description (oceanic set-up):
TFLUX   	total heat flux (match heat-content variations), >0 increases theta
SFLUX   	total salt flux (match salt-content variations), >0 increases salt

-----------------------------------------------
III) New diagnostics have been added:

name:		description (oceanic set-up):
oceQnet		net surface heat flux into the ocean (+=down), >0 increases theta
oceFWflx	net surface Fresh-Water flux into the ocean (+=down), >0 decreases salinity
oceSflux 	net surface Salt flux into the ocean (+=down), >0 increases salinity
surForcT	model surface forcing for Temperature, >0 increases theta
surForcS	model surface forcing for Salinity, >0 increases salinity

--------------------------------------------------------------------------
Relation between surForcT,surForcS and others surface forcing diagnostics:
 [x] = average of model variable "x" over the diagnostic time interval

a) if useRealFreshWaterFlux=F or (nonlinFreeSurf=0 & usingZCoords=T)
 surForcT = oceQnet  + TRELAX - oceQsw
 surForcS = oceSflux + SRELAX - [PmEpR*So]
  (with So = local Sea-Surface Salinity (SSS) if convertFW2Salt=-1
    and So = convertFW2Salt otherwise) 
 oceFWflx = [PmEpR]
 TFLUX = surForcT + oceQsw + oceFreez
 SFLUX = surForcS

b) if useRealFreshWaterFlux=T & (nonlinFreeSurf>0 or usingPCoords=T):
with the default value: temp_EvPrRn=UNSET_RL & salt_EvPrRn=0.:
 surForcT = oceQnet  + TRELAX - oceQsw
 surForcS = oceSflux + SRELAX - [PmEpR*SSS]
 oceFWflx = [PmEpR]
 TFLUX = surForcT + oceQsw + oceFreez + [PmEpR*SST]
 SFLUX = surForcS + [PmEpR*SSS]

Notes:
1) here PmEpR is assumed to be the fresh-water mass flux 
 per surface area [units: kg/m^2/s] whereas the model variable
 EmPmR is still a volume flux per surface area [units: m/s].
2) with Linear Free surface (nonlinFreeSurf=0), the term 
corresponding to w_surface*SST,SSS is missing in TFLUX,SFLUX 
[might be added later ?] to match exactly the Heat and Salt 
budget evolution.
--------------------------------------------------------------------------