Skip to content

Commit

Permalink
Merge pull request #119 from gsketefian/community_develop
Browse files Browse the repository at this point in the history 
Community develop
  • Loading branch information
@gsketefian
gsketefian authored Oct 3, 2019
2 parents e62e7a0 + cce234c commit 7e4b67a
Show file tree
Hide file tree
Showing 2 changed files with 5 additions and 88 deletions.
13 changes: 5 additions & 8 deletions ush/NCL/make_FV3_RAP_domain_plots.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -36,14 +36,11 @@ fields='"PWATclm", "HGTsfc"' # Fields to plot
nlev="50" # Vertical index to plot for 3D fields
fcst_index="1,2" # Time index '2' or indices '1, 2' of forecast to plot or '"all"'
#fcst_index='"all"' # Time index '2' or indices '1, 2' of forecast to plot or '"all"'
#grid_dir="/scratch3/BMC/det/beck/FV3-CAM/work.C384r0p7n3_regional_RAP/INPUT"
#grid_dir="/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/work_FV3_regional_C96_2018032900/INPUT"
#grid_dir="/scratch3/BMC/fim/Julie.Schramm/regional_FV3_EMC_visit_20180509/work_FV3_regional_C96_2018032900/INPUT"
grid_dir="/scratch4/BMC/gmtb/Julie.Schramm/C96fv3gfs2016092900/INPUT"
#grid_dir="/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/work_FV3_regional_C96_2018032900/INPUT"

#RAP_grid_fn="/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/geo_em.d01.nc"
RAP_grid_fn="/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/geo_em.d01.RAP.nc"

grid_dir="/scratch2/BMC/gmtb/Julie.Schramm/C96fv3gfs2016092900/INPUT"

#RAP_grid_fn="/scratch2/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/geo_em.d01.nc"
RAP_grid_fn="/scratch2/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/geo_em.d01.RAP.nc"

if [ 1 = 1 ]; then
#
Expand Down
80 changes: 0 additions & 80 deletions ush/NCL/plot_FV3SAR_field_native.ncl
View file
Original file line number Diff line number Diff line change
Expand Up @@ -18,33 +18,6 @@

;help = True

;run_dir = "ABCD"
;run_dir = "/home/Gerard/fv3_regional/fv3sar_workflow/ush/NCL"
run_dir = "/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/run_dirs/C768_S1p5_RR3_EMCCONUS_UUUU"
;run_dir = "/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/run_dirs/C384_S1p65_RR5_HRRR_new_chgres_fv3sar12"
;run_dir = "/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/fv3sar_workflow/sorc/regional_grid"
run_dir = "/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/run_dirs/C384_S0p63_RR3_test_new_grid"
run_dir = "/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/run_dirs/C384_S0p63_RR3_test_new_grid05"
run_dir = "/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/run_dirs/C384_S0p63_RR3_test_new_grid05"
run_dir = "/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/run_dirs/C384_S0p63_RR3_test_new_grid06"
;run_dir = "/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/run_dirs/NX960_NY960_A0p21423_Kmns0p23209_HRRR_crap01"
;run_dir = "/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/run_dirs/NX1800_NY1120_A0p21423_Kmns0p23209_HRRR_crap01"
;run_dir = "/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/run_dirs/NX1800_NY1120_A0p21423_Kmns0p23209_HRRR_works"
run_dir = "/scratch3/BMC/fim/Gerard.Ketefian/regional_FV3_EMC_visit_20180509/run_dirs/C384_S0p63_RR3_test_new_grid05"
run_dir = "/scratch3/BMC/det/Gerard.Ketefian/UFS_CAM_test_instructions/expt_dirs/NX1700_NY1056_A0p21423_Kmns0p23209_HRRR_test_HRRRgrid"
run_dir = "/scratch3/BMC/det/Gerard.Ketefian/UFS_CAM_test_instructions/expt_dirs/NX1700_NY1056_A0p21423_Kmns0p23209_HRRR_test_HRRRgrid02"
run_dir = "/scratch3/BMC/det/Gerard.Ketefian/UFS_CAM_test_instructions/expt_dirs/NX1768_NY1032_A0p21423_Kmns0p23209_HRRR_test_HRRRgrid02"
run_dir = "/scratch3/BMC/det/Gerard.Ketefian/UFS_CAM_test_instructions/expt_dirs/NX1734_NY1008_A0p21423_Kmns0p23209_HRRR_test_HRRRgrid02"
;run_dir = "/scratch3/BMC/det/Gerard.Ketefian/UFS_CAM_test_chgres_cube/expt_dirs/test2"

run_dir = "/scratch4/BMC/gmtb/Linlin.Pan/ice/Linlin.Pan/rtgen.112531/tmp/fv3_regional_control"
run_dir = "/scratch3/BMC/det/beck/FV3-CAM/run_dirs/real_time/2019062300"

run_dir = "/scratch3/BMC/det/Gerard.Ketefian/FV3SAR_demo_PASPBL_20190725/expt_dirs/FFFFF"
run_dir = "/scratch3/BMC/det/beck/CAM_chgres/run_dirs/test_HRRR25km_FV3GFSICBCs_GFSphysics"
run_dir = "/scratch3/BMC/det/Gerard.Ketefian/FV3SAR_demo_PASPBL_20190725/expt_dirs/test_GSD_HRRR13km"
run_dir = "/scratch2/BMC/det/Gerard.Ketefian/UFS_CAM/FV3SAR_demo_PASPBL_20190725/expt_dirs/test_get_RAP_arcv_01"

run_dir = "/scratch2/BMC/det/Gerard.Ketefian/UFS_CAM/FV3SAR_demo_PASPBL_20190725/expt_dirs/test_GSD_HRRR25km"
run_dir = "/scratch2/BMC/det/Gerard.Ketefian/UFS_CAM/FV3SAR_demo_PASPBL_20190725/expt_dirs/test_GSD_HRRR13km"
run_dir = "/scratch2/BMC/det/Gerard.Ketefian/UFS_CAM/FV3SAR_demo_PASPBL_20190725/expt_dirs/test_GSD_HRRR3km"
Expand Down Expand Up @@ -672,59 +645,6 @@ end if

wrtcmp_config_fn := rd + "/model_configure"
wrtcmp_config_tmpl_fn = wrtcmp_config_fn + ".tmpl"

if (True) then
;if (False) then

sys_output = systemfunc("cp " + wrtcmp_config_tmpl_fn + " " + wrtcmp_config_fn)

nn = 0
lon_grid_cntr_native = lon_tile_cntr_tiles_to_plot(nn)
lat_grid_cntr_native = lat_tile_cntr_tiles_to_plot(nn)
lon_tile_corners_face_midpts_native \
:= lon_tile_corners_face_midpts_tiles_to_plot(nn,:)
lat_tile_corners_face_midpts_native \
:= lat_tile_corners_face_midpts_tiles_to_plot(nn,:)

;
; Read in the native grid parameters from the variable defintions file
; in the experiment directory.
;
var_defns_fp = run_dir + "/var_defns.sh"
param_names = (/ "delx", "dely" /)
num_params = dimsizes(param_names)
do np=0, num_params-1

param_name = param_names(np)
regex_search = "^\s*" + param_name + "=(" + char_dq + "([^" \
+ char_dq + "]+)" + char_dq + "|([^ " + char_dq + "]+))$"
regex_print = "\2\3"
sed_cmd = "sed -r -n -e " + char_sq + "s/" + regex_search + "/" \
+ regex_print + "/p" + char_sq + " " + var_defns_fp
sed_output = systemfunc(sed_cmd)

if (strcmp_exact(param_name, "delx")) then
dx_native = todouble(sed_output)
else if (strcmp_exact(param_name, "dely")) then
dy_native = todouble(sed_output)
end if
end if

end do

angle_units = "deg"

out := calc_wrtcmp_grid_params_lambert_cnfrml( \
wrtcmp_config_fn, \
lon_grid_cntr_native, lat_grid_cntr_native, \
lon_tile_corners_face_midpts_native, \
lat_tile_corners_face_midpts_native, \
dx_native, dy_native, \
angle_units)

print("DONE calculating write-component grid parameters.")
pause
end if
;
; **********************************************************************
;
Expand Down

0 comments on commit 7e4b67a

Please sign in to comment.