Skip to content

Commit

Permalink
Added S2S LIS forecast files.
Browse files Browse the repository at this point in the history
  • Loading branch information
emkemp committed Sep 30, 2021
1 parent bc4d8d0 commit 272bc1b
Show file tree
Hide file tree
Showing 25 changed files with 3,343 additions and 0 deletions.
Original file line number Diff line number Diff line change
@@ -0,0 +1,9 @@
#nfields
3
#varmin varmax
Incident Shortwave Radiation Level 001
0.0 1600.0
Incident Longwave Radiation Level 001
-500.0 1000.0
Rainfall Rate Level 001
0.0 1.0
Original file line number Diff line number Diff line change
@@ -0,0 +1,7 @@
#ptype std std_max zeromean tcorr xcorr ycorr ccorr
Incident Shortwave Radiation Level 001
1 0.30 2.5 1 86400 0 0 1.0 -0.5 -0.8
Incident Longwave Radiation Level 001
0 50.0 2.5 1 86400 0 0 -0.5 1.0 0.5
Rainfall Rate Level 001
1 0.50 2.5 1 86400 0 0 -0.8 0.5 1.0
Original file line number Diff line number Diff line change
@@ -0,0 +1,11 @@
#nfields
4
#name varmin varmax
Soil Moisture Layer 1
0.01 0.55
Soil Moisture Layer 2
0.01 0.55
Soil Moisture Layer 3
0.01 0.55
Soil Moisture Layer 4
0.01 0.55
Original file line number Diff line number Diff line change
@@ -0,0 +1,9 @@
#perttype std std_max zeromean tcorr xcorr ycorr ccorr
Soil Moisture Layer 1
0 0.004 2.5 1 43200 0 0 1.0 0.0 0.0 0.0
Soil Moisture Layer 2
0 0.00 2.5 1 43200 0 0 0.0 1.0 0.0 0.0
Soil Moisture Layer 3
0 0.00 2.5 1 43200 0 0 0.0 0.0 1.0 0.0
Soil Moisture Layer 4
0 0.00 2.5 1 43200 0 0 0.0 0.0 0.0 1.0
Original file line number Diff line number Diff line change
@@ -0,0 +1,7 @@
# number of state variables
2
#varmin varmax
SWE
0.0 10000.0
Snowdepth
0.0 100.0
Original file line number Diff line number Diff line change
@@ -0,0 +1,6 @@
#perttype std std_max zeromean tcorr xcorr ycorr ccorr
SWE
1 0.01 2.5 1 10800 0 0 1.0 0.9
Snowdepth
1 0.01 2.5 1 10800 0 0 0.9 1.0

Original file line number Diff line number Diff line change
@@ -0,0 +1,5 @@
#nfields
1
#name varmin varmax
SMAP soil moisture
0.01 1.0
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
#perttype std std_max zeromean tcorr xcorr ycorr ccorr
SMAP soil moisture
0 0.02 2.5 1 43200 0 0 1.0
Original file line number Diff line number Diff line change
@@ -0,0 +1,5 @@
#nfields
1
#name varmin varmax
SMOPS soil moisture
0.01 1.0
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
#perttype std std_max zeromean tcorr xcorr ycorr ccorr
SMOPS soil moisture
0 0.02 2.5 1 43200 0 0 1.0
Original file line number Diff line number Diff line change
@@ -0,0 +1,5 @@
#nfields
1
#varmin varmax
SNODEP obs
0. 10000.0
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
#perttype std std_max zeromean tcorr xcorr ycorr ccorr
SNODEP obs
0 0.05 2.5 1 10800 0 0 1
Original file line number Diff line number Diff line change
@@ -0,0 +1,5 @@
#nfields
1
#varmin varmax
USAFSI obs
0. 10000.0
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
#perttype std std_max zeromean tcorr xcorr ycorr ccorr
USAFSI obs
0 0.05 2.5 1 10800 0 0 1

Large diffs are not rendered by default.

Original file line number Diff line number Diff line change
@@ -0,0 +1,249 @@
#Overall driver options
Running mode: "retrospective"
Map projection of the LIS domain: latlon
Number of nests: 1
Number of surface model types: 2
Surface model types: "LSM" "Openwater"
Land surface model: "Noah-MP.4.0.1"
Open water model: "template open water"
Surface model output interval: 1da

Blending method for forcings: overlay
Number of met forcing sources: 2
# Generic Ensemble Forecast readers:
Met forcing sources: "GenEnsFcst" "PPTEnsFcst"
Temporal interpolation method (met forcing): "linear" "linear"
Spatial interpolation method (met forcing): "neighbor" "neighbor" # 25KM_res grid
Spatial upscaling method (met forcing): "none" "none"
Enable spatial downscaling of precipitation: 0 0

Topographic correction method (met forcing): "slope-aspect" "none"

# Number of hindcast ensemble members:
Number of ensembles per tile: 12 # To be modified, depending on forecast model

# Initial Condition (IC) mode for the LSM:
Start mode: restart

#Runtime options
Forcing variables list file: ./forcing_variables.txt
Output methodology: "2d ensemble gridspace"
Output model restart files: 1
Output data format: netcdf
Output naming style: "3 level hierarchy"
Enable output statistics: .false.

Starting year: 2019
Starting month: 10 # Start of 9-month hindcast
Starting day: 1
Starting hour: 0
Starting minute: 0
Starting second: 0
Ending year: 2020
Ending month: 7 # End of 9-month hindcast
Ending day: 1
Ending hour: 0
Ending minute: 0
Ending second: 0
Undefined value: -9999
Output directory: ./Oct1/
Diagnostic output file: ./Oct1/201910/logs/lislog

#The following options are used for subgrid tiling based on vegetation
Maximum number of surface type tiles per grid: 1
Minimum cutoff percentage (surface type tiles): 0.05
Maximum number of soil texture tiles per grid: 1
Minimum cutoff percentage (soil texture tiles): 0.05
Maximum number of soil fraction tiles per grid: 1
Minimum cutoff percentage (soil fraction tiles): 0.05
Maximum number of elevation bands per grid: 1
Minimum cutoff percentage (elevation bands): 0.05
Maximum number of slope bands per grid: 1
Minimum cutoff percentage (slope bands): 0.05
Maximum number of aspect bands per grid: 1
Minimum cutoff percentage (aspect bands): 0.05

#Processor layout
# AFRICOM/NHyFAS specs (please do not change):
Number of processors along x: 1
Number of processors along y: 28
Halo size along x: 0
Halo size along y: 0

#Sub-models
Radiative transfer model: none
Number of application models: 0

#---------------------DATA ASSIMILATION ----------------------------------

Number of data assimilation instances: 0
Data assimilation algorithm: none
Data assimilation set: none
Data assimilation exclude analysis increments: 0
Data assimilation output interval for diagnostics: 1da
Data assimilation number of observation types: 0
Data assimilation output ensemble members: 0
Data assimilation output processed observations: 0
Data assimilation output innovations: 0
Data assimilation use a trained forward model: 0
Data assimilation trained forward model output file: none

#Bias estimation options
Bias estimation algorithm: none
Bias estimation attributes file: none
Bias estimation restart output frequency: 1da
Bias estimation start mode: none
Bias estimation restart file: none

#Perturbation options
Perturbations start mode: coldstart
Perturbations restart output interval: 1da
Perturbations restart filename: none
Forcing perturbation algorithm: none
Forcing perturbation frequency: 1hr
Forcing attributes file: none
Forcing perturbation attributes file: none
State perturbation algorithm: none
State perturbation frequency: 1hr
State attributes file: none
State perturbation attributes file: none
Observation perturbation algorithm: none
Observation perturbation frequency: 1hr
Observation attributes file: none
Observation perturbation attributes file: none

#------------------------DOMAIN SPECIFICATION--------------------------
#The following options list the choice of parameter maps to be used

LIS domain and parameter data file: ./lis_input.s2s_africom.noahmp401_hymap2.25km.nc

Landmask data source: LDT
Landcover data source: LDT
Soil texture data source: LDT
Soil fraction data source: none
Soil color data source: none
Elevation data source: LDT
Slope data source: LDT
Aspect data source: LDT
Curvature data source: none
LAI data source: none
SAI data source: none
Albedo data source: LDT
Max snow albedo data source: LDT
Greenness data source: LDT
Roughness data source: none
Porosity data source: none
Ksat data source: none
B parameter data source: none
Quartz data source: none
Emissivity data source: none

TBOT lag skin temperature update option: 0
TBOT skin temperature lag days: 0

#--------------------------------FORCINGS----------------------------------

# Ensemble Forecast (user-generated):
Generic ensemble forecast directory: ./CFSv2/6-hourly/
Generic ensemble forecast number of ensemble members: 12

Precipitation ensemble forecast directory: ./CFSv2/6-hourly
Precipitation ensemble forecast number of ensemble members: 12

# Example:
#Generic ensemble forecast directory: ./BCSD/6-hourly/1982/may01/ens#/GEOS5.all_forc_198205.nc4

#-----------------------LAND SURFACE MODELS--------------------------
# Noah-MP.4.0.1 LSM:
Noah-MP.4.0.1 model timestep: 15mn

# Restart ICs for hindcast/forecast runs:
Noah-MP.4.0.1 restart file: ./DA_ICs/CFSv2/LIS_RST_NOAHMP401_201909302345.d01.nc
Noah-MP.4.0.1 restart output interval: 1mo
Noah-MP.4.0.1 restart file format: netcdf

Noah-MP.4.0.1 soil parameter table: ./noahmp401_parms/SOILPARM.TBL
Noah-MP.4.0.1 general parameter table: ./noahmp401_parms/GENPARM.TBL
Noah-MP.4.0.1 MP parameter table: ./noahmp401_parms/MPTABLE.TBL
Noah-MP.4.0.1 number of soil layers: 4
Noah-MP.4.0.1 thickness of soil layers: 0.1 0.3 0.6 1.0
Noah-MP.4.0.1 dynamic vegetation option: 4 # 9 options available
Noah-MP.4.0.1 canopy stomatal resistance option: 1 # 1=Ball-Berry; 2=Jarvis
Noah-MP.4.0.1 soil moisture factor for stomatal resistance: 1 # 1=Noah; 2=CLM; 3=SSiB
Noah-MP.4.0.1 runoff and groundwater option: 1 # 1=SIMGM; 2=SIMTOP; 3=Schaake96; 4=BATS
Noah-MP.4.0.1 surface layer drag coefficient option: 1 # 1=M-O; 2=Chen97
Noah-MP.4.0.1 supercooled liquid water option: 1 # 1=NY06; 2=Koren99
Noah-MP.4.0.1 frozen soil permeability option: 1 # 1=NY06; 2=Koren99
Noah-MP.4.0.1 radiation transfer option: 3 # 1=gap=F(3D;cosz); 2=gap=0; 3=gap=1-Fveg
Noah-MP.4.0.1 snow surface albedo option: 2 # 1=BATS; 2=CLASS
Noah-MP.4.0.1 rainfall & snowfall option: 1 # 1=Jordan91; 2=BATS; 3=Noah
Noah-MP.4.0.1 lower boundary of soil temperature option: 2 # 1=zero-flux; 2=Noah
Noah-MP.4.0.1 snow&soil temperature time scheme option: 1 # 1=semi-implicit; 2=fully implicit
Noah-MP.4.0.1 glacier option: 1 # 1=include phase change; 2=slab ice (Noah)
Noah-MP.4.0.1 surface resistance option: 1 # 1=Sakaguchi and Zeng 2009; 2=Sellers (1992); 3=adjusted Sellers; 4=option1 for non-snow and rsurf_snow for snow
Noah-MP.4.0.1 soil configuration option: 1 # 1=input dominant soil texture; 2=input soil texture varies that varies with depth; 3=soil composition and pedotransfer functions
Noah-MP.4.0.1 soil pedotransfer function option: 1 # 1=Saxton and Rawls (2006) (used when soil_opt=3)
Noah-MP.4.0.1 crop model option: 0 # 0=No crop model; 1=Liu et al. 2016; 2=Gecros
Noah-MP.4.0.1 urban physics option: 0 # 0=No; 1=Single-layer; 2=Multi-layer BEP scheme; 3=Multi-layer BEM scheme
Noah-MP.4.0.1 initial surface skin temperature: 288.0
Noah-MP.4.0.1 initial soil temperatures: 288.0 288.0 288.0 288.0
Noah-MP.4.0.1 initial total soil moistures: 0.20 0.20 0.20 0.20
Noah-MP.4.0.1 initial snow water equivalent: 0.0
Noah-MP.4.0.1 initial snow depth: 0.0
Noah-MP.4.0.1 initial total canopy surface water: 0.0
Noah-MP.4.0.1 initial leaf area index: 0.5
Noah-MP.4.0.1 initial water table depth: 2.5
Noah-MP.4.0.1 initial water in the aquifer: 4900.0
Noah-MP.4.0.1 initial water in aquifer and saturated soil: 4900.0
Noah-MP.4.0.1 reference height of temperature and humidity: 10.0

Template open water timestep: 15mn


#----------------------- ROUTING --------------------------
# ROUTING MODEL: HYMAP
Routing model: "HYMAP2 router"

# Restart ICs for hindcast/forecast runs:
HYMAP2 routing model start mode: restart
HYMAP2 routing model restart file: ./DA_ICs/CFSv2/LIS_RST_HYMAP2_router_201909302345.d01.nc
HYMAP2 routing model restart interval: "1mo"

HYMAP2 routing model time step: "15mn"
HYMAP2 routing model output interval: "1da"
HYMAP2 enable 2-way coupling: 0
#HYMAP2 run in ensemble mode: 1 # Run HYMAP in ensemble mode
HYMAP2 run in ensemble mode: 0 # Run HYMAP as single member

HYMAP2 routing method: "kinematic" #"local inertia"
HYMAP2 routing model time step method: "constant" # "adaptive"
HYMAP2 routing model adaptive time step alfa coefficient: 1.0
HYMAP2 floodplain dynamics: 1
HYMAP2 routing model dwi flag: 0

HYMAP2 routing model evaporation option: none # "penman"
HYMAP2 routing model linear reservoir flag: 0
HYMAP2 reservoir operation option: 0
HYMAP2 number of reservoirs: 1
HYMAP2 reservoir operation input time series size: 2
HYMAP2 reservoir operation input directory: ./
HYMAP2 reservoir operation header filename: ./header_test.txt
HYMAP2 reservoir operation input data type: "water level" # "streamflow" #

TEMPLATE model timestep: "15mn"


#---------------------------MODEL OUTPUT CONFIGURATION-----------------------
#Specify the list of ALMA variables that need to be featured in the
#LSM model output

Model output attributes file: ./tables/MODEL_OUTPUT_LIST.TBL.noahmp401_hymap2

# Start Date/Time when output is to be started:
Output start year:
Output start month:
Output start day:
Output start hour:
Output start minutes:
Output start seconds:

29 changes: 29 additions & 0 deletions lis/configs/557WW-7.4-FOC/S2S_AFRICOM/forcing_variables.txt
Original file line number Diff line number Diff line change
@@ -0,0 +1,29 @@
#ALMA Name select vlevels units
Tair: 1 1 K # Near surface air temperature
Qair: 1 1 kg/kg # Near surface specific humidity
SWdown: 1 1 W/m2 # Incident shortwave radiation (total)
SWdirect: 0 1 W/m2 # Incident shortwave radiation (direct)
SWdiffuse: 0 1 W/m2 # Incident shortwave radiation (diffuse)
LWdown: 1 1 W/m2 # Incident longwave radiation
Wind_E: 1 1 m/s # Eastward wind
Wind_N: 1 1 m/s # Northward wind
Psurf: 1 1 Pa # Surface pressure
Rainf: 1 1 kg/m2s # Rainfall rate
Snowf: 0 1 kg/m2s # Snowfall rate
CRainf: 0 1 kg/m2s # Convective rainfall rate
Forc_Hgt: 0 1 m # Height of forcing variables
Ch: 0 1 - # Surface exchange coefficient for heat
Cm: 0 1 - # Surface exchange coefficient for momentum
Q2sat: 0 1 - # Saturated mixing ratio
Emiss: 0 1 - # Surface emissivity
Cosz: 0 1 - # Cosine of zenith angle
Albedo: 0 1 - # Surface albedo
PARDR: 0 1 - # Photosynthetically Active Radiation (direct)
PARDF: 0 1 - # Photosynthetically Active Radiation (diffuse)
SWnet: 0 1 - # Net shortwave radiation at the surface
PET: 0 1 kg/m2 # Potential ET
RefET: 0 1 kg/m2s # Reference ET
CAPE: 0 1 J/kg # Convective Available Potential Energy
LPressure: 0 1 Pa # Level pressure
O3: 0 1 - # Ozone concentration
Xice: 0 1 - # Sea ice mask
Loading

0 comments on commit 272bc1b

Please sign in to comment.