ccdfheatcm : computes the heat content in the specified area (mono-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 [T-file]
  • a dataset with 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

>>> # How to get required files for cdftools cdfheatc binary
>>> fixed_fields('ccdfheatcm',
 ... ('',    '/data/climaf/${project}/${model}/'),
 ... ('','/data/climaf/${project}/${model}/'),
 ... ('','/data/climaf/${project}/${model}/'))
>>> dso=ds(simulation="PRE6CPLCr2alb", variable="so", period="199807", realm="O") # dataset with salinity
>>> dtho=ds(simulation="PRE6CPLCr2alb", variable="thetao", period="199807", realm="O") # dataset with temperature
>>> my_cdfheatc=ccdfheatcm(dso,dtho,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.