Issues of Creating landuse.timeseries file by mksurfdata.pl

[Dong Peng]

Hello everyone. I was trying to use mksurfdata.pl to create landuse.timeseries file with the option -no_surfdata. Which was not working. Here is the history of the command line:
perl mksurfdata.pl -r 360x720cru -no_surfdata
ERROR: Tried to run mksurfdata_map without creating either a surface dataset or a landuse.timeseries file at mksurfdata.pl line 789.

But I can successfully create the surfdata without the option -no_surfdata. What should I do if I want a single landuse.timeseries file.

Thanks.

 

[oleson]

If you want a landuse timeseries file then you would need to request it, e.g.,

./mksurfdata.pl -r 360x720cru -no_surfdata -years 1850-2000

 

[Dong Peng]

If you want a landuse timeseries file then you would need to request it, e.g.,

./mksurfdata.pl -r 360x720cru -no_surfdata -years 1850-2000

Thanks, Keith. I got another question. I want to add some parameter to pftcon.f90, it requires reading in some parameters from the surfdata pft variables. How can I add the parameter to the surfdata by using mksurfdata.pl? The parameter file is located at https://web.lcrc.anl.gov/public/e3sm/inputdata/lnd/clm2/paramdata/clm_params.ctc.erosion.c220120.nc

 

[oleson]

I don't think I understand your question. The parameter file (and pftcon.f90 which reads in the parameters on the parameter file) are meant for spatially-invariant and time-invariant pft-dependent parameters. The surface dataset is for parameters that are spatially-dependent and/or time-dependent.
What exactly are you trying to do?

 

[Dong Peng]

I don't think I understand your question. The parameter file (and pftcon.f90 which reads in the parameters on the parameter file) are meant for spatially-invariant and time-invariant pft-dependent parameters. The surface dataset is for parameters that are spatially-dependent and/or time-dependent.
What exactly are you trying to do?

Hi Keith, sorry for the confusion. I try to make up the pft-specific parameters. Here I try to read in 4 parameters from pft data in pftcon.f90. The code is:

call ncd_io('gcbc_p',this%gcbc_p, 'read', ncid, readvar=readv, posNOTonfile=.true.)
if ( .not. readv ) this%gcbc_p(:) = 0._r8
call ncd_io('gcbc_q',this%gcbc_q, 'read', ncid, readvar=readv, posNOTonfile=.true.)
if ( .not. readv ) this%gcbc_q(:) = 0._r8
call ncd_io('gcbr_p',this%gcbr_p, 'read', ncid, readvar=readv, posNOTonfile=.true.)
if ( .not. readv ) this%gcbr_p(:) = 0._r8
call ncd_io('gcbr_q',this%gcbr_q, 'read', ncid, readvar=readv, posNOTonfile=.true.)
if ( .not. readv ) this%gcbr_q(:) = 0._r8

The parameters are required to read in PFT physiological data. The question is how can I add the parameter files to the PFT physiological data? It seems the PFT physiological data are built in surfdata. Can I build the pft-specific parameters by mksurfdata.pl? Here I attached the pftcon.f90 file for you. Hope that will help you understand my question.

 

[oleson]

If you are looking for a method to add the parameter variable name and numerical values to the parameter file (e.g., clm_params.ctc.erosion.c220120.nc), then one method is shown here:

Moving a parameter to the params file

Community Terrestrial Systems Model (includes the Community Land Model of CESM) - ESCOMP/CTSM

github.com

If you want to use that parameter file in a model run, then you would point to the parameter file in your user_nl_clm in your case directory, e.g.,

paramfile = 'YOUR_PARAM_FILE'

None of this has anything to do with mksurfdata.pl Does that answer your question?

 

[Dong Peng]

If you are looking for a method to add the parameter variable name and numerical values to the parameter file (e.g., clm_params.ctc.erosion.c220120.nc), then one method is shown here:

Moving a parameter to the params file

Community Terrestrial Systems Model (includes the Community Land Model of CESM) - ESCOMP/CTSM

github.com

If you want to use that parameter file in a model run, then you would point to the parameter file in your user_nl_clm in your case directory, e.g.,

paramfile = 'YOUR_PARAM_FILE'

None of this has anything to do with mksurfdata.pl Does that answer your question?

Thanks， Keith. That's the answer I am seeking. Really appreciate!