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

SP mode with custom LAI

Status
Not open for further replies.

czarakas@uw_edu

New Member
Hi -

I'm interested in running some experiments with CLM in SP mode, but where I modify the satellite phenology to be consistent with a BGC run (i.e. instead of LAI being prescribed to match satellite-derived LAI, LAI is prescribed to be the LAI from a different model experiment). Could you point me to the file for the default SP driver file is, so that I could model an alternative SP file based on that? And should I be able to construct alternative SP driver files from saved monthly output for PFT-specific LAI, stem area index, and vegetation height, or will I need to save any other variables (or these variables at a different time frequency)?

Thanks in advance for your help!
Claire
 

oleson

Keith Oleson
CSEG and Liaisons
Staff member
Yes, you should be able to do that. The LAI, SAI, etc data is provided by the surface dataset. The surface dataset is identified by the namelist parameter "fsurdat" which you can find in your lnd_in file after you setup a case.
For example, for CESM2.1.3, for a 2 degree year 2000 SP case, the surface dataset located on cheyenne is:

/glade/p/cesmdata/cseg/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/surfdata_1.9x2.5_hist_16pfts_Irrig_CMIP6_simyr2000_c190304.nc

The fields you mention are:

double MONTHLY_LAI(time, lsmpft, lsmlat, lsmlon) ;
MONTHLY_LAI:long_name = "monthly leaf area index" ;
MONTHLY_LAI:units = "unitless" ;
double MONTHLY_SAI(time, lsmpft, lsmlat, lsmlon) ;
MONTHLY_SAI:long_name = "monthly stem area index" ;
MONTHLY_SAI:units = "unitless" ;
double MONTHLY_HEIGHT_TOP(time, lsmpft, lsmlat, lsmlon) ;
MONTHLY_HEIGHT_TOP:long_name = "monthly height top" ;
MONTHLY_HEIGHT_TOP:units = "meters" ;
double MONTHLY_HEIGHT_BOT(time, lsmpft, lsmlat, lsmlon) ;
MONTHLY_HEIGHT_BOT:long_name = "monthly height bottom" ;
MONTHLY_HEIGHT_BOT:units = "meters" ;

These are dimensioned by time (monthly average), lsmpft (plant functional type), and latitude/longitude.
 
Status
Not open for further replies.
Top