12 use utilities,
only : error_handler, to_lower
156 character(len=*),
intent(in),
optional :: filename
157 character(len=250),
allocatable :: filename_to_use
159 integer :: is, ie, ierr
198 print*,
"- READ SETUP NAMELIST" 200 if (
present(filename))
then 201 filename_to_use = filename
203 filename_to_use =
"./fort.41" 206 open(41, file=trim(filename_to_use), iostat=ierr)
207 if (ierr /= 0)
call error_handler(
"OPENING SETUP NAMELIST.", ierr)
208 read(41, nml=config, iostat=ierr)
209 if (ierr /= 0)
call error_handler(
"READING SETUP NAMELIST.", ierr)
224 if (is == 0 .or. ie == 0)
then 225 call error_handler(
"CANT DETERMINE CRES FROM MOSAIC FILE.", 1)
231 call error_handler(
"MUST CONVERT EITHER AN ATM OR SFC FILE.", 1)
243 print*,
"- PROCESSING A REGIONAL NEST WITH A BOUNDARY HALO OF ",
halo_bndy 244 print*,
"- PROCESSING A REGIONAL NEST WITH A BLENDING HALO OF ",
halo_blend 252 if (trim(
tracers(is)) ==
"NULL")
exit 254 print*,
"- TRACER NAME IN OUTPUT FILE ", trim(
tracers(is))
261 print*,
"- TRACER NAME IN INPUT FILE ", trim(
tracers_input(is))
271 if(trim(
tracers(is)) ==
"spo" ) ierr = ierr - 1
272 if(trim(
tracers(is)) ==
"spo2" ) ierr = ierr - 1
273 if(trim(
tracers(is)) ==
"spo3" ) ierr = ierr - 1
276 print*,
"-ERROR: spo, spo2, and spo3 should be in tracers namelist" 277 call error_handler(
"WAM TRACER NAMELIST.", ierr)
279 print*,
"- WAM COLDSTART OPTION IS TURNED ON." 288 print*,
'- INPUT DATA FROM FV3 TILED RESTART FILES.' 290 print*,
'- INPUT DATA FROM FV3 TILED HISTORY FILES.' 291 case (
"gaussian_nemsio")
292 print*,
'- INPUT DATA FROM FV3 GAUSSIAN NEMSIO FILE.' 293 case (
"gfs_gaussian_nemsio")
294 print*,
'- INPUT DATA FROM SPECTRAL GFS GAUSSIAN NEMSIO FILE.' 296 print*,
'- INPUT DATA FROM SPECTRAL GFS SIGIO/SFCIO FILE.' 297 case (
"gaussian_netcdf")
298 print*,
'- INPUT DATA FROM FV3 GAUSSIAN NETCDF FILE.' 300 print*,
'- INPUT DATA FROM A GRIB2 FILE' 302 call error_handler(
"UNRECOGNIZED INPUT DATA TYPE.", 1)
311 call error_handler(
"FOR GRIB2 DATA, PLEASE PROVIDE GRIB2_FILE_INPUT_GRID", 1)
320 if (.not. any((/
character(4)::
"GFS",
"NAM",
"RAP",
"HRRR"/)==trim(
external_model)))
then 321 call error_handler(
"KNOWN SUPPORTED external_model INPUTS ARE GFS, NAM, RAP, AND HRRR. " // &
322 "IF YOU WISH TO PROCESS GRIB2 DATA FROM ANOTHER MODEL, YOU MAY ATTEMPT TO DO SO AT YOUR OWN RISK. " // &
323 "ONE WAY TO DO THIS IS PROVIDE NAM FOR external_model AS IT IS A RELATIVELY STRAIGHT-" // &
324 "FORWARD REGIONAL GRIB2 FILE. YOU MAY ALSO COMMENT OUT THIS ERROR MESSAGE IN " // &
325 "program_setup.f90 LINE 389. NO GUARANTEE IS PROVIDED THAT THE CODE WILL WORK OR "// &
326 "THAT THE RESULTING DATA WILL BE CORRECT OR WORK WITH THE ATMOSPHERIC MODEL.", 1)
337 print*,
"HRRR DATA DOES NOT CONTAIN SOIL TYPE INFORMATION. WITHOUT" 338 print*,
"GEOGRID_FILE_INPUT_GRID SPECIFIED, SOIL MOISTURE INTERPOLATION MAY BE LESS ACCURATE." 379 integer :: istat, k, nvars
380 character(len=500) :: line
381 character(len=20),
allocatable :: var_type(:)
385 print*,
"OPEN VARIABLE MAPPING FILE: ", trim(
varmap_file)
386 open(14, file=trim(
varmap_file), form=
'formatted', iostat=istat)
388 call error_handler(
"OPENING VARIABLE MAPPING FILE", istat)
396 read(14,
'(A)', iostat=istat) line
399 if ( trim(line) .eq.
'' ) cycle
402 if ( nvars == 0)
call error_handler(
"VARMAP FILE IS EMPTY.", -1)
409 allocate(var_type(nvars))
416 if (istat /= 0)
call error_handler(
"READING VARIABLE MAPPING FILE", istat)
417 if(trim(var_type(k))==
'T')
then 422 call error_handler(
"VARMAP TABLE CONTAINS TRACER ENTRIES FOR THOMPSON AEROSOLS liq_aero or "// &
423 "ice_aero. REMOVE THESE ENTRIES OR REMOVE THE NAMELIST ENTRY FOR "// &
424 "thomp_mp_climo_file AND TRY AGAIN.",1)
445 subroutine get_var_cond(var_name,this_miss_var_method,this_miss_var_value, &
446 this_field_var_name, loc)
449 character(len=20),
intent(in) :: var_name
451 character(len=20),
optional,
intent(out) :: this_miss_var_method, &
453 real(esmf_kind_r4),
optional,
intent(out):: this_miss_var_value
455 integer,
optional,
intent(out) :: loc
466 if (maxval(tmp).eq.0)
then 467 print*,
"WARNING: NO ENTRY FOR ", trim(var_name),
" IN VARMAP TABLE. WILL SKIP " // &
468 "VARIABLE IF NOT FOUND IN EXTERNAL MODEL FILE" 470 if(
present(this_miss_var_method)) this_miss_var_method =
"skip" 471 if(
present(this_miss_var_value)) this_miss_var_value = -9999.9_esmf_kind_r4
472 if(
present(this_field_var_name)) this_field_var_name =
"NULL" 473 if(
present(loc)) loc = 9999
477 if(
present(this_field_var_name)) this_field_var_name =
field_var_names(i)
478 if(
present(loc)) loc = i
520 integer,
intent(in) :: localpet
522 integer,
parameter :: num_statsgo = 16
523 real,
parameter :: smlow_statsgo = 0.5
524 real,
parameter :: smhigh_statsgo = 6.0
527 integer,
parameter :: num_zobler = 9
528 real,
parameter :: smlow_zobler = 0.5
529 real,
parameter :: smhigh_zobler = 6.0
531 integer :: num_soil_cats
533 real :: bb_statsgo(num_statsgo)
534 real :: maxsmc_statsgo(num_statsgo)
535 real :: satdk_statsgo(num_statsgo)
536 real :: satpsi_statsgo(num_statsgo)
538 real :: bb_zobler(num_zobler)
539 real :: maxsmc_zobler(num_zobler)
540 real :: satdk_zobler(num_zobler)
541 real :: satpsi_zobler(num_zobler)
543 real,
allocatable :: bb(:)
544 real :: smlow, smhigh
545 real,
allocatable :: satdk(:)
546 real,
allocatable :: satpsi(:)
547 real,
allocatable :: satdw(:)
550 data bb_statsgo /4.05, 4.26, 4.74, 5.33, 5.33, 5.25, &
551 6.77, 8.72, 8.17, 10.73, 10.39, 11.55, &
552 5.25, -9.99, 4.05, 4.26/
554 data maxsmc_statsgo /0.395, 0.421, 0.434, 0.476, 0.476, 0.439, &
555 0.404, 0.464, 0.465, 0.406, 0.468, 0.457, &
556 0.464, -9.99, 0.200, 0.421/
558 data satdk_statsgo /1.7600e-4, 1.4078e-5, 5.2304e-6, 2.8089e-6, 2.8089e-6, &
559 3.3770e-6, 4.4518e-6, 2.0348e-6, 2.4464e-6, 7.2199e-6, &
560 1.3444e-6, 9.7384e-7, 3.3770e-6, -9.99, 1.4078e-5, &
563 data satpsi_statsgo /0.0350, 0.0363, 0.1413, 0.7586, 0.7586, 0.3548, &
564 0.1349, 0.6166, 0.2630, 0.0977, 0.3236, 0.4677, &
565 0.3548, -9.99, 0.0350, 0.0363/
567 data bb_zobler /4.26, 8.72, 11.55, 4.74, 10.73, 8.17, &
570 data maxsmc_zobler /0.421, 0.464, 0.468, 0.434, 0.406, 0.465, &
573 data satdk_zobler /1.41e-5, 0.20e-5, 0.10e-5, 0.52e-5, 0.72e-5, &
574 0.25e-5, 0.45e-5, 0.34e-5, 1.41e-5/
576 data satpsi_zobler /0.040, 0.620, 0.470, 0.140, 0.100, 0.260, &
585 print*,
'- INPUT GRID USED ZOBLER SOIL TYPES.' 586 num_soil_cats = num_zobler
588 print*,
'- INPUT GRID USED STATSGO SOIL TYPES.' 589 num_soil_cats = num_statsgo
596 allocate(bb(num_soil_cats))
597 allocate(satdk(num_soil_cats))
598 allocate(satpsi(num_soil_cats))
599 allocate(satdw(num_soil_cats))
604 smhigh = smhigh_zobler
608 satpsi = satpsi_zobler
610 smlow = smlow_statsgo
611 smhigh = smhigh_statsgo
614 satdk = satdk_statsgo
615 satpsi = satpsi_statsgo
621 deallocate(bb, satdk, satpsi, satdw)
623 if (localpet == 0) print*,
'maxsmc input grid ',
maxsmc_input 624 if (localpet == 0) print*,
'wltsmc input grid ',
wltsmc_input 630 print*,
'- TARGET GRID USEING STATSGO SOIL TYPES.' 632 num_soil_cats = num_statsgo
640 allocate(satdk(num_soil_cats))
641 allocate(satdw(num_soil_cats))
643 smlow = smlow_statsgo
644 smhigh = smhigh_statsgo
647 satdk = satdk_statsgo
653 deallocate(satdk, satdw)
655 if (localpet == 0) print*,
'maxsmc target grid ',
maxsmc_target 677 maxsmc, bb, satpsi, satdw, refsmc, drysmc, wltsmc)
681 integer,
intent(in) :: num_soil_cats
683 real,
intent(in) :: smlow, smhigh
684 real,
intent(in) :: bb(num_soil_cats)
685 real,
intent(in) :: maxsmc(num_soil_cats)
686 real,
intent(in) :: satdk(num_soil_cats)
687 real,
intent(in) :: satpsi(num_soil_cats)
689 real,
intent(out) :: satdw(num_soil_cats)
690 real,
intent(out) :: refsmc(num_soil_cats)
691 real,
intent(out) :: drysmc(num_soil_cats)
692 real,
intent(out) :: wltsmc(num_soil_cats)
704 do i = 1, num_soil_cats
706 if (maxsmc(i) > 0.0)
then 708 satdw(i) = bb(i)*satdk(i)*(satpsi(i)/maxsmc(i))
709 refsmc1 = maxsmc(i)*(5.79e-9/satdk(i)) **(1.0/(2.0*bb(i)+3.0))
710 refsmc(i) = refsmc1 + (maxsmc(i)-refsmc1) / smhigh
711 wltsmc1 = maxsmc(i) * (200.0/satpsi(i))**(-1.0/bb(i))
712 wltsmc(i) = wltsmc1 - smlow * wltsmc1
719 drysmc(i) = wltsmc(i)
character(len=500), public fix_dir_target_grid
Directory containing target grid pre-computed fixed data (ex: soil type).
This module contains code to read the setup namelist file, handle the varmap file for GRIB2 data...
character(len=500), dimension(6), public atm_tracer_files_input_grid
File names of input atmospheric restart tracer files.
integer, public num_tracers_input
Number of atmospheric tracers in input file.
character(len=500), public nst_files_input_grid
File name of input nst data.
character(len=500), public varmap_file
REQUIRED.
real, dimension(:), allocatable, public wltsmc_input
Plant wilting point soil moisture content input grid.
integer, public cycle_mon
Cycle month.
logical, public lai_from_climo
If false, interpolate leaf area index from the input data to the target grid instead of using data fr...
character(len=6), public cres_target_grid
Target grid resolution, i.e., C768.
integer, public regional
For regional target grids.
character(len=20), dimension(max_tracers), public tracers
Name of each atmos tracer to be processed.
real, dimension(:), allocatable, public satpsi_target
Saturated soil potential, target grid.
real, dimension(:), allocatable, public maxsmc_input
Maximum soil moisture content input grid.
integer, public cycle_year
Cycle year.
integer, public nsoill_out
Number of soil levels desired in the output data.
character(len=20), dimension(max_tracers), public tracers_input
Name of each atmos tracer record in the input file.
real, dimension(:), allocatable, public bb_target
Soil 'b' parameter, target grid.
subroutine, public read_varmap
Reads the variable mapping table, which is required for initializing with GRIB2 data.
integer, public halo_bndy
Number of row/cols of lateral halo, where pure lateral bndy conditions are applied (regional target g...
subroutine, public get_var_cond(var_name, this_miss_var_method, this_miss_var_value, this_field_var_name, loc)
Search the variable mapping table to find conditions for handling missing variables.
logical, public sotyp_from_climo
If false, interpolate soil type from the input data to the target grid instead of using data from sta...
integer, public cycle_day
Cycle day.
logical, public minmax_vgfrc_from_climo
If false, interpolate min/max vegetation fraction from the input data to the target grid instead of u...
character(len=500), public thomp_mp_climo_file
Path/name to the Thompson MP climatology file.
character(len=20), dimension(:), allocatable, public missing_var_methods
Method to replace missing GRIB2 input records.
logical, public convert_sfc
Convert sfc data when true.
character(len=500), public orog_dir_target_grid
Directory containing the target grid orography files.
character(len=500), dimension(6), public orog_files_input_grid
Input grid orography files.
character(len=20), dimension(:), allocatable, public chgres_var_names
Varmap table variable name as recognized by this program.
character(len=20), public external_model
The model that the input data is derived from.
character(len=500), public orog_dir_input_grid
Directory containing the input grid orography files.
logical, public wam_cold_start
When true, cold start for whole atmosphere model.
subroutine calc_soil_params(num_soil_cats, smlow, smhigh, satdk, maxsmc, bb, satpsi, satdw, refsmc, drysmc, wltsmc)
Compute soil parameters.
character(len=500), public mosaic_file_target_grid
Target grid mosaic file.
character(len=500), public mosaic_file_input_grid
Input grid mosaic file.
logical, public convert_nst
Convert nst data when true.
character(len=500), dimension(6), public atm_files_input_grid
File names of input atmospheric data.
character(len=500), public vcoord_file_target_grid
Vertical coordinate definition file.
character(len=500), dimension(6), public orog_files_target_grid
Target grid orography files.
logical, public tg3_from_soil
If false, use lowest level soil temperature for the base soil temperature instead of using data from ...
character(len=20), dimension(:), allocatable, public field_var_names
The GRIB2 variable name in the varmap table.
logical, public convert_atm
Convert atmospheric data when true.
integer, parameter, public max_tracers
Maximum number of atmospheric tracers processed.
real, dimension(:), allocatable, public drysmc_input
Air dry soil moisture content input grid.
character(len=500), public geogrid_file_input_grid
Name of "geogrid" file, which contains static surface fields on the input grid.
real, dimension(:), allocatable, public wltsmc_target
Plant wilting point soil moisture content target grid.
character(len=500), public atm_weight_file
File containing pre-computed weights to horizontally interpolate atmospheric fields.
logical, public vgfrc_from_climo
If false, interpolate vegetation fraction from the input data to the target grid instead of using dat...
subroutine, public read_setup_namelist(filename)
Reads program configuration namelist.
subroutine, public calc_soil_params_driver(localpet)
Driver routine to compute soil parameters for each soil type.
character(len=500), public data_dir_input_grid
Directory containing input atm or sfc files.
character(len=25), public input_type
Input data type:
character(len=500), dimension(7), public atm_core_files_input_grid
File names of input atmospheric restart core files.
integer, public num_tracers
Number of atmospheric tracers to be processed.
real, dimension(:), allocatable, public refsmc_target
Reference soil moisture content target grid (onset of soil moisture stress).
character(len=500), public wam_parm_file
Full path to msis21.parm for WAM initialization.
character(len=500), public grib2_file_input_grid
REQUIRED.
logical, public use_thomp_mp_climo
When true, read and process Thompson MP climatological tracers.
logical, public vgtyp_from_climo
If false, interpolate vegetation type from the input data to the target grid instead of using data fr...
integer, public cycle_hour
Cycle hour.
logical, dimension(:), allocatable, public read_from_input
When false, variable was not read from GRIB2 input file.
real(kind=esmf_kind_r4), dimension(:), allocatable, public missing_var_values
If input GRIB2 record is missing, the variable is set to this value.
integer, public halo_blend
Number of row/cols of blending halo, where model tendencies and lateral boundary tendencies are appli...
character(len=500), dimension(6), public sfc_files_input_grid
File names containing input surface data.
real, dimension(:), allocatable, public refsmc_input
Reference soil moisture content input grid (onset of soil moisture stress).
real, dimension(:), allocatable, public maxsmc_target
Maximum soil moisture content target grid.
real, dimension(:), allocatable, public drysmc_target
Air dry soil moisture content target grid.