diff --git a/GEOS_Util/plots/nettoa/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc similarity index 100% rename from GEOS_Util/plots/nettoa/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc rename to GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc diff --git a/GEOS_Util/plots/nettoac/VERIFICATION.CERES_EBAF.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.CERES_EBAF.rc similarity index 100% rename from GEOS_Util/plots/nettoac/VERIFICATION.CERES_EBAF.rc rename to GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.CERES_EBAF.rc diff --git a/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.COADS.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.COADS.rc new file mode 100644 index 00000000..c207dacc --- /dev/null +++ b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.COADS.rc @@ -0,0 +1,48 @@ + +OBSNAM: COADS +OBSDSC: COADS_CLIMATOLOGY_(1945-1989) + + filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' + filetype: 'BIN' + fields: 'EMP' , 'SURFACE' , 'EMINUSP' , 8 + :: + + filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' + filetype: 'BIN' + fields: 'TPREC' , 'MOIST' , 'PRECIP6' , 8 + :: + + filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' + filetype: 'BIN' + fields: 'EVAP' , 'SURFACE' , 'EVAPRATE' , 8 + :: + + filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' + filetype: 'BIN' + fields: 'VENT', 'SURFACE' , 'W3' , 1 + :: + + filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' + filetype: 'BIN' + fields: 'SH' , 'SURFACE' , 'SENSIB3' , 1 + :: + + filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' + filetype: 'BIN' + fields: 'LHFX' , 'SURFACE' , 'LATENT3' , 1 + :: + + filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' + filetype: 'BIN' + fields: 'CLDTT' , 'SOLAR' , 'CLOUD' , 1 , + :: + + filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' + filetype: 'BIN' + fields: 'TAUX' , 'SURFACE' , 'TAUX3' , 100 + :: + + filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' + filetype: 'BIN' + fields: 'TAUY' , 'SURFACE' , 'TAUY3' , 100 + :: diff --git a/GEOS_Util/plots/hcmp/VERIFICATION.ECANA.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ECANA.rc similarity index 100% rename from GEOS_Util/plots/hcmp/VERIFICATION.ECANA.rc rename to GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ECANA.rc diff --git a/GEOS_Util/plots/lcmp/VERIFICATION.ECINTERIM.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ECINTERIM.rc similarity index 80% rename from GEOS_Util/plots/lcmp/VERIFICATION.ECINTERIM.rc rename to GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ECINTERIM.rc index 78330b32..86b66017 100644 --- a/GEOS_Util/plots/lcmp/VERIFICATION.ECINTERIM.rc +++ b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ECINTERIM.rc @@ -17,3 +17,7 @@ OBSDSC: ECMWF_INTERIM_ANA fields: 'SLP' , 'DYN' , 'MSL' , 1 :: + filename: 'VERIFICATION/ECMWF_Interim_FullRes/Interim_monthly_FLX.ctl' + filetype: 'BIN' + fields: 'TPREC' , 'MOIST' , 'tp' , 1000 + :: diff --git a/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ERA5.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ERA5.rc new file mode 100644 index 00000000..387f8f22 --- /dev/null +++ b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ERA5.rc @@ -0,0 +1,49 @@ +OBSNAM: ERA5 +OBSDSC: ECMWF_ERA5_ANA + + filename: 'VERIFICATION//ERA5_Monthly/era5_pl.tabl' + filetype: 'HDF' + fields: 'ZLE' , 'DYN' , 'h' , 1 + 'T' , 'DYN' , 't' , 1 + 'U' , 'DYN' , 'u' , 1 + 'V' , 'DYN' , 'v' , 1 + 'OMEGA' , 'DYN' , 'w' , 1 + 'Q' , 'MOIST' , 'q' , 1 + 'RH2' , 'MOIST' , 'r' , 1/100 + 'QITOT' , 'AGCM' , 'ciwc' , 1 + 'QLTOT' , 'AGCM' , 'clwc' , 1 + 'O3' , 'CHEMISTRY' , 'o3' , 1 + :: + + filename: 'VERIFICATION//ERA5_Monthly/era5_sl-atm.tabl' + filetype: 'HDF' + fields: 'SLP' , 'DYN' , 'msl' , 1 + 'PS' , 'DYN' , 'sp' , 1 + 'TS' , 'SURFACE' , 'skt' , 1 + 'T2M' , 'SURFACE' , 't2' , 1 + 'U10M' , 'SURFACE' , 'u10' , 1 + 'V10M' , 'SURFACE' , 'v10' , 1 + 'U10N' , 'SURFACE' , 'u10n' , 1 + 'V10N' , 'SURFACE' , 'v10n' , 1 + 'CLDTT' , 'SOLAR' , 'tcc' , 1 + 'CLDLO' , 'SOLAR' , 'lcc' , 1 + 'CLDMD' , 'SOLAR' , 'mcc' , 1 + 'CLDHI' , 'SOLAR' , 'hcc' , 1 + 'ZPBL' , 'TURBULENCE' , 'blh' , 1 + 'CN_PRCP' , 'MOIST' , 'cp' , 1000 + 'LS_PRCP' , 'MOIST' , 'lsp' , 1000 + 'TPREC' , 'MOIST' , 'tp' , 1000 + 'TQV' , 'AGCM' , 'tcwv' , 1 + 'TQI' , 'AGCM' , 'tciw' , 1 + 'TQL' , 'AGCM' , 'tclw' , 1 + 'EVAP' , 'SURFACE' , 'e' , -1/100 + 'MSNLWRF' , 'ERA5' , 'msnlwrf' , -1 + 'MSNLWRFCS' , 'ERA5' , 'msnlwrfcs' , -1 + 'MSDWLWRF' , 'ERA5' , 'msdwlwrf' , 1 + 'MSDWLWRFCS' , 'ERA5' , 'msdwlwrfcs' , 1 + 'RSRS' , 'SOLAR' , 'msnswrf' , 1 + 'RSCS' , 'SOLAR' , 'msnswrfcs' , 1 + 'SLRSF' , 'SOLAR' , 'msdwswrf' , 1 + 'SLRSFC' , 'SOLAR' , 'msdwswrfcs' , 1 + :: + diff --git a/GEOS_Util/plots/nettoac/VERIFICATION.ERBE.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ERBE.rc similarity index 99% rename from GEOS_Util/plots/nettoac/VERIFICATION.ERBE.rc rename to GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ERBE.rc index af1cad7b..802eae40 100644 --- a/GEOS_Util/plots/nettoac/VERIFICATION.ERBE.rc +++ b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ERBE.rc @@ -12,4 +12,3 @@ OBSDSC: ERBE_Monthly_Data 'RSR' , 'SOLAR' , 'NSR' , 1 'RSC' , 'SOLAR' , 'NSRCLR' , 1 :: - diff --git a/GEOS_Util/plots/hcmp/VERIFICATION.GDAS1.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.GDAS1.rc similarity index 100% rename from GEOS_Util/plots/hcmp/VERIFICATION.GDAS1.rc rename to GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.GDAS1.rc diff --git a/GEOS_Util/plots/precip/VERIFICATION.GPCP.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.GPCP.rc similarity index 100% rename from GEOS_Util/plots/precip/VERIFICATION.GPCP.rc rename to GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.GPCP.rc diff --git a/GEOS_Util/plots/eflux/VERIFICATION.GSSTF.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.GSSTF.rc similarity index 100% rename from GEOS_Util/plots/eflux/VERIFICATION.GSSTF.rc rename to GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.GSSTF.rc diff --git a/GEOS_Util/plots/cldhi/VERIFICATION.ISCCP.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ISCCP.rc similarity index 100% rename from GEOS_Util/plots/cldhi/VERIFICATION.ISCCP.rc rename to GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.ISCCP.rc diff --git a/GEOS_Util/plots/hcmp/VERIFICATION.NCEP.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.NCEP.rc similarity index 100% rename from GEOS_Util/plots/hcmp/VERIFICATION.NCEP.rc rename to GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.NCEP.rc diff --git a/GEOS_Util/plots/radnetgc/VERIFICATION.SRB.rc b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.SRB.rc similarity index 70% rename from GEOS_Util/plots/radnetgc/VERIFICATION.SRB.rc rename to GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.SRB.rc index 11803bcc..ec2deddf 100644 --- a/GEOS_Util/plots/radnetgc/VERIFICATION.SRB.rc +++ b/GEOS_Util/plots/VERIFICATION_RCS/VERIFICATION.SRB.rc @@ -1,35 +1,40 @@ OBSNAM: SRB -OBSDSC: SRB_Net_Sfc_Clear_Sky_Radiation_(198307-200506) +OBSDSC: SRB_Rel2.5_Monthly_Data_(198307-200506) filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 + fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 , :: filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 + fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 , :: filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 + fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 , :: filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 + fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 , :: filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 + fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 , + :: + + filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' + filetype: 'BIN' + fields: 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 , :: filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 + fields: 'OLR' , 'IRRAD' , 'OLR' , 1 , :: filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' @@ -52,8 +57,22 @@ OBSDSC: SRB_Net_Sfc_Clear_Sky_Radiation_(198307-200506) fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , :: + filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' + filetype: 'BIN' + fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , + :: + filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' filetype: 'BIN' fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , :: + filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' + filetype: 'BIN' + fields: 'TPW' , 'MOIST' , 'CWV' , 1 , + :: + + filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' + filetype: 'BIN' + fields: 'TQV' , 'AGCM' , 'CWV' , 1 , + :: diff --git a/GEOS_Util/plots/albedo/VERIFICATION.SRB.rc b/GEOS_Util/plots/albedo/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/albedo/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/cldlo/VERIFICATION.ISCCP.rc b/GEOS_Util/plots/cldlo/VERIFICATION.ISCCP.rc deleted file mode 100644 index 20d56c32..00000000 --- a/GEOS_Util/plots/cldlo/VERIFICATION.ISCCP.rc +++ /dev/null @@ -1,11 +0,0 @@ - -OBSNAM: ISCCP -OBSDSC: ISCCP_N36_Monthly_Data - - filename: 'VERIFICATION/Clouds_radiation/isccp/isccp.tabl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'CLDTOT' , 0.01 , - 'CLDLO' , 'SOLAR' , 'CLDLOW' , 0.01 , - 'CLDMD' , 'SOLAR' , 'CLDMID' , 0.01 , - 'CLDHI' , 'SOLAR' , 'CLDHI' , 0.01 , - :: diff --git a/GEOS_Util/plots/cldmd/VERIFICATION.ISCCP.rc b/GEOS_Util/plots/cldmd/VERIFICATION.ISCCP.rc deleted file mode 100644 index 20d56c32..00000000 --- a/GEOS_Util/plots/cldmd/VERIFICATION.ISCCP.rc +++ /dev/null @@ -1,11 +0,0 @@ - -OBSNAM: ISCCP -OBSDSC: ISCCP_N36_Monthly_Data - - filename: 'VERIFICATION/Clouds_radiation/isccp/isccp.tabl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'CLDTOT' , 0.01 , - 'CLDLO' , 'SOLAR' , 'CLDLOW' , 0.01 , - 'CLDMD' , 'SOLAR' , 'CLDMID' , 0.01 , - 'CLDHI' , 'SOLAR' , 'CLDHI' , 0.01 , - :: diff --git a/GEOS_Util/plots/cldtt/VERIFICATION.COADS.rc b/GEOS_Util/plots/cldtt/VERIFICATION.COADS.rc deleted file mode 100644 index cf81c05f..00000000 --- a/GEOS_Util/plots/cldtt/VERIFICATION.COADS.rc +++ /dev/null @@ -1,8 +0,0 @@ - -OBSNAM: COADS -OBSDSC: COADS_CLIMATOLOGY_(1945-1989) - - filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'CLOUD' , 1 , - :: diff --git a/GEOS_Util/plots/cldtt/VERIFICATION.ISCCP.rc b/GEOS_Util/plots/cldtt/VERIFICATION.ISCCP.rc deleted file mode 100644 index 20d56c32..00000000 --- a/GEOS_Util/plots/cldtt/VERIFICATION.ISCCP.rc +++ /dev/null @@ -1,11 +0,0 @@ - -OBSNAM: ISCCP -OBSDSC: ISCCP_N36_Monthly_Data - - filename: 'VERIFICATION/Clouds_radiation/isccp/isccp.tabl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'CLDTOT' , 0.01 , - 'CLDLO' , 'SOLAR' , 'CLDLOW' , 0.01 , - 'CLDMD' , 'SOLAR' , 'CLDMID' , 0.01 , - 'CLDHI' , 'SOLAR' , 'CLDHI' , 0.01 , - :: diff --git a/GEOS_Util/plots/cldtt/VERIFICATION.SRB.rc b/GEOS_Util/plots/cldtt/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/cldtt/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/configure b/GEOS_Util/plots/configure index 33528560..f4032e5a 100755 --- a/GEOS_Util/plots/configure +++ b/GEOS_Util/plots/configure @@ -82,19 +82,132 @@ s?ARCH?$arch?g EOF sed -f sedfile $GASCRP/udft.template > $GAUDFT -# Link Radiation Directories for Cloud-Cleared Diagnostics +# Link VERIFICATION Directories +# ----------------------------------------------------------------------------------- + +# VERIFICATION.SRB +# --------------------- +set dirs = `echo olrc olrcf radnetg radnetgc radswt tpw cldtt albedo \ + lwgnet lwgnetc lwgnetcf lwgup lwgupc lwgupcf lwgdwn lwgdwnc lwgdwncf \ + swgnet swgnetc swgnetcf swgup swgupc swgupcf swgdwn swgdwnc swgdwncf ` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.SRB.rc . +end + +# VERIFICATION.ERBE +# ---------------------- +set dirs = `echo osr osrc osrcf olr olrc olrcf nettoa nettoac toacf` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.ERBE.rc . +end + +# VERIFICATION.CERES_EBAF-TOA_Ed4.0 +# -------------------------------------- +set dirs = `echo toacf tseries osrcf osrc osr radswt olrc olr olrcf nettoac nettoa` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc . +end + +# VERIFICATION.CERES_EBAF +# ----------------------- +set dirs = `echo toacf olrc olrcf nettoac tseries nettoa olr osrcf osrc osr radswt` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.CERES_EBAF.rc . +end + +# VERIFICATION.COADS +# ----------------------- +set dirs = `echo emp speed eflux hflux cldtt taux tauy` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.COADS.rc . +end + +# VERIFICATION.ISCCP +# ----------------------- +set dirs = `echo cldhi cldlo cldmd cldtt` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.ISCCP.rc . +end + +# VERIFICATION.GSSTF +# ----------------------- +set dirs = `echo eflux hflux speed taux tauy tpw` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.GSSTF.rc . +end + +# VERIFICATION.ECINTERIM +# --------------------------- +set dirs = `echo hcmp lcmp zcmp zcmp1 tseries precip` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.ECINTERIM.rc . +end + +# VERIFICATION.ECANA +# ----------------------- +set dirs = `echo hcmp lcmp zcmp` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.ECANA.rc . +end + +# VERIFICATION.ERA5 +# ---------------------- +set dirs = `echo lwgdwnc lwgdwn lwgnetc lwgnet lwgupc lwgup \ + swgdwnc swgdwn swgnetc swgnet swgupc swgup \ + precipcn precipls precip ` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.ERA5.rc . +end + +# VERIFICATION.GDAS1 +# ----------------------- +set dirs = `echo hcmp lcmp zcmp` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.GDAS1.rc . +end + +# VERIFICATION.NCEP +# ---------------------- +set dirs = `echo hcmp lcmp zcmp` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.NCEP.rc . +end + +# VERIFICATION.GPCP +# ---------------------- +set dirs = `echo precip tseries` +foreach dir ($dirs) +cd $GEOSUTIL/plots/$dir +/bin/ln -sf ../VERIFICATION_RCS/VERIFICATION.GPCP.rc . +end + + +# Link Radiation Directories for 5% Cloud-Cleared (cc5) Diagnostics # ----------------------------------------------------------------------------------- cd $GEOSUTIL/plots -/bin/ln -sf lwgupcf lwgupcf_cc5 -/bin/ln -sf lwgupc lwgupc_cc5 -/bin/ln -sf lwgnetcf lwgnetcf_cc5 -/bin/ln -sf lwgnetc lwgnetc_cc5 -/bin/ln -sf lwgdwncf lwgdwncf_cc5 -/bin/ln -sf lwgdwnc lwgdwnc_cc5 -/bin/ln -sf olrcf olrcf_cc5 -/bin/ln -sf olrc olrc_cc5 -/bin/ln -sf nettoac nettoac_cc5 -/bin/ln -sf toacf toacf_cc5 +if(! -e lwgupcf_cc5 ) /bin/ln -sf lwgupcf lwgupcf_cc5 +if(! -e lwgupc_cc5 ) /bin/ln -sf lwgupc lwgupc_cc5 +if(! -e lwgnetcf_cc5 ) /bin/ln -sf lwgnetcf lwgnetcf_cc5 +if(! -e lwgnetc_cc5 ) /bin/ln -sf lwgnetc lwgnetc_cc5 +if(! -e lwgdwncf_cc5 ) /bin/ln -sf lwgdwncf lwgdwncf_cc5 +if(! -e lwgdwnc_cc5 ) /bin/ln -sf lwgdwnc lwgdwnc_cc5 +if(! -e olrcf_cc5 ) /bin/ln -sf olrcf olrcf_cc5 +if(! -e olrc_cc5 ) /bin/ln -sf olrc olrc_cc5 +if(! -e nettoac_cc5 ) /bin/ln -sf nettoac nettoac_cc5 +if(! -e toacf_cc5 ) /bin/ln -sf toacf toacf_cc5 + # Link Cubed-Sphere FRAC Dataset for dc.gs Utility and CONST_2d Dataset for TopoShade # ----------------------------------------------------------------------------------- diff --git a/GEOS_Util/plots/eflux/VERIFICATION.COADS.rc b/GEOS_Util/plots/eflux/VERIFICATION.COADS.rc deleted file mode 100644 index 2f02b562..00000000 --- a/GEOS_Util/plots/eflux/VERIFICATION.COADS.rc +++ /dev/null @@ -1,8 +0,0 @@ - -OBSNAM: COADS -OBSDSC: COADS_CLIMATOLOGY_(1945-1989) - - filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' - filetype: 'BIN' - fields: 'LHFX' , 'SURFACE' , 'LATENT3' , 1 - :: diff --git a/GEOS_Util/plots/emp/VERIFICATION.COADS.rc b/GEOS_Util/plots/emp/VERIFICATION.COADS.rc deleted file mode 100644 index eba78bca..00000000 --- a/GEOS_Util/plots/emp/VERIFICATION.COADS.rc +++ /dev/null @@ -1,18 +0,0 @@ - -OBSNAM: COADS -OBSDSC: COADS_CLIMATOLOGY_(1945-1989) - - filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' - filetype: 'BIN' - fields: 'EMP' , 'SURFACE' , 'EMINUSP' , 8 - :: - - filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' - filetype: 'BIN' - fields: 'TPREC' , 'MOIST' , 'PRECIP6' , 8 - :: - - filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' - filetype: 'BIN' - fields: 'EVAP' , 'SURFACE' , 'EVAPRATE' , 8 - :: diff --git a/GEOS_Util/plots/fcst/gif_stats b/GEOS_Util/plots/fcst/gif_stats index 021404e3..9d7ec3d7 100755 --- a/GEOS_Util/plots/fcst/gif_stats +++ b/GEOS_Util/plots/fcst/gif_stats @@ -3,7 +3,7 @@ set plots = "hght tmpu uwnd vwnd" set types = "tcor cor rms" set seasons = "JAN" -set regions = "EUR GLO NAM NEQ NHE NWQ SEQ SHE SWQ TRO NPO SPO" +set regions = "EUR GLO NAM NEQ NHE NWQ SEQ SHE SWQ TRO NPO SPO XPO" set levs = "1000 850 700 500 400 300 250 200 150 100 z" # Create Postscript from Grads diff --git a/GEOS_Util/plots/fcst/pltcor.gs b/GEOS_Util/plots/fcst/pltcor.gs index a781cb74..073a5807 100644 --- a/GEOS_Util/plots/fcst/pltcor.gs +++ b/GEOS_Util/plots/fcst/pltcor.gs @@ -47,6 +47,8 @@ endif 'set dfile 1' 'getinfo nvars' nvars = result +'getinfo xdim' + xdim = result n=1 while(n<=nvars) @@ -70,7 +72,7 @@ endif z = 1 while ( z<=zmax ) x = 1 -while ( x<=12 ) +while ( x<=xdim ) 'set x 'x 'set z 'z 'cor 'root' 'DESC diff --git a/GEOS_Util/plots/fcst/pltcorcmp.gs b/GEOS_Util/plots/fcst/pltcorcmp.gs index 7a4730b0..9e9a1dea 100644 --- a/GEOS_Util/plots/fcst/pltcorcmp.gs +++ b/GEOS_Util/plots/fcst/pltcorcmp.gs @@ -112,6 +112,8 @@ endwhile say 'ctlinfo 'result 'getinfo nvars' nvars = result +'getinfo xdim' + xdim = result 'getinfo zdim' zdim = result m=0 @@ -163,7 +165,7 @@ if( field = subwrd(fields,j) ) while ( z<=numlevs ) if( level.z = 100 | level.z = 150 | level.z = 200 | level.z = 250 | level.z = 300 | level.z = 400 | level.z = 500 | level.z = 600 | level.z = 700 | level.z = 750 | level.z = 800 | level.z = 850 | level.z = 900 | level.z = 925 | level.z = 950 | level.z = 975 | level.z = 1000 ) x = 1 - while ( x<=12 ) + while ( x<=xdim ) 'run 'geosutil'/plots/grads_util/corcmp -x 'x' -lev 'level.z' -field 'field' -rc 'rcfile' -rms 'rms' -desc 'DESC' -debug FALSE' 'c' x = x + 1 diff --git a/GEOS_Util/plots/fcst/pltrms.gs b/GEOS_Util/plots/fcst/pltrms.gs index 49729e99..fac22c62 100644 --- a/GEOS_Util/plots/fcst/pltrms.gs +++ b/GEOS_Util/plots/fcst/pltrms.gs @@ -47,6 +47,8 @@ endif 'set dfile 1' 'getinfo nvars' nvars = result +'getinfo xdim' + xdim = result n=1 while(n<=nvars) @@ -70,7 +72,7 @@ endif z = 1 while ( z<=zmax ) x = 1 -while ( x<=12 ) +while ( x<=xdim ) 'set x 'x 'set z 'z 'rms 'root' 'DESC diff --git a/GEOS_Util/plots/fcst/pltrmsz.gs b/GEOS_Util/plots/fcst/pltrmsz.gs index b0e8076e..51366b6b 100644 --- a/GEOS_Util/plots/fcst/pltrmsz.gs +++ b/GEOS_Util/plots/fcst/pltrmsz.gs @@ -41,9 +41,11 @@ endif 'set display color white' 'rgbset' 'c' +'getinfo xdim' + xdim = result x = 1 -while ( x<=12 ) +while ( x<=xdim ) 'set x 'x 'rmsz h 'DESC 'c' diff --git a/GEOS_Util/plots/fcst/plttcor.gs b/GEOS_Util/plots/fcst/plttcor.gs index 2553ddab..dfbcf6a3 100644 --- a/GEOS_Util/plots/fcst/plttcor.gs +++ b/GEOS_Util/plots/fcst/plttcor.gs @@ -41,10 +41,12 @@ endif 'set display color white' 'rgbset' 'c' +'getinfo xdim' + xdim = result z = 1 x = 1 -while ( x<=12 ) +while ( x<=xdim ) 'set x 'x 'set z 'z n = 1 @@ -59,7 +61,7 @@ endwhile z = 1 while ( z<=10 ) x = 1 -while ( x<=10 ) +while ( x<=xdim ) 'set x 'x 'set z 'z n = 1 diff --git a/GEOS_Util/plots/get_exports b/GEOS_Util/plots/get_exports index 7c494a46..71841f73 100755 --- a/GEOS_Util/plots/get_exports +++ b/GEOS_Util/plots/get_exports @@ -1,7 +1,7 @@ #!/bin/csh -f set datetime = `date +%Y%m%d_%H%M%S` -cat $1 | sed -e 's/,/ , /g' > HISTORYRC.$datetime +cat $1 | sed -e 's/,/ , /g' | sed -e 's/*/@/g' | sed -e 's/?/ /g' > HISTORYRC.$datetime if( .$2 == . ) then set PLOT_SCRIPTS = `echo $GEOSUTIL/plots/landscape.script $GEOSUTIL/plots/portrait.script` diff --git a/GEOS_Util/plots/grads_util/cor.gs b/GEOS_Util/plots/grads_util/cor.gs index ec07191d..519ab735 100644 --- a/GEOS_Util/plots/grads_util/cor.gs +++ b/GEOS_Util/plots/grads_util/cor.gs @@ -92,6 +92,7 @@ if( xpos = 9 ) ; region = "North America (Lons:-140,-60 Lats: 20,60)" ; reg = if( xpos = 10 ) ; region = "Europe (Lons:-10,30 Lats: 30,60)" ; reg = "EUR" ; endif if( xpos = 11 ) ; region = "N.Polar (Lats: 60,90)" ; reg = "NPO" ; endif if( xpos = 12 ) ; region = "S.Polar (Lats: -60,-90)" ; reg = "SPO" ; endif +if( xpos = 13 ) ; region = "X.Polar (Lats: -60,60)" ; reg = "XPO" ; endif 'getinfo tdim' tdim = result diff --git a/GEOS_Util/plots/grads_util/corave.gs b/GEOS_Util/plots/grads_util/corave.gs index 2af440d1..68988997 100644 --- a/GEOS_Util/plots/grads_util/corave.gs +++ b/GEOS_Util/plots/grads_util/corave.gs @@ -133,6 +133,7 @@ if( xpos = 9 ) ; region = "North America (Lons:-140,-60 Lats: 20,60)" ; reg = if( xpos = 10 ) ; region = "Europe (Lons:-10,30 Lats: 30,60)" ; reg = "EUR" ; endif if( xpos = 11 ) ; region = "N.Polar (Lats: 60,90)" ; reg = "NPO" ; endif if( xpos = 12 ) ; region = "S.Polar (Lats: -60,-90)" ; reg = "SPO" ; endif +if( xpos = 13 ) ; region = "X.Polar (Lats: -60,60)" ; reg = "XPO" ; endif * Initialize Plot diff --git a/GEOS_Util/plots/grads_util/corcmp_plot.gs b/GEOS_Util/plots/grads_util/corcmp_plot.gs index d86b35c1..7a3316fb 100644 --- a/GEOS_Util/plots/grads_util/corcmp_plot.gs +++ b/GEOS_Util/plots/grads_util/corcmp_plot.gs @@ -122,7 +122,8 @@ if( xpos = 8 ) ; region = "S.E. Quadrant (Lons: 0,180 Lats: 0,-90)" ; reg = if( xpos = 9 ) ; region = "North America (Lons:-140,-60 Lats: 20,60)" ; reg = "NAM" ; endif if( xpos = 10 ) ; region = "Europe (Lons:-10,30 Lats: 30,60)" ; reg = "EUR" ; endif if( xpos = 11 ) ; region = "N.Polar (Lats: 60,90)" ; reg = "NPO" ; endif -if( xpos = 12 ) ; region = "S.Polar (Lats: -60,-90)" ; reg = "SPO" ; endif +if( xpos = 12 ) ; region = "S.Polar (Lats:-90,-60)" ; reg = "SPO" ; endif +if( xpos = 13 ) ; region = "X.Polar (Lats: -60,60)" ; reg = "XPO" ; endif filebeg = 1 diff --git a/GEOS_Util/plots/grads_util/corcmpz_plot.gs b/GEOS_Util/plots/grads_util/corcmpz_plot.gs index 461b5274..4b3b0067 100644 --- a/GEOS_Util/plots/grads_util/corcmpz_plot.gs +++ b/GEOS_Util/plots/grads_util/corcmpz_plot.gs @@ -119,6 +119,7 @@ if( xpos = 9 ) ; region = "North America (Lons:-140,-60 Lats: 20,60)" ; reg = if( xpos = 10 ) ; region = "Europe (Lons:-10,30 Lats: 30,60)" ; reg = "EUR" ; endif if( xpos = 11 ) ; region = "N.Polar (Lats: 60,90)" ; reg = "NPO" ; endif if( xpos = 12 ) ; region = "S.Polar (Lats: -60,-90)" ; reg = "SPO" ; endif +if( xpos = 13 ) ; region = "X.Polar (Lats: -60,60)" ; reg = "XPO" ; endif filebeg = 1 @@ -653,14 +654,44 @@ while( flag = '' ) 'set gxout shaded' 'rgbset' +* Define New Variables for Montage Plots +* ------------------------------------- +'set dfile 'ddif.m +'set t 'tbeg.m' 'tdif.m + +* For ravediff > 0 +* ---------------- +'define sigdiffp90 = 1000 * ( ravediff-rUp90diff )' +'define sigdiffp95 = 1000 * ( ravediff-rUp95diff )' +'define sigdiffp99 = 1000 * ( ravediff-rUp99diff )' + +* For ravediff < 0 +* ---------------- +'define sigdiffm90 = 1000 * ( ravediff+rUp90diff )' +'define sigdiffm95 = 1000 * ( ravediff+rUp95diff )' +'define sigdiffm99 = 1000 * ( ravediff+rUp99diff )' + +'define maskm90 = ( sigdiffm90 - abs(sigdiffm90) )/2' +'define maskm95 = ( sigdiffm95 - abs(sigdiffm95) )/2' +'define maskm99 = ( sigdiffm99 - abs(sigdiffm99) )/2' + +'define maskp90 = ( sigdiffp90 + abs(sigdiffp90) )/2' +'define maskp95 = ( sigdiffp95 + abs(sigdiffp95) )/2' +'define maskp99 = ( sigdiffp99 + abs(sigdiffp99) )/2' + +'define sigdiff90 = maskm90 + maskp90' +'define sigdiff95 = maskm95 + maskp95' +'define sigdiff99 = maskm99 + maskp99' + * Find maximum value of rUp90diff across all levels (at end of forecast period) * ----------------------------------------------------------------------------- -'set t 'tdif.m -'minmax rUp90diff' - dcint = subwrd(result,1) * 1000 / 6 + 'set t 'tdif.m + 'minmax rUp90diff' + rUp90diffmax = result +dcint = subwrd(rUp90diffmax,1) * 1000 / 6 -* Save dcint base on Global Region -* -------------------------------- +* Save dcint based on Global Region +* --------------------------------- if( xpos = 1 ) 'run setenv corcmp_'field'_DCINT 'dcint else @@ -668,12 +699,19 @@ else dcint = result endif -* WMP -* Fixed Color Scales to better judge magnitude of change... -* dcint = 7.5 - 'set dfile 'ddif.m 'set t 'tbeg.m' 'tdif.m + +' minmax sigdiff90' + qmax = subwrd(result,1) + qmin = subwrd(result,2) + qmax = math_abs(qmax) + qmin = math_abs(qmin) +if( qmin > qmax ) + qmax = qmin +endif + dcint = qmax / 9 + 'set csmooth on' 'set datawarn off' @@ -681,54 +719,41 @@ dcintx = dcint * 100 'getint 'dcintx dcintx = result / 100 -* Shade ravediff that is > 90% confidence diffs (grey shaded) -* ----------------------------------------------------------- -'set gxout shaded' -'shades 1.0' -'set csmooth on' -'d ravediff/rUp90diff' -'cbarn -xmid 6 -snum 0.80 -ndot 1' - -* Draw ravediff that is > 90% confidence diffs (thin black outline contours) -* -------------------------------------------------------------------------- -'shades 1' -'set gxout contour' -'run getenv SHADES_CLEVS' - clevs = result - ccols = "1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1" -say 'CCOLS = 'ccols -say 'CLEVS = 'clevs -'set ccols 'ccols -'set clevs 'clevs -'set clab off' -'set cthick 1' -'set csmooth on' -'d ravediff/rUp90diff' - -* Draw ravediff (grey contours without labels) -* -------------------------------------------- -'set gxout contour' -'set clab off' -'set cthick 1' -'set ccolor 99' -'set cint 'dcint -'set csmooth on' - -'black 'dcint - -'d ravediff*1000' -'q contours' -say 'Contours = 'result - -'set cthick 6' -'set clab off' -'set ccolor 99' -'set clevs 0' -'set csmooth on' -'d ravediff*1000' - -* Draw ravediff that is = 99% confidence diffs (black contour without label) -* -------------------------------------------------------------------------- +* Shade where sigdiff > 90% confidence error bar (color shaded) +* ------------------------------------------------------------- + ' rgbset' + ' set gxout grfill ' + ' set csmooth off' + ' shades 'dcint + ' run getenv SHADES_CLEVS' + SHADES_CLEVS = result + + clevs = subwrd( SHADES_CLEVS,1 ) + ii = 2 + while( ii <= 9 ) + clev = subwrd( SHADES_CLEVS,ii ) + clevs = clevs' 'clev + ii = ii + 1 + endwhile + clevs = clevs' -0.0001 0.0001 ' + ii = 10 + while( ii <= 18 ) + clev = subwrd( SHADES_CLEVS,ii ) + clevs = clevs' 'clev + ii = ii + 1 + endwhile + + 'set clevs 'clevs + 'set ccols 59 57 55 47 44 37 36 34 33 31 0 20 21 22 23 24 25 26 27 28 29' + +*' set gxout shaded ' + ' d sigdiff90 ' + ' cbarn -xmid 6 -snum 0.80 -ndot 1' + + + +* Contour sigdiff that is = 90, 95, & 99% confidence diffs (black lines without label) +* ------------------------------------------------------------------------------------ 'set gxout contour' 'set csmooth on' 'set clab off' @@ -736,83 +761,67 @@ say 'Contours = 'result 'set cstyle 2' 'set cthick 8' 'set ccolor 1' -'set clevs 1' -'d ravediff/rUp99diff' +'set clevs 0' +'d sigdiffp99' 'set cstyle 2' 'set cthick 8' 'set ccolor 1' -'set clevs -1' -'d ravediff/rUp99diff' +'set clevs 0' +'d sigdiffm99' 'set cstyle 6' 'set cthick 4' 'set ccolor 1' -'set clevs 1' -'d ravediff/rUp95diff' +'set clevs 0' +'d sigdiffp95' 'set cstyle 6' 'set cthick 4' 'set ccolor 1' -'set clevs -1' -'d ravediff/rUp95diff' +'set clevs 0' +'d sigdiffm95' 'set cstyle 1' 'set cthick 1' 'set ccolor 1' -'set clevs 1' -'d ravediff/rUp90diff' +'set clevs 0' +'d sigdiffp90' 'set cstyle 1' 'set cthick 1' 'set ccolor 1' -'set clevs -1' -'d ravediff/rUp90diff' +'set clevs 0' +'d sigdiffm90' * reset some background values +* ---------------------------- 'set datawarn on' 'set clab on' -* Re-Draw Individual Plots -* ------------------------ -'set gxout contour' -'set clab on' -'set xaxis 0 'nday' .5' -'set cmark 0' -'set cthick 6' -'set ccolor rainbow' -'set rbrange 'rbrange -'set cint 'cint -'set cstyle 2' -*'d ravem' -'set cmark 0' -'set cthick 6' -'set ccolor rainbow' -'set rbrange 'rbrange -'set cint 'cint -'set cstyle 3' -*'d raveo' - +* Draw Labels +* ----------- 'draw ylab Pressure (hPa)' 'set string 1 c 6 0' 'set strsiz .132' -'draw string 6.0 8.15 'expdsc.m' - 'expdsc.0' ('numfiles') 'region -'draw string 6.0 7.90 'name' Anomaly Correlation Difference (x10`a-3`n) CINT: 'dcintx +'draw string 6.0 8.15 'expdsc.m' - 'expdsc.0' ('numfiles') 'name' 'region +'draw string 6.0 7.90 Anomaly Correlation Difference (x10`a-3`n)' 'set strsiz .125' 'draw string 6.0 7.65 'desc 'set strsiz .12' 'draw string 6.0 0.72 Forecast Day' -'set string 1 l 8 0' -'set strsiz .32' 'run uppercase 'field FIELD = result + +'set string 1 l 8 0' +'set strsiz .32' 'draw string 0.70 7.38 'reg 'draw string 0.70 7.01 'FIELD 'set string 1 l 3 0' 'set strsiz .087' -'draw string 0.25 1.50 Color Shaded (>90%)' -'draw string 0.25 1.35 Dot -Dash Line (>95%)' -'draw string 0.25 1.20 Long-Dash Line (>99%)' +'draw string 0.23 1.50 Thin Solid Line (>90%)' +'draw string 0.23 1.35 Dot -Dash Line (>95%)' +'draw string 0.23 1.20 Long-Dash Line (>99%)' 'set string 1 c 6 90' 'set strsiz .18' diff --git a/GEOS_Util/plots/grads_util/corz.gs b/GEOS_Util/plots/grads_util/corz.gs index 80448ea0..ec0f0c96 100644 --- a/GEOS_Util/plots/grads_util/corz.gs +++ b/GEOS_Util/plots/grads_util/corz.gs @@ -79,6 +79,7 @@ if( xpos = 9 ) ; region = "North America (Lons:-140,-60 Lats: 20,60)" ; reg = if( xpos = 10 ) ; region = "Europe (Lons:-10,30 Lats: 30,60)" ; reg = "EUR" ; endif if( xpos = 11 ) ; region = "N.Polar (Lats: 60,90)" ; reg = "NPO" ; endif if( xpos = 12 ) ; region = "S.Polar (Lats: -60,-90)" ; reg = "SPO" ; endif +if( xpos = 13 ) ; region = "X.Polar (Lats: -60,60)" ; reg = "XPO" ; endif 'setz' 'getinfo tdim' diff --git a/GEOS_Util/plots/grads_util/genplot.gs b/GEOS_Util/plots/grads_util/genplot.gs index 91cf4031..6ec88147 100644 --- a/GEOS_Util/plots/grads_util/genplot.gs +++ b/GEOS_Util/plots/grads_util/genplot.gs @@ -739,6 +739,108 @@ while( k <= numrc ) if(RC != 'NULL') ; RCFILE = geosutil'/plots/'NAME'/VERIFICATION.'RC'.rc' ; endif if( RCFILE = rcfile ) + say 'k = 'k' RCFILE = 'RCFILE + +* Check for VERIFICATION Formula File +* ----------------------------------- +'!remove OBSFORM.txt' +'!echo `basename 'RCFILE' | cut -d. -f2 > OBSFORM.txt`' +'run getenv OBSFORM' + OBSNAME = result + OBSFORM = OBSNAME'form.gs' + +say 'Checking for OBSFORM: 'geosutil'/plots/'NAME'/'OBSFORM + +'!remove CHECKFILE.txt' +'!chckfile 'geosutil'/plots/'NAME'/'OBSFORM + 'run getenv CHECKFILE' + CHECKFILE = result + +say 'CHECKFILE = 'CHECKFILE +if( CHECKFILE != 'NULL' ) + + LOBSFORM = 'TRUE' + 'run 'geosutil'/plots/'NAME'/'OBSFORM' -OBS' + OBS = result + say 'OBSFORM OBS: 'OBS + 'numargs 'OBS + numobs = result + say ' NUMOBS: 'numobs + +* Read Verification OBS:OBSGC +* --------------------------- + n = 1 + OBSNAME.n = subwrd(OBS,n ) + next = subwrd(OBS,n+1 ) + bit = checkbit(next) + say 'OBSNAME.'n': 'OBSNAME.n + while( bit != '' ) + n = n + 1 + OBSNAME.n = subwrd(OBS,n ) + next = subwrd(OBS,n+1 ) + bit = checkbit(next) + say 'OBSNAME.'n': 'OBSNAME.n + endwhile + +* Construct OBS GCs from OBS EXPORTS +* ---------------------------------- + m = 0 + n = 1 +while ( n <= numobs ) + EX = '' + j = 1 + bit = substr(OBSNAME.n,j,1) + while(bit != ':' & bit != '') + EX = EX''bit + j = j + 1 + bit = substr(OBSNAME.n,j,1) + endwhile + if( EX != OBSNAME.n ) + m = m + 1 + j = j + 1 + OBSNAMEGC.m = '' + bit = substr(OBSNAME.n,j,1) + while(bit != '') + OBSNAMEGC.m = OBSNAMEGC.m''bit + j = j + 1 + bit = substr(OBSNAME.n,j,1) + endwhile + OBSNAME.n = EX + endif +n = n + 1 +endwhile + +* Get Verification Variables +* -------------------------- +FOUND = TRUE + n = 1 +while ( n <= numobs & FOUND = TRUE ) +'run getobs 'OBSNAME.n' 'OBSNAMEGC.n' 'rcfile + oname.n = subwrd(result,1) + ofile.n = subwrd(result,2) + oscal.n = subwrd(result,3) + odesc.n = subwrd(result,4) + otag.n = subwrd(result,5) + +say 'VERIFICATION_EXPORT_name: 'oname.n +say ' description: 'odesc.n +say ' tag: ' otag.n +say ' file: 'ofile.n +say ' scaling: 'oscal.n +say '' + if( oname.n = 'NULL' ) ; FOUND = FALSE ; endif + n = n + 1 +endwhile + +* Save Original numn for other RCs +* -------------------------------- + numnorig = numn + numn = numobs + +else + + LOBSFORM = 'FALSE' + OBSFORM = 'obsform.gs' * Get Verification Variables * -------------------------- @@ -762,6 +864,8 @@ say '' n = n + 1 endwhile +endif + if( FOUND = TRUE ) 'setlons' 'setlats' @@ -798,7 +902,7 @@ endif if( numn = 1 ) 'define qobs = 'oname.1'.'ofile.1'*'oscal.1 else - filename = geosutil'/plots/'NAME'/obsform.gs' + filename = geosutil'/plots/'NAME'/'OBSFORM ioflag = sublin( read(filename),1 ) if(ioflag = 0) close = close(filename) @@ -808,7 +912,7 @@ else ostring = ostring' 'oname.n'.'ofile.n'*'oscal.n n = n + 1 endwhile - 'run 'geosutil'/plots/'NAME'/obsform 'ostring + 'run 'geosutil'/plots/'NAME'/'OBSFORM' 'ostring else ostring = NAME n = 1 @@ -996,6 +1100,12 @@ endif * ------------------------------ k = k + 1 +* Restore Original numn for other RCs +* ----------------------------------- +if( LOBSFORM = 'TRUE' ) + numn = numnorig +endif + * End Verification Loop * --------------------- endwhile diff --git a/GEOS_Util/plots/grads_util/rms.gs b/GEOS_Util/plots/grads_util/rms.gs index f065c0c6..e1e71f5b 100644 --- a/GEOS_Util/plots/grads_util/rms.gs +++ b/GEOS_Util/plots/grads_util/rms.gs @@ -101,6 +101,7 @@ if( xpos = 9 ) ; region = "North America (Lons:-140,-60 Lats: 20,60)" ; reg = if( xpos = 10 ) ; region = "Europe (Lons:-10,30 Lats: 30,60)" ; reg = "EUR" ; endif if( xpos = 11 ) ; region = "N.Polar (Lats: 60,90)" ; reg = "NPO" ; endif if( xpos = 12 ) ; region = "S.Polar (Lats: -60,-90)" ; reg = "SPO" ; endif +if( xpos = 13 ) ; region = "X.Polar (Lats: -60,60)" ; reg = "XPO" ; endif 'getinfo tdim' tdim = result diff --git a/GEOS_Util/plots/grads_util/rmscmp_plot.gs b/GEOS_Util/plots/grads_util/rmscmp_plot.gs index fcb9b10b..761827cd 100644 --- a/GEOS_Util/plots/grads_util/rmscmp_plot.gs +++ b/GEOS_Util/plots/grads_util/rmscmp_plot.gs @@ -154,7 +154,8 @@ if( xpos = 8 ) ; region = "S.E. Quadrant (Lons: 0,180 Lats: 0,-90)" ; reg = if( xpos = 9 ) ; region = "North America (Lons:-140,-60 Lats: 20,60)" ; reg = "NAM" ; endif if( xpos = 10 ) ; region = "Europe (Lons:-10,30 Lats: 30,60)" ; reg = "EUR" ; endif if( xpos = 11 ) ; region = "N.Polar (Lats: 60,90)" ; reg = "NPO" ; endif -if( xpos = 12 ) ; region = "S.Polar (Lats: -60,-90)" ; reg = "SPO" ; endif +if( xpos = 12 ) ; region = "S.Polar (Lats:-90,-60)" ; reg = "SPO" ; endif +if( xpos = 13 ) ; region = "X.Polar (Lats: -60,60)" ; reg = "XPO" ; endif filebeg = 1 diff --git a/GEOS_Util/plots/grads_util/rmscmpz_plot.gs b/GEOS_Util/plots/grads_util/rmscmpz_plot.gs index 6982e666..5e08688c 100644 --- a/GEOS_Util/plots/grads_util/rmscmpz_plot.gs +++ b/GEOS_Util/plots/grads_util/rmscmpz_plot.gs @@ -136,6 +136,7 @@ if( xpos = 9 ) ; region = "North America (Lons:-140,-60 Lats: 20,60)" ; reg = if( xpos = 10 ) ; region = "Europe (Lons:-10,30 Lats: 30,60)" ; reg = "EUR" ; endif if( xpos = 11 ) ; region = "N.Polar (Lats: 60,90)" ; reg = "NPO" ; endif if( xpos = 12 ) ; region = "S.Polar (Lats: -60,-90)" ; reg = "SPO" ; endif +if( xpos = 13 ) ; region = "X.Polar (Lats: -60,60)" ; reg = "XPO" ; endif filebeg = 1 @@ -740,14 +741,44 @@ while( flag = '' ) 'set gxout shaded' 'rgbset' +* Define New Variables for Montage Plots +* ------------------------------------- +'set dfile 'ddif.m +'set t 'tbeg.m' 'tdif.m + +* For ravediff > 0 +* ---------------- +'define sigdiffp90 = 1000 * ( ravediff-rUp90diff )' +'define sigdiffp95 = 1000 * ( ravediff-rUp95diff )' +'define sigdiffp99 = 1000 * ( ravediff-rUp99diff )' + +* For ravediff < 0 +* ---------------- +'define sigdiffm90 = 1000 * ( ravediff+rUp90diff )' +'define sigdiffm95 = 1000 * ( ravediff+rUp95diff )' +'define sigdiffm99 = 1000 * ( ravediff+rUp99diff )' + +'define maskm90 = ( sigdiffm90 - abs(sigdiffm90) )/2' +'define maskm95 = ( sigdiffm95 - abs(sigdiffm95) )/2' +'define maskm99 = ( sigdiffm99 - abs(sigdiffm99) )/2' + +'define maskp90 = ( sigdiffp90 + abs(sigdiffp90) )/2' +'define maskp95 = ( sigdiffp95 + abs(sigdiffp95) )/2' +'define maskp99 = ( sigdiffp99 + abs(sigdiffp99) )/2' + +'define sigdiff90 = maskm90 + maskp90' +'define sigdiff95 = maskm95 + maskp95' +'define sigdiff99 = maskm99 + maskp99' + * Find maximum value of rUp90diff across all levels (at end of forecast period) * ----------------------------------------------------------------------------- -'set t 'tdif.m -'minmax rUp90diff' - dcint = subwrd(result,1) * 1000 / 6 + 'set t 'tdif.m + 'minmax rUp90diff' + rUp90diffmax = result +dcint = subwrd(rUp90diffmax,1) * 1000 / 6 -* Save dcint base on Global Region -* -------------------------------- +* Save dcint based on Global Region +* --------------------------------- if( xpos = 1 ) 'run setenv rmscmp_'field'_DCINT 'dcint else @@ -755,73 +786,56 @@ else dcint = result endif -* WMP -* Fixed Color Scales to better judge magnitude of change... -'run uppercase 'field - FIELD = result -*if (FIELD='H') ; dcint = 750.000 ; endif -*if (FIELD='U') ; dcint = 90.000 ; endif -*if (FIELD='V') ; dcint = 90.000 ; endif -*if (FIELD='T') ; dcint = 30.000 ; endif -*if (FIELD='Q') ; dcint = 0.012 ; endif - 'set dfile 'ddif.m 'set t 'tbeg.m' 'tdif.m -'set csmooth on' -'set datawarn off' - -dcintx = dcint * 100 -'getint 'dcintx - dcintx = result / 100 - -* Shade ravediff that is > 90% confidence diffs (grey shaded) -* ----------------------------------------------------------- -'set gxout shaded' -'shades 1.0' -'set csmooth on' -'d ravediff/rUp90diff' -'cbarn -xmid 6 -snum 0.80 -ndot 1' -* Draw ravediff that is > 90% confidence diffs (thin black outline contours) -* -------------------------------------------------------------------------- -'shades 1' -'set gxout contour' -'run getenv SHADES_CLEVS' - clevs = result - ccols = "1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1" -say 'CCOLS = 'ccols -say 'CLEVS = 'clevs -'set ccols 'ccols -'set clevs 'clevs -'set clab off' -'set cthick 1' -'set csmooth on' -'d ravediff/rUp90diff' +' minmax sigdiff90' + qmax = subwrd(result,1) + qmin = subwrd(result,2) + qmax = math_abs(qmax) + qmin = math_abs(qmin) +if( qmin > qmax ) + qmax = qmin +endif + dcint = qmax / 9 -* Draw ravediff (grey contours without labels) -* -------------------------------------------- -'set gxout contour' -'set clab off' -'set cthick 1' -'set ccolor 99' -'set cint 'dcint 'set csmooth on' +'set datawarn off' -'black 'dcint - -'d ravediff*1000' -'q contours' -say 'Contours = 'result - -'set cthick 6' -'set clab off' -'set ccolor 99' -'set clevs 0' -'set csmooth on' -'d ravediff*1000' -* Draw ravediff that is = 99% confidence diffs (black contour without label) -* -------------------------------------------------------------------------- +* Shade where sigdiff > 90% confidence error bar (color shaded) +* ------------------------------------------------------------- + ' rgbset' + ' set gxout grfill ' + ' set csmooth off' + ' shades 'dcint + ' run getenv SHADES_CLEVS' + SHADES_CLEVS = result + + clevs = subwrd( SHADES_CLEVS,1 ) + ii = 2 + while( ii <= 9 ) + clev = subwrd( SHADES_CLEVS,ii ) + clevs = clevs' 'clev + ii = ii + 1 + endwhile + clevs = clevs' -0.0001 0.0001 ' + ii = 10 + while( ii <= 18 ) + clev = subwrd( SHADES_CLEVS,ii ) + clevs = clevs' 'clev + ii = ii + 1 + endwhile + + 'set clevs 'clevs + 'set ccols 59 57 55 47 44 37 36 34 33 31 0 20 21 22 23 24 25 26 27 28 29' + +*' set gxout shaded ' + ' d sigdiff90 ' + ' cbarn -xmid 6 -snum 0.80 -ndot 1' + +* Contour sigdiff that is = 90, 95, & 99% confidence diffs (black lines without label) +* ------------------------------------------------------------------------------------ 'set gxout contour' 'set csmooth on' 'set clab off' @@ -829,60 +843,43 @@ say 'Contours = 'result 'set cstyle 2' 'set cthick 8' 'set ccolor 1' -'set clevs 1' -'d ravediff/rUp99diff' +'set clevs 0' +'d sigdiffp99' 'set cstyle 2' 'set cthick 8' 'set ccolor 1' -'set clevs -1' -'d ravediff/rUp99diff' +'set clevs 0' +'d sigdiffm99' 'set cstyle 6' 'set cthick 4' 'set ccolor 1' -'set clevs 1' -'d ravediff/rUp95diff' +'set clevs 0' +'d sigdiffp95' 'set cstyle 6' 'set cthick 4' 'set ccolor 1' -'set clevs -1' -'d ravediff/rUp95diff' +'set clevs 0' +'d sigdiffm95' 'set cstyle 1' 'set cthick 1' 'set ccolor 1' -'set clevs 1' -'d ravediff/rUp90diff' +'set clevs 0' +'d sigdiffp90' 'set cstyle 1' 'set cthick 1' 'set ccolor 1' -'set clevs -1' -'d ravediff/rUp90diff' +'set clevs 0' +'d sigdiffm90' * reset some background values +* ---------------------------- 'set datawarn on' 'set clab on' -* Re-Draw Individual Plots -* ------------------------ -'set gxout contour' -'set clab on' -'set xaxis 0 'nday' .5' -'set cmark 0' -'set cthick 6' -'set ccolor rainbow' -'set rbrange 0 'qmax -'set cint 'qcnt -'set cstyle 2' -*'d ravem' -'set cmark 0' -'set cthick 6' -'set ccolor rainbow' -'set rbrange 0 'qmax -'set cint 'qcnt -'set cstyle 3' -*'d raveo' - +* Draw Labels +* ----------- 'draw ylab Pressure (hPa)' 'set string 1 c 6 0' @@ -892,14 +889,21 @@ say 'Contours = 'result if( rms = 3 ) ; rms_label = 'RMS for AMPLITUDE Error Difference' ; endif if( rms = 4 ) ; rms_label = 'RMS for PHASE Error Difference' ; endif +dcintx = dcint * 100 +'getint 'dcintx + dcintx = result / 100 + 'set strsiz .132' -'draw string 6.0 8.15 'expdsc.m' - 'expdsc.0' ('numfiles') 'region -'draw string 6.0 7.90 'name' 'rms_label' (x10`a-3`n) CINT: 'dcintx +'draw string 6.0 8.15 'expdsc.m' - 'expdsc.0' ('numfiles') 'name' 'region +'draw string 6.0 7.90 'rms_label' (x10`a-3`n)' 'set strsiz .125' 'draw string 6.0 7.65 'desc 'set strsiz .12' 'draw string 6.0 0.72 Forecast Day' +'run uppercase 'field + FIELD = result + 'set string 1 l 8 0' 'set strsiz .32' 'draw string 0.70 7.38 'reg @@ -907,9 +911,9 @@ say 'Contours = 'result 'set string 1 l 3 0' 'set strsiz .087' -'draw string 0.25 1.50 Color Shaded (>90%)' -'draw string 0.25 1.35 Dot -Dash Line (>95%)' -'draw string 0.25 1.20 Long-Dash Line (>99%)' +'draw string 0.23 1.50 Thin Solid Line (>90%)' +'draw string 0.23 1.35 Dot -Dash Line (>95%)' +'draw string 0.23 1.20 Long-Dash Line (>99%)' 'set string 1 c 6 90' 'set strsiz .18' diff --git a/GEOS_Util/plots/grads_util/rmsz.gs b/GEOS_Util/plots/grads_util/rmsz.gs index 4546b5bb..9b95eff6 100644 --- a/GEOS_Util/plots/grads_util/rmsz.gs +++ b/GEOS_Util/plots/grads_util/rmsz.gs @@ -79,6 +79,7 @@ if( xpos = 9 ) ; region = "North America (Lons:-140,-60 Lats: 20,60)" ; reg = if( xpos = 10 ) ; region = "Europe (Lons:-10,30 Lats: 30,60)" ; reg = "EUR" ; endif if( xpos = 11 ) ; region = "N.Polar (Lats: 60,90)" ; reg = "NPO" ; endif if( xpos = 12 ) ; region = "S.Polar (Lats: -60,-90)" ; reg = "SPO" ; endif +if( xpos = 13 ) ; region = "X.Polar (Lats: -60,60)" ; reg = "XPO" ; endif 'setz' 'getinfo tdim' diff --git a/GEOS_Util/plots/grads_util/shades.gs b/GEOS_Util/plots/grads_util/shades.gs index 2ef814f8..4473d762 100644 --- a/GEOS_Util/plots/grads_util/shades.gs +++ b/GEOS_Util/plots/grads_util/shades.gs @@ -186,10 +186,11 @@ if( fixed = '' | fixed = '-quad' ) 'set ccols 59 57 55 47 44 37 36 34 33 0 21 22 23 24 25 26 27 28 29' say 'Inside Shades1, Quad CLEVS: 'clevs cint = cint2 + endif + 'run setenv SHADES_CLEVS 'clevs 'run setenv SHADES_CCOLS 59 57 55 47 44 37 36 34 33 0 21 22 23 24 25 26 27 28 29' 'run setenv SHADES_CINT 'cint - endif return cint diff --git a/GEOS_Util/plots/grads_util/stats_montage b/GEOS_Util/plots/grads_util/stats_montage index e9a9e6d5..799f958d 100755 --- a/GEOS_Util/plots/grads_util/stats_montage +++ b/GEOS_Util/plots/grads_util/stats_montage @@ -23,7 +23,7 @@ echo EXPS: $EXPS set STATSDIR = corcmp cd $STATSDIR -set files = `/bin/ls -1 *_*_stats_uwnd_corcmp_GLO_z_*.gif` +set files = `/bin/ls -1 *_*_stats_*_*cmp_GLO_z_*.gif` if( $#files >= 1 ) then @@ -31,7 +31,7 @@ if( $#files >= 1 ) then # ----------------------- @ n = 0 set node = '' -while( .$node != '.corcmp' ) +while( (.$node != '.corcmp') & (.$node != '.rmscmp') ) @ n = $n + 1 set node = `echo $files[1] | cut -d_ -f$n` end @@ -73,10 +73,10 @@ foreach EXP ($EXPS) foreach field ('hght' 'uwnd' 'vwnd' 'tmpu' 'sphu' ) montage ${DSC0}_${EXP}_stats_${field}_${VAR}_SHE_z_${SEASON}.gif \ - ${DSC0}_${EXP}_stats_${field}_${VAR}_GLO_z_${SEASON}.gif \ + ${DSC0}_${EXP}_stats_${field}_${VAR}_TRO_z_${SEASON}.gif \ ${DSC0}_${EXP}_stats_${field}_${VAR}_NHE_z_${SEASON}.gif \ ${DSC0}_${EXP}_stats_${field}_${VAR}_SPO_z_${SEASON}.gif \ - ${DSC0}_${EXP}_stats_${field}_${VAR}_TRO_z_${SEASON}.gif \ + ${DSC0}_${EXP}_stats_${field}_${VAR}_XPO_z_${SEASON}.gif \ ${DSC0}_${EXP}_stats_${field}_${VAR}_NPO_z_${SEASON}.gif \ -tile 3x2 -geometry +0+0 ${DSC0}_${EXP}_${field}_${VAR}_${SEASON}_montage.gif & end diff --git a/GEOS_Util/plots/grads_util/stats_rms_montage b/GEOS_Util/plots/grads_util/stats_rms_montage index 54feeb5e..204abe2f 100755 --- a/GEOS_Util/plots/grads_util/stats_rms_montage +++ b/GEOS_Util/plots/grads_util/stats_rms_montage @@ -4,7 +4,7 @@ set STATSDIR = corcmp cd $STATSDIR set fields = `echo slp uwnd vwnd tmpu hght sphu chi psi` -set regions = `echo GLO NHE TRO SHE NWQ NEQ SWQ SEQ NAM EUR NPO SPO` +set regions = `echo GLO NHE TRO SHE NWQ NEQ SWQ SEQ NAM EUR NPO SPO XPO` foreach field ($fields) foreach region ($regions) diff --git a/GEOS_Util/plots/grads_util/tcor.gs b/GEOS_Util/plots/grads_util/tcor.gs index 162aa68e..bdedb2e6 100644 --- a/GEOS_Util/plots/grads_util/tcor.gs +++ b/GEOS_Util/plots/grads_util/tcor.gs @@ -157,6 +157,7 @@ if( xpos = 9 ) ; region = "North America (Lons:-140,-60 Lats: 20,60)" ; reg = if( xpos = 10 ) ; region = "Europe (Lons:-10,30 Lats: 30,60)" ; reg = "EUR" ; endif if( xpos = 11 ) ; region = "N.Polar (Lats: 60,90)" ; reg = "NPO" ; endif if( xpos = 12 ) ; region = "S.Polar (Lats: -60,-90)" ; reg = "SPO" ; endif +if( xpos = 13 ) ; region = "X.Polar (Lats: -60,60)" ; reg = "XPO" ; endif if( fday = 5 ) diff --git a/GEOS_Util/plots/hcmp/VERIFICATION.ECINTERIM.rc b/GEOS_Util/plots/hcmp/VERIFICATION.ECINTERIM.rc deleted file mode 100644 index 78330b32..00000000 --- a/GEOS_Util/plots/hcmp/VERIFICATION.ECINTERIM.rc +++ /dev/null @@ -1,19 +0,0 @@ -OBSNAM: ecinterim -OBSDSC: ECMWF_INTERIM_ANA - - filename: 'VERIFICATION/ECMWF_Interim_FullRes/ecint_prs.ctl' - filetype: 'BIN' - fields: 'ZLE' , 'DYN' , 'Z' , 1/9.80665 - 'T' , 'DYN' , 'T' , 1 - 'U' , 'DYN' , 'U' , 1 - 'V' , 'DYN' , 'V' , 1 - 'RH2' , 'MOIST' , 'R' , 1/100 - 'Q' , 'MOIST' , 'Q' , 1 - 'OMEGA' , 'DYN' , 'W' , 1 - :: - - filename: 'VERIFICATION/ECMWF_Interim_FullRes/surface/MSL_197901_201111.nc' - filetype: 'SDF' - fields: 'SLP' , 'DYN' , 'MSL' , 1 - :: - diff --git a/GEOS_Util/plots/hflux/VERIFICATION.COADS.rc b/GEOS_Util/plots/hflux/VERIFICATION.COADS.rc deleted file mode 100644 index 0319883e..00000000 --- a/GEOS_Util/plots/hflux/VERIFICATION.COADS.rc +++ /dev/null @@ -1,8 +0,0 @@ - -OBSNAM: COADS -OBSDSC: COADS_CLIMATOLOGY_(1945-1989) - - filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' - filetype: 'BIN' - fields: 'SH' , 'SURFACE' , 'SENSIB3' , 1 - :: diff --git a/GEOS_Util/plots/hflux/VERIFICATION.GSSTF.rc b/GEOS_Util/plots/hflux/VERIFICATION.GSSTF.rc deleted file mode 100644 index 2506871e..00000000 --- a/GEOS_Util/plots/hflux/VERIFICATION.GSSTF.rc +++ /dev/null @@ -1,14 +0,0 @@ - -OBSNAM: GSSTF -OBSDSC: Goddard_V2_Sat-based_Surf_Turb_Fluxes - - filename: 'VERIFICATION/GSSTF_V2_Surf_Flux/testv2b.ctl' - filetype: 'BIN' - fields: 'VENT' , 'SURFACE' , 'S10' , 1 , - 'LHFX' , 'SURFACE' , 'EFLUX' , 1 , - 'SH' , 'SURFACE' , 'HFLUX' , 1 , - 'TAUX' , 'SURFACE' , 'UFLUX' , 100 , - 'TAUY' , 'SURFACE' , 'VFLUX' , 100 , - 'TQV' , 'AGCM' , 'TPWA' , 1 , - 'Q10M' , 'SURFACE' , 'Q10M' , 0.001 , - :: diff --git a/GEOS_Util/plots/landscape.script b/GEOS_Util/plots/landscape.script index 44518408..57ecc6ca 100755 --- a/GEOS_Util/plots/landscape.script +++ b/GEOS_Util/plots/landscape.script @@ -367,6 +367,7 @@ end if( $plot == $plotname.$n ) $grads $batch -l -c "run template -EXPID $expid -EXPORT RSRS:SOLAR RSRSNA:SOLAR -NAME SWGNETAF -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -l -c "run template -EXPID $expid -EXPORT RSCS:SOLAR RSCSNA:SOLAR -NAME SWGNETCAF -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -l -c "run template -EXPID $expid -EXPORT ALBEDO:SOLAR -SEASON $seasons -OUTPUT $output -DEBUG $debug " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -l -c "run template -EXPID $expid -EXPORT ALBVR:SOLAR -SEASON $seasons -OUTPUT $output -DEBUG $debug " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -l -c "run template -EXPID $expid -EXPORT ALBVF:SOLAR -SEASON $seasons -OUTPUT $output -DEBUG $debug " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -l -c "run template -EXPID $expid -EXPORT ALBNR:SOLAR -SEASON $seasons -OUTPUT $output -DEBUG $debug " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 diff --git a/GEOS_Util/plots/lcmp/VERIFICATION.ECANA.rc b/GEOS_Util/plots/lcmp/VERIFICATION.ECANA.rc deleted file mode 100644 index ae4df4db..00000000 --- a/GEOS_Util/plots/lcmp/VERIFICATION.ECANA.rc +++ /dev/null @@ -1,13 +0,0 @@ - -OBSNAM: ecops -OBSDSC: ECMWF_OPNL_ANA - - filename: 'VERIFICATION/ECMWF_OPS.NCDF3/xdf.tabl' - filetype: 'HDF' - fields: 'ZLE' , 'DYN' , 'HGHT' , 1 - 'SLP' , 'DYN' , 'SLP' , 1 - 'T' , 'DYN' , 'T' , 1 - 'U' , 'DYN' , 'U' , 1 - 'V' , 'DYN' , 'V' , 1 - 'RH2' , 'MOIST' , 'RH' , 1/100 - :: diff --git a/GEOS_Util/plots/lcmp/VERIFICATION.GDAS1.rc b/GEOS_Util/plots/lcmp/VERIFICATION.GDAS1.rc deleted file mode 100644 index 2b0c7545..00000000 --- a/GEOS_Util/plots/lcmp/VERIFICATION.GDAS1.rc +++ /dev/null @@ -1,13 +0,0 @@ - -OBSNAM: gdas1 -OBSDSC: NCEP_OPNL_gdas1 - - filename: 'VERIFICATION/NCEP_GDAS-1.NCDF3/xdf.tabl' - filetype: 'HDF' - fields: 'ZLE' , 'DYN' , 'HGHT' , 1 - 'SLP' , 'DYN' , 'SLP' , 1 - 'T' , 'DYN' , 'TMPU' , 1 - 'U' , 'DYN' , 'UWND' , 1 - 'V' , 'DYN' , 'VWND' , 1 - 'RH2' , 'MOIST' , 'RH' , 1/100 - :: diff --git a/GEOS_Util/plots/lcmp/VERIFICATION.NCEP.rc b/GEOS_Util/plots/lcmp/VERIFICATION.NCEP.rc deleted file mode 100644 index 1a015ade..00000000 --- a/GEOS_Util/plots/lcmp/VERIFICATION.NCEP.rc +++ /dev/null @@ -1,20 +0,0 @@ - -OBSNAM: ncep -OBSDSC: NCEP_REANAL_Ver-2_Jan_1982-Dec_1998 - - filename: 'VERIFICATION/NCEP-R2/pgb.anl.197901.avrg.grib.ctl' - filetype: 'BIN' - fields: 'ZLE' , 'DYN' , 'HGHT' , 1 - 'SLP' , 'DYN' , 'SLP' , 1 - 'T' , 'DYN' , 'TMPU' , 1 - 'U' , 'DYN' , 'UWND' , 1 - 'V' , 'DYN' , 'VWND' , 1 - 'OMEGA' , 'DYN' , 'OMEGA' , 1 - :: - - filename: 'VERIFICATION/NCEP_REANAL-2/ncep_sphu.ctl' - filetype: 'BIN' - fields: 'RH2' , 'MOIST' , 'RH' , 1/100 - 'Q' , 'MOIST' , 'SPHU' , 1 - :: - diff --git a/GEOS_Util/plots/lwgdwn/ERA5form.gs b/GEOS_Util/plots/lwgdwn/ERA5form.gs new file mode 100644 index 00000000..b7b46e38 --- /dev/null +++ b/GEOS_Util/plots/lwgdwn/ERA5form.gs @@ -0,0 +1,23 @@ +function formula (args) + +'numargs 'args + numargs = result + +* OBS EXPORTS & GCS +* ----------------- +if( subwrd(args,1) = '-OBS' ) + return "MSDWLWRF:ERA5" +endif + +* Formula: lwgdwn = lwgdwn +* --------------------------- + n = 1 +while( n <= numargs ) + oname.n = subwrd(args,n) + n = n + 1 +endwhile + +say 'define qobs = 'oname.1 + 'define qobs = 'oname.1 + +return diff --git a/GEOS_Util/plots/lwgdwn/VERIFICATION.SRB.rc b/GEOS_Util/plots/lwgdwn/VERIFICATION.SRB.rc deleted file mode 100644 index c7b3cfbe..00000000 --- a/GEOS_Util/plots/lwgdwn/VERIFICATION.SRB.rc +++ /dev/null @@ -1,33 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_Longwave_Monthly_(198307-200506) - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - :: diff --git a/GEOS_Util/plots/lwgdwnc/ERA5form.gs b/GEOS_Util/plots/lwgdwnc/ERA5form.gs new file mode 100644 index 00000000..ca40f05d --- /dev/null +++ b/GEOS_Util/plots/lwgdwnc/ERA5form.gs @@ -0,0 +1,23 @@ +function formula (args) + +'numargs 'args + numargs = result + +* OBS EXPORTS & GCS +* ----------------- +if( subwrd(args,1) = '-OBS' ) + return "MSDWLWRFCS:ERA5" +endif + +* Formula: lwgdwnc = lwgdwnc +* --------------------------- + n = 1 +while( n <= numargs ) + oname.n = subwrd(args,n) + n = n + 1 +endwhile + +say 'define qobs = 'oname.1 + 'define qobs = 'oname.1 + +return diff --git a/GEOS_Util/plots/lwgdwnc/VERIFICATION.SRB.rc b/GEOS_Util/plots/lwgdwnc/VERIFICATION.SRB.rc deleted file mode 100644 index 97ca40e0..00000000 --- a/GEOS_Util/plots/lwgdwnc/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_Longwave_Monthly_(198307-200506) - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - :: diff --git a/GEOS_Util/plots/lwgdwncf/VERIFICATION.SRB.rc b/GEOS_Util/plots/lwgdwncf/VERIFICATION.SRB.rc deleted file mode 100644 index 97ca40e0..00000000 --- a/GEOS_Util/plots/lwgdwncf/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_Longwave_Monthly_(198307-200506) - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - :: diff --git a/GEOS_Util/plots/lwgnet/ERA5form.gs b/GEOS_Util/plots/lwgnet/ERA5form.gs new file mode 100644 index 00000000..3c0ac5a4 --- /dev/null +++ b/GEOS_Util/plots/lwgnet/ERA5form.gs @@ -0,0 +1,23 @@ +function formula (args) + +'numargs 'args + numargs = result + +* OBS EXPORTS & GCS +* ----------------- +if( subwrd(args,1) = '-OBS' ) + return "MSNLWRF:ERA5" +endif + +* Formula: lwgnet = lwgnet +* --------------------------- + n = 1 +while( n <= numargs ) + oname.n = subwrd(args,n) + n = n + 1 +endwhile + +say 'define qobs = 'oname.1 + 'define qobs = 'oname.1 + +return diff --git a/GEOS_Util/plots/lwgnet/VERIFICATION.SRB.rc b/GEOS_Util/plots/lwgnet/VERIFICATION.SRB.rc deleted file mode 100644 index c7b3cfbe..00000000 --- a/GEOS_Util/plots/lwgnet/VERIFICATION.SRB.rc +++ /dev/null @@ -1,33 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_Longwave_Monthly_(198307-200506) - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - :: diff --git a/GEOS_Util/plots/lwgnetc/ERA5form.gs b/GEOS_Util/plots/lwgnetc/ERA5form.gs new file mode 100644 index 00000000..3b84f55a --- /dev/null +++ b/GEOS_Util/plots/lwgnetc/ERA5form.gs @@ -0,0 +1,23 @@ +function formula (args) + +'numargs 'args + numargs = result + +* OBS EXPORTS & GCS +* ----------------- +if( subwrd(args,1) = '-OBS' ) + return "MSNLWRFCS:ERA5" +endif + +* Formula: lwgnetc = lwgnetc +* --------------------------- + n = 1 +while( n <= numargs ) + oname.n = subwrd(args,n) + n = n + 1 +endwhile + +say 'define qobs = 'oname.1 + 'define qobs = 'oname.1 + +return diff --git a/GEOS_Util/plots/lwgnetc/VERIFICATION.SRB.rc b/GEOS_Util/plots/lwgnetc/VERIFICATION.SRB.rc deleted file mode 100644 index 97ca40e0..00000000 --- a/GEOS_Util/plots/lwgnetc/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_Longwave_Monthly_(198307-200506) - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - :: diff --git a/GEOS_Util/plots/lwgnetcf/VERIFICATION.SRB.rc b/GEOS_Util/plots/lwgnetcf/VERIFICATION.SRB.rc deleted file mode 100644 index 97ca40e0..00000000 --- a/GEOS_Util/plots/lwgnetcf/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_Longwave_Monthly_(198307-200506) - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - :: diff --git a/GEOS_Util/plots/lwgup/ERA5form.gs b/GEOS_Util/plots/lwgup/ERA5form.gs new file mode 100644 index 00000000..9e0ad497 --- /dev/null +++ b/GEOS_Util/plots/lwgup/ERA5form.gs @@ -0,0 +1,24 @@ +function formula (args) + +'numargs 'args + numargs = result + +* OBS EXPORTS & GCS +* ----------------- +if( subwrd(args,1) = '-OBS' ) + return " MSDWLWRF:ERA5 MSNLWRF:ERA5" +endif + +* Formula: lwgup = msdwlwrf - (-1 * msnlwrf) +* Import of msnlwrf is scaled by -1 in VERIFICATION.ERA5.rc for comparison with lwgnet +* ------------------------------------ + n = 1 +while( n <= numargs ) + oname.n = subwrd(args,n) + n = n + 1 +endwhile + +say 'define qobs = 'oname.1' + 'oname.2 + 'define qobs = 'oname.1' + 'oname.2 + +return diff --git a/GEOS_Util/plots/lwgup/VERIFICATION.SRB.rc b/GEOS_Util/plots/lwgup/VERIFICATION.SRB.rc deleted file mode 100644 index c7b3cfbe..00000000 --- a/GEOS_Util/plots/lwgup/VERIFICATION.SRB.rc +++ /dev/null @@ -1,33 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_Longwave_Monthly_(198307-200506) - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - :: diff --git a/GEOS_Util/plots/lwgupc/ERA5form.gs b/GEOS_Util/plots/lwgupc/ERA5form.gs new file mode 100644 index 00000000..c474a658 --- /dev/null +++ b/GEOS_Util/plots/lwgupc/ERA5form.gs @@ -0,0 +1,24 @@ +function formula (args) + +'numargs 'args + numargs = result + +* OBS EXPORTS & GCS +* ----------------- +if( subwrd(args,1) = '-OBS' ) + return "MSDWLWRFCS:ERA5 MSNLWRFCS:ERA5" +endif + +* Formula: lwgup = msdwlwrfcs - (-1 * msnlwrfcs) +* Import of msnlwrfcs is scaled by -1 in VERIFICATION.ERA5.rc for comparison with lwgnetc +* ------------------------------------ + n = 1 +while( n <= numargs ) + oname.n = subwrd(args,n) + n = n + 1 +endwhile + +say 'define qobs = 'oname.1' + 'oname.2 + 'define qobs = 'oname.1' + 'oname.2 + +return diff --git a/GEOS_Util/plots/lwgupc/VERIFICATION.SRB.rc b/GEOS_Util/plots/lwgupc/VERIFICATION.SRB.rc deleted file mode 100644 index 97ca40e0..00000000 --- a/GEOS_Util/plots/lwgupc/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_Longwave_Monthly_(198307-200506) - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - :: diff --git a/GEOS_Util/plots/lwgupcf/VERIFICATION.SRB.rc b/GEOS_Util/plots/lwgupcf/VERIFICATION.SRB.rc deleted file mode 100644 index 97ca40e0..00000000 --- a/GEOS_Util/plots/lwgupcf/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_Longwave_Monthly_(198307-200506) - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - :: diff --git a/GEOS_Util/plots/nettoa/VERIFICATION.CERES_EBAF.rc b/GEOS_Util/plots/nettoa/VERIFICATION.CERES_EBAF.rc deleted file mode 100644 index a3072c1e..00000000 --- a/GEOS_Util/plots/nettoa/VERIFICATION.CERES_EBAF.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF -OBSDSC: CERES_EBAF_L4_Ed2-6r - - filename: 'VERIFICATION/CERES-EBAF/CERES-EBAF_L4_Ed2-6r.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSC' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/nettoa/VERIFICATION.ERBE.rc b/GEOS_Util/plots/nettoa/VERIFICATION.ERBE.rc deleted file mode 100644 index d2c6d694..00000000 --- a/GEOS_Util/plots/nettoa/VERIFICATION.ERBE.rc +++ /dev/null @@ -1,14 +0,0 @@ - -OBSNAM: ERBE -OBSDSC: ERBE_Monthly_Data - - filename: 'VERIFICATION/Clouds_radiation/erbe/erbe2x25.tabl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - 'OSR' , 'SOLAR' , 'OSR' , 1 - 'OSRCLR' , 'SOLAR' , 'OSRCLR' , 1 - 'RSR' , 'SOLAR' , 'NSR' , 1 - 'RSC' , 'SOLAR' , 'NSRCLR' , 1 - :: - diff --git a/GEOS_Util/plots/nettoac/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/nettoac/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc deleted file mode 100644 index 821f87ef..00000000 --- a/GEOS_Util/plots/nettoac/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc +++ /dev/null @@ -1,15 +0,0 @@ -OBSNAM: CERES_EBAF-4 -OBSDSC: CERES_EBAF-TOA_Ed4.0_Subset_200003-201701 - - filename: 'VERIFICATION/CERES_EBAF-TOA_Ed4.0/CERES_EBAF-TOA_Ed4.0_Subset_200003-201701.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OLCC5' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/olr/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/olr/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc deleted file mode 100644 index 821f87ef..00000000 --- a/GEOS_Util/plots/olr/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc +++ /dev/null @@ -1,15 +0,0 @@ -OBSNAM: CERES_EBAF-4 -OBSDSC: CERES_EBAF-TOA_Ed4.0_Subset_200003-201701 - - filename: 'VERIFICATION/CERES_EBAF-TOA_Ed4.0/CERES_EBAF-TOA_Ed4.0_Subset_200003-201701.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OLCC5' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/olr/VERIFICATION.CERES_EBAF.rc b/GEOS_Util/plots/olr/VERIFICATION.CERES_EBAF.rc deleted file mode 100644 index a3072c1e..00000000 --- a/GEOS_Util/plots/olr/VERIFICATION.CERES_EBAF.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF -OBSDSC: CERES_EBAF_L4_Ed2-6r - - filename: 'VERIFICATION/CERES-EBAF/CERES-EBAF_L4_Ed2-6r.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSC' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/olr/VERIFICATION.ERBE.rc b/GEOS_Util/plots/olr/VERIFICATION.ERBE.rc deleted file mode 100644 index d2c6d694..00000000 --- a/GEOS_Util/plots/olr/VERIFICATION.ERBE.rc +++ /dev/null @@ -1,14 +0,0 @@ - -OBSNAM: ERBE -OBSDSC: ERBE_Monthly_Data - - filename: 'VERIFICATION/Clouds_radiation/erbe/erbe2x25.tabl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - 'OSR' , 'SOLAR' , 'OSR' , 1 - 'OSRCLR' , 'SOLAR' , 'OSRCLR' , 1 - 'RSR' , 'SOLAR' , 'NSR' , 1 - 'RSC' , 'SOLAR' , 'NSRCLR' , 1 - :: - diff --git a/GEOS_Util/plots/olrc/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/olrc/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc deleted file mode 100644 index 821f87ef..00000000 --- a/GEOS_Util/plots/olrc/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc +++ /dev/null @@ -1,15 +0,0 @@ -OBSNAM: CERES_EBAF-4 -OBSDSC: CERES_EBAF-TOA_Ed4.0_Subset_200003-201701 - - filename: 'VERIFICATION/CERES_EBAF-TOA_Ed4.0/CERES_EBAF-TOA_Ed4.0_Subset_200003-201701.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OLCC5' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/olrc/VERIFICATION.CERES_EBAF.rc b/GEOS_Util/plots/olrc/VERIFICATION.CERES_EBAF.rc deleted file mode 100644 index 83bbcaca..00000000 --- a/GEOS_Util/plots/olrc/VERIFICATION.CERES_EBAF.rc +++ /dev/null @@ -1,15 +0,0 @@ -OBSNAM: CERES_EBAF -OBSDSC: CERES_EBAF_L4_Ed2-6r - - filename: 'VERIFICATION/CERES-EBAF/CERES-EBAF_L4_Ed2-6r.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OLCC5' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/olrc/VERIFICATION.ERBE.rc b/GEOS_Util/plots/olrc/VERIFICATION.ERBE.rc deleted file mode 100644 index b00f4acb..00000000 --- a/GEOS_Util/plots/olrc/VERIFICATION.ERBE.rc +++ /dev/null @@ -1,11 +0,0 @@ - -OBSNAM: ERBE -OBSDSC: ERBE_Monthly_Data - - filename: 'VERIFICATION/Clouds_radiation/erbe/erbe2x25.tabl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 - :: - diff --git a/GEOS_Util/plots/olrc/VERIFICATION.SRB_LW.rc b/GEOS_Util/plots/olrc/VERIFICATION.SRB_LW.rc deleted file mode 100644 index 3d5c42b0..00000000 --- a/GEOS_Util/plots/olrc/VERIFICATION.SRB_LW.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_LW_Monthly_Data - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 , - :: diff --git a/GEOS_Util/plots/olrcc5/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/olrcc5/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc deleted file mode 100644 index e0f7dcc2..00000000 --- a/GEOS_Util/plots/olrcc5/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc +++ /dev/null @@ -1,15 +0,0 @@ -OBSNAM: CERES_EBAF-4 -OBSDSC: CERES_EBAF-TOA_Ed4.0_Subset_200003-201701 - - filename: 'VERIFICATION/CERES_EBAF-TOA_Ed4.0/CERES_EBAF-TOA_Ed4.0_Subset_200003-201701.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OLCC5' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/olrcf/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/olrcf/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc deleted file mode 100644 index 821f87ef..00000000 --- a/GEOS_Util/plots/olrcf/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc +++ /dev/null @@ -1,15 +0,0 @@ -OBSNAM: CERES_EBAF-4 -OBSDSC: CERES_EBAF-TOA_Ed4.0_Subset_200003-201701 - - filename: 'VERIFICATION/CERES_EBAF-TOA_Ed4.0/CERES_EBAF-TOA_Ed4.0_Subset_200003-201701.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OLCC5' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/olrcf/VERIFICATION.CERES_EBAF.rc b/GEOS_Util/plots/olrcf/VERIFICATION.CERES_EBAF.rc deleted file mode 100644 index 83bbcaca..00000000 --- a/GEOS_Util/plots/olrcf/VERIFICATION.CERES_EBAF.rc +++ /dev/null @@ -1,15 +0,0 @@ -OBSNAM: CERES_EBAF -OBSDSC: CERES_EBAF_L4_Ed2-6r - - filename: 'VERIFICATION/CERES-EBAF/CERES-EBAF_L4_Ed2-6r.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OLCC5' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/olrcf/VERIFICATION.ERBE.rc b/GEOS_Util/plots/olrcf/VERIFICATION.ERBE.rc deleted file mode 100644 index b00f4acb..00000000 --- a/GEOS_Util/plots/olrcf/VERIFICATION.ERBE.rc +++ /dev/null @@ -1,11 +0,0 @@ - -OBSNAM: ERBE -OBSDSC: ERBE_Monthly_Data - - filename: 'VERIFICATION/Clouds_radiation/erbe/erbe2x25.tabl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 - :: - diff --git a/GEOS_Util/plots/olrcf/VERIFICATION.SRB_LW.rc b/GEOS_Util/plots/olrcf/VERIFICATION.SRB_LW.rc deleted file mode 100644 index 3d5c42b0..00000000 --- a/GEOS_Util/plots/olrcf/VERIFICATION.SRB_LW.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_LW_Monthly_Data - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 , - :: diff --git a/GEOS_Util/plots/osr/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/osr/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc deleted file mode 100644 index 2f56bfff..00000000 --- a/GEOS_Util/plots/osr/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF-4 -OBSDSC: CERES_EBAF-TOA_Ed4.0_Subset_200003-201701 - - filename: 'VERIFICATION/CERES_EBAF-TOA_Ed4.0/CERES_EBAF-TOA_Ed4.0_Subset_200003-201701.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/osr/VERIFICATION.CERES_EBAF.rc b/GEOS_Util/plots/osr/VERIFICATION.CERES_EBAF.rc deleted file mode 100644 index a3072c1e..00000000 --- a/GEOS_Util/plots/osr/VERIFICATION.CERES_EBAF.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF -OBSDSC: CERES_EBAF_L4_Ed2-6r - - filename: 'VERIFICATION/CERES-EBAF/CERES-EBAF_L4_Ed2-6r.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSC' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/osr/VERIFICATION.ERBE.rc b/GEOS_Util/plots/osr/VERIFICATION.ERBE.rc deleted file mode 100644 index d2c6d694..00000000 --- a/GEOS_Util/plots/osr/VERIFICATION.ERBE.rc +++ /dev/null @@ -1,14 +0,0 @@ - -OBSNAM: ERBE -OBSDSC: ERBE_Monthly_Data - - filename: 'VERIFICATION/Clouds_radiation/erbe/erbe2x25.tabl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - 'OSR' , 'SOLAR' , 'OSR' , 1 - 'OSRCLR' , 'SOLAR' , 'OSRCLR' , 1 - 'RSR' , 'SOLAR' , 'NSR' , 1 - 'RSC' , 'SOLAR' , 'NSRCLR' , 1 - :: - diff --git a/GEOS_Util/plots/osrc/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/osrc/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc deleted file mode 100644 index 2f56bfff..00000000 --- a/GEOS_Util/plots/osrc/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF-4 -OBSDSC: CERES_EBAF-TOA_Ed4.0_Subset_200003-201701 - - filename: 'VERIFICATION/CERES_EBAF-TOA_Ed4.0/CERES_EBAF-TOA_Ed4.0_Subset_200003-201701.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/osrc/VERIFICATION.CERES_EBAF.rc b/GEOS_Util/plots/osrc/VERIFICATION.CERES_EBAF.rc deleted file mode 100644 index a3072c1e..00000000 --- a/GEOS_Util/plots/osrc/VERIFICATION.CERES_EBAF.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF -OBSDSC: CERES_EBAF_L4_Ed2-6r - - filename: 'VERIFICATION/CERES-EBAF/CERES-EBAF_L4_Ed2-6r.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSC' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/osrc/VERIFICATION.ERBE.rc b/GEOS_Util/plots/osrc/VERIFICATION.ERBE.rc deleted file mode 100644 index d2c6d694..00000000 --- a/GEOS_Util/plots/osrc/VERIFICATION.ERBE.rc +++ /dev/null @@ -1,14 +0,0 @@ - -OBSNAM: ERBE -OBSDSC: ERBE_Monthly_Data - - filename: 'VERIFICATION/Clouds_radiation/erbe/erbe2x25.tabl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - 'OSR' , 'SOLAR' , 'OSR' , 1 - 'OSRCLR' , 'SOLAR' , 'OSRCLR' , 1 - 'RSR' , 'SOLAR' , 'NSR' , 1 - 'RSC' , 'SOLAR' , 'NSRCLR' , 1 - :: - diff --git a/GEOS_Util/plots/osrcf/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/osrcf/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc deleted file mode 100644 index 2f56bfff..00000000 --- a/GEOS_Util/plots/osrcf/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF-4 -OBSDSC: CERES_EBAF-TOA_Ed4.0_Subset_200003-201701 - - filename: 'VERIFICATION/CERES_EBAF-TOA_Ed4.0/CERES_EBAF-TOA_Ed4.0_Subset_200003-201701.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/osrcf/VERIFICATION.CERES_EBAF.rc b/GEOS_Util/plots/osrcf/VERIFICATION.CERES_EBAF.rc deleted file mode 100644 index a3072c1e..00000000 --- a/GEOS_Util/plots/osrcf/VERIFICATION.CERES_EBAF.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF -OBSDSC: CERES_EBAF_L4_Ed2-6r - - filename: 'VERIFICATION/CERES-EBAF/CERES-EBAF_L4_Ed2-6r.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSC' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/osrcf/VERIFICATION.ERBE.rc b/GEOS_Util/plots/osrcf/VERIFICATION.ERBE.rc deleted file mode 100644 index d2c6d694..00000000 --- a/GEOS_Util/plots/osrcf/VERIFICATION.ERBE.rc +++ /dev/null @@ -1,14 +0,0 @@ - -OBSNAM: ERBE -OBSDSC: ERBE_Monthly_Data - - filename: 'VERIFICATION/Clouds_radiation/erbe/erbe2x25.tabl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - 'OSR' , 'SOLAR' , 'OSR' , 1 - 'OSRCLR' , 'SOLAR' , 'OSRCLR' , 1 - 'RSR' , 'SOLAR' , 'NSR' , 1 - 'RSC' , 'SOLAR' , 'NSRCLR' , 1 - :: - diff --git a/GEOS_Util/plots/portrait.script b/GEOS_Util/plots/portrait.script index 58ffa19b..23b43173 100755 --- a/GEOS_Util/plots/portrait.script +++ b/GEOS_Util/plots/portrait.script @@ -283,6 +283,9 @@ endif if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT $EXPORTS -OBS $OBS -NAME precip -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 86400 -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT $EXPORTS -OBS $OBS -NAME precip -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 86400 -PREFIX LAND -LAND -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT CN_PRCP:MOIST AN_PRCP:MOIST -OBS CN_PRCP:MOIST -NAME precipcn -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 86400 -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LS_PRCP:MOIST -OBS LS_PRCP:MOIST -NAME precipls -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 86400 -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run $GEOSUTIL/plots/precip/precip $expid $output $debug $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run $GEOSUTIL/plots/olr/olr $expid $output $debug $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run $GEOSUTIL/plots/olr/olrcc5 $expid $output $debug $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 @@ -302,6 +305,10 @@ endif if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT TAUX:SURFACE -NAME taux -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 100 -OCEAN -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT TAUY:SURFACE -NAME tauy -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 100 -OCEAN -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT ALBEDO:SOLAR -NAME albedo -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT ALBNF:SOLAR -NAME albnf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT ALBNR:SOLAR -NAME albnr -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT ALBVF:SOLAR -NAME albvf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT ALBVR:SOLAR -NAME albvr -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT TBISCCP:SATSIM -NAME tb -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT RADSRF:RADIATION LHFX:SURFACE SH:SURFACE -OBS QNET:RADIATION -NAME NETHEATG -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 @@ -320,11 +327,11 @@ endif # Radiation: Shortwave # -------------------- - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME radswt -EXPID $expid -EXPORT SLRTP:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgnet -EXPID $expid -EXPORT RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME radswt -EXPID $expid -EXPORT SLRTP:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgnetc -EXPID $expid -EXPORT RSCS:SOLAR -OBS SLRSFC:SOLAR ALBEDO:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgnetcf -EXPID $expid -EXPORT RSCS:SOLAR RSRS:SOLAR -OBS SLRSFC:SOLAR ALBEDO:SOLAR RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgnet -EXPID $expid -EXPORT RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgnetc -EXPID $expid -EXPORT RSCS:SOLAR -OBS SLRSFC:SOLAR ALBEDO:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgnetcf -EXPID $expid -EXPORT RSCS:SOLAR RSRS:SOLAR -OBS SLRSFC:SOLAR ALBEDO:SOLAR RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgdwn -EXPID $expid -EXPORT SLRSF:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgdwnc -EXPID $expid -EXPORT SLRSFC:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 @@ -333,10 +340,10 @@ endif if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME osr -EXPID $expid -EXPORT RSR:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME osrc -EXPID $expid -EXPORT RSC:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME osrcf -EXPID $expid -EXPORT RSC:SOLAR RSR:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgup -EXPID $expid -EXPORT SLRSF:SOLAR RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgupc -EXPID $expid -EXPORT SLRSFC:SOLAR RSCS:SOLAR -OBS ALBEDO:SOLAR SLRSFC:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgupcf -EXPID $expid -EXPORT SLRSFC:SOLAR RSCS:SOLAR SLRSF:SOLAR RSRS:SOLAR -OBS ALBEDO:SOLAR SLRSFC:SOLAR SLRSF:SOLAR RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgup -EXPID $expid -EXPORT SLRSF:SOLAR RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgupc -EXPID $expid -EXPORT SLRSFC:SOLAR RSCS:SOLAR -OBS ALBEDO:SOLAR SLRSFC:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -NAME swgupcf -EXPID $expid -EXPORT SLRSFC:SOLAR RSCS:SOLAR SLRSF:SOLAR RSRS:SOLAR -OBS ALBEDO:SOLAR SLRSFC:SOLAR SLRSF:SOLAR RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 # Radiation: Longwave @@ -348,25 +355,28 @@ endif if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT OLC:IRRAD -NAME olrc -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT OLCC5:IRRAD -NAME olrc_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT OLC:IRRAD OLR:IRRAD -NAME olrcf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT OLCC5:IRRAD OLR:IRRAD -NAME olrcf_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD EMIS:SURFACE -OBS SFCDWN:IRRAD -NAME lwgdwn -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LCS:IRRAD EMIS:SURFACE -OBS SFCDWNC:IRRAD -NAME lwgdwnc -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LCSC5:IRRAD EMIS:SURFACE -OBS SFCDWNC:IRRAD -NAME lwgdwnc_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCS:IRRAD EMIS:SURFACE -OBS SFCDWN:IRRAD SFCDWNC:IRRAD -NAME lwgdwncf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCSC5:IRRAD EMIS:SURFACE -OBS SFCDWN:IRRAD SFCDWNC:IRRAD -NAME lwgdwncf_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD LWS:IRRAD -OBS SFCUP:IRRAD SFCDWN:IRRAD -NAME lwgnet -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD LCS:IRRAD -OBS SFCUPC:IRRAD SFCDWNC:IRRAD -NAME lwgnetc -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD LCSC5:IRRAD -OBS SFCUPC:IRRAD SFCDWNC:IRRAD -NAME lwgnetc_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT OLC:IRRAD OLR:IRRAD -NAME olrcf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT OLCC5:IRRAD OLR:IRRAD -NAME olrcf_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD EMIS:SURFACE -OBS SFCDWN:IRRAD -NAME lwgdwn -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LCS:IRRAD EMIS:SURFACE -OBS SFCDWNC:IRRAD -NAME lwgdwnc -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LCSC5:IRRAD EMIS:SURFACE -OBS SFCDWNC:IRRAD -NAME lwgdwnc_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCS:IRRAD EMIS:SURFACE -OBS SFCDWN:IRRAD SFCDWNC:IRRAD -NAME lwgdwncf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCSC5:IRRAD EMIS:SURFACE -OBS SFCDWN:IRRAD SFCDWNC:IRRAD -NAME lwgdwncf_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD LWS:IRRAD -OBS SFCUP:IRRAD SFCDWN:IRRAD -NAME lwgnet -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD LCS:IRRAD -OBS SFCUPC:IRRAD SFCDWNC:IRRAD -NAME lwgnetc -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD LCSC5:IRRAD -OBS SFCUPC:IRRAD SFCDWNC:IRRAD -NAME lwgnetc_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCS:IRRAD -OBS SFCUPC:IRRAD SFCDWNC:IRRAD SFCUP:IRRAD SFCDWN:IRRAD -NAME lwgnetcf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCSC5:IRRAD -OBS SFCUPC:IRRAD SFCDWNC:IRRAD SFCUP:IRRAD SFCDWN:IRRAD -NAME lwgnetcf_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD EMIS:SURFACE LWS:IRRAD -OBS SFCUP:IRRAD -NAME lwgup -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD EMIS:SURFACE LCS:IRRAD -OBS SFCUPC:IRRAD -NAME lwgupc -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD EMIS:SURFACE LCSC5:IRRAD -OBS SFCUPC:IRRAD -NAME lwgupc_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCS:IRRAD EMIS:SURFACE -OBS SFCUP:IRRAD SFCUPC:IRRAD -NAME lwgupcf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 - if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCSC5:IRRAD EMIS:SURFACE -OBS SFCUP:IRRAD SFCUPC:IRRAD -NAME lwgupcf_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD EMIS:SURFACE LWS:IRRAD -OBS SFCUP:IRRAD -NAME lwgup -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD EMIS:SURFACE LCS:IRRAD -OBS SFCUPC:IRRAD -NAME lwgupc -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD EMIS:SURFACE LCSC5:IRRAD -OBS SFCUPC:IRRAD -NAME lwgupc_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCS:IRRAD EMIS:SURFACE -OBS SFCUP:IRRAD SFCUPC:IRRAD -NAME lwgupcf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n ) $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCSC5:IRRAD EMIS:SURFACE -OBS SFCUP:IRRAD SFCUPC:IRRAD -NAME lwgupcf_cc5 -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 # Radiation: Net @@ -577,6 +587,11 @@ if( $MAKE_LIST == TRUE | "$CMPEXP" != "NULL" ) then if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT EVAP:SURFACE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT T2M:SURFACE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT Q2M:SURFACE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT QS:SURFACE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT U10M:SURFACE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT V10M:SURFACE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT U10N:SURFACE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT V10N:SURFACE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT THAT:SURFACE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT QHAT:SURFACE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT WET1:SURFACE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 @@ -586,6 +601,10 @@ if( $MAKE_LIST == TRUE | "$CMPEXP" != "NULL" ) then if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT CLDTT:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT ZPBL:TURBULENCE -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT DMDT:PHYSICS -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT TAULO:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT TAUMD:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT TAUHI:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 + if( $plot == $plotname.$n) $grads $batch -p -c "run gencmp -EXPID $expid -EXPORT TAUTT:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " ; set plots = `echo $plots $plotname.$n` ; @ n = $n + 1 # Dynamics Diagnostics # -------------------- diff --git a/GEOS_Util/plots/precip/VERIFICATION.ECINTERIM.rc b/GEOS_Util/plots/precip/VERIFICATION.ECINTERIM.rc deleted file mode 100644 index c19f6682..00000000 --- a/GEOS_Util/plots/precip/VERIFICATION.ECINTERIM.rc +++ /dev/null @@ -1,8 +0,0 @@ - -OBSNAM: ecinterim -OBSDSC: ECMWF_INTERIM - - filename: 'VERIFICATION/ECMWF_Interim_FullRes/Interim_monthly_FLX.ctl' - filetype: 'BIN' - fields: 'TPREC' , 'MOIST' , 'tp' , 1000 - :: diff --git a/GEOS_Util/plots/precipcn/modform.gs b/GEOS_Util/plots/precipcn/modform.gs new file mode 100644 index 00000000..3dcee52d --- /dev/null +++ b/GEOS_Util/plots/precipcn/modform.gs @@ -0,0 +1,17 @@ +function formula (args) + +'numargs 'args + numargs = result + + n = 1 +while( n <= numargs ) + qname.n = subwrd(args,n) + n = n + 1 +endwhile + +* Formula: precipconv = cnprcp + anprcp +* ---------------------------------------------------------------- +say 'define qmod = 'qname.1' + 'qname.2 + 'define qmod = 'qname.1' + 'qname.2 + +return diff --git a/GEOS_Util/plots/precipcn/plot.rc b/GEOS_Util/plots/precipcn/plot.rc new file mode 100644 index 00000000..f71573af --- /dev/null +++ b/GEOS_Util/plots/precipcn/plot.rc @@ -0,0 +1,15 @@ + CCOLS: 0 50 42 43 44 45 46 47 48 49 39 38 37 36 35 34 33 32 31 21 22 23 24 25 26 27 28 29 + CLEVS: .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 2.3 2.6 3.0 3.3 3.6 4.0 4.5 5.0 5.5 6.0 7.0 8.0 9.0 10. 11. 12. 15. + + DCOLS: 55 49 47 45 44 36 34 33 32 0 21 22 23 24 25 26 27 28 69 + DLEVS: -8 -7 -6 -5 -4 -3 -2 -1 -0.5 0.5 1 2 3 4 5 6 7 8 + +LAND_CLEVS: .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 5.0 6.0 8.0 9.0 10 12.0 15 +LAND_DLEVS: -8 -6 -4 -2 -1 -.8 -.6 -.4 -.2 .2 .4 .6 .8 1 2 4 6 8 + + TITLE: Convective Precipitation (mm/day) +LAND_TITLE: Convective Precipitation Over Land (mm/day) + +GRID: on +YLAB: %.0f +AXLIM: 0 9 diff --git a/GEOS_Util/plots/precipls/plot.rc b/GEOS_Util/plots/precipls/plot.rc new file mode 100644 index 00000000..8fce369f --- /dev/null +++ b/GEOS_Util/plots/precipls/plot.rc @@ -0,0 +1,15 @@ + CCOLS: 0 50 42 43 44 45 46 47 48 49 39 38 37 36 35 34 33 32 31 21 22 23 24 25 26 27 28 29 + CLEVS: .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 2.3 2.6 3.0 3.3 3.6 4.0 4.5 5.0 5.5 6.0 7.0 8.0 9.0 10. 11. 12. 15. + + DCOLS: 55 49 47 45 44 36 34 33 32 0 21 22 23 24 25 26 27 28 69 + DLEVS: -8 -7 -6 -5 -4 -3 -2 -1 -0.5 0.5 1 2 3 4 5 6 7 8 + +LAND_CLEVS: .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 5.0 6.0 8.0 9.0 10 12.0 15 +LAND_DLEVS: -8 -6 -4 -2 -1 -.8 -.6 -.4 -.2 .2 .4 .6 .8 1 2 4 6 8 + + TITLE: Large-Scale Precipitation (mm/day) +LAND_TITLE: Large-Scale Precipitation Over Land (mm/day) + +GRID: on +YLAB: %.0f +AXLIM: 0 9 diff --git a/GEOS_Util/plots/quickplot b/GEOS_Util/plots/quickplot index 97a8dc4e..6e1344fd 100755 --- a/GEOS_Util/plots/quickplot +++ b/GEOS_Util/plots/quickplot @@ -352,8 +352,9 @@ source $plots_dir/$plot.$datetime.$type/.quickplotrc # Set Default OPS Comparison Experiments # -------------------------------------- -set cmpops = "$VERIFICATION/MERRA_MEANS:A \ - $VERIFICATION/MERRA2_MEANS:A" +set cmpops = "$VERIFICATION/MERRA_MEANS:A \ + $VERIFICATION/MERRA2_MEANS:A \ + $VERIFICATION/ERA5_Monthly:V " set cmpops = `echo $cmpops` # Update CMPEXP diff --git a/GEOS_Util/plots/quickstat b/GEOS_Util/plots/quickstat index 98d8687d..f2cf933c 100755 --- a/GEOS_Util/plots/quickstat +++ b/GEOS_Util/plots/quickstat @@ -434,12 +434,14 @@ if( "$plots" =~ *corcmp* ) then set pltdir = $SOURCE/scratch_corcmp$rmstype/$field.`date +%Y%m%d_%H%M%S` mkdir -p $pltdir cd $pltdir ; set tmpdir = `echo $tmpdir $pltdir` + if( $debug == 'TRUE' ) then $GRADS -l -c "run $GEOSUTIL/plots/fcst/pltcorcmp -source $SOURCE -rc $statsrc -fields $field -rms $rmstype" else $GRADS -l -b -c "run $GEOSUTIL/plots/fcst/pltcorcmp -source $SOURCE -rc $statsrc -fields $field -rms $rmstype" & endif sleep 5 + end end @@ -471,8 +473,8 @@ if( "$plots" =~ *corcmp* ) then else $GRADS -l -b -c "run $GEOSUTIL/plots/fcst/pltcorcmpz -source $SOURCE -rc $pltdir/corcmp_$exp.rc -fields $field" & endif - sleep 5 + endif end @@ -502,8 +504,8 @@ if( "$plots" =~ *corcmp* ) then else $GRADS -l -b -c "run $GEOSUTIL/plots/fcst/pltrmscmpz -source $SOURCE -rc $pltdir/corcmp_$exp.rc -fields $field -rms $rmstype" & endif - sleep 5 + end endif end diff --git a/GEOS_Util/plots/radnetg/VERIFICATION.SRB.rc b/GEOS_Util/plots/radnetg/VERIFICATION.SRB.rc deleted file mode 100644 index 40bbf49a..00000000 --- a/GEOS_Util/plots/radnetg/VERIFICATION.SRB.rc +++ /dev/null @@ -1,54 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Sfc_Net_Radiation_(198307-200506) - - filename: 'VERIFICATION/SRB/LW/clr_sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWNC' , 'IRRAD' , 'SFCDWNC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_down_monthly.ctl' - filetype: 'BIN' - fields: 'SFCDWN' , 'IRRAD' , 'SFCDWN' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUP' , 'IRRAD' , 'SFCUP' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_sfc_up_monthly.ctl' - filetype: 'BIN' - fields: 'SFCUPC' , 'IRRAD' , 'SFCUPC' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/clr_toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - :: - - filename: 'VERIFICATION/SRB/LW/toa_up_monthly.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - :: - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - diff --git a/GEOS_Util/plots/radswt/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/radswt/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc deleted file mode 100644 index 2f56bfff..00000000 --- a/GEOS_Util/plots/radswt/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF-4 -OBSDSC: CERES_EBAF-TOA_Ed4.0_Subset_200003-201701 - - filename: 'VERIFICATION/CERES_EBAF-TOA_Ed4.0/CERES_EBAF-TOA_Ed4.0_Subset_200003-201701.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/radswt/VERIFICATION.CERES_EBAF.rc b/GEOS_Util/plots/radswt/VERIFICATION.CERES_EBAF.rc deleted file mode 100644 index a3072c1e..00000000 --- a/GEOS_Util/plots/radswt/VERIFICATION.CERES_EBAF.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF -OBSDSC: CERES_EBAF_L4_Ed2-6r - - filename: 'VERIFICATION/CERES-EBAF/CERES-EBAF_L4_Ed2-6r.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSC' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/radswt/VERIFICATION.SRB.rc b/GEOS_Util/plots/radswt/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/radswt/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/speed/VERIFICATION.COADS.rc b/GEOS_Util/plots/speed/VERIFICATION.COADS.rc deleted file mode 100644 index 14958537..00000000 --- a/GEOS_Util/plots/speed/VERIFICATION.COADS.rc +++ /dev/null @@ -1,8 +0,0 @@ - -OBSNAM: COADS -OBSDSC: COADS_CLIMATOLOGY_(1945-1989) - - filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' - filetype: 'BIN' - fields: 'VENT', 'SURFACE' , 'W3' , 1 - :: diff --git a/GEOS_Util/plots/speed/VERIFICATION.GSSTF.rc b/GEOS_Util/plots/speed/VERIFICATION.GSSTF.rc deleted file mode 100644 index 2506871e..00000000 --- a/GEOS_Util/plots/speed/VERIFICATION.GSSTF.rc +++ /dev/null @@ -1,14 +0,0 @@ - -OBSNAM: GSSTF -OBSDSC: Goddard_V2_Sat-based_Surf_Turb_Fluxes - - filename: 'VERIFICATION/GSSTF_V2_Surf_Flux/testv2b.ctl' - filetype: 'BIN' - fields: 'VENT' , 'SURFACE' , 'S10' , 1 , - 'LHFX' , 'SURFACE' , 'EFLUX' , 1 , - 'SH' , 'SURFACE' , 'HFLUX' , 1 , - 'TAUX' , 'SURFACE' , 'UFLUX' , 100 , - 'TAUY' , 'SURFACE' , 'VFLUX' , 100 , - 'TQV' , 'AGCM' , 'TPWA' , 1 , - 'Q10M' , 'SURFACE' , 'Q10M' , 0.001 , - :: diff --git a/GEOS_Util/plots/swgdwn/VERIFICATION.SRB.rc b/GEOS_Util/plots/swgdwn/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/swgdwn/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/swgdwnc/VERIFICATION.SRB.rc b/GEOS_Util/plots/swgdwnc/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/swgdwnc/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/swgdwncf/VERIFICATION.SRB.rc b/GEOS_Util/plots/swgdwncf/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/swgdwncf/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/swgnet/VERIFICATION.SRB.rc b/GEOS_Util/plots/swgnet/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/swgnet/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/swgnetc/ERA5form.gs b/GEOS_Util/plots/swgnetc/ERA5form.gs new file mode 100644 index 00000000..ddb23282 --- /dev/null +++ b/GEOS_Util/plots/swgnetc/ERA5form.gs @@ -0,0 +1,23 @@ +function formula (args) + +'numargs 'args + numargs = result + +* OBS EXPORTS & GCS +* ----------------- +if( subwrd(args,1) = '-OBS' ) + return "RSCS:SOLAR" +endif + +* Formula: swgnetc = swgnetc +* --------------------------- + n = 1 +while( n <= numargs ) + oname.n = subwrd(args,n) + n = n + 1 +endwhile + +say 'define qobs = 'oname.1 + 'define qobs = 'oname.1 + +return diff --git a/GEOS_Util/plots/swgnetc/VERIFICATION.SRB.rc b/GEOS_Util/plots/swgnetc/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/swgnetc/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/swgnetcf/VERIFICATION.SRB.rc b/GEOS_Util/plots/swgnetcf/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/swgnetcf/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/swgup/VERIFICATION.SRB.rc b/GEOS_Util/plots/swgup/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/swgup/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/swgupc/ERA5form.gs b/GEOS_Util/plots/swgupc/ERA5form.gs new file mode 100644 index 00000000..af3596e2 --- /dev/null +++ b/GEOS_Util/plots/swgupc/ERA5form.gs @@ -0,0 +1,23 @@ +function formula (args) + +'numargs 'args + numargs = result + +* OBS EXPORTS & GCS +* ----------------- +if( subwrd(args,1) = '-OBS' ) + return "SLRSFC:SOLAR RSCS:SOLAR" +endif + +* Formula: swgupc = msdwswrf - msnswrf +* --------------------------------- + n = 1 +while( n <= numargs ) + oname.n = subwrd(args,n) + n = n + 1 +endwhile + +say 'define qobs = 'oname.1' - 'oname.2 + 'define qobs = 'oname.1' - 'oname.2 + +return diff --git a/GEOS_Util/plots/swgupc/VERIFICATION.SRB.rc b/GEOS_Util/plots/swgupc/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/swgupc/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/swgupcf/VERIFICATION.SRB.rc b/GEOS_Util/plots/swgupcf/VERIFICATION.SRB.rc deleted file mode 100644 index 7f643649..00000000 --- a/GEOS_Util/plots/swgupcf/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TPW' , 'MOIST' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/taux/VERIFICATION.COADS.rc b/GEOS_Util/plots/taux/VERIFICATION.COADS.rc deleted file mode 100644 index 7066c545..00000000 --- a/GEOS_Util/plots/taux/VERIFICATION.COADS.rc +++ /dev/null @@ -1,8 +0,0 @@ - -OBSNAM: COADS -OBSDSC: COADS_CLIMATOLOGY_(1945-1989) - - filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' - filetype: 'BIN' - fields: 'TAUX' , 'SURFACE' , 'TAUX3' , 100 - :: diff --git a/GEOS_Util/plots/taux/VERIFICATION.GSSTF.rc b/GEOS_Util/plots/taux/VERIFICATION.GSSTF.rc deleted file mode 100644 index 2506871e..00000000 --- a/GEOS_Util/plots/taux/VERIFICATION.GSSTF.rc +++ /dev/null @@ -1,14 +0,0 @@ - -OBSNAM: GSSTF -OBSDSC: Goddard_V2_Sat-based_Surf_Turb_Fluxes - - filename: 'VERIFICATION/GSSTF_V2_Surf_Flux/testv2b.ctl' - filetype: 'BIN' - fields: 'VENT' , 'SURFACE' , 'S10' , 1 , - 'LHFX' , 'SURFACE' , 'EFLUX' , 1 , - 'SH' , 'SURFACE' , 'HFLUX' , 1 , - 'TAUX' , 'SURFACE' , 'UFLUX' , 100 , - 'TAUY' , 'SURFACE' , 'VFLUX' , 100 , - 'TQV' , 'AGCM' , 'TPWA' , 1 , - 'Q10M' , 'SURFACE' , 'Q10M' , 0.001 , - :: diff --git a/GEOS_Util/plots/tauy/VERIFICATION.COADS.rc b/GEOS_Util/plots/tauy/VERIFICATION.COADS.rc deleted file mode 100644 index 2249fcba..00000000 --- a/GEOS_Util/plots/tauy/VERIFICATION.COADS.rc +++ /dev/null @@ -1,8 +0,0 @@ - -OBSNAM: COADS -OBSDSC: COADS_CLIMATOLOGY_(1945-1989) - - filename: 'VERIFICATION/surface_fields/coads/coadsgrads2.tabl' - filetype: 'BIN' - fields: 'TAUY' , 'SURFACE' , 'TAUY3' , 100 - :: diff --git a/GEOS_Util/plots/tauy/VERIFICATION.GSSTF.rc b/GEOS_Util/plots/tauy/VERIFICATION.GSSTF.rc deleted file mode 100644 index 2506871e..00000000 --- a/GEOS_Util/plots/tauy/VERIFICATION.GSSTF.rc +++ /dev/null @@ -1,14 +0,0 @@ - -OBSNAM: GSSTF -OBSDSC: Goddard_V2_Sat-based_Surf_Turb_Fluxes - - filename: 'VERIFICATION/GSSTF_V2_Surf_Flux/testv2b.ctl' - filetype: 'BIN' - fields: 'VENT' , 'SURFACE' , 'S10' , 1 , - 'LHFX' , 'SURFACE' , 'EFLUX' , 1 , - 'SH' , 'SURFACE' , 'HFLUX' , 1 , - 'TAUX' , 'SURFACE' , 'UFLUX' , 100 , - 'TAUY' , 'SURFACE' , 'VFLUX' , 100 , - 'TQV' , 'AGCM' , 'TPWA' , 1 , - 'Q10M' , 'SURFACE' , 'Q10M' , 0.001 , - :: diff --git a/GEOS_Util/plots/toacf/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/toacf/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc deleted file mode 100644 index 821f87ef..00000000 --- a/GEOS_Util/plots/toacf/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc +++ /dev/null @@ -1,15 +0,0 @@ -OBSNAM: CERES_EBAF-4 -OBSDSC: CERES_EBAF-TOA_Ed4.0_Subset_200003-201701 - - filename: 'VERIFICATION/CERES_EBAF-TOA_Ed4.0/CERES_EBAF-TOA_Ed4.0_Subset_200003-201701.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OLCC5' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/toacf/VERIFICATION.CERES_EBAF.rc b/GEOS_Util/plots/toacf/VERIFICATION.CERES_EBAF.rc deleted file mode 100644 index 83bbcaca..00000000 --- a/GEOS_Util/plots/toacf/VERIFICATION.CERES_EBAF.rc +++ /dev/null @@ -1,15 +0,0 @@ -OBSNAM: CERES_EBAF -OBSDSC: CERES_EBAF_L4_Ed2-6r - - filename: 'VERIFICATION/CERES-EBAF/CERES-EBAF_L4_Ed2-6r.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OLCC5' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/toacf/VERIFICATION.ERBE.rc b/GEOS_Util/plots/toacf/VERIFICATION.ERBE.rc deleted file mode 100644 index af1cad7b..00000000 --- a/GEOS_Util/plots/toacf/VERIFICATION.ERBE.rc +++ /dev/null @@ -1,15 +0,0 @@ - -OBSNAM: ERBE -OBSDSC: ERBE_Monthly_Data - - filename: 'VERIFICATION/Clouds_radiation/erbe/erbe2x25.tabl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'OLR' , 1 - 'OLC' , 'IRRAD' , 'OLRCLR' , 1 - 'OLCC5' , 'IRRAD' , 'OLRCLR' , 1 - 'OSR' , 'SOLAR' , 'OSR' , 1 - 'OSRCLR' , 'SOLAR' , 'OSRCLR' , 1 - 'RSR' , 'SOLAR' , 'NSR' , 1 - 'RSC' , 'SOLAR' , 'NSRCLR' , 1 - :: - diff --git a/GEOS_Util/plots/tpw/VERIFICATION.GSSTF.rc b/GEOS_Util/plots/tpw/VERIFICATION.GSSTF.rc deleted file mode 100644 index 2506871e..00000000 --- a/GEOS_Util/plots/tpw/VERIFICATION.GSSTF.rc +++ /dev/null @@ -1,14 +0,0 @@ - -OBSNAM: GSSTF -OBSDSC: Goddard_V2_Sat-based_Surf_Turb_Fluxes - - filename: 'VERIFICATION/GSSTF_V2_Surf_Flux/testv2b.ctl' - filetype: 'BIN' - fields: 'VENT' , 'SURFACE' , 'S10' , 1 , - 'LHFX' , 'SURFACE' , 'EFLUX' , 1 , - 'SH' , 'SURFACE' , 'HFLUX' , 1 , - 'TAUX' , 'SURFACE' , 'UFLUX' , 100 , - 'TAUY' , 'SURFACE' , 'VFLUX' , 100 , - 'TQV' , 'AGCM' , 'TPWA' , 1 , - 'Q10M' , 'SURFACE' , 'Q10M' , 0.001 , - :: diff --git a/GEOS_Util/plots/tpw/VERIFICATION.SRB.rc b/GEOS_Util/plots/tpw/VERIFICATION.SRB.rc deleted file mode 100644 index 9837e997..00000000 --- a/GEOS_Util/plots/tpw/VERIFICATION.SRB.rc +++ /dev/null @@ -1,38 +0,0 @@ - -OBSNAM: SRB -OBSDSC: SRB_Rel2.5_SW_Monthly_Data - - filename: 'VERIFICATION/SRB/SW/FTOA_monthly.ctl' - filetype: 'BIN' - fields: 'SLRTP' , 'SOLAR' , 'FTOA' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FCLR_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSFC' , 'SOLAR' , 'FCLR' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FABS_monthly.ctl' - filetype: 'BIN' - fields: 'RSRS' , 'SOLAR' , 'FABS' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/FALL_monthly.ctl' - filetype: 'BIN' - fields: 'SLRSF' , 'SOLAR' , 'FALL' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/TCLD_monthly.ctl' - filetype: 'BIN' - fields: 'CLDTT' , 'SOLAR' , 'TCLD' , 0.01 , - :: - - filename: 'VERIFICATION/SRB/SW/SALB_monthly.ctl' - filetype: 'BIN' - fields: 'ALBEDO' , 'SOLAR' , 'SALB' , 1 , - :: - - filename: 'VERIFICATION/SRB/SW/CWV_monthly.ctl' - filetype: 'BIN' - fields: 'TQV' , 'AGCM' , 'CWV' , 1 , - :: diff --git a/GEOS_Util/plots/tseries/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc b/GEOS_Util/plots/tseries/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc deleted file mode 100644 index 2f56bfff..00000000 --- a/GEOS_Util/plots/tseries/VERIFICATION.CERES_EBAF-TOA_Ed4.0.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF-4 -OBSDSC: CERES_EBAF-TOA_Ed4.0_Subset_200003-201701 - - filename: 'VERIFICATION/CERES_EBAF-TOA_Ed4.0/CERES_EBAF-TOA_Ed4.0_Subset_200003-201701.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSRCLR' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/tseries/VERIFICATION.CERES_EBAF.rc b/GEOS_Util/plots/tseries/VERIFICATION.CERES_EBAF.rc deleted file mode 100644 index a3072c1e..00000000 --- a/GEOS_Util/plots/tseries/VERIFICATION.CERES_EBAF.rc +++ /dev/null @@ -1,14 +0,0 @@ -OBSNAM: CERES_EBAF -OBSDSC: CERES_EBAF_L4_Ed2-6r - - filename: 'VERIFICATION/CERES-EBAF/CERES-EBAF_L4_Ed2-6r.ctl' - filetype: 'BIN' - fields: 'OLR' , 'IRRAD' , 'RLUT' , 1 - 'OLC' , 'IRRAD' , 'RLUTCS' , 1 - 'OSR' , 'SOLAR' , 'RSUT' , 1 - 'OSC' , 'SOLAR' , 'RSUTCS' , 1 - 'SLRTP' , 'SOLAR' , 'RSDT' , 1 - 'RSR' , 'SOLAR' , 'RSNUT' , 1 - 'RSC' , 'SOLAR' , 'RSNUTCS' , 1 - :: - diff --git a/GEOS_Util/plots/tseries/VERIFICATION.ECINTERIM.rc b/GEOS_Util/plots/tseries/VERIFICATION.ECINTERIM.rc deleted file mode 100644 index c19f6682..00000000 --- a/GEOS_Util/plots/tseries/VERIFICATION.ECINTERIM.rc +++ /dev/null @@ -1,8 +0,0 @@ - -OBSNAM: ecinterim -OBSDSC: ECMWF_INTERIM - - filename: 'VERIFICATION/ECMWF_Interim_FullRes/Interim_monthly_FLX.ctl' - filetype: 'BIN' - fields: 'TPREC' , 'MOIST' , 'tp' , 1000 - :: diff --git a/GEOS_Util/plots/tseries/VERIFICATION.GPCP.rc b/GEOS_Util/plots/tseries/VERIFICATION.GPCP.rc deleted file mode 100644 index 14a24728..00000000 --- a/GEOS_Util/plots/tseries/VERIFICATION.GPCP.rc +++ /dev/null @@ -1,8 +0,0 @@ - -OBSNAM: GPCPv2.1 -OBSDSC: GPCP_Version-2.1_satellite/gauge - - filename: 'VERIFICATION/GPCPv2.1/gpcp_v2.1_psg.ctl' - filetype: 'BIN' - fields: 'TPREC' , 'MOIST' , 'prec' , 1 - :: diff --git a/GEOS_Util/plots/zcmp/VERIFICATION.ECANA.rc b/GEOS_Util/plots/zcmp/VERIFICATION.ECANA.rc deleted file mode 100644 index ae4df4db..00000000 --- a/GEOS_Util/plots/zcmp/VERIFICATION.ECANA.rc +++ /dev/null @@ -1,13 +0,0 @@ - -OBSNAM: ecops -OBSDSC: ECMWF_OPNL_ANA - - filename: 'VERIFICATION/ECMWF_OPS.NCDF3/xdf.tabl' - filetype: 'HDF' - fields: 'ZLE' , 'DYN' , 'HGHT' , 1 - 'SLP' , 'DYN' , 'SLP' , 1 - 'T' , 'DYN' , 'T' , 1 - 'U' , 'DYN' , 'U' , 1 - 'V' , 'DYN' , 'V' , 1 - 'RH2' , 'MOIST' , 'RH' , 1/100 - :: diff --git a/GEOS_Util/plots/zcmp/VERIFICATION.ECINTERIM.rc b/GEOS_Util/plots/zcmp/VERIFICATION.ECINTERIM.rc deleted file mode 100644 index 78330b32..00000000 --- a/GEOS_Util/plots/zcmp/VERIFICATION.ECINTERIM.rc +++ /dev/null @@ -1,19 +0,0 @@ -OBSNAM: ecinterim -OBSDSC: ECMWF_INTERIM_ANA - - filename: 'VERIFICATION/ECMWF_Interim_FullRes/ecint_prs.ctl' - filetype: 'BIN' - fields: 'ZLE' , 'DYN' , 'Z' , 1/9.80665 - 'T' , 'DYN' , 'T' , 1 - 'U' , 'DYN' , 'U' , 1 - 'V' , 'DYN' , 'V' , 1 - 'RH2' , 'MOIST' , 'R' , 1/100 - 'Q' , 'MOIST' , 'Q' , 1 - 'OMEGA' , 'DYN' , 'W' , 1 - :: - - filename: 'VERIFICATION/ECMWF_Interim_FullRes/surface/MSL_197901_201111.nc' - filetype: 'SDF' - fields: 'SLP' , 'DYN' , 'MSL' , 1 - :: - diff --git a/GEOS_Util/plots/zcmp/VERIFICATION.GDAS1.rc b/GEOS_Util/plots/zcmp/VERIFICATION.GDAS1.rc deleted file mode 100644 index 2b0c7545..00000000 --- a/GEOS_Util/plots/zcmp/VERIFICATION.GDAS1.rc +++ /dev/null @@ -1,13 +0,0 @@ - -OBSNAM: gdas1 -OBSDSC: NCEP_OPNL_gdas1 - - filename: 'VERIFICATION/NCEP_GDAS-1.NCDF3/xdf.tabl' - filetype: 'HDF' - fields: 'ZLE' , 'DYN' , 'HGHT' , 1 - 'SLP' , 'DYN' , 'SLP' , 1 - 'T' , 'DYN' , 'TMPU' , 1 - 'U' , 'DYN' , 'UWND' , 1 - 'V' , 'DYN' , 'VWND' , 1 - 'RH2' , 'MOIST' , 'RH' , 1/100 - :: diff --git a/GEOS_Util/plots/zcmp/VERIFICATION.NCEP.rc b/GEOS_Util/plots/zcmp/VERIFICATION.NCEP.rc deleted file mode 100644 index 43c69413..00000000 --- a/GEOS_Util/plots/zcmp/VERIFICATION.NCEP.rc +++ /dev/null @@ -1,19 +0,0 @@ - -OBSNAM: ncep -OBSDSC: NCEP_REANAL_Ver-2 - - filename: 'VERIFICATION/NCEP-R2/pgb.anl.197901.avrg.grib.ctl' - filetype: 'BIN' - fields: 'ZLE' , 'DYN' , 'HGHT' , 1 - 'SLP' , 'DYN' , 'SLP' , 1 - 'T' , 'DYN' , 'TMPU' , 1 - 'U' , 'DYN' , 'UWND' , 1 - 'V' , 'DYN' , 'VWND' , 1 - 'OMEGA' , 'DYN' , 'OMEGA' , 1 - :: - - filename: 'VERIFICATION/NCEP_REANAL-2/ncep_sphu.ctl' - filetype: 'BIN' - fields: 'RH2' , 'MOIST' , 'RH' , 1/100 - 'Q' , 'MOIST' , 'SPHU' , 1 - :: diff --git a/GEOS_Util/plots/zcmp1/VERIFICATION.ECINTERIM.rc b/GEOS_Util/plots/zcmp1/VERIFICATION.ECINTERIM.rc deleted file mode 100644 index 78330b32..00000000 --- a/GEOS_Util/plots/zcmp1/VERIFICATION.ECINTERIM.rc +++ /dev/null @@ -1,19 +0,0 @@ -OBSNAM: ecinterim -OBSDSC: ECMWF_INTERIM_ANA - - filename: 'VERIFICATION/ECMWF_Interim_FullRes/ecint_prs.ctl' - filetype: 'BIN' - fields: 'ZLE' , 'DYN' , 'Z' , 1/9.80665 - 'T' , 'DYN' , 'T' , 1 - 'U' , 'DYN' , 'U' , 1 - 'V' , 'DYN' , 'V' , 1 - 'RH2' , 'MOIST' , 'R' , 1/100 - 'Q' , 'MOIST' , 'Q' , 1 - 'OMEGA' , 'DYN' , 'W' , 1 - :: - - filename: 'VERIFICATION/ECMWF_Interim_FullRes/surface/MSL_197901_201111.nc' - filetype: 'SDF' - fields: 'SLP' , 'DYN' , 'MSL' , 1 - :: - diff --git a/GEOS_Util/post/regrid.pl b/GEOS_Util/post/regrid.pl index 1b9da229..7bdb23a2 100755 --- a/GEOS_Util/post/regrid.pl +++ b/GEOS_Util/post/regrid.pl @@ -22,20 +22,24 @@ # global variables #----------------- -my ($bcsHEAD, $bcsHEAD_ops, $bcsTagIN, $bcsTagOUT, $bkgFLG, $bkg_regrid_FLG); -my ($capture, $debug, $dbHash, $dyn2dynX, $c2cX, $drymassFLG); -my ($ESMABIN, $ESMATAG, $expid, $gcmFLG, $g5modules, $getinput, $getsponsorX); -my ($grIN, $grINocean, $grINocean_, $grOUT, $grOUTocean, $grpID, $grouplist); -my ($hr, $interactive, $interp_restartsX, $landIceDT, $lblFLG, $lcvFLG); -my ($levsIN, $levsOUT, $logfile, $merra, $mk_catch, $mk_catchcn, $mk_route); -my ($mk_RestartsX, $mkdrstdateX, $month, $newid, $node, $noprompt, $outdir); -my ($regridj, $qos, $rs_hinterpX, $rs_scaleX, $rsFLG, $rstdir, $rst_tarfile); -my ($rst_template, $rst_templateB, $rstIN_template, $rstIN_templateB); -my ($slurmjob, $scale_catchX, $scale_catchcnX, $surfFLG, $surflay, $surflayIN); -my ($tagIN, $tagOUT, $upairFLG, $verbose, $workdir, $year, $ymd, $zoom, $zoom_); -my (%IN, %iceIN, %OUT, %CS, %CSo, %atmLevs, %input_restarts, %hgrd); -my (%im, %im4, %imo, %imo4, %jm, %jm4, %jm5, %jmo, %jmo4); -my (%SURFACE, %UPPERAIR_OPT, %UPPERAIR_REQ, @anafiles, @warnings); +my ($ESMABIN, $ESMATAG, $bcsHEAD, $bcsHEAD_ops, $bcsTagIN, $bcsTagOUT); +my ($bkgFLG, $bkg_regrid_FLG, $c2cX, $capture, $coupled_model_dir); +my ($dbHash, $debug, $drymassFLG, $dyn2dynX, $expid); +my ($g5modules, $gcmFLG, $getinput, $grIN, $grINocean, $grINocean_); +my ($grOUT, $grOUTocean, $grouplist, $grpID, $hr, $interactive); +my ($interp_restartsX, $landIceDT, $lblFLG, $lcvFLG, $levsIN, $levsOUT); +my ($logfile, $merra, $mk_RestartsX, $mk_catch, $mk_catchcn, $mk_route); +my ($mkdrstdateX, $month, $newid, $node, $noprompt, $outdir, $outdir_save); +my ($qos, $regridj, $rsFLG, $rstdir, $rstTAR, $rs_hinterpX, $rs_scaleX); +my ($rstIN_template, $rstIN_templateB, $rst_tarfile); +my ($rst_template, $rst_templateB, $scale_catchX, $scale_catchcnX); +my ($slurmjob, $surfFLG, $surflay, $surflayIN, $tagIN, $tagOUT); +my ($tarFLG, $upairFLG, $verbose, $workdir, $year, $ymd, $zoom, $zoom_); +my (%CS, %CSo, %IN, %OUT, %SURFACE, %UPPERAIR_OPT, %UPPERAIR_REQ); +my (%atmLevs, %coupledFLG, %coupled_model_tile, %hgrd, %iceIN); +my (%im, %im4, %imo, %imo4, %input_restarts); +my (%jm, %jm4, %jm5, %jmo, %jmo4); +my (@anafiles, @warnings); # global tag variables #--------------------- @@ -78,6 +82,35 @@ $jmo{"C$_"} = 6*$_; } +# coupled ocean grids +#-------------------- +%coupledFLG = ( "cc" => 1, "dd" => 1, "ee" => 1 ); +$imo{"cc"} = "360"; $jmo{"cc"} = "200"; +$imo{"dd"} = "720"; $jmo{"dd"} = "410"; +$imo{"ee"} = "1440"; $jmo{"ee"} = "1080"; + +# until an official location for coupled-model tiles is created +#-------------------------------------------------------------- +$coupled_model_dir = "/discover/nobackup/yvikhlia/coupled/Forcings"; + +$coupled_model_tile{"CF0048x6C_TM0360xTM0200-Pfafstetter.til"} = + "${coupled_model_dir}/Ganymed/a48x288_o360x200/"; + +$coupled_model_tile{"CF0090x6C_TM0360xTM0200-Pfafstetter.til"} = + "${coupled_model_dir}/a90x540_o360x200"; + +$coupled_model_tile{"CF0090x6C_TM0720xTM0410-Pfafstetter.til"} = + "${coupled_model_dir}/Ganymed/a90x540_o720x410"; + +$coupled_model_tile{"CF0090x6C_TM1440xTM1080-Pfafstetter.til"} = + "${coupled_model_dir}/Ganymed/a90x540_o1440x1080"; + +$coupled_model_tile{"CF0180x6C_TM0720xTM0410-Pfafstetter.til"} = + "${coupled_model_dir}/a180x1080_o720x410"; + +$coupled_model_tile{"CF0180x6C_TM1440xTM1080-Pfafstetter.til"} = + "${coupled_model_dir}/Ganymed/a180x1080_o1440x1080"; + # atmosphere cubed-sphere grids #------------------------------ foreach (qw/ 12 24 48 90 180 360 500 720 1000 1440 2000 2880 5760 /) { @@ -227,6 +260,7 @@ sub init { "newid=s" => \$newid, "tagin=s" => \$tagIN, "tagout=s" => \$tagOUT, + "tar" => \$tarFLG, "rs=i" => \$rsFLG, "catchcn" => \$mk_catchcn, "route" => \$mk_route, @@ -396,9 +430,9 @@ sub init_tag_arrays_and_hashes { GEOSadas-5_22_0_p1 GEOSadas-5_22_0_p2 GEOSadas-5_23_0 GEOSadas-5_23_0_p1 GEOSadas-5_24_0 GEOSadas-5_24_0_p1 ); - # BCS Tags: Icarus-NLv2 (New Land Parameters) + # BCS Tags: Icarus-NLv3 (New Land Parameters) #--------------------------------------------------------------------------- - @INL = qw( INL Icarus-NL Icarus-NLv2 ); + @INL = qw( INL Icarus-NL Icarus-NLv3 ); foreach (@F14) { $landIceVER{$_} = 1; $bcsTAG{$_} = "Fortuna-1_4" } foreach (@F20) { $landIceVER{$_} = 1; $bcsTAG{$_} = "Fortuna-2_0" } @@ -409,7 +443,7 @@ sub init_tag_arrays_and_hashes { foreach (@G30) { $landIceVER{$_} = 2; $bcsTAG{$_} = "Ganymed-1_0_Reynolds" } foreach (@G40) { $landIceVER{$_} = 2; $bcsTAG{$_} = "Ganymed-4_0_Reynolds" } foreach (@ICA) { $landIceVER{$_} = 2; $bcsTAG{$_} = "Icarus_Reynolds" } - foreach (@INL) { $landIceVER{$_} = 2; $bcsTAG{$_} = "Icarus-NLv2_Reynolds" } + foreach (@INL) { $landIceVER{$_} = 2; $bcsTAG{$_} = "Icarus-NLv3_Reynolds" } foreach (@D214) { $landIceVER{$_} = 1; $bcsTAG{$_} = "Fortuna-1_4" } foreach (@D540) { $landIceVER{$_} = 1; $bcsTAG{$_} = "Fortuna-1_4" } @@ -426,9 +460,9 @@ sub init_tag_arrays_and_hashes { # rank of BCS tags #----------------- - %rank = ( "Icarus-NLv2_Ostia" => 20, - "Icarus-NLv2_MERRA-2" => 19, - "Icarus-NLv2_Reynolds" => 18, + %rank = ( "Icarus-NLv3_Ostia" => 20, + "Icarus-NLv3_MERRA-2" => 19, + "Icarus-NLv3_Reynolds" => 18, "Icarus_Ostia" => 17, "Icarus_MERRA-2" => 16, "Icarus_Reynolds" => 15, @@ -449,8 +483,8 @@ sub init_tag_arrays_and_hashes { # minimum tag for catchcn and route options #------------------------------------------ - $rank_1_catchcn = $rank{"Icarus-NLv2_Reynolds"}; - $rank_1_route = $rank{"Icarus-NLv2_Reynolds"}; + $rank_1_catchcn = $rank{"Icarus-NLv3_Reynolds"}; + $rank_1_route = $rank{"Icarus-NLv3_Reynolds"}; $rank_saltwater_split = $rank{"Icarus_Reynolds"}; } @@ -460,10 +494,15 @@ sub init_tag_arrays_and_hashes { # not supplied (if interactive mode) #======================================================================= sub check_inputs { - my ($arcdir, $fvrst, $ans, $prompt, $len, $msg, $warnFLG); - my ($grINocean_dflt, $grOUTocean_, $grOUTocean_dflt, $levsOUTdflt); - my ($fname, $rstlcvIN, $newid_dflt, $bkg_dflt, $lcv_dflt, $lbl_dflt); - my ($dflt, $landIceVERin, $landIceVERout); + my ($ans, $arcdir, $bkg_dflt, $dflt, $fname, $fvrst); + my ($grINocean_dflt, $grOUTocean_dflt); + my ($label, $landIceVERin, $landIceVERout, $lbl_dflt); + my ($lcv_dflt, $len, $levsOUTdflt, $msg, $newid_dflt); + my ($prompt, $rstlcvIN, $warnFLG); + + # check for input tarfile + #------------------------ + get_inputs_from_tarname() if $rstdir; # use MERRA input restarts? #-------------------------- @@ -491,13 +530,14 @@ sub check_inputs { . "Inputs from archive directories are not available.\n" . "==================================================\n"; } - $rstdir = query("Enter INPUT restart directory:"); + $rstdir = query("Enter INPUT tarfile or restart directory:"); print "\nCannot find restart dir: $rstdir\n" unless -d $rstdir; } if ($rstdir) { $rstdir =~ s/\/*$//; # remove trailing '/'s $rstdir = realpath_($rstdir); } + get_inputs_from_tarname($rstdir); # check $outdir #-------------- @@ -509,6 +549,17 @@ sub check_inputs { $outdir = realpath_($outdir); chdir_($outdir, $verbose); + # check $tarFLG + #-------------- + $dflt = "n"; + if ($interactive) { + until (defined($tarFLG)) { + $ans = query("Do you want to tar the output restarts? ", $dflt); + $tarFLG = 0 if lc($ans) eq "n"; + $tarFLG = 1 if lc($ans) eq "y"; + } + } else { $tarFLG = 0 unless defined($tarFLG) } + # check $ymd and $hr #------------------- until ($ymd and $ymd =~ /^\d{8}$/) { @@ -525,10 +576,18 @@ sub check_inputs { # create $workdir #---------------- - $workdir = "$outdir/${ymd}_$$"; + chomp($label = `date +%d%H%M%S`); + $workdir = "$outdir/$label"; unlink_($workdir, $verbose) if -d $workdir; mkpath_($workdir, $verbose); + $outdir_save = $outdir; + if ($tarFLG) { + $outdir = "$workdir/outdir"; + mkpath_($outdir); + chdir_($outdir, $verbose); + } + # set values for MERRA restarts #------------------------------ if ($merra == 1) { @@ -625,34 +684,34 @@ sub check_inputs { # input ocean grid: $grINocean #----------------------------- - # Two variables are used to represent the input ocean grid + # The value, "CSi", is used to represent the input cubed-sphere ocean grid + # whereas, "CS", is used to represent the output cubed-sphere ocean grid. + # The values are different, because one agrees with the input atmosphere + # grid while the other agrees with the output atmosphere grid. # - # 1. If input ocean grid is cubed-sphere ("CS"), then - # --------------------- - # * $grINocean = "CSi" - # --------------------- - # Indicates that input ocean grid agrees with input atmosphere cubed - # grid. This distinguishes input ocean grid from $grOUTocean = "CS" - # which indicates that output ocean grid agrees with output atmosphere - # cubed grid. - # --------------------- - # * $grINocean_ = "CS" - # --------------------- - # This value is used for display purposes and for writing the *.CMD file + # The variable, $grINocean_ (with underscore), is for display purposes + # only and will show the input cubed-sphere ocean grid as "CS" even though + # internally, the value is actually "CSi". # - # 2. If input ocean grid it not cubed sphere, then $grINocean_ = $grINocean + # For non-cubed-sphere grids, $grINocean_ = $grINocean #--------------------------------------------------------------------------- $grINocean_dflt = "c"; $grOUTocean_dflt = "c"; - unless ($grINocean and $imo{$grINocean} and + unless ($grINocean and ($imo{$grINocean}) and $grOUTocean and $imo{$grOUTocean}) { print "\nOcean Grids\n" . "-----------\n" . "c = 360x180 (Reynolds)\n" . "e = 1440x720 (MERRA-2)\n" . "f = 2880x1440 (OSTIA)\n" - . "CS = same as atmosphere (OSTIA cubed-sphere)\n\n"; + . "CS = same as atmosphere (OSTIA cubed-sphere)\n" + . "\n" + . "Coupled Ocean Grids\n" + . "-------------------\n" + . "cc = 360x200\n" + . "dd = 720x410\n" + . "ee = 1440x1080\n\n"; } until ($grINocean and $imo{$grINocean}) { $grINocean = query("Enter INPUT ocean grid:", $grINocean_dflt); @@ -814,20 +873,15 @@ sub check_inputs { } # check $newid #------------- - if ($grOUTocean eq "CS") { $grOUTocean_ = "-$grOUTocean" } - else { $grOUTocean_ = $grOUTocean } - - if ($expid) { $newid_dflt = "$grOUT${grOUTocean_}_$expid" } - else { $newid_dflt = "$grOUT$grOUTocean_" } + if ($expid) { $newid_dflt = "$grOUT${grOUTocean}_$expid" } + else { $newid_dflt = "$grOUT$grOUTocean" } if ($interactive) { until ($newid) { $newid = $newid_dflt unless $newid; $newid = query("Enter Experiment ID for OUTPUT restarts", $newid); } - } else { - $newid = $newid_dflt unless $newid; - } + } else { $newid = $newid_dflt unless $newid } # check $bkgFLG, $lblFLG, and $lcvFLG #------------------------------------ @@ -949,6 +1003,38 @@ sub check_inputs { } } +#======================================================================= +# name - get_inputs_from_tarname +# purpose - if $rstdir is an input tarfile, then get $rstdir, $expid, +# $ymd, and $hr from its name +#======================================================================= +sub get_inputs_from_tarname { + my ($base, $dir, @match); + my ($expidT, $ymdT, $hrT); + return unless -f $rstdir and $rstdir =~ m/\.tar$/; + + $base = basename($rstdir); + $dir = dirname($rstdir); + + @match = ($base =~ m/(\S+)\.rst\.(\d{8})_(\d{2})z\.tar/); + unless (scalar(@match) == 3) { + die "Error. Cannot extract info from tarfile: $rstdir;" } + ($expidT, $ymdT, $hrT) = @match; + + if ($expid and ($expid ne $expidT)) { + die "Error. expid mismatch, $expid != $expidT;" } + if ($ymd and ($ymd ne $ymdT)) { + die "Error. ymd mismatch, $ymd != $ymdT;" } + if ($hr and ($hr ne $hrT)) { + die "Error. hr mismatch, $hr != $hrT;" } + + $rstTAR = $rstdir; + $rstdir = $dir; + $expid = $expidT; + $ymd = $ymdT; + $hr = $hrT; +} + #======================================================================= # name - get_fvrst # purpose - find fvcore restart file in INPUT restart directory and @@ -968,7 +1054,9 @@ sub get_fvrst { # check for rst tar files #------------------------ unless (@fvrsts) { - $pattern = "*\.rst\.*\.tar"; + if ($expid) { $pattern = "$expid\.rst\.*\.tar" } + else { $pattern = "*\.rst\.*\.tar" } + @fvrsts = <$rstdir/$pattern>; @fvrsts = filter_datetimes(@fvrsts); } @@ -1313,11 +1401,9 @@ sub check_programs { } if ($CS{$grOUT}) { $interp_restartsX = "$ESMABIN/interp_restarts.x"; - $getsponsorX = "$ESMABIN/getsponsor.pl"; $g5modules = "$ESMABIN/g5_modules"; $c2cX = "$ESMABIN/c2c.x"; die "Error. $interp_restartsX not found;" unless -e $interp_restartsX; - die "Error. $getsponsorX not found;" unless -e $getsponsorX; die "Error. $g5modules not found;" unless -e $g5modules; die "Error. $c2cX not found;" unless -x $c2cX; } @@ -1687,9 +1773,11 @@ sub getLandIceInput { # purpose - set values in the %IN and %OUT hashes #======================================================================= sub set_IN_OUT { - my ($HH, $bcsTAG, $agrid, $ogrid, $hgrid, $bcsdir, $tile, $topo, $val); - my ($oceanID1, $oceanID2, $atmosID1, $atmosID2, $atmosID3, $atmosID4); - my ($gridID, $im, $jm, $im4, $jm4, $jm5); + my ($HH, $agrid, $atmosID1, $atmosID2, $atmosID3, $atmosID4); + my ($bcsTAG, $bcsdir, $gridID, $gridID_tile, $hgrid); + my ($im, $im4, $imo, $imo4, $jm, $jm4, $jm5, $jmo, $jmo4); + my ($oceanID1, $oceanID2, $ogrid, $ogrid_); + my ($tile, $topo, $val); # tag values #----------- @@ -1719,6 +1807,9 @@ sub set_IN_OUT { $IN{"ogrid"} = $grINocean; $OUT{"ogrid"} = $grOUTocean; + $IN{"ogrid_"} = $grINocean_; # used for display purposes + $OUT{"ogrid_"} = $grOUTocean; + # bkg_eta grid value #------------------- if ($bkg_regrid_FLG) { @@ -1773,41 +1864,63 @@ sub set_IN_OUT { $bcsTAG = $$HH{"bcsTAG"}; $agrid = $$HH{"agrid"}; $ogrid = $$HH{"ogrid"}; - - $oceanID1 = "$imo{$ogrid}x$jmo{$ogrid}"; - $oceanID2 = "DE$imo4{$ogrid}xPE$jmo4{$ogrid}"; - + $ogrid_ = $$HH{"ogrid_"}; + + $im = $im{$agrid}; + $im4 = $im4{$agrid}; + $jm = $jm{$agrid}; + $jm4 = $jm4{$agrid}; + $jm5 = $jm5{$agrid}; + + $imo = $imo{$ogrid}; + $imo4 = $imo4{$ogrid}; + $jmo = $jmo{$ogrid}; + $jmo4 = $jmo4{$ogrid}; + + # atmosphere IDs + #--------------- + # * atmosID1: used in bcs directory name + # * atmosID2: used in topological and vegdyn file names + # * atmosID3: used for display and in jobscript file names + # * atmosID4: used as extension of GCM restart file names + #---------------------------------------------------------- if ($CS{$agrid}) { - $im = $im{$agrid}; - $jm = $jm{$agrid}; - $im4 = $im4{$agrid}; - $jm4 = $jm4{$agrid}; - $jm5 = $jm5{$agrid}; - $atmosID1 = "c$im"; $atmosID2 = "${im}x${jm}"; $atmosID3 = $atmosID2; $atmosID4 = "CF${im4}x6C"; - - if ($ogrid eq "CS" or $ogrid eq "CSi") { - $oceanID1 = $atmosID3; - $oceanID2 = $atmosID4; - } - - if ($bcsTAG eq "Ganymed-1_0") { - $gridID = "${atmosID4}${jm5}_${oceanID2}" } - else { $gridID = "${atmosID4}_${oceanID2}" } } else { - $im = $im{$agrid}; - $jm = $jm{$agrid}; - $im4 = $im4{$agrid}; - $jm4 = $jm4{$agrid}; - $atmosID1 = "${im}x${jm}"; $atmosID2 = "${im}x${jm}_DC"; $atmosID3 = $atmosID1; $atmosID4 = "${im4}x${jm4}"; + } + + # ocean IDs + #---------- + # * oceanID1: used in tile file name for non-CS ocean grids + # * oceanID2: used in grid ID + #---------------------------------------------------------- + $oceanID1 = "${imo}x${jmo}"; + + if ($coupledFLG{$ogrid}) { $oceanID2 = "TM${imo4}xTM${jmo4}" } + elsif ($ogrid_ eq "CS") { $oceanID2 = $atmosID4 } + else { $oceanID2 = "DE${imo4}xPE${jmo4}" } + + # grid IDs + #--------- + if ($coupledFLG{$ogrid}) { + $gridID = "${atmosID4}"; + $gridID_tile = "${atmosID4}_${oceanID2}"; + } + elsif ($CS{$agrid}) { + $gridID = "${atmosID4}_${oceanID2}"; + if ($bcsTAG eq "Ganymed-1_0") { + $gridID = "${atmosID4}${jm5}_${oceanID2}" + } + } + else { $gridID = "DC${im4}xPC${jm4}_${oceanID2}"; } @@ -1819,14 +1932,15 @@ sub set_IN_OUT { $$HH{"atmos4"} = $atmosID4; $$HH{"ocean" } = $oceanID1; $$HH{"gridID"} = $gridID; + $$HH{"gridID"} = $gridID_tile if $coupledFLG{$ogrid}; # bcs directory #-------------- if ($atmosID2 eq "5760x34560") { $bcsdir = "/discover/nobackup/projects/gmao/osse2/stage/BCS_FILES/C5760"; } - elsif ($rank{$bcsTAG} >= $rank{"Icarus-NLv2_Reynolds"}) { - $bcsdir = "$bcsHEAD/Icarus-NLv2/$bcsTAG/$gridID"; + elsif ($rank{$bcsTAG} >= $rank{"Icarus-NLv3_Reynolds"}) { + $bcsdir = "$bcsHEAD/Icarus-NLv3/$bcsTAG/$gridID"; } elsif ($rank{$bcsTAG} >= $rank{"Icarus_Reynolds"}) { if ($bcsHEAD eq $bcsHEAD_ops) { @@ -1847,6 +1961,10 @@ sub set_IN_OUT { else { $bcsdir = "$bcsHEAD/$bcsTAG/$gridID"; } + if ($coupledFLG{$ogrid}) { + $bcsdir = (<${bcsdir}*>)[0]; + } + die "Error; Cannot find bcs directory: $bcsdir;" unless -d $bcsdir; $$HH{"bcsdir"} = $bcsdir; @@ -1859,13 +1977,19 @@ sub set_IN_OUT { # tile file #---------- if ($rank{$bcsTAG} <= $rank{"Ganymed-1_0"}) { - if ($CS{$agrid}) { $tile = "${gridID}_Pfafstetter.til"} + if ($CS{$agrid}) { $tile = "${gridID}_Pfafstetter.til" } else { $tile = "FV_${atmosID1}_DC_${oceanID1}_DE.til" } } else { $tile = "${gridID}-Pfafstetter.til"; } $tile = "$bcsdir/$tile"; + + if ($coupledFLG{$ogrid}) { + $tile = "${gridID_tile}-Pfafstetter.til"; + $tile = "$coupled_model_tile{$tile}/$tile"; + } + die "Error. Cannot find tile file: $tile" unless -f $tile; $$HH{"tile"} = "$tile"; } @@ -1947,7 +2071,9 @@ sub confirm_inputs { . ". hour: $hr\n" . ". atmos grid: $IN{atmos3} ($IN{agrid})\n" . ". atmos levs: $IN{levs}\n" - . ". ocean grid: $IN{ocean} ($IN{ogrid})\n" + . ". ocean grid: $IN{ocean} ($IN{ogrid_})\n" + . ". bcsdir: $IN{bcsdir}\n" + . ". tile file: $IN{tile}\n" . ". BCS tag: $IN{bcsTAG}\n" . ". surflay: $surflayIN\n" . ". rstdir: " .display($rstdir) ."\n"); @@ -1970,12 +2096,14 @@ sub confirm_inputs { . ". hour: $hr\n" . ". atmos grid: $OUT{atmos3} ($OUT{agrid})\n" . ". atmos levs: $OUT{levs}\n" - . ". ocean grid: $OUT{ocean} ($OUT{ogrid})\n"); + . ". ocean grid: $OUT{ocean} ($OUT{ogrid})\n" + . ". bcsdir: $OUT{bcsdir}\n" + . ". tile file: $OUT{tile}\n"); print_( ". bkg_eta grid: $OUT{bkg_regrid}\n") if $OUT{"bkg_regrid"}; print_( ". BCS tag: $OUT{bcsTAG}\n" . ". surflay: $surflay\n" . ". rescale: $rescale_catch\n" - . ". outdir: " .display($outdir) ."\n" + . ". outdir: " .display($outdir_save) ."\n" . ". workdir: " .display($workdir) ."\n\n"); confirm("y"); @@ -2293,7 +2421,7 @@ sub regrid_upperair_rsts_CS { # ------------------------------------- $qsublog =~ s/.o%j//g; - system_("$regridj 1>$qsublog 2>&1"); + system_("$regridj 1>$qsublog 2>&1") && die "Error with $regridj;"; } else { print_("\nThe CS regridding is MPI based; submitting job to PBS\n"); system_("qsub -W block=true -o $qsublog $regridj"); @@ -2342,7 +2470,7 @@ sub regrid_upperair_rsts_LATLON { . " -topo_old " .$IN{"topo"} . " -topo_new " .$OUT{"topo"} . " -im $im{$grOUT} -jm $jm{$grOUT}"; - system_($cmd); + system_($cmd) && die "Error with $rs_hinterpX;"; } #======================================================================= @@ -2350,7 +2478,7 @@ sub regrid_upperair_rsts_LATLON { # purpose - rename the new upperair restarts #======================================================================= sub rename_upperair_rsts { - my ($label, $tagID, $gridID); + my ($tagID, $gridID); my ($arrAddr, $vFLG, $type, $newrst, $newname); $tagID = $OUT{"bcsTAG"}; @@ -2458,7 +2586,7 @@ sub set_dry_mass { $cmd = "$rs_scaleX $fvrst $moist"; $cmd .= " $AREA_bin" if $CS{$grOUT}; - $status = system_($cmd); + $status = system_($cmd) && die "Error with $rs_scaleX;"; die "Error. Setting dry mass of atmosphere;" if $status; # rename outputs @@ -2482,7 +2610,7 @@ sub regrid_surface_rsts { my ($hash1Addr, $hash2Addr, %H1, %H2); my (@SFC, $flags, $InData_dir, $OutData_dir); my ($rstdir1, $expid1, $template1, $tile1, $type, $src, $alt); - my ($rstdir2, $expid2, $template2, $tile2, $cmd, $status); + my ($rstdir2, $expid2, $template2, $tile2, $cmd, $status, $clsm); my ($catchName, $catchIN, $catch, $catch_regrid, $catch_scaled); my ($catchcnName, $catchcnIN, $catchcn, $catchcn_regrid, $catchcn_scaled); my ($label, $tagID, $gridID); @@ -2615,7 +2743,7 @@ sub regrid_surface_rsts { # run mk_Restarts program #------------------------ chdir_($workdir, $verbose); - system_("\n$mk_RestartsX $flags"); + system_("\n$mk_RestartsX $flags") && die "Error with $mk_RestartsX;"; # split saltwater into openwater and seaicethermo, if necessary #-------------------------------------------------------------- @@ -2673,7 +2801,9 @@ sub regrid_surface_rsts { # link clsm directory to OutData #------------------------------- - symlinkinput("$H2{bcsdir}/clsm", $OutData_dir); + $clsm = dirname($tile2) ."/clsm"; + $clsm = "$H2{bcsdir}/clsm" unless -d $clsm; + symlinkinput($clsm, $OutData_dir); # catch and/or catchcn restart only during this pass #--------------------------------------------------- @@ -2696,7 +2826,7 @@ sub regrid_surface_rsts { # run mk_Restarts program #------------------------ chdir_($workdir, $verbose); - system_("\n$mk_RestartsX $flags"); + system_("\n$mk_RestartsX $flags") && die "Error with mk_RestartsX;"; if ($mk_catchcn) { move_("\n$mk_catchcn_j", "$outdir/$mk_catchcn_j.2"); @@ -2728,6 +2858,9 @@ sub regrid_surface_rsts { # check vegdyn_internal_rst #-------------------------- + # commented by S.Mahanama + # at some point in the future this code should be removed + #-------------------------------------------------------- # check_vegdyn(\%H2, $OutData_dir) if $H1{"vegdyn_check"}; } @@ -2927,6 +3060,7 @@ sub get_anafiles { foreach (qw(bkg ana_satb)) { @ana_tar = (`tar tf $rst_tarfile | grep $_`); foreach $ana (@ana_tar) { + next if $ana =~ m/\.orig$/i; chomp($ana); $dbase = basename($ana); $dbase =~ s/$expid/$newid/ if $expid ne $newid; @@ -2972,6 +3106,7 @@ sub satbias_edit { open NEW, "> $satbias" or die "Error. Opening new satbias: $satbias;"; foreach (@biasNEW) { print NEW $_ } close NEW; + unlink_($sbORIG) unless $debug; print_("( airs281SUBSET_aqua -> airs281_aqua )\n"); } } @@ -2998,7 +3133,8 @@ sub regrid_bkg_eta { #-------------------- $flags = "-g5 -res $hgrd{$grOUT} -nlevs $levsOUT"; $cmd = "$dyn2dynX $flags -o $bkg_eta $bkg_eta_orig"; - system_($cmd); + system_($cmd) && die "Error with $dyn2dynX;"; + unlink_($bkg_eta_orig) unless $debug; } #======================================================================= @@ -3009,7 +3145,22 @@ sub write_rst_lcv { my ($rstlcvOUT, $hhmmss); $rstlcvOUT = rstname($newid, "rst.lcv", "$outdir/$rst_templateB"); $hhmmss = $hr."0000"; - system_("\n$mkdrstdateX $ymd $hhmmss $rstlcvOUT"); + system_("\n$mkdrstdateX $ymd $hhmmss $rstlcvOUT") + && die "Error with $mkdrstdateX"; +} + +#======================================================================= +# name - tar_outputs +# purpose - tar all outputs into a single file +#======================================================================= +sub tar_outputs { + my ($tarfile); + return unless $tarFLG; + chdir_($outdir); + $tarfile = "${newid}.rst.${ymd}_${hr}z.tar"; + system_("tar --remove-files -cf $tarfile \*") && die "Error with tar;"; + move_($tarfile, $outdir_save); + $outdir = $outdir_save; } #======================================================================= @@ -3022,8 +3173,10 @@ sub write_CMD_file { # update captured command inputs #------------------------------- $capture =~ s/\s+\-i\b//; - $capture .= " -ymd $ymd" if $capture !~ m/\s+\-ymd\b/; - $capture .= " -hr $hr" if $capture !~ m/\s+\-hr\b/; + unless ($rstTAR) { + $capture .= " -ymd $ymd" if $capture !~ m/\s+\-ymd\b/; + $capture .= " -hr $hr" if $capture !~ m/\s+\-hr\b/; + } $capture .= " -grout $grOUT" if $capture !~ m/\s+\-grout\b/ and $capture !~ m/\s+\-gridout\b/; $capture .= " -levsout $levsOUT" if $capture !~ m/\s+\-levsout\b/ @@ -3038,8 +3191,13 @@ sub write_CMD_file { or $capture !~ m/\s+\-merra\b/; } else { - $capture .= " -d $rstdir" if $capture !~ m/\s+\-d\b/; - $capture .= " -expid $expid" if $capture !~ m/\s+\-expid\b/ and $expid; + if ($rstTAR) { + $capture .= " -d $rstTAR" if $capture !~ m/\s+\-d\b/; + } + else { + $capture .= " -d $rstdir" if $capture !~ m/\s+\-d\b/; + $capture .= " -expid $expid" if $capture !~ m/\s+\-expid\b/ and $expid; + } $capture .= " -tagin $tagIN" if $capture !~ m/\s+\-tagin\b/; $capture .= " -oceanin $grINocean_" if $capture !~ m/\s+\-oceanin\b/; } @@ -3106,6 +3264,15 @@ sub write_CMD_file { #======================================================================= sub cleanup { + # print warnings before closing logfile + #-------------------------------------- + foreach (@warnings) { print_("WARNING: $_\n") } + closeLOG(); + + # tar outputs if requested + #------------------------- + tar_outputs() if $tarFLG; + # remove work directory #---------------------- unless ($debug) { @@ -3114,11 +3281,6 @@ sub cleanup { unlink_($workdir, $verbose); } } - - # print warnings before closing logfile - #-------------------------------------- - foreach (@warnings) { print_("WARNING: $_\n") } - closeLOG(); exit; } @@ -3404,6 +3566,20 @@ sub pause { } } +#======================================================================= +# name - pauseDB +# purpose - print message and then pause the interactive session until +# user enters input; the $noprompt flag is ignored +#======================================================================= +sub pauseDB { + my ($msg); + $msg = shift(@_); + $msg = "" unless defined($msg); + print "$msg\n" if $msg; + print "Hit to continue ... "; + my $dummy = ; +} + #======================================================================= # name - printlabel # purpose - print a label or msg sandwiched between upper and lower lines @@ -3567,7 +3743,7 @@ sub usage { -merra2 get input restarts from OPS MERRA-2 data archives REQUIRED OPTIONS FOR NON-MERRA INPUTS - -d rstdir location of input restart files + -d rstdir pathname for input tarfile or restart directory -expid expid experiment ID of input restart files INTERACTIVE OPTION @@ -3687,6 +3863,6 @@ sub usage { close INFO; select $FH; - system("echo \"$info\" | more"); + system("echo \"$info\" | less"); exit; } diff --git a/GEOS_Util/post/stats.F90 b/GEOS_Util/post/stats.F90 index 4082bd06..cd1cbde2 100644 --- a/GEOS_Util/post/stats.F90 +++ b/GEOS_Util/post/stats.F90 @@ -38,7 +38,7 @@ program stats parameter ( im = 144 ) ! x-dimension for unified comparison parameter ( jm = 91 ) ! y-dimension for unified comparison - parameter ( nr = 12 ) ! r-dimension (number of geographic regions ) + parameter ( nr = 13 ) ! r-dimension (number of geographic regions ) integer,parameter :: luout = 90 ! output file unit for formatted datafile w/ results @@ -69,15 +69,16 @@ program stats real zlon1(nr), lon1 real zlon2(nr), lon2 - data zlat1 / -90., 20., -20., -80., 0., 0., -90.,-90., 20., 30., 60., -90. / - data zlat2 / 90., 80., 20., -20., 90., 90., 0., 0., 60., 60., 90., -60. / - data zlon1 /-180.,-180.,-180.,-180.,-180., 0.,-180., 0.,-140., -10., -180., -180. / - data zlon2 / 180., 180., 180., 180., 0.,180., 0.,180.,- 60., 30., 180., 180. / + data zlat1 / -90., 20., -20., -80., 0., 0., -90.,-90., 20., 30., 60., -90., -60. / + data zlat2 / 90., 80., 20., -20., 90., 90., 0., 0., 60., 60., 90., -60., 60. / + data zlon1 /-180.,-180.,-180.,-180.,-180., 0.,-180., 0.,-140., -10., -180., -180., -180. / + data zlon2 / 180., 180., 180., 180., 0.,180., 0.,180.,- 60., 30., 180., 180., 180. / character(len=7) region(nr) data region / 'global ', 'n.hem ','tropics','s.hem ', & 'nw.quad', 'ne.quad','sw.quad','se.quad', & - 'america', 'europe ','npolar ','spolar ' / + 'america', 'europe ','npolar ','spolar ', & + 'xpolar '/ character*1 char character*256 name, suffix diff --git a/GEOS_Util/post/stats.rc b/GEOS_Util/post/stats.rc index caaf6bcb..60ee6865 100644 --- a/GEOS_Util/post/stats.rc +++ b/GEOS_Util/post/stats.rc @@ -34,7 +34,7 @@ # --------------------------------------------------------------------------------------------------- TYPE_p: scalar - ALIAS_p: slp slprs mean_sea_level_ Mean_sea_level_pressure + ALIAS_p: slp slprs mean_sea_level_ Mean_sea_level_pressure msl TYPE_u10m: vector ALIAS_u10m: u10m diff --git a/GMAO_etc/pyrob b/GMAO_etc/pyrob index 5fc5638a..e8628305 100755 --- a/GMAO_etc/pyrob +++ b/GMAO_etc/pyrob @@ -174,7 +174,7 @@ class Writer(object): # Descr = Coll.LongName.split()[4:] # descr = ' '.join(Descr) - descr = descr.replace('Forecast,','').replace('Assimilation,','').replace('Analysis,','') + descr = descr.replace('Forecast,','').replace('Assimilation,','').replace('Analysis,','').replace('Model,','') # Section Header # -------------- @@ -249,12 +249,12 @@ class stdoutWriter(Writer): def table(self,choice='header'): if choice == 'header': print - print " -----------|------|------------------------------------------------------------|------------" - print " Name | Dims | Description | Units " - print " -----------|------|------------------------------------------------------------|------------" + print " -------------|------|--------------------------------------------------------------------------------|------------" + print " Name | Dims | Description | Units " + print " -------------|------|--------------------------------------------------------------------------------|------------" return def row(self,vname,dim,descr,units): - print ' %-10s | %4s | %-58s | %-10s '%(vname,dim,descr,units) + print ' %-12s | %4s | %-78s | %-10s '%(vname,dim,descr,units) def close(self): pass @@ -282,12 +282,12 @@ class txtWriter(Writer): def table(self,choice='header'): if choice == 'header': self.doc.write("\n") - self.doc.write( " -----------|------|------------------------------------------------------------|------------\n" ) - self.doc.write( " Name | Dims | Description | Units \n" ) - self.doc.write( " -----------|------|------------------------------------------------------------|------------\n" ) + self.doc.write( " -------------|------|--------------------------------------------------------------------------------|------------\n" ) + self.doc.write( " Name | Dims | Description | Units \n" ) + self.doc.write( " -------------|------|--------------------------------------------------------------------------------|------------\n" ) return def row(self,vname,dim,descr,units): - self.doc.write( ' %-10s | %4s | %-58s | %-10s \n'%(vname,dim,descr,units) ) + self.doc.write( ' %-12s | %4s | %-78s | %-10s \n'%(vname,dim,descr,units) ) def close(self): self.doc.close()