ccdfheatc : computes the heat content in the specified area (multi-variable input file)

Computes the heat content in the specified area (Joules). A sub-domain can be specified in option. This is the wrapping around the native cdfheatc operator assuming its usage is:

cdfheatc  T-file ...
... **[imin imax jmin jmax kmin kmax]** **[-full]**

CliMAF optional arguments are the ones surrounded with ‘**’.

References :

Provider / contact : climaf at meteo dot for the wrapping

Inputs (in the order of CliMAF call):

  • a dataset with salinity and temperature [T-file]

Mandatory arguments: None

Optional arguments:

  • imin, imax, jmin, jmax, kmin, kmax : limit of a sub domain where the heat content will be calculated
  • opt : may be used to pass key -full (assume full step model output instead of default partial steps)

Required files: Files,, must be in the current directory (use function fixed_fields() for that; see example below).


  • main output : a netcdf file (variables : heatc_2D, heatc_3D) and standard output

Climaf call example:: For more examples which are systematically tested, see

>>> # Use "data_CNRM" project, where netcdf files with values of so and thetao are defined
>>> cdef("project","data_CNRM")
>>> # How to get required files for cdftools cdfheatc binary
>>> fixed_fields('ccdfheatc',
 ... ('',    '/data/climaf/${project}/${model}/'),
 ... ('','/data/climaf/${project}/${model}/'),
 ... ('','/data/climaf/${project}/${model}/'))
>>> dT=ds(simulation="PRE6CPLCr2alb", variable="so,thetao", period="199807", realm="O")
>>> my_cdfheatc=ccdfheatc(dT,imin=100,imax=102,jmin=117,jmax=118,kmin=1,kmax=2,opt='-full')
>>> cfile(my_cdfheatc) # to compute the heat content in the specified area and get the result in a netcdf file

Implementation: The operator is implemented as a binary using cdftools cdfheatc operator.