Merge pull request #904 from emkemp/fix/s2s_scripts

Fix/s2s scripts

ldt/configs/557WW-7.4-FOC/S2S_AFRICOM/README

@@ -0,0 +1,67 @@
+
+Noah-MP-4.0.1 S2S-AFRICOM LDT ensemble restart file generation script
+-------------------------------------------------------------------------
+
+The following script can be used to generate the six NMME 
+model-based ensemble restart files required to initialize 
+the same six NMME model-based forecasts using the LDT
+software for the 557WW S2S AFRICOM ~25KM domain for the 
+Noah-MP-4.0.1 LSM. Also, this script is designed to run as an
+automated script to run in part of the near-realtime (NRT)
+fashion of generating the S2S forecasts.
+
+HyMAP2 routing model restart files are single-member at this time, so
+they are copied from the LIS DA run directories during operations.
+
+Note: Initial conditions (ICs) is interchangeable with "restart file".
+
+Here is the script:
+
+  generate_ldtconfig_files_ensrst_nrt.sh
+
+and the steps to run it and process the IC files include:
+
+ 1. Need the following directory:  ./template_files
+    and the TEMPLATE-based config file within it:
+
+     template_files/ldt.config_noahmp401_nmme_TEMPLATE
+
+ 2. Inside the script, set your working directory and any 
+    other paths to where your input restart files are found,
+    for NoahMP401 and HYMAP2. 
+
+ 3. Copy your LDT executable to the local directory where you
+    will run this script.
+
+ 4. Run the script as:
+
+    sh generate_ldtconfig_files_ensrst_nrt.sh
+
+ 5. The following config and output restart directories/files 
+    are generated using the LDT executable, along with 
+    the automated date strings required to run this in NRT.
+
+   * The six NMME-ensemble based LDT configuration files:
+
+    ldt.config_noahmp401_nmme_ccm4
+    ldt.config_noahmp401_nmme_ccsm4
+    ldt.config_noahmp401_nmme_cfsv2
+    ldt.config_noahmp401_nmme_geosv2
+    ldt.config_noahmp401_nmme_gfdl
+    ldt.config_noahmp401_nmme_gnemo
+
+    written to a directory named:  ldt.config_files
+
+   * The six different LDT runs should generate six different output directories,
+    one for each NMME model, and the respective output restart files will be 
+    contained in the following list output directories, with varying numbers
+    of ensemble members based on the NMME model:
+
+     CCM4
+     CCSM4
+     CFSv2
+     GEOSv2
+     GFDL
+     GNEMO
+
+

ldt/configs/557WW-7.4-FOC/S2S_AFRICOM/for_discover/ldt.config.s2s_africom.hymap2.ensemble

@@ -0,0 +1,209 @@
+# Overall driver options
+LDT running mode:                       "Ensemble restart processing"
+Processed LSM parameter filename:       ./lis_input.s2s_africom.noahmp401_hymap2.25km.nc
+
+LIS number of nests:                    1
+Number of surface model types:          2
+Surface model types:                    "LSM"  "Openwater"
+Land surface model:                     "Noah-MP.4.0.1"
+Routing model:                          "HYMAP2"
+Lake model:                             none
+Water fraction cutoff value:            0.5
+Incorporate crop information:           .false.
+Number of met forcing sources:          0
+Met forcing sources:                    none
+LDT diagnostic file:                    ldtlog.s2safricom.hymap2.ensemble
+Mask-parameter fill diagnostic file:    MaskParamFill.s2safricom.hymap2.ensemble.log
+LDT output directory:                   OUTPUT
+Undefined value:                        -9999.0
+
+#AFRICOM 25KM domain:
+Map projection of the LIS domain:       latlon
+Run domain lower left lat:             -39.8750
+Run domain upper right lat:             39.875
+Run domain lower left lon:             -19.8750
+Run domain upper right lon:             59.875
+Run domain resolution (dx):              0.25
+Run domain resolution (dy):              0.25
+
+#Runtime options
+Number of ensembles per tile:           1
+
+#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
+
+# Ensemble restart options
+LIS restart source:                               "Routing"
+Ensemble restart generation mode:                 "upscale"
+Number of ensembles per tile (input restart):     1
+Number of ensembles per tile (output restart):    12
+
+Input restart filename:                        ./Run-D/noahmp401_2dgrid/ROUTING/202010/LIS_RST_HYMAP2_router_202010312345.d01.nc
+Output restart filename:                       ./LDT_RST_202010312345_HYMAP2_ENS12.nc
+
+
+# ________________
+
+
+#Landcover parameter inputs
+Landcover data source:                  "MODIS_Native"
+Landcover classification:               "IGBPNCEP"
+Landcover file:                         ./input/LS_PARAMETERS/noah_2dparms/igbp.bin
+Landcover spatial transform:            tile
+Landcover map projection:               latlon
+Landcover fill option:                  neighbor
+Landcover fill radius:                  5
+Landcover fill value:                   10
+
+#Landmask parameter inputs
+Create or readin landmask:              "readin"
+Landmask data source:                   "UKMO_CAP_Netcdf"
+Landmask file:                          ./data/cap2ldt_ps41.nc
+Landmask spatial transform:             mode    # Go from ~10 KM to 25KM domain ...
+Landmask map projection:                latlon
+Landmask lower left lat:               -89.9531250
+Landmask lower left lon:              -179.9296875
+Landmask upper right lat:               89.9531250
+Landmask upper right lon:              179.9296875
+Landmask resolution (dx):                0.1406250
+Landmask resolution (dy):                0.0937500
+
+#Soil parameter inputs
+Soil fraction data source:              none
+Soils spatial transform:                none
+Soils map projection:                   latlon
+Soils fill option:                      none
+Porosity data source:                   none
+Porosity map:                           none
+
+#Soil texture map:
+Soil texture data source:    ISRIC
+Soil texture map:        ./input/LS_PARAMETERS/soil_parms/ISRIC/v2017/TEXMHT_M_sl1_250m.tif     # v2017 file
+Soil texture map projection:     latlon
+Soil texture spatial transform:   mode                  # none | mode | neighbor | tile
+Soil texture fill option:       neighbor                # none | neighbor
+Soil texture fill radius:         5                     # Number of pixels to search for neighbor
+Soil texture fill value:          6                     # Static value to fill where missing
+Soil texture fill value for Antarctica:  4              # 4 -- USDA value for silty-loam
+Soil texture force exclusion of water points during fill: true
+
+#Topography parameter inputs
+Elevation data source:                  "MERIT_1K"
+Elevation number of bands:              1
+Elevation map:                          ./input/LS_PARAMETERS/topo_parms/MERIT
+Elevation fill option:                  average
+Elevation fill radius:                  5
+Elevation fill value:                   300.
+
+Slope data source:                      "MERIT_1K"
+Slope number of bands:                  1
+Slope map:                              ./input/LS_PARAMETERS/topo_parms/MERIT
+Slope fill option:                      average
+Slope fill radius:                      5
+Slope fill value:                       0
+
+Aspect data source:                     "MERIT_1K"
+Aspect number of bands:                 1
+Aspect map:                             ./input/LS_PARAMETERS/topo_parms/MERIT
+Aspect fill option:                     average
+Aspect fill radius:                     5
+Aspect fill value:                      3.14159
+
+Topography spatial transform:           average
+Topography map projection:              latlon
+
+#Albedo inputs
+Albedo data source:                     "NCEP_Native"
+Albedo map:                             ./input/LS_PARAMETERS/noah_2dparms/albedo
+Albedo climatology interval:            monthly
+Albedo spatial transform:               "budget-bilinear"
+Albedo map projection:                  latlon
+Albedo fill option:                     neighbor
+Albedo fill radius:                     5
+Albedo fill value:                      0.15
+
+#Maximum snow albedo inputs
+Max snow albedo data source:            "Barlage_Native"
+Max snow albedo map:                    ./input/LS_PARAMETERS/noah_2dparms/maximum_snow_albedo.hdf
+Max snow albedo spatial transform:      average
+Max snow albedo map projection:         latlon
+Max snow albedo fill option:            neighbor
+Max snow albedo fill radius:            5
+Max snow albedo fill value:             0.3
+
+#Greenness inputs
+Greenness data source:                  "NCEP_Native"
+Greenness fraction map:                 ./input/LS_PARAMETERS/noah_2dparms/gfrac
+Greenness climatology interval:         monthly
+Calculate min-max greenness fraction:   .false.
+Greenness maximum map:                  ./input/LS_PARAMETERS/noah_2dparms/gfrac_max.asc
+Greenness minimum map:                  ./input/LS_PARAMETERS/noah_2dparms/gfrac_min.asc
+Greenness spatial transform:            "budget-bilinear"
+Greenness map projection:               latlon
+Greenness fill option:                  neighbor
+Greenness fill radius:                  5
+Greenness fill value:                   0.3
+Greenness maximum fill value:           1.0
+Greenness minimum fill value:           0.0
+
+#Slope type inputs
+Slope type data source:         "none"
+
+#Bottom temperature inputs
+Bottom temperature data source:         "ISLSCP1"
+Bottom temperature map:                 ./input/LS_PARAMETERS/noah_2dparms/SOILTEMP.60
+Bottom temperature spatial transform:   "budget-bilinear"
+Bottom temperature map projection:      latlon
+Bottom temperature fill option:         average
+Bottom temperature fill radius:         5
+Bottom temperature fill value:          287.0
+Bottom temperature topographic downscaling:  "lapse-rate"
+
+Noah-MP PBL Height Value:               900.
+
+
+# 25KM GLDAS/Global Runs; MERIT1KM DEM
+HYMAP params spatial transform:       none
+HYMAP params map projection:         latlon
+HYMAP params lower left lat:        -59.875
+HYMAP params lower left lon:       -179.875
+HYMAP params upper right lat:        89.875
+HYMAP params upper right lon:       179.875
+HYMAP params resolution (dx):         0.25
+HYMAP params resolution (dy):         0.25
+HYMAP river width map:                ./input/LS_PARAMETERS/hymap_025/lis_rivwth.bin
+HYMAP river height map:               ./input/LS_PARAMETERS/hymap_025/lis_rivhgt.bin
+HYMAP river roughness map:            ./input/LS_PARAMETERS/hymap_025/lis_rivman.bin
+HYMAP floodplain roughness map:       ./input/LS_PARAMETERS/hymap_025/lis_fldman.bin
+HYMAP river length map:               ./input/LS_PARAMETERS/hymap_025/lis_rivlen.bin
+HYMAP floodplain height map:          ./input/LS_PARAMETERS/hymap_025/lis_fldhgt.bin
+HYMAP floodplain height levels:        10
+HYMAP flow direction x map:           ./input/LS_PARAMETERS/hymap_025/lis_nextx.bin
+HYMAP flow direction y map:           ./input/LS_PARAMETERS/hymap_025/lis_nexty.bin
+HYMAP grid elevation map:             ./input/LS_PARAMETERS/hymap_025/lis_elevtn.bin
+HYMAP grid distance map:              ./input/LS_PARAMETERS/hymap_025/lis_nxtdst.bin
+HYMAP grid area map:                  ./input/LS_PARAMETERS/hymap_025/lis_grarea.bin
+HYMAP runoff time delay map:          ./input/LS_PARAMETERS/hymap_025/lis_getirana_paiva.bin
+HYMAP runoff time delay multiplier map:  ./input/LS_PARAMETERS/hymap_025/lis_trunoff.bin
+HYMAP baseflow time delay map:        ./input/LS_PARAMETERS/hymap_025/lis_tbasflw.bin
+HYMAP reference discharge map:        ./input/LS_PARAMETERS/hymap_025/lis_qrefer.bin
+HYMAP basin mask map:                 ./input/LS_PARAMETERS/hymap_025/lis_mask.bin
+HYMAP drainage area map:              ./input/LS_PARAMETERS/hymap_025/lis_uparea.bin
+HYMAP basin map:                      ./input/LS_PARAMETERS/hymap_025/lis_basin.bin
+HYMAP river flow type map:            ./input/LS_PARAMETERS/hymap_025/lis_trunoff.bin
+HYMAP baseflow dwi ratio map:         ./input/LS_PARAMETERS/hymap_025/lis_trunoff.bin
+HYMAP runoff dwi ratio map:           ./input/LS_PARAMETERS/hymap_025/lis_trunoff.bin
+
+
+