CCLM2 Forcing height is below canopy height for pft index – in #14: CCLM-CLM

in #14: CCLM-CLM

<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>

  @merjatölle in #c7a90e5

<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

MT wrote

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

View in channel
<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>

  @redc_migration in #7a28028

<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>
re replied

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

View in channel
<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>

  @merjatölle in #99bfea9

<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>
MT replied

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

View in channel
<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>

  @redc_migration in #66b0b88

<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>
re replied

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

View in channel
<p> Thanks, Wim, I will try that. What value did you set for vcflat in your case? <br/> Cheers, <br/> Merja </p>

  @merjatölle in #79e0c95

<p> Thanks, Wim, I will try that. What value did you set for vcflat in your case? <br/> Cheers, <br/> Merja </p>
MT replied

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

View in channel