namelist_enkf.sh

#!/bin/bash
. $CONFIG_FILE
domain_id=$1
dx=`echo ${DX[$domain_id-1]}/1000 |bc -l`


##This if statement swiths the radar rv data off for parent domains
##  the radar data is only assimilated for d03
if [[ $domain_id != 3 ]]; then USE_RADAR_RV=false; fi

buffer=4 #buffer=0 if update bc; buffer=spec_bdy_width-1 if bc is fixed (perfect model)

cat << EOF
&enkf_parameter
numbers_en   = $NUM_ENS, 
expername    = '$EXP_NAME',  
enkfvar      = 'U         ', 'V         ', 'W         ', 'T         ', 'QVAPOR    ', 'QCLOUD    ', 'QRAIN     ', 'PH        ', 'MU        ', 'PSFC      ', 'P         ', 'PHB       ', 'PB        ', 'MUB       ',
updatevar    = 'U         ', 'V         ', 'W         ', 'T         ', 'QVAPOR    ', 'QCLOUD    ', 'QRAIN     ', 'PH        ', 'MU        ', 'PSFC      ', 'P         ',
update_is    = `echo 1+$buffer |bc`,
update_ie    = `echo ${E_WE[$domain_id-1]}-1-$buffer |bc`,
update_js    = `echo 1+$buffer |bc`,
update_je    = `echo ${E_SN[$domain_id-1]}-1-$buffer |bc`,
update_ks    = 1,
update_ke    = `echo ${E_VERT[$domain_id-1]}-1 |bc`,
inflate      = $INFLATION_COEF,
relax_opt    = $RELAX_OPT,
relax_adaptive = .$RELAX_ADAPTIVE.,
mixing       = $RELAXATION_COEF,
random_order = .false.,
print_detail = 0,
/

&parallel
manual_parallel = .true.,
nmcpu  = $NMCPU,
nicpu  = $NICPU,
njcpu  = $NJCPU,
/

&osse
use_ideal_obs    = .false.,
gridobs_is   = 25,
gridobs_ie   = `echo ${E_WE[$domain_id-1]}-20 |bc`,
gridobs_js   = 25,
gridobs_je   = `echo ${E_SN[$domain_id-1]}-20 |bc`,
gridobs_ks   = 1,
gridobs_ke   = `echo ${E_VERT[$domain_id-1]}-1 |bc`,
gridobs_int_x= 40,
gridobs_int_k= 1,
use_simulated= .false.,
/

&hurricane_PI 
use_hurricane_PI  = .false.,
hroi_hurricane_PI = 60,
vroi_hurricane_PI = 35,
/

&surface_obs
use_surface      = .$USE_SURFOBS.,
datathin_surface = $THIN_SURFACE,
hroi_surface     = $(printf %.0f `echo $HROI_SFC/$dx |bc -l`),
vroi_surface     = $VROI,
/

&sounding_obs
use_sounding      = .$USE_SOUNDOBS.,
datathin_sounding = $THIN_SOUNDING,
hroi_sounding     = $(printf %.0f `echo $HROI_UPPER/$dx |bc -l`),
vroi_sounding     = $VROI,
/

&profiler_obs
use_profiler      = .$USE_PROFILEROBS.,
datathin_profiler = $THIN_PROFILER,
hroi_profiler     = $(printf %.0f `echo $HROI_UPPER/$dx |bc -l`),
vroi_profiler     = $VROI,
/

&aircft_obs
use_aircft      = .$USE_AIREPOBS.,
datathin_aircft = $THIN_AIRCFT,
hroi_aircft     = $(printf %.0f `echo $HROI_UPPER/$dx |bc -l`),
vroi_aircft     = $VROI,
/

&metar_obs
use_metar      = .$USE_METAROBS.,
datathin_metar = $THIN_METAR,
hroi_metar     = $(printf %.0f `echo $HROI_SFC/$dx |bc -l`),
vroi_metar     = $VROI,
/

&sfcshp_obs
use_sfcshp      = .$USE_SHIPSOBS.,
datathin_sfcshp = $THIN_SFCSHP,
hroi_sfcshp     = $(printf %.0f `echo $HROI_SFC/$dx |bc -l`),
vroi_sfcshp     = $VROI,
/

&spssmi_obs
use_spssmi      = .$USE_SSMIOBS.,
datathin_spssmi = $THIN_SPSSMI,
hroi_spssmi     = $(printf %.0f `echo $HROI_UPPER/$dx |bc -l`),
vroi_spssmi     = $VROI,
/

&atovs_obs
use_atovs      = .$USE_ATOVS.,
datathin_atovs = $THIN_ATOVS,
hroi_atovs     = $(printf %.0f `echo $HROI_UPPER/$dx |bc -l`),
vroi_atovs     = $VROI,
/

&satwnd_obs
use_satwnd      = .$USE_GEOAMVOBS.,
datathin_satwnd = $THIN_SATWND,
hroi_satwnd     = $(printf %.0f `echo $HROI_UPPER/$dx |bc -l`),
vroi_satwnd     = $VROI,
/

&gpspw_obs
use_gpspw      = .$USE_GPSPWOBS.,
datathin_gpspw = $THIN_GPSPW,
hroi_gpspw     = $(printf %.0f `echo $HROI_SFC/$dx |bc -l`),
vroi_gpspw     = $VROI,
/

&radar_obs
radar_number   = 1,
use_radar_rf   = .$USE_RADAR_RF.,
use_radar_rv   = .$USE_RADAR_RV., 
datathin_radar = $THIN_RADAR,
hroi_radar     = $(printf %.0f `echo $HROI_RADAR/$dx |bc -l`),
vroi_radar     = $VROI_RADAR,
/

&airborne_radar   
use_airborne_rf   = .$USE_AIRBORNE_RF.,
use_airborne_rv   = .$USE_AIRBORNE_RV.,
datathin_airborne = $THIN_RADAR,
hroi_airborne     = $(printf %.0f `echo $HROI_RADAR/$dx |bc -l`),
vroi_airborne     = $VROI_RADAR,
/
EOF