Opened 8 years ago

Closed 8 years ago

#103 closed defect (fixed)

Error in units of an output variable with the ALMA convention

Reported by: jpolcher Owned by: aducharne
Priority: major Milestone:
Component: Anthropogenic processes Version:
Keywords: Cc:

Description

When the ALMA convention is activated soil moisture over the 11 layers of the CWRR parametrisation is written into the history files.

The units of this variable are declared as 'kg/m2' when in fact it should be 'm³/m³'.

Attached a script is proposed to convert the m³/m3 of the file into kg/m² integrated over a given thickness :

Attachments (1)

ExtractSM (2.4 KB) - added by jpolcher 8 years ago.
Script to generate a file containing the thickness of each layer and integarting then the soil moisture over a number of layers

Download all attachments as: .zip

Change History (5)

Changed 8 years ago by jpolcher

Script to generate a file containing the thickness of each layer and integarting then the soil moisture over a number of layers

comment:1 Changed 8 years ago by peylin

  • Owner changed from somebody to aducharne
  • Status changed from new to assigned

comment:2 Changed 8 years ago by aducharne

  • Status changed from assigned to accepted

comment:3 Changed 8 years ago by aducharne

In hydrol, I checked all the soil moisture variables in the outputs and found the following bugs/questions:

NOT almaoutput:

  • 'moist_c1', 'moist_c2', 'moist_c3' hold mc(:,:,jst) and should be in m3/m3 instead of %
  • 'humrel', which holds the variable humrel: the unit is '1', shouldn't it rather be '-' (as for 'SWI')

almaoutput:

  • For 'SoilMoist?', the ALMA standard requires the mean soil moisture (over the three soiltiles) in the 11 layer, in mm. The variable soilmoist in hydrol.f90 has the good dimensions, but holds the spatial mean of mc (in m3/m3) => it is the way soilmoist is computed that needs to be corrected, not the unit of SoilMoist? in intersurf.f90
  • 'SoilWet?', holding the variable soilwet calculated in hydrol_alma.f90, is declared in intersurf.f90 in kg/m2, but the actual unit is [-], as required by the ALMA standards

comment:4 Changed 8 years ago by aducharne

  • Resolution set to fixed
  • Status changed from accepted to closed

The above 4 points have been corrected in rev [1943] including for XML.

Note: See TracTickets for help on using tickets.