Project

General

Profile

CCLM2 Forcing height is below canopy height for pft index

Added by Merja Tölle over 6 years ago

Hi,

by attempting to run CCLM2 over SE Asia an error occurred that the forcing height is below the canopy height for several pft indices. Could you please suggest a way to solve this problem with me?
The whole error file is attached to this email.

Best regards,
Merja


Replies (4)

RE: CCLM2 Forcing height is below canopy height for pft index - Added by Edouard Davin over 6 years ago

Hi Merja,

It looks like you have a problem when creating the regional input file for CLM3.5. This doesn’t look good:

-----------
0-127: A file or directory in the path name does not exist. /scratch/g/g260039/cc2cERAIntInd/clmglobal_surfacedata/mksrf_lanwat.060929.nc
0-127: X/Y Dimensions of global input data (soil texture): 4320 2160
0-127: X/Y Dimensions of global input data (all other) : 720 360
16: NetCDF: Not a valid ID
1-5 7-13 15 17-24 27-37 39-51 53-63 65-67 69-70 72-84 86-87 89 92-98 101-107 109-110 112-116 118-120 124 127: NetCDF: Not a valid ID
0 6 14 25-26 38 52 64 68 71 85 88 90-91 99-100 108 111 117 121-123 125-126: NetCDF: Not a valid ID
94: Creating CLM3 type NetCDF file from SiB boundary data
0-93 95-127: Creating CLM3 type NetCDF file from SiB boundary data
----------

It might be that you are not pointing to the right global raw datasets. Or maybe you have changed something to write_surfacedata.f90? I suggest to check carefully the input file (clm3_surfacedata.nc) that is produced out of this.

cheers
Edouard

RE: CCLM2 Forcing height is below canopy height for pft index - Added by Merja Tölle over 6 years ago

Hi Edouard,
thanks that helped a step but still not solved the problem. As you said the model did not find mksrf_lanwat.060929.nc and that is why the PFTDATA_MASK in clm3_surfacedata.nc was over Africa and not Indonesia. Now, the program finds mksrf_lanwat.060929.nc and still I am getting the error that the forcing height is below the canopy height. Looking at clm3_surfacedata.nc I find that the Landmask and PFTDATA_MASK are both over SE Asia but the PFTDATA_MASK is blurred and the land points are much bigger than in Landmask. If in Landmask it is a sea point maybe that is the reason for the error then? How to solve it?
Cheers,
Merja

RE: CCLM2 Forcing height is below canopy height for pft index - Added by Wim Thiery over 6 years ago

Hi Merja,

When setting up CCLM² for East-Africa I also encountered this error. I was able to resolve it by increasing the vertical coordinate of the lowest level (vcoord_d) to 43 m or higher in the Int2lm namelist. I now use

vcoord_d =  30000.00,  29000.00,  28010.00,  27030.00,  26060.00,
25100.00, 24160.00, 23220.00, 22300.00, 21390.00,
20500.00, 19620.00, 18750.00, 17890.00, 17050.00,
16230.00, 15420.00, 14630.00, 13850.00, 13090.00,
12350.00, 11630.00, 10920.00, 10230.00, 9560.00,
8910.00, 8280.00, 7670.00, 7080.00, 6500.00,
5960.00, 5430.00, 4920.00, 4440.00, 3980.00,
3540.00, 3120.00, 2730.00, 2370.00, 2020.00,
1710.00, 1420.00, 1150.00, 910.00, 700.00,
510.00, 350.00, 220.00, 120.00, 45.00,
0.00,

These values are calculated with a small fortran program that I got from Hans-Jürgen (CLM_Schichtdicke.f). So in that case you would have to re-run int2lm.

Cheers,
Wim

RE: CCLM2 Forcing height is below canopy height for pft index - Added by Merja Tölle over 6 years ago

Thanks, Wim, I will try that. What value did you set for vcflat in your case?
Cheers,
Merja

    (1-4/4)