From 7222f84041f2a2f221e85d6543657f229408328c Mon Sep 17 00:00:00 2001 From: "fanglin.yang" Date: Thu, 21 May 2020 14:09:36 +0000 Subject: [PATCH 1/9] modified: checkout.sh to check out modle tag GFS.v16.0.4. changes include: 1. Remove constraints on mixing length and background diffusitivity over inversion layers on land 2. Enhance mass flux for deep convection, hence to increase subsidzing warming to reduce cold bias in the lower tropospehre 3. Fix a RRTMg solar radiaiton bug which has impact in SW abosrption in the UV region in the upper atmospehre. --- sorc/checkout.sh | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/sorc/checkout.sh b/sorc/checkout.sh index c672f195f4..90670de550 100755 --- a/sorc/checkout.sh +++ b/sorc/checkout.sh @@ -9,7 +9,7 @@ if [[ ! -d fv3gfs.fd ]] ; then rm -f ${topdir}/checkout-fv3gfs.log git clone https://github.com/ufs-community/ufs-weather-model fv3gfs.fd >> ${topdir}/checkout-fv3gfs.log 2>&1 cd fv3gfs.fd - git checkout GFS.v16.0.3 + git checkout GFS.v16.0.4 git submodule update --init --recursive cd ${topdir} else From 96460e6e5f4b09b9fd10d39550bf022dc62c37ab Mon Sep 17 00:00:00 2001 From: "russ.treadon" Date: Fri, 22 May 2020 13:02:06 +0000 Subject: [PATCH 2/9] Issue #1: replace "nawips" with "gempak" in hpssarch_gen.sh path to gfs sfc and snd files --- ush/hpssarch_gen.sh | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/ush/hpssarch_gen.sh b/ush/hpssarch_gen.sh index e7af4f122d..8dd812453a 100755 --- a/ush/hpssarch_gen.sh +++ b/ush/hpssarch_gen.sh @@ -70,8 +70,8 @@ if [ $type = "gfs" ]; then echo "${dirname}trak.gfso.atcfunix.altg.${PDY}${cyc} " >>gfsa.txt echo "${dirname}storms.gfso.atcf_gen.${PDY}${cyc} " >>gfsa.txt echo "${dirname}storms.gfso.atcf_gen.altg.${PDY}${cyc} " >>gfsa.txt - echo "${dirname}nawips/gfs_${PDY}${cyc}.sfc " >>gfsa.txt - echo "${dirname}nawips/gfs_${PDY}${cyc}.snd " >>gfsa.txt + echo "${dirname}gempak/gfs_${PDY}${cyc}.sfc " >>gfsa.txt + echo "${dirname}gempak/gfs_${PDY}${cyc}.snd " >>gfsa.txt echo "${dirname}bufr.t${cyc}z " >>gfsa.txt echo "./logs/${CDATE}/gfs*.log " >>gfsa.txt From e248236b233b6bcf2dc008006a82ccd459d7e3f3 Mon Sep 17 00:00:00 2001 From: "fanglin.yang" Date: Fri, 22 May 2020 15:23:59 +0000 Subject: [PATCH 3/9] modified: config.fcst to set if [ $LEVS = "128" -a $CDUMP = "gdas" ]; then ... lheatstrg=".false." --- parm/config/config.fcst | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/parm/config/config.fcst b/parm/config/config.fcst index e5d8e5eb64..1425676dc2 100755 --- a/parm/config/config.fcst +++ b/parm/config/config.fcst @@ -59,7 +59,7 @@ export do_tofd=".true." export launch_level=$(echo "$LEVS/2.35" |bc) # Sponge layer settings for L127 -if [ $LEVS = "128" ]; then +if [ $LEVS = "128" -a $CDUMP = "gdas" ]; then export tau=5.0 export rf_cutoff=1.0e3 export d2_bg_k1=0.20 @@ -79,7 +79,7 @@ if [ $lsm -eq 2 ]; then export lheatstrg=".false." export landice=".false." else - export lheatstrg=".true." + export lheatstrg=".false." export landice=".true." fi From 0dff61f103bab78f9d2426599f976dc318eff4e0 Mon Sep 17 00:00:00 2001 From: "fanglin.yang" Date: Mon, 25 May 2020 03:08:27 +0000 Subject: [PATCH 4/9] modified: config.fv3 to reduce tasks assigned to the WAVE component. 70 tasks at C768 is adequate for wave. --- parm/config/config.fv3 | 28 ++++++++++++++-------------- 1 file changed, 14 insertions(+), 14 deletions(-) diff --git a/parm/config/config.fv3 b/parm/config/config.fv3 index f8fbbfce16..8aaaee3d8e 100755 --- a/parm/config/config.fv3 +++ b/parm/config/config.fv3 @@ -41,8 +41,8 @@ case $case_in in export layout_y=4 export layout_x_gfs=2 export layout_y_gfs=4 - export npe_wav=70 - export npe_wav_gfs=70 + export npe_wav=14 + export npe_wav_gfs=14 export nth_fv3=1 export cdmbgwd="0.071,2.1,1.0,1.0" # mountain blocking, ogwd, cgwd, cgwd src scaling export WRITE_GROUP=1 @@ -57,8 +57,8 @@ case $case_in in export layout_y=4 export layout_x_gfs=4 export layout_y_gfs=4 - export npe_wav=70 - export npe_wav_gfs=70 + export npe_wav=14 + export npe_wav_gfs=14 export nth_fv3=1 export cdmbgwd="0.14,1.8,1.0,1.0" # mountain blocking, ogwd, cgwd, cgwd src scaling export WRITE_GROUP=1 @@ -73,8 +73,8 @@ case $case_in in export layout_y=6 export layout_x_gfs=4 export layout_y_gfs=6 - export npe_wav=70 - export npe_wav_gfs=70 + export npe_wav=21 + export npe_wav_gfs=21 export nth_fv3=2 export cdmbgwd="0.23,1.5,1.0,1.0" # mountain blocking, ogwd, cgwd, cgwd src scaling export WRITE_GROUP=1 @@ -89,8 +89,8 @@ case $case_in in export layout_y=8 export layout_x_gfs=6 export layout_y_gfs=6 - export npe_wav=270 - export npe_wav_gfs=270 + export npe_wav=35 + export npe_wav_gfs=35 export nth_fv3=1 export cdmbgwd="1.1,0.72,1.0,1.0" # mountain blocking, ogwd, cgwd, cgwd src scaling export WRITE_GROUP=1 @@ -105,8 +105,8 @@ case $case_in in export layout_y=8 export layout_x_gfs=16 export layout_y_gfs=12 - export npe_wav=280 - export npe_wav_gfs=280 + export npe_wav=70 + export npe_wav_gfs=70 export nth_fv3=4 export cdmbgwd="4.0,0.15,1.0,1.0" # mountain blocking, ogwd, cgwd, cgwd src scaling export WRITE_GROUP=2 @@ -121,8 +121,8 @@ case $case_in in export layout_y=16 export layout_x_gfs=8 export layout_y_gfs=16 - export npe_wav=270 - export npe_wav_gfs=270 + export npe_wav=140 + export npe_wav_gfs=140 export nth_fv3=4 export cdmbgwd="4.0,0.10,1.0,1.0" # mountain blocking, ogwd, cgwd, cgwd src scaling export WRITE_GROUP=4 @@ -137,8 +137,8 @@ case $case_in in export layout_y=32 export layout_x_gfs=16 export layout_y_gfs=32 - export npe_wav=270 - export npe_wav_gfs=270 + export npe_wav=140 + export npe_wav_gfs=140 export nth_fv3=4 export cdmbgwd="4.0,0.05,1.0,1.0" # mountain blocking, ogwd, cgwd, cgwd src scaling export WRITE_GROUP=4 From 38cd82133b64af289648ca5e88fde1ac3dcaeab0 Mon Sep 17 00:00:00 2001 From: "russ.treadon" Date: Tue, 26 May 2020 18:04:59 +0000 Subject: [PATCH 5/9] Issue #1: update config files * parm/config/config.anal - add logic to use correct global_convinfo.txt prior to GFS v15.3 implementation (2020052612) * parm/config/config.awips - set NAWIPSGRP to equal NPOSTGRP (config.post) * parm/config/config.fcst - add double quotes around CDUMP on levs test to prevent setup_workflow.py runtime error --- parm/config/config.anal | 18 ++++++++++-------- parm/config/config.awips | 2 +- parm/config/config.fcst | 2 +- 3 files changed, 12 insertions(+), 10 deletions(-) diff --git a/parm/config/config.anal b/parm/config/config.anal index adc893f266..7120c8a3a7 100755 --- a/parm/config/config.anal +++ b/parm/config/config.anal @@ -82,16 +82,18 @@ if [[ "$CDATE" -ge "2019110706" && "$CDATE" -lt "2020040718" ]]; then export OBERROR=$FIXgsi/gfsv16_historical/prepobs_errtable.global.2019110706 fi +# Assimilate 135 (T) & 235 (uv) Canadian AMDAR observations +if [[ "$CDATE" -ge "2020040718" && "$CDATE" -lt "2020052612" ]]; then_break + export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2020040718 +fi + # NOTE: -# As of 2020040718, gfsv16_historical/global_convinfo.txt.2020040718 is +# As of 2020052612, gfsv16_historical/global_convinfo.txt.2020052612 is # identical to ../global_convinfo.txt. Thus, the logic below is not -# needed at this time. A similar comment applies to ../prepobs_errtable.global -# and gfsv16_historical/prepobs_errtable.global.2019110706 - -# Assimilate 135 (T) & 235 (uv) Canadian AMDAR observations -##if [[ "$CDATE" -ge "2020040718" && "$CDATE" -lt "YYYYMMDDHH" ]]; then_break -## export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2020040718 -## export OBERROR=$FIXgsi/gfsv16_historical/prepobs_errtable.global.2020040718 +# needed at this time. +# Assimilate COSMIC-2 GPS +##if [[ "$CDATE" -ge "2020052612" && "$CDATE" -lt "YYYYMMDDHH" ]]; then_break +## export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2020052612 ##fi diff --git a/parm/config/config.awips b/parm/config/config.awips index c5a8d45caa..76a9c21536 100755 --- a/parm/config/config.awips +++ b/parm/config/config.awips @@ -12,6 +12,6 @@ export AWIPS20SH=$HOMEgfs/jobs/JGFS_AWIPS_20KM_1P0DEG export AWIPSG2SH=$HOMEgfs/jobs/JGFS_AWIPS_G2 # No. of concurrent awips jobs -export NAWIPSGRP=21 +export NAWIPSGRP=42 echo "END: config.awips" diff --git a/parm/config/config.fcst b/parm/config/config.fcst index 1425676dc2..a26729bd08 100755 --- a/parm/config/config.fcst +++ b/parm/config/config.fcst @@ -59,7 +59,7 @@ export do_tofd=".true." export launch_level=$(echo "$LEVS/2.35" |bc) # Sponge layer settings for L127 -if [ $LEVS = "128" -a $CDUMP = "gdas" ]; then +if [ $LEVS = "128" -a "$CDUMP" = "gdas" ]; then export tau=5.0 export rf_cutoff=1.0e3 export d2_bg_k1=0.20 From 2e12e63b13533a8a1b741bf55cbbaba49beb14cc Mon Sep 17 00:00:00 2001 From: "russ.treadon" Date: Tue, 26 May 2020 19:18:57 +0000 Subject: [PATCH 6/9] Issue #1: update config.anal logic to point ABIBF at the correct GDA directory --- parm/config/config.anal | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/parm/config/config.anal b/parm/config/config.anal index 7120c8a3a7..70331e94e7 100755 --- a/parm/config/config.anal +++ b/parm/config/config.anal @@ -48,7 +48,7 @@ export AMSR2BF=${AMSR2BF:-/dev/null} export ABIBF="/dev/null" if [[ "$CDATE" -ge "2019022800" ]] ; then export ABIBF="$DMPDIR/${CDUMP}x.${PDY}/${cyc}/${CDUMP}.t${cyc}z.gsrcsr.tm00.bufr_d" - if [[ "$CDATE" -ge "2019111000" ]]; then + if [[ "$CDATE" -ge "2019111000" && "$CDATE" -le "2020052612" ]]; then export ABIBF="$DMPDIR/${CDUMP}y.${PDY}/${cyc}/${CDUMP}.t${cyc}z.gsrcsr.tm00.bufr_d" fi fi From 914cb8dff96e45e690c8e1a2932c8d846afb1655 Mon Sep 17 00:00:00 2001 From: "russ.treadon" Date: Tue, 26 May 2020 20:37:27 +0000 Subject: [PATCH 7/9] Issue #1: correct typo in parm/config/config.anal --- parm/config/config.anal | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/parm/config/config.anal b/parm/config/config.anal index 70331e94e7..60ce611678 100755 --- a/parm/config/config.anal +++ b/parm/config/config.anal @@ -83,7 +83,7 @@ if [[ "$CDATE" -ge "2019110706" && "$CDATE" -lt "2020040718" ]]; then fi # Assimilate 135 (T) & 235 (uv) Canadian AMDAR observations -if [[ "$CDATE" -ge "2020040718" && "$CDATE" -lt "2020052612" ]]; then_break +if [[ "$CDATE" -ge "2020040718" && "$CDATE" -lt "2020052612" ]]; then export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2020040718 fi @@ -92,7 +92,7 @@ fi # identical to ../global_convinfo.txt. Thus, the logic below is not # needed at this time. # Assimilate COSMIC-2 GPS -##if [[ "$CDATE" -ge "2020052612" && "$CDATE" -lt "YYYYMMDDHH" ]]; then_break +##if [[ "$CDATE" -ge "2020052612" && "$CDATE" -lt "YYYYMMDDHH" ]]; then ## export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2020052612 ##fi From 70d71310132813c1066a9b0f881ef686602a23c4 Mon Sep 17 00:00:00 2001 From: "russ.treadon" Date: Tue, 26 May 2020 20:39:33 +0000 Subject: [PATCH 8/9] Issue #1: remove "_break" from commented out lines in config.anal and config.prep. "_break" will cause failue if line active --- parm/config/config.anal | 2 +- parm/config/config.prep | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/parm/config/config.anal b/parm/config/config.anal index 60ce611678..5b045a7d56 100755 --- a/parm/config/config.anal +++ b/parm/config/config.anal @@ -119,7 +119,7 @@ fi # needed at this time # # Turn off assmilation of all Metop-A MHS -## if [[ "$CDATE" -ge "2020022012" && "$CDATE" -lt "YYYYMMDDHH" ]]; then_break +## if [[ "$CDATE" -ge "2020022012" && "$CDATE" -lt "YYYYMMDDHH" ]]; then ## export SATINFO=$FIXgsi/gfsv16_historical/global_satinfo.txt.2020022012 ## fi diff --git a/parm/config/config.prep b/parm/config/config.prep index 15c58c6121..71479e9b9d 100755 --- a/parm/config/config.prep +++ b/parm/config/config.prep @@ -48,7 +48,7 @@ fi # needed at this time # Set observation errors for type 135 (T) & 235 (uv) Canadian AMDAR observations -##if [[ "$CDATE" -ge "2020040718" && "$CDATE" -lt "YYYMMDDHH" ]]; then_break +##if [[ "$CDATE" -ge "2020040718" && "$CDATE" -lt "YYYMMDDHH" ]]; then ## export PRVT=$EXPDIR/prepobs_errtable.global ##fi From e87b5a18f9dbad2230b9d3324f0c48a7b25d9b62 Mon Sep 17 00:00:00 2001 From: wx20ha Date: Wed, 27 May 2020 02:01:57 +0000 Subject: [PATCH 9/9] JWAVE_PREP updated to use CDO_ROOT defined in config config.waveprep updated to default to WCOSS rtofs operational cdo if no module found wave_prnc_cur.sh bug fixed now provides proper fhr in temp file names exwave_prep.sh adjusted for WCOSS and Hera. --- jobs/JWAVE_PREP | 3 +-- parm/config/config.waveprep | 3 +++ scripts/exwave_prep.sh | 45 ++++++++++++++++++++++++------------- ush/wave_prnc_cur.sh | 21 +++++++++-------- 4 files changed, 46 insertions(+), 26 deletions(-) diff --git a/jobs/JWAVE_PREP b/jobs/JWAVE_PREP index 5505b0d52a..6e37281010 100755 --- a/jobs/JWAVE_PREP +++ b/jobs/JWAVE_PREP @@ -60,8 +60,7 @@ export pgmout=OUTPUT.$$ export MP_PULSE=0 # CDO required for processing RTOFS currents -# export CDO=${COMROOTp2}/nwprod/rtofs_glo.v1.2.0/bin/cdo -export CDO=/gpfs/dell2/emc/verification/noscrub/Todd.Spindler/CDO/bin/cdo +export CDO=${CDO_ROOT}/bin/cdo # Path to HOME Directory export FIXwave=${FIXwave:-$HOMEwave/fix/fix_wave_${NET}} diff --git a/parm/config/config.waveprep b/parm/config/config.waveprep index ad57a78b1a..19353c2440 100755 --- a/parm/config/config.waveprep +++ b/parm/config/config.waveprep @@ -36,6 +36,9 @@ export WAV_CUR_HF_DT=${WAV_CUR_HF_DT:-3} export WAV_CUR_HF_FH=${WAV_CUR_HF_FH:-72} export WAV_CUR_CDO_SMOOTH="NO" +# Location of CDO module +export CDO_ROOT=${CDO_ROOT:-/gpfs/dell1/nco/ops/nwprod/rtofs_shared/rtofs_cdo.v1.6.0} + if [ "${WW3ICEINP}" = "YES" ]; then export WAVICEFILE=${CDUMP}.t${cyc}z.seaice.5min.grib2 fi diff --git a/scripts/exwave_prep.sh b/scripts/exwave_prep.sh index 97957e9f91..65979a55f0 100755 --- a/scripts/exwave_prep.sh +++ b/scripts/exwave_prep.sh @@ -311,9 +311,15 @@ # 3.a Gather and pre-process grib2 files ymdh=$ymdh_beg + if [ ${CFP_MP:-"NO"} = "YES" ]; then nm=0 ; fi # Counter for MP CFP while [ "$ymdh" -le "$ymdh_end" ] do - echo "$USHwave/wave_g2ges.sh $ymdh > grb_$ymdh.out 2>&1" >> cmdfile + if [ ${CFP_MP:-"NO"} = "YES" ]; then + echo "$nm $USHwave/wave_g2ges.sh $ymdh > grb_$ymdh.out 2>&1" >> cmdfile + nm=`expr $nm + 1` + else + echo "$USHwave/wave_g2ges.sh $ymdh > grb_$ymdh.out 2>&1" >> cmdfile + fi ymdh=`$NDATE $WAV_WND_HOUR_INC $ymdh` done @@ -334,7 +340,11 @@ if [ "$wavenproc" -gt '1' ] then - ${wavempexec} ${wavenproc} ${wave_mpmd} cmdfile + if [ ${CFP_MP:-"NO"} = "YES" ]; then + ${wavempexec} -n ${wavenproc} ${wave_mpmd} cmdfile + else + ${wavempexec} ${wavenproc} ${wave_mpmd} cmdfile + fi exit=$? else ./cmdfile @@ -626,30 +636,24 @@ # Prepare files for cfp process rm -f cmdfile touch cmdfile - chmod 744 cmfile + chmod 744 cmdfile ymdh_rtofs=${PDY}00 # RTOFS runs once daily ymdh_end=`$NDATE ${FHMAX_WAV_CUR} ${ymdh_rtofs}` NDATE_DT=${WAV_CUR_HF_DT} FLGHF='T' + if [ ${CFP_MP:-"NO"} = "YES" ]; then nm=0 ; fi # Counter for MP CFP while [ "$ymdh_rtofs" -le "$ymdh_end" ] do # Timing has to be made relative to the single 00z RTOFS cycle for that PDY fhr_rtofs=`${NHOUR} ${ymdh_rtofs} ${PDY}00` fext='f' - if [ ${fhr_rtofs} -lt 0 ] - then -# Data from nowcast phase - fhr_rtofs=`expr 48 + ${fhr_rtofs}` - fext='n' - fi + fh3_rtofs=`printf "%03d" "${fhr_rtofs#0}"` - fhr_rtofs=`printf "%03d\n" ${fhr_rtofs}` - - curfile1h=${COMIN_WAV_CUR}/rtofs_glo_2ds_${fext}${fhr_rtofs}_1hrly_prog.nc - curfile3h=${COMIN_WAV_CUR}/rtofs_glo_2ds_${fext}${fhr_rtofs}_3hrly_prog.nc + curfile1h=${COMIN_WAV_CUR}/rtofs_glo_2ds_${fext}${fh3_rtofs}_1hrly_prog.nc + curfile3h=${COMIN_WAV_CUR}/rtofs_glo_2ds_${fext}${fh3_rtofs}_3hrly_prog.nc if [ -s ${curfile1h} ] && [ "${FLGHF}" = "T" ] ; then curfile=${curfile1h} @@ -671,7 +675,13 @@ echo ' ' fi - echo "$USHwave/wave_prnc_cur.sh $ymdh_rtofs $curfile > cur_$ymdh_rtofs.out 2>&1" >> cmdfile + if [ ${CFP_MP:-"NO"} = "YES" ]; then + echo "$nm $USHwave/wave_prnc_cur.sh $ymdh_rtofs $curfile $fhr_rtofs > cur_$ymdh_rtofs.out 2>&1" >> cmdfile + nm=`expr $nm + 1` + else + echo "$USHwave/wave_prnc_cur.sh $ymdh_rtofs $curfile $fhr_rtofs > cur_$ymdh_rtofs.out 2>&1" >> cmdfile + fi + if [ $fhr_rtofs -ge ${WAV_CUR_HF_FH} ] ; then NDATE_DT=${WAV_CUR_DT} fi @@ -691,7 +701,11 @@ if [ $wavenproc -gt '1' ] then - ${wavempexec} ${wavenproc} ${wave_mpmd} cmdfile + if [ ${CFP_MP:-"NO"} = "YES" ]; then + ${wavempexec} -n ${wavenproc} ${wave_mpmd} cmdfile + else + ${wavempexec} ${wavenproc} ${wave_mpmd} cmdfile + fi exit=$? else chmod 744 ./cmdfile @@ -731,6 +745,7 @@ for file in $files do + echo $file cat $file >> cur.${WAVECUR_FID} rm -f $file done diff --git a/ush/wave_prnc_cur.sh b/ush/wave_prnc_cur.sh index d94b78aff1..d05255bc6a 100755 --- a/ush/wave_prnc_cur.sh +++ b/ush/wave_prnc_cur.sh @@ -23,38 +23,41 @@ set -x ymdh_rtofs=$1 curfile=$2 +fhr=$3 +fh3=`printf "%03d" "${fhr#0}"` # Timing has to be made relative to the single 00z RTOFS cycle for that PDY mkdir -p rtofs_${ymdh_rtofs} cd rtofs_${ymdh_rtofs} -ncks -x -v sst,sss,layer_density $curfile cur_uv_${PDY}_${fext}${fhr}.nc -ncks -O -a -h -x -v Layer cur_uv_${PDY}_${fext}${fhr}.nc cur_temp1.nc +ncks -x -v sst,sss,layer_density $curfile cur_uv_${PDY}_${fext}${fh3}.nc +ncks -O -a -h -x -v Layer cur_uv_${PDY}_${fext}${fh3}.nc cur_temp1.nc ncwa -h -O -a Layer cur_temp1.nc cur_temp2.nc ncrename -h -O -v MT,time cur_temp2.nc ncrename -h -O -d MT,time cur_temp2.nc ncks -v u_velocity,v_velocity cur_temp2.nc cur_temp3.nc -mv -f cur_temp3.nc cur_uv_${PDY}_${fext}${fhr}_flat.nc +mv -f cur_temp3.nc cur_uv_${PDY}_${fext}${fh3}_flat.nc # Convert to regular lat lon file - +# If weights need to be regenerated due to CDO ver change, use: +# $CDO genbil,r4320x2160 rtofs_glo_2ds_f000_3hrly_prog.nc weights.nc cp ${FIXwave}/weights_rtofs_to_r4320x2160.nc ./weights.nc # Interpolate to regular 5 min grid -$CDO remap,r4320x2160,weights.nc cur_uv_${PDY}_${fext}${fhr}_flat.nc cur_5min_01.nc +$CDO remap,r4320x2160,weights.nc cur_uv_${PDY}_${fext}${fh3}_flat.nc cur_5min_01.nc # Perform 9-point smoothing twice to make RTOFS data less noisy when # interpolating from 1/12 deg RTOFS grid to 1/6 deg wave grid if [ "WAV_CUR_CDO_SMOOTH" = "YES" ]; then $CDO -f nc -smooth9 cur_5min_01.nc cur_5min_02.nc - $CDO -f nc -smooth9 cur_5min_02.nc cur_glo_uv_${PDY}_${fext}${fhr}_5min.nc + $CDO -f nc -smooth9 cur_5min_02.nc cur_glo_uv_${PDY}_${fext}${fh3}_5min.nc else - mv cur_5min_01.nc cur_glo_uv_${PDY}_${fext}${fhr}_5min.nc + mv cur_5min_01.nc cur_glo_uv_${PDY}_${fext}${fh3}_5min.nc fi # Cleanup -rm -f cur_temp[123].nc cur_5min_??.nc cur_glo_uv_${PDY}_${fext}${fhr}.nc weights.nc +rm -f cur_temp[123].nc cur_5min_??.nc cur_glo_uv_${PDY}_${fext}${fh3}.nc weights.nc if [ ${fhr} -gt ${WAVHINDH} ] then @@ -64,7 +67,7 @@ else fi rm -f cur.nc -ln -s cur_glo_uv_${PDY}_${fext}${fhr}_5min.nc cur.nc +ln -s cur_glo_uv_${PDY}_${fext}${fh3}_5min.nc cur.nc ln -s ${DATA}/mod_def.${WAVECUR_FID} ./mod_def.ww3 $EXECwave/ww3_prnc 1> prnc_${WAVECUR_FID}_${ymdh_rtofs}.out 2>&1