From 756afe6d29a7e16ccffb09bbfddb921559aae3f6 Mon Sep 17 00:00:00 2001 From: AliciaBentley-NOAA Date: Wed, 30 Sep 2020 21:53:19 +0000 Subject: [PATCH] Final versions of GFSv16 retrospective plotting scripts (9/30/20). --- check_for_data.ncl | 12 ++-- create_plot_lowlevel.sh | 12 ++-- create_plot_precip.sh | 12 ++-- create_plot_refl.sh | 12 ++-- create_plot_snod.sh | 12 ++-- create_plot_snow.sh | 12 ++-- create_plot_standard.sh | 12 ++-- drive_download_data.sh | 8 +-- drive_plot_images.sh | 10 ++-- htar_retros_prod.ncl | 101 +++++++++++++++++++++------------ htar_st4_anl.ncl | 15 +++++ plot_gfs_lambert_highcloud.ncl | 6 +- plot_gfs_lambert_lowcloud.ncl | 6 +- plot_gfs_lambert_midcloud.ncl | 2 + rename_image_dir.ncl | 6 +- 15 files changed, 143 insertions(+), 95 deletions(-) diff --git a/check_for_data.ncl b/check_for_data.ncl index d01e565..7e94012 100644 --- a/check_for_data.ncl +++ b/check_for_data.ncl @@ -11,17 +11,17 @@ print("======================") print("Script has started!") print("======================") -data_path = "/scratch2/NCEPDEV/stmp3/Alicia.Bentley/scripts/gfsv16/Karen2019/" -beginyyyymmddhh = "2019091600" -endyyyymmddhh = "2019092600" +data_path = "/scratch2/NCEPDEV/stmp3/Alicia.Bentley/scripts/gfsv16/Dorian2019/" +beginyyyymmddhh = "2019082700" +endyyyymmddhh = "2019090600" int_inc = 12 fcst_inc = 6 beginyyyymmddhh_anl = beginyyyymmddhh -endyyyymmddhh_anl = "2019100600" ;Add 10 days to endyyyymmddhh +endyyyymmddhh_anl = "2019091600" ;Add 10 days to endyyyymmddhh -data_to_check = (/"gfs","rap","st4","retros"/) -;data_to_check = (/"rap"/) +data_to_check = (/"rap","st4","gfs","retros"/) +;data_to_check = (/"retros"/) ;================================================================== ;================================================================== diff --git a/create_plot_lowlevel.sh b/create_plot_lowlevel.sh index e2cf0e2..a4008bd 100755 --- a/create_plot_lowlevel.sh +++ b/create_plot_lowlevel.sh @@ -16,11 +16,11 @@ cd ${DATA_PATH}/${CASE}/ if [[ $DOMAIN = "new" ]]; then -REGIONNAME="NC" -minlat=29.0 -maxlat=37.5 -minlon=-86.5 -maxlon=-72.5 +REGIONNAME="tn" +minlat=31.0 +maxlat=40.0 +minlon=-93.5 +maxlon=-79.5 echo "submitting ${REGIONNAME} domain script for ${CYCLE}" @@ -46,7 +46,7 @@ ncl 'scriptyyyymmddhh="${CYCLE}"' 'eventname="${CASE}"' 'scriptregion="${DOMAIN} ncl 'scriptyyyymmddhh="${CYCLE}"' 'eventname="${CASE}"' 'scriptregion="${DOMAIN}"' 'regionname="${REGIONNAME}"' 'minlat="${minlat}"' 'maxlat="${maxlat}"' 'minlon="${minlon}"' 'maxlon="${maxlon}"' 'fhr_inc="${FHR_INC}"' 'fhr_start="${FHR_START}"' 'fhr_end="${FHR_END}"' 'images_path="${IMAGES_PATH}"' 'cycle_number="${CYCLE_NUMBER}"' ${DATA_PATH}/${CASE}/${CYCLE}/plot_gfs_lambert_lowcloud.ncl -touch ../plot_${CYCLE}_${REGIONNAME}done +touch ../plot_${CYCLE}_${REGIONNAME}_done exit diff --git a/create_plot_precip.sh b/create_plot_precip.sh index 225c947..566bd8e 100755 --- a/create_plot_precip.sh +++ b/create_plot_precip.sh @@ -16,11 +16,11 @@ cd ${DATA_PATH}/${CASE}/ if [[ $DOMAIN = "new" ]]; then -REGIONNAME="NC" -minlat=29.0 -maxlat=37.5 -minlon=-86.5 -maxlon=-72.5 +REGIONNAME="tn" +minlat=31.0 +maxlat=40.0 +minlon=-93.5 +maxlon=-79.5 echo "submitting ${REGIONNAME} domain script for ${CYCLE}" @@ -43,7 +43,7 @@ ncl 'scriptyyyymmddhh="${CYCLE}"' 'eventname="${CASE}"' 'scriptregion="${DOMAIN} ncl 'scriptyyyymmddhh="${CYCLE}"' 'eventname="${CASE}"' 'scriptregion="${DOMAIN}"' 'regionname="${REGIONNAME}"' 'minlat="${minlat}"' 'maxlat="${maxlat}"' 'minlon="${minlon}"' 'maxlon="${maxlon}"' 'fhr_inc="${FHR_INC}"' 'fhr_start="${FHR_START}"' 'fhr_end="${FHR_END}"' 'images_path="${IMAGES_PATH}"' 'cycle_number="${CYCLE_NUMBER}"' ${DATA_PATH}/${CASE}/${CYCLE}/plot_gfs_lambert_highcloud.ncl -touch ../plot_${CYCLE}_${REGIONNAME}done +touch ../plot_${CYCLE}_${REGIONNAME}_done exit diff --git a/create_plot_refl.sh b/create_plot_refl.sh index b2b9a23..1c538df 100755 --- a/create_plot_refl.sh +++ b/create_plot_refl.sh @@ -16,11 +16,11 @@ cd ${DATA_PATH}/${CASE}/ if [[ $DOMAIN = "new" ]]; then -REGIONNAME="NC" -minlat=29.0 -maxlat=37.5 -minlon=-86.5 -maxlon=-72.5 +REGIONNAME="tn" +minlat=31.0 +maxlat=40.0 +minlon=-93.5 +maxlon=-79.5 echo "submitting ${REGIONNAME} domain script for ${CYCLE}" @@ -40,7 +40,7 @@ cd ${DATA_PATH}/${CASE}/${CYCLE}/ ncl 'scriptyyyymmddhh="${CYCLE}"' 'eventname="${CASE}"' 'scriptregion="${DOMAIN}"' 'regionname="${REGIONNAME}"' 'minlat="${minlat}"' 'maxlat="${maxlat}"' 'minlon="${minlon}"' 'maxlon="${maxlon}"' 'fhr_inc="${FHR_INC}"' 'fhr_start="${FHR_START}"' 'fhr_end="${FHR_END}"' 'images_path="${IMAGES_PATH}"' 'cycle_number="${CYCLE_NUMBER}"' ${DATA_PATH}/${CASE}/${CYCLE}/plot_gfs_lambert_refl.ncl -touch ../plot_${CYCLE}_${REGIONNAME}done +touch ../plot_${CYCLE}_${REGIONNAME}_done exit diff --git a/create_plot_snod.sh b/create_plot_snod.sh index ed6e796..f8d681f 100755 --- a/create_plot_snod.sh +++ b/create_plot_snod.sh @@ -16,11 +16,11 @@ cd ${DATA_PATH}/${CASE}/ if [[ $DOMAIN = "new" ]]; then -REGIONNAME="NC" -minlat=29.0 -maxlat=37.5 -minlon=-86.5 -maxlon=-72.5 +REGIONNAME="tn" +minlat=31.0 +maxlat=40.0 +minlon=-93.5 +maxlon=-79.5 echo "submitting ${REGIONNAME} domain script for ${CYCLE}" @@ -42,7 +42,7 @@ ncl 'scriptyyyymmddhh="${CYCLE}"' 'eventname="${CASE}"' 'scriptregion="${DOMAIN} ncl 'scriptyyyymmddhh="${CYCLE}"' 'eventname="${CASE}"' 'scriptregion="${DOMAIN}"' 'regionname="${REGIONNAME}"' 'minlat="${minlat}"' 'maxlat="${maxlat}"' 'minlon="${minlon}"' 'maxlon="${maxlon}"' 'fhr_inc="${FHR_INC}"' 'fhr_start="${FHR_START}"' 'fhr_end="${FHR_END}"' 'images_path="${IMAGES_PATH}"' 'cycle_number="${CYCLE_NUMBER}"' ${DATA_PATH}/${CASE}/${CYCLE}/plot_gfs_lambert_instsnod.ncl -touch ../plot_${CYCLE}_${REGIONNAME}done +touch ../plot_${CYCLE}_${REGIONNAME}_done exit diff --git a/create_plot_snow.sh b/create_plot_snow.sh index 8f18f71..3496d4f 100755 --- a/create_plot_snow.sh +++ b/create_plot_snow.sh @@ -16,11 +16,11 @@ cd ${DATA_PATH}/${CASE}/ if [[ $DOMAIN = "new" ]]; then -REGIONNAME="NC" -minlat=29.0 -maxlat=37.5 -minlon=-86.5 -maxlon=-72.5 +REGIONNAME="tn" +minlat=31.0 +maxlat=40.0 +minlon=-93.5 +maxlon=-79.5 echo "submitting ${REGIONNAME} domain script for ${CYCLE}" @@ -41,7 +41,7 @@ ncl 'scriptyyyymmddhh="${CYCLE}"' 'eventname="${CASE}"' 'scriptregion="${DOMAIN} ncl 'scriptyyyymmddhh="${CYCLE}"' 'eventname="${CASE}"' 'scriptregion="${DOMAIN}"' 'regionname="${REGIONNAME}"' 'minlat="${minlat}"' 'maxlat="${maxlat}"' 'minlon="${minlon}"' 'maxlon="${maxlon}"' 'fhr_inc="${FHR_INC}"' 'fhr_start="${FHR_START}"' 'fhr_end="${FHR_END}"' 'images_path="${IMAGES_PATH}"' 'cycle_number="${CYCLE_NUMBER}"' ${DATA_PATH}/${CASE}/${CYCLE}/plot_gfs_lambert_24hweasd.ncl -touch ../plot_${CYCLE}_${REGIONNAME}done +touch ../plot_${CYCLE}_${REGIONNAME}_done exit diff --git a/create_plot_standard.sh b/create_plot_standard.sh index d7a567e..719d393 100755 --- a/create_plot_standard.sh +++ b/create_plot_standard.sh @@ -16,11 +16,11 @@ cd ${DATA_PATH}/${CASE}/ if [[ $DOMAIN = "new" ]]; then -REGIONNAME="NC" -minlat=29.0 -maxlat=37.5 -minlon=-86.5 -maxlon=-72.5 +REGIONNAME="tn" +minlat=31.0 +maxlat=40.0 +minlon=-93.5 +maxlon=-79.5 echo "submitting ${REGIONNAME} domain script for ${CYCLE}" @@ -46,7 +46,7 @@ ncl 'scriptyyyymmddhh="${CYCLE}"' 'eventname="${CASE}"' 'scriptregion="${DOMAIN} ncl 'scriptyyyymmddhh="${CYCLE}"' 'eventname="${CASE}"' 'scriptregion="${DOMAIN}"' 'regionname="${REGIONNAME}"' 'minlat="${minlat}"' 'maxlat="${maxlat}"' 'minlon="${minlon}"' 'maxlon="${maxlon}"' 'fhr_inc="${FHR_INC}"' 'fhr_start="${FHR_START}"' 'fhr_end="${FHR_END}"' 'images_path="${IMAGES_PATH}"' 'cycle_number="${CYCLE_NUMBER}"' ${DATA_PATH}/${CASE}/${CYCLE}/plot_gfs_lambert_ptype.ncl -touch ../plot_${CYCLE}_${REGIONNAME}done +touch ../plot_${CYCLE}_${REGIONNAME}_done exit diff --git a/drive_download_data.sh b/drive_download_data.sh index 6fc1183..7d21227 100755 --- a/drive_download_data.sh +++ b/drive_download_data.sh @@ -1,6 +1,6 @@ #!/bin/bash # Author: Alicia Bentley -# Modified: 4 June 2020 +# Modified: 4 June 2019 # ################################################### # Script to get download necessary data from HPSS # @@ -11,7 +11,7 @@ export SCRIPTS_PATH='/scratch2/NCEPDEV/stmp3/Alicia.Bentley/scripts/gfsv16' export DATA_PATH='/scratch2/NCEPDEV/stmp3/Alicia.Bentley/scripts/gfsv16' -export CASE='Karen2019' +export CASE='Dorian2019' export FHR_START=0 export FHR_END=240 export FHR_INC=6 @@ -21,8 +21,8 @@ export GET_GFS=true export GET_RAP=true export GET_ST4=true -#for cycle in 2019092900 -for cycle in 2019091600 2019091612 2019091700 2019091712 2019091800 2019091812 2019091900 2019091912 2019092000 2019092012 2019092100 2019092112 2019092200 2019092212 2019092300 2019092312 2019092400 2019092412 2019092500 2019092512 2019092600 +for cycle in 2019082700 2019082712 2019082800 2019082812 2019082900 2019082912 2019083000 2019083012 2019083100 2019083112 2019090100 2019090112 2019090200 2019090212 2019090300 2019090312 2019090400 2019090412 2019090500 2019090512 2019090600 +#for cycle in 2020072918 do export CYCLE=${cycle} diff --git a/drive_plot_images.sh b/drive_plot_images.sh index c1d8a41..de9f840 100755 --- a/drive_plot_images.sh +++ b/drive_plot_images.sh @@ -11,7 +11,7 @@ export SCRIPTS_PATH='/scratch2/NCEPDEV/stmp3/Alicia.Bentley/scripts/gfsv16' export DATA_PATH='/scratch2/NCEPDEV/stmp3/Alicia.Bentley/scripts/gfsv16' export IMAGES_PATH='/scratch2/NCEPDEV/stmp3/Alicia.Bentley/scripts/gfsv16' -export CASE='Barbara2019' +export CASE='Dorian2019' export FHR_START=0 export FHR_END=240 export FHR_INC=6 @@ -19,15 +19,15 @@ export FHR_INC=6 export PLOT_FCST=true export PLOT_LLVL=true export PLOT_PPCP=true -export PLOT_REFL=true +export PLOT_REFL=true export PLOT_SNOD=false # leave false unless you have asked for NOHRSC snowfall analyses from Alicia export PLOT_WEASD=false # leave false unless you have asked for NOHRSC snowfall analyses from Alicia -#for cycle in 2019121000 -for cycle in 2019062712 2019062800 2019062812 2019062900 2019062912 2019063000 2019063012 2019070100 2019070112 2019070200 2019070212 2019070300 2019070312 2019070400 2019070412 2019070500 2019070512 2019070600 2019070612 2019070700 2019070712 +#for cycle in 2020020212 +for cycle in 2019082700 2019082712 2019082800 2019082812 2019082900 2019082912 2019083000 2019083012 2019083100 2019083112 2019090100 2019090112 2019090200 2019090212 2019090300 2019090312 2019090400 2019090412 2019090500 2019090512 2019090600 do -for domain in us hi epactc +for domain in pr se nwatl us do export CYCLE=${cycle} diff --git a/htar_retros_prod.ncl b/htar_retros_prod.ncl index c34d67d..7e11ba2 100644 --- a/htar_retros_prod.ncl +++ b/htar_retros_prod.ncl @@ -152,9 +152,9 @@ print("======================") do h=0,dimsizes(fdates)-1 ;print((/fhours(h)/)) - thisdate=fdates(h) - thisdate@units=timeunits - modelstrings=tostring(ut_calendar(thisdate,-3)) + intdate=datestart + intdate@units=timeunits + modelstrings=tostring(ut_calendar(intdate,-3)) ms=tostring(ut_calendar(datemodel,-3)) ;model date yearm=str_get_cols(ms,0,3) @@ -178,21 +178,46 @@ print("======================") print("Int: YYYYMMDDHH (FXXX): "+ymdm+hourm+" (F"+sprinti("%0.3i",fhours(h))+")") ; ========================================================================== -;Obtaining GEFSv16 forecast files +;Obtaining GFSv16 forecast files print("Downloading GFSv16 forecast file from HPSS...") + retro_cdate6 = cd_inv_calendar(2020,8,20,12,00,00,timeunits,0) ; Specify date formating changed + retro_cdate5e = cd_inv_calendar(2018,10,13,00,00,00,timeunits,0) ; Specify date formating changed retro_cdate5s = cd_inv_calendar(2018,8,31,00,00,00,timeunits,0) ; Specify date formating changed retro_cdate4 = cd_inv_calendar(2020,5,19,00,00,00,timeunits,0) ; Specify date formating changed + + retro_cdate3d = cd_inv_calendar(2020,3,5,00,00,00,timeunits,0) ; Specify date formating changed + retro_cdate3c = cd_inv_calendar(2020,2,23,00,00,00,timeunits,0) ; Specify date formating changed + + retro_cdate3b = cd_inv_calendar(2020,2,8,00,00,00,timeunits,0) ; Specify date formating changed + retro_cdate3a = cd_inv_calendar(2020,1,27,00,00,00,timeunits,0) ; Specify date formating changed + retro_cdate3 = cd_inv_calendar(2019,12,1,00,00,00,timeunits,0) ; Specify date formating changed retro_cdate2 = cd_inv_calendar(2019,9,1,00,00,00,timeunits,0) ; Specify date formating changed retro_cdate1 = cd_inv_calendar(2019,6,1,00,00,00,timeunits,0) ; Specify date formating changed - print("thisdate: "+modelstrings) + print("intdate: "+modelstrings) + print("Getting gfsv16!") + + if (intdate .ge. retro_cdate6) then + print("Initialized in v16rt2 stream [realtime (after /atmos)]!") + + cmd = "htar -xvf /NCEPDEV/emc-global/5year/emc.glopara/WCOSS_D/gfsv16/v16rt2/"+ymdm+hourm+"/gfsa.tar" + \ + " ./gfs."+ymdm+"/"+hourm+"/gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) + print(cmd) + system(cmd) + + cmd = "mv gfs."+ymdm+"/"+hourm+"/atmos/gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) + \ + " ./gfsv16."+ymdm+".t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h))+".grib2" + system(cmd) + +; cmd = "rm -fR gfs."+ymdm+"/" +; system(cmd) - if (thisdate .ge. retro_cdate4) then - print("Initialized in v16rt2 stream (realtime)!") + else if (intdate .ge. retro_cdate4) .and. (intdate .lt. retro_cdate6) then + print("Initialized in v16rt2 stream [realtime (but before /atmos)]!") cmd = "htar -xvf /NCEPDEV/emc-global/5year/emc.glopara/WCOSS_D/gfsv16/v16rt2/"+ymdm+hourm+"/gfsa.tar" + \ " ./gfs."+ymdm+"/"+hourm+"/gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) @@ -206,7 +231,7 @@ print("======================") ; cmd = "rm -fR gfs."+ymdm+"/" ; system(cmd) - else if (thisdate .ge. retro_cdate5s) .and. (thisdate .lt. retro_cdate5e) then + else if (intdate .ge. retro_cdate5s) .and. (intdate .lt. retro_cdate5e) then print("Initialized in v16retro5e!") cmd = "htar -xvf /NCEPDEV/emc-global/5year/emc.glopara/WCOSS_D/gfsv16/v16retro5e/"+ymdm+hourm+"/gfsa.tar" + \ @@ -221,22 +246,22 @@ print("======================") ; cmd = "rm -fR gfs."+ymdm+"/" ; system(cmd) - else if (thisdate .ge. retro_cdate3) .and. (thisdate .lt. retro_cdate4) then + else if (intdate .ge. retro_cdate3) .and. (intdate .lt. retro_cdate3a) .or. (intdate .ge. retro_cdate3b) .and. (intdate .lt. retro_cdate3c) .or. (intdate .ge. retro_cdate3d) .and. (intdate .lt. retro_cdate4) then print("Initialized in v16retro3e!") - + cmd = "htar -xvf /NCEPDEV/emc-global/5year/glopara/HERA/gfsv16/v16retro3e/"+ymdm+hourm+"/gfsa.tar" + \ " ./gfs."+ymdm+"/"+hourm+"/gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) print(cmd) system(cmd) - + cmd = "mv gfs."+ymdm+"/"+hourm+"/gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) + \ " ./gfsv16."+ymdm+".t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h))+".grib2" system(cmd) - + ; cmd = "rm -fR gfs."+ymdm+"/" ; system(cmd) - else if (thisdate .ge. retro_cdate2) .and. (thisdate .lt. retro_cdate3) then + else if (intdate .ge. retro_cdate2) .and. (intdate .lt. retro_cdate3) then print("Initialized in v16retro2e!") cmd = "htar -xvf /NCEPDEV/emc-global/5year/emc.glopara/WCOSS_D/gfsv16/v16retro2e/"+ymdm+hourm+"/gfsa.tar" + \ @@ -251,7 +276,7 @@ print("======================") ; cmd = "rm -fR gfs."+ymdm+"/" ; system(cmd) - else if (thisdate .ge. retro_cdate1) .and. (thisdate .lt. retro_cdate2) then + else if (intdate .ge. retro_cdate1) .and. (intdate .lt. retro_cdate2) then print("Initialized in v16retro1e!") cmd = "htar -xvf /NCEPDEV/emc-global/5year/emc.glopara/WCOSS_D/gfsv16/v16retro1e/"+ymdm+hourm+"/gfsa.tar" + \ @@ -266,10 +291,10 @@ print("======================") ; cmd = "rm -fR gfs."+ymdm+"/" ; system(cmd) - else if (thisdate .lt. retro_cdate1) then - print("Initialized in v16retro0e!") + else if (intdate .lt. retro_cdate1) .or. (intdate .ge. retro_cdate3a) .and. (intdate .lt. retro_cdate3b) .or. (intdate .ge. retro_cdate3c) .and. (intdate .lt. retro_cdate3d) then + print("Initialized in v16retro0ex!") - cmd = "htar -xvf /NCEPDEV/emc-global/5year/emc.glopara/WCOSS_D/gfsv16/v16retro0e/"+ymdm+hourm+"/gfsa.tar" + \ + cmd = "htar -xvf /NCEPDEV/emc-global/5year/glopara/HERA/gfsv16/v16retro0ex/"+ymdm+hourm+"/gfsa.tar" + \ " ./gfs."+ymdm+"/"+hourm+"/gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) print(cmd) system(cmd) @@ -287,6 +312,7 @@ print("======================") end if end if end if + end if ;----------------------------------------------------------------------------------- @@ -294,12 +320,13 @@ print("======================") print("Downloading GFSv15 forecast file from HPSS...") gfs_cdate3 = cd_inv_calendar(2020,2,26,00,00,00,timeunits,0) ; Specify date formating changed - gfs_cdate2 = cd_inv_calendar(2019,6,12,00,00,00,timeunits,0) ; Specify date formating changed - gfs_cdate1 = cd_inv_calendar(2017,7,20,00,00,00,timeunits,0) ; Specify date formating changed + gfs_cdate2 = cd_inv_calendar(2019,6,12,12,00,00,timeunits,0) ; Specify date formating changed + gfs_cdate1 = cd_inv_calendar(2018,12,15,00,00,00,timeunits,0) ; Specify date formating changed - print("thisdate: "+modelstrings) + print("intdate: "+modelstrings) + print("Getting gfsv15!") - if (thisdate .ge. gfs_cdate3) then + if (intdate .ge. gfs_cdate3) then print("Valid date is after GFS naming convention was FINALLY corrected!") cmd = "htar -xvf /NCEPPROD/hpssprod/runhistory/rh"+yearm+"/"+str_get_cols(ymdm,0,5)+"/"+ymdm+"/com_gfs_prod_gfs."+ymdm+"_"+hourm+".gfs_pgrb2.tar" + \ @@ -314,7 +341,7 @@ print("======================") ; cmd = "rm -fR gfs."+ymdm+"/" ; system(cmd) - else if (thisdate .ge. gfs_cdate2 .and. thisdate .lt. gfs_cdate3) then + else if (intdate .ge. gfs_cdate2 .and. intdate .lt. gfs_cdate3) then print("Valid date is after GFS moved to Dell") cmd = "htar -xvf /NCEPPROD/hpssprod/runhistory/rh"+yearm+"/"+str_get_cols(ymdm,0,5)+"/"+ymdm+"/gpfs_dell1_nco_ops_com_gfs_prod_gfs."+ymdm+"_"+hourm+".gfs_pgrb2.tar" + \ @@ -329,34 +356,34 @@ print("======================") ; cmd = "rm -fR gfs."+ymdm+"/" ; system(cmd) - else if (thisdate .ge. gfs_cdate1 .and. thisdate .lt. gfs_cdate2) then - print("Valid date is before GFS moved to Dell") + else if (intdate .ge. gfs_cdate1 .and. intdate .lt. gfs_cdate2) then + print("Getting GFSv15 from its real-time archive (wasn't ops yet)") - cmd = "htar -xvf /NCEPPROD/hpssprod/runhistory/rh"+yearm+"/"+str_get_cols(ymdm,0,5)+"/"+ymdm+"/gpfs_hps_nco_ops_com_gfs_prod_gfs."+ymdm+hourm+".pgrb2_0p25.tar" + \ - " ./gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) - ;print(cmd) + cmd = "htar -xvf /NCEPDEV/emc-global/5year/emc.glopara/WCOSS_C/Q2FY19/prfv3rt3/"+ymdm+hourm+"/gfsa.tar" + \ + " ./gfs."+ymdm+"/"+hourm+"/gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) + print(cmd) system(cmd) - - cmd = "mv gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) + \ + + cmd = "mv gfs."+ymdm+"/"+hourm+"/gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) + \ " ./gfsv15."+ymdm+".t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h))+".grib2" system(cmd) -; cmd = "rm -fR gfs."+ymdm+hourm+"/" +; cmd = "rm -fR gfs."+ymdm+"/" ; system(cmd) - else if (thisdate .lt. gfs_cdate1) then - print("Valid date is before GFS name change at 0000 UTC 20 July 2017") + else if (intdate .lt. gfs_cdate1) then + print("Getting GFSv15 from it's retrospective run (wasn't ops yet)") - cmd = "htar -xvf /NCEPPROD/hpssprod/runhistory/rh"+yearm+"/"+str_get_cols(ymdm,0,5)+"/"+ymdm+"/com2_gfs_prod_gfs."+ymdm+hourm+".pgrb2_0p25.tar" + \ - " ./gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) + cmd = "htar -xvf /NCEPDEV/emc-global/5year/emc.glopara/WCOSS_C/Q2FY19/prfv3rt3s/"+ymdm+hourm+"/gfsa.tar" + \ + " ./gfs."+ymdm+"/"+hourm+"/gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) print(cmd) system(cmd) - cmd = "mv gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) + \ - " ./gfs15."+ymdm+".t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h))+".grib2" + cmd = "mv gfs."+ymdm+"/"+hourm+"/gfs.t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h)) + \ + " ./gfsv15."+ymdm+".t"+hourm+"z.pgrb2.0p25.f"+sprinti("%0.3i",fhours(h))+".grib2" system(cmd) -; cmd = "rm -fR gfs."+ymdm+hourm+"/" +; cmd = "rm -fR gfs."+ymdm+"/" ; system(cmd) end if diff --git a/htar_st4_anl.ncl b/htar_st4_anl.ncl index 2230085..5b21335 100644 --- a/htar_st4_anl.ncl +++ b/htar_st4_anl.ncl @@ -139,6 +139,7 @@ print("======================") ;;;;; If not found, pull file from HPSS. ;;; define start and end dates for changing ST4 tarfile naming conventions on HPSS + st4_date4 = cd_inv_calendar(2020,7,19,00,00,00,timeunits,0) ; Specify date formating changed st4_date3 = cd_inv_calendar(2020,2,19,00,00,00,timeunits,0) ; Specify date formating changed st4_date2 = cd_inv_calendar(2019,8,13,00,00,00,timeunits,0) ; Specify date formating changed st4_date1 = cd_inv_calendar(2017,4,27,00,00,00,timeunits,0) ; Specify date when HPSS changed @@ -158,6 +159,8 @@ print("======================") print("Downloading Stage-IV file analysis file from HPSS") + + if (thisdate .lt. st4_date4) then cmd = "htar -xvf /NCEPPROD/hpssprod/runhistory/rh"+year+"/"+str_get_cols(ymd,0,5)+"/"+ymd+"/"+st4tar+"."+ymd+".tar" + \ " ./ST4."+modelstrings+".06h.gz" system(cmd) @@ -168,6 +171,18 @@ print("======================") cmd = "mv ST4."+modelstrings+".06h ./ST4."+modelstrings+".06h.grb" system(cmd) + else + cmd = "htar -xvf /NCEPPROD/hpssprod/runhistory/rh"+year+"/"+str_get_cols(ymd,0,5)+"/"+ymd+"/"+st4tar+"."+ymd+".tar" + \ + " ./st4_conus."+modelstrings+".06h.grb2" + system(cmd) + + cmd = "mv st4_conus."+modelstrings+".06h.grb2 ./ST4."+modelstrings+".06h.grb" + system(cmd) + + end if + + + end do diff --git a/plot_gfs_lambert_highcloud.ncl b/plot_gfs_lambert_highcloud.ncl index 88aa816..c17c0ae 100644 --- a/plot_gfs_lambert_highcloud.ncl +++ b/plot_gfs_lambert_highcloud.ncl @@ -183,10 +183,12 @@ if (fhours(h).eq.0) then slppara=slppara*0. else if (fhours(h).eq.6) then slpprod = prod_gfs->TCDC_P8_L234_GLL0_avg - slppara = para_gfs->TCDC_P8_L234_GLL0_avg +; slppara = para_gfs->TCDC_P8_L234_GLL0_avg + slppara = para_gfs->HCDC_P8_L234_GLL0_avg else if (fhours(h).ge.12) then slpprod = prod_gfs->TCDC_P8_L234_GLL0_avg6h - slppara = para_gfs->TCDC_P8_L234_GLL0_avg6h +; slppara = para_gfs->TCDC_P8_L234_GLL0_avg6h + slppara = para_gfs->HCDC_P8_L234_GLL0_avg6h end if end if end if diff --git a/plot_gfs_lambert_lowcloud.ncl b/plot_gfs_lambert_lowcloud.ncl index 9a97841..a5072c8 100644 --- a/plot_gfs_lambert_lowcloud.ncl +++ b/plot_gfs_lambert_lowcloud.ncl @@ -183,10 +183,12 @@ if (fhours(h).eq.0) then slppara=slppara*0. else if (fhours(h).eq.6) then slpprod = prod_gfs->TCDC_P8_L214_GLL0_avg - slppara = para_gfs->TCDC_P8_L214_GLL0_avg +; slppara = para_gfs->TCDC_P8_L214_GLL0_avg + slppara = para_gfs->LCDC_P8_L214_GLL0_avg else if (fhours(h).ge.12) then slpprod = prod_gfs->TCDC_P8_L214_GLL0_avg6h - slppara = para_gfs->TCDC_P8_L214_GLL0_avg6h +; slppara = para_gfs->TCDC_P8_L214_GLL0_avg6h + slppara = para_gfs->LCDC_P8_L214_GLL0_avg6h end if end if end if diff --git a/plot_gfs_lambert_midcloud.ncl b/plot_gfs_lambert_midcloud.ncl index 7957e8f..72768ff 100644 --- a/plot_gfs_lambert_midcloud.ncl +++ b/plot_gfs_lambert_midcloud.ncl @@ -184,9 +184,11 @@ if (fhours(h).eq.0) then else if (fhours(h).eq.6) then slpprod = prod_gfs->TCDC_P8_L224_GLL0_avg slppara = para_gfs->TCDC_P8_L224_GLL0_avg +; slppara = para_gfs->MCDC_P8_L224_GLL0_avg else if (fhours(h).ge.12) then slpprod = prod_gfs->TCDC_P8_L224_GLL0_avg6h slppara = para_gfs->TCDC_P8_L224_GLL0_avg6h +; slppara = para_gfs->MCDC_P8_L224_GLL0_avg6h end if end if end if diff --git a/rename_image_dir.ncl b/rename_image_dir.ncl index 542e6db..4fedd3c 100644 --- a/rename_image_dir.ncl +++ b/rename_image_dir.ncl @@ -11,9 +11,9 @@ print("======================") print("Script has started!") print("======================") -image_path = "/scratch2/NCEPDEV/stmp3/Alicia.Bentley/scripts/gfsv16/Karen2019/images/" -beginyyyymmddhh = "2019091600" -endyyyymmddhh = "2019092600" +image_path = "/scratch2/NCEPDEV/stmp3/Alicia.Bentley/scripts/gfsv16/Dorian2019/images/" +beginyyyymmddhh = "2019082700" +endyyyymmddhh = "2019090600" int_inc = 12 ;==================================================================