CCLM2 Forcing height is below canopy height for pft index – in #14: CCLM-CLM
in #14: CCLM-CLM
Cookies disclaimer
Our site saves small pieces of text information (cookies) on your
device in order to verify your login. These cookies are essential
to provide access to resources on this website and it will not
work properly without.
Learn more
<p>
Hi,
</p>
<p>
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?
<br/>
The whole error file is attached to this email.
</p>
<p>
Best regards,
<br/>
Merja
</p>
<p>
Hi,
</p>
<p>
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?
<br/>
The whole error file is attached to this email.
</p>
<p>
Best regards,
<br/>
Merja
</p>
CCLM2 Forcing height is below canopy height for pft index
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.
<p>
Hi Merja,
</p>
<p>
It looks like you have a problem when creating the regional input file for CLM3.5. This doesn’t look good:
</p>
<p>
—————-
<br/>
0-127: A file or directory in the path name does not exist. /scratch/g/g260039/cc2cERAIntInd/clmglobal_surfacedata/mksrf_lanwat.060929.nc
<br/>
0-127: X/Y Dimensions of global input data (soil texture): 4320 2160
<br/>
0-127: X/Y Dimensions of global input data (all other) : 720 360
<br/>
16: NetCDF: Not a valid ID
<br/>
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
<br/>
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
<br/>
94: Creating CLM3 type NetCDF file from SiB boundary data
<br/>
0-93 95-127: Creating CLM3 type NetCDF file from SiB boundary data
<br/>
—————
</p>
<p>
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.
</p>
<p>
cheers
<br/>
Edouard
</p>
<p>
Hi Merja,
</p>
<p>
It looks like you have a problem when creating the regional input file for CLM3.5. This doesn’t look good:
</p>
<p>
—————-
<br/>
0-127: A file or directory in the path name does not exist. /scratch/g/g260039/cc2cERAIntInd/clmglobal_surfacedata/mksrf_lanwat.060929.nc
<br/>
0-127: X/Y Dimensions of global input data (soil texture): 4320 2160
<br/>
0-127: X/Y Dimensions of global input data (all other) : 720 360
<br/>
16: NetCDF: Not a valid ID
<br/>
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
<br/>
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
<br/>
94: Creating CLM3 type NetCDF file from SiB boundary data
<br/>
0-93 95-127: Creating CLM3 type NetCDF file from SiB boundary data
<br/>
—————
</p>
<p>
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.
</p>
<p>
cheers
<br/>
Edouard
</p>
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.
<p>
Hi Edouard,
<br/>
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
<span class="caps">
PFTDATA
</span>
_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
<span class="caps">
PFTDATA
</span>
_MASK are both over SE Asia but the
<span class="caps">
PFTDATA
</span>
_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?
<br/>
Cheers,
<br/>
Merja
</p>
<p>
Hi Edouard,
<br/>
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
<span class="caps">
PFTDATA
</span>
_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
<span class="caps">
PFTDATA
</span>
_MASK are both over SE Asia but the
<span class="caps">
PFTDATA
</span>
_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?
<br/>
Cheers,
<br/>
Merja
</p>
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
<p>
Hi Merja,
</p>
<p>
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
</p>
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,
<p>
These values are calculated with a small fortran program that I got from Hans-Jürgen (
<span class="caps">
CLM
</span>
_Schichtdicke.f). So in that case you would have to re-run int2lm.
</p>
<p>
Cheers,
<br/>
Wim
</p>
<p>
Hi Merja,
</p>
<p>
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
</p>
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,
<p>
These values are calculated with a small fortran program that I got from Hans-Jürgen (
<span class="caps">
CLM
</span>
_Schichtdicke.f). So in that case you would have to re-run int2lm.
</p>
<p>
Cheers,
<br/>
Wim
</p>
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
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.
CCLM2 Forcing height is below canopy height for pft index
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
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
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
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
Thanks, Wim, I will try that. What value did you set for vcflat in your case?
Cheers,
Merja