Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

Single-column CLM with dynamic vegetation and crop model

jaedynw

Jaedyn Williams
New Member
What I am trying to do:
I am looking at single column CAM+CLM (SCAM) and sensitivity to surface cover. I will be running the model for a year at a time. I want the crop model to be active with predictive phenology. This link shows the forcing data that I am using. The ARM data is already processed into forcing files for SCAM.

Next steps I want to do:
  1. Get dynamic vegetation (predictive phenology) working with the single column model to have BGC running to predict GPP, LAI, etc. From the CESM forum and wiki page, "CLM_BLDNML_OPTS" has a value of "-bgc bgc" in the env_run.xml to turn on the dynamic vegetation.​
  2. Look at land use (monocultures versus small family-owned farms) with different crops which require having the crop model running. Are the steps in this Wiki page how I should activate the crop model? It says I can choose a compset with "crop" in the name or by adding -crop to CLM_BLDNML_OPTS.
What version of the code are you using?
cesm2_3_beta17

Describe every step you took leading up to the problem:
This is part of my build script for running SCAM with a single IOP:
set CASETITLE=dyn_crop_test_case
set CASESET=FSCAM
set CASERES=T42_T42
set COMPILER=intel
set CESMROOT=/glade/work/$USER/cesm2_3_beta17
set IOP = arm97
set IOPNAME = scam_$IOP

## create new case
cd $CESMROOT/cime/scripts

#Create full casename
set CASENAME=${CASESET}.${CASERES}.${CASETITLE}

./create_newcase --case /glade/work/$USER/cases/$CASENAME --res $CASERES --compset $CASESET --compiler $COMPILER --driver nuopc --project $PROJNAME --user-mods-dir ${MODSDIR}/${IOPNAME} --run-unsupported

cd /glade/work/$USER/cases/$CASENAME

# Various XMLchange options to control settings of SCAM
./xmlchange ATM_NCPL=288 # Timesteps per day (given 5m dt)
./xmlchange STOP_OPTION=ndays
./xmlchange STOP_N=30
./xmlchange RUN_STARTDATE="2012-05-01"
./xmlchange JOB_WALLCLOCK_TIME=01:00:00
./xmlchange START_TOD=18060
./xmlchange CLM_BLDNML_OPTS="-bgc bgc -crop -ignore_ic_date -ignore_warnings"

./case.setup

I have the compset set to FSCAM, but do I need to use a specific compset to turn on the crop model?

This shows the use_nl_clm from my build:

Errors that show up after submitting the case:


When I check med.log.4869360.desched1.260131-114459 there aren't any errors that show up.

Describe your problem or question:
- Does my setup look correct in order to have the crop model running and dynamic vegetation?
- Do I need to modify assumptions about local conditions on the surface dataset and parameter files?
- Do I need to have datm files?
 

jaedynw

Jaedyn Williams
New Member
Sorry, the images got removed.
user_nl_clm from my build:
1770070901122.png

Errors that show up after submitting the case:
1770070956050.png
 

slevis

Moderator
Staff member
@jaedynw your steps seem logical and correct to me, and I am not familiar with the specific error messages, nor with running SCAM. As sanity checks, I would suggest starting a few cases where you change nothing, to see that you can get some simple simulations started. For example, an "out-of-the-box" SCAM simulation with no changes, to see if that works for you. Also an out-of-the-box BgcCrop I-case (land only) to see that you can also run that, and maybe an out-of-the-box F2000. Then you can use the "diff" command to difference the contents of the case directories that worked and the ones that didn't, in case this helps with the troubleshooting.
 

slevis

Moderator
Staff member
If the above approach does not get you up and running eventually, then I would suggest posting in the CAM forum, where you are likely to find people more experienced with SCAM.
 

oleson

Keith Oleson
CSEG and Liaisons
Staff member
I see this in your cesm log:

dec2027.hsn.de.hpc.ucar.edu 0: ERROR: hycoef_read: ERROR: file lev does not match model.
dec2027.hsn.de.hpc.ucar.edu 0: Image PC Routine Line Source
dec2027.hsn.de.hpc.ucar.edu 0: cesm.exe 0000000002C0862D shr_abort_mod_mp_ 114 shr_abort_mod.F90
dec2027.hsn.de.hpc.ucar.edu 0: cesm.exe 00000000006597D9 hycoef_mp_hycoef_ 334 hycoef.F90
dec2027.hsn.de.hpc.ucar.edu 0: cesm.exe 000000000064ECA6 dyn_grid_mp_dyn_g 125 dyn_grid.F90
dec2027.hsn.de.hpc.ucar.edu 0: cesm.exe 00000000005A622B cam_comp_mp_cam_i 165 cam_comp.F90
dec2027.hsn.de.hpc.ucar.edu 0: cesm.exe 000000000058ECAA atm_comp_nuopc_mp 635 atm_comp_nuopc.F90

which evidently has something to do with CAM.

And this in the atm log:

hycoef_read: ERROR: file lev does not match model. lev (file, model):
64 32
ERROR: hycoef_read: ERROR: file lev does not match model.

I interpret this as meaning a file you are reading in has a different number of levels (64) than the model expects (32).
As @slevis says, you will likely get some better help by posting on the CAM Forum.
 
Top