From 60cba4feeabda1ca117cfdf8bcc1e3ab009cb681 Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Thu, 12 Mar 2020 14:12:17 -0400 Subject: [PATCH 01/16] conv_gefs2datm --- utilities/conv_gefs2datm/0readme | 1 + 1 file changed, 1 insertion(+) create mode 100644 utilities/conv_gefs2datm/0readme diff --git a/utilities/conv_gefs2datm/0readme b/utilities/conv_gefs2datm/0readme new file mode 100644 index 0000000..cef4e24 --- /dev/null +++ b/utilities/conv_gefs2datm/0readme @@ -0,0 +1 @@ +conv_gefs2datm From 047b1c5eb3370d2c7da0a45afd3115569a9ed133 Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Thu, 12 Mar 2020 14:13:35 -0400 Subject: [PATCH 02/16] Add files via upload --- .../conv_gefs2datm/0Read_me_missing_GEFS | 233 +++++ utilities/conv_gefs2datm/comp_f77_code.csh | 6 + utilities/conv_gefs2datm/conv_gefs2datm.f | 888 ++++++++++++++++++ .../conv_gefs2datm/conv_gefs2datm.fort.csh | 52 + 4 files changed, 1179 insertions(+) create mode 100644 utilities/conv_gefs2datm/0Read_me_missing_GEFS create mode 100644 utilities/conv_gefs2datm/comp_f77_code.csh create mode 100644 utilities/conv_gefs2datm/conv_gefs2datm.f create mode 100644 utilities/conv_gefs2datm/conv_gefs2datm.fort.csh diff --git a/utilities/conv_gefs2datm/0Read_me_missing_GEFS b/utilities/conv_gefs2datm/0Read_me_missing_GEFS new file mode 100644 index 0000000..60e47ef --- /dev/null +++ b/utilities/conv_gefs2datm/0Read_me_missing_GEFS @@ -0,0 +1,233 @@ +#---- Followings are missing or failed in opening GEFS NEMSio data from HPSS archive +#---- hyun-chul.lee@noaa.gov + +Followings were all recovered from errors, by gary.bates@noaa.gov + + +2000103100 --> recovered (01/22) +2000103118 --> recovered (01/23) +2001121118 --> recovered (01/29) +2004062712 --> recovered (01/29) +2004062718 --> recovered (01/29) +2004100500 --> recovered (01/29) +2004100518 --> recovered (01/29) +2005080512 --> recovered (01/29) +2005082000 --> recovered (01/29) +2005082200 --> recovered (01/29) +2009080612 --> recovered (02/04) +2009080618 --> recovered (02/04) +2009080700 --> recovered (02/04) +2010121512 --> recovered (02/04) +2012120212 --> recovered (01/23) +2013051800 --> recovered (02/04) +2015022712 --> recovered (02/04) +2015083018 --> recovered (02/04) +2016043000 --> recovered (02/04) +2017031512 --> recovered (02/04) +2017031618 --> recovered (02/04) +2017031812 --> recovered (02/04) +2018041700 --> recovered (01/22) +2019070806 --> recovered (02/04) + + +#---- massages in error + + +2000103100 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2000103100.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2000103100.tar 2000103100/bfg_2000103100_fhr00_control + + +2000103118 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2000103118.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2000103118.tar 2000103118/bfg_2000103118_fhr00_control + + +2001121118 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2001121118.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2001121118.tar 2001121118/bfg_2001121118_fhr00_control + + +2004062712 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004062712.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004062712.tar 2004062712/bfg_2004062712_fhr00_control + + +2004062718 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004062718.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004062718.tar 2004062718/bfg_2004062718_fhr00_control + + +2004100500 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004100500.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004100500.tar 2004100500/bfg_2004100500_fhr00_control + + +2004100518 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004100518.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004100518.tar 2004100518/bfg_2004100518_fhr00_control + + +2005080512 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005080512.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005080512.tar 2005080512/bfg_2005080512_fhr00_control + + +2005082000 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005082000.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005082000.tar 2005082000/bfg_2005082000_fhr00_control + + +2005082200 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005082200.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005082200.tar 2005082200/bfg_2005082200_fhr00_control + + +2009080612 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080612.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080612.tar 2009080612/bfg_2009080612_fhr00_control + + +2009080618 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080618.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080618.tar 2009080618/bfg_2009080618_fhr00_control + + +2009080700 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080700.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080700.tar 2009080700/bfg_2009080700_fhr00_control + + +2010121512 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2010121512.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2010121512.tar 2010121512/bfg_2010121512_fhr00_control + + +2012120212 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2012120212.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2012120212.tar 2012120212/bfg_2012120212_fhr00_control + + +2013051800 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2013051800.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2013051800.tar 2013051800/bfg_2013051800_fhr00_control + + +2015022712 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] Archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2015022712.tar size = 23691833856 +ERROR: is inconsistent with size of 23691907584 from Index File /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2015022712.tar.idx +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2015022712.tar 2015022712/bfg_2015022712_fhr00_control + + +2015083018 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2015083018.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2015083018.tar 2015083018/bfg_2015083018_fhr00_control + + +2016043000 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2016043000.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2016043000.tar 2016043000/bfg_2016043000_fhr00_control + + +2017031512 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031512.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031512.tar 2017031512/bfg_2017031512_fhr00_control + + +2017031618 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +Opening archive file +ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031618.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031618.tar 2017031618/bfg_2017031618_fhr00_control + + +2017031812 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031812.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031812.tar 2017031812/bfg_2017031812_fhr00_control + + +2018041700 +-rw-r--r-- 1 Hyun-Chul.Lee climate 262144 Dec 26 11:41 bfg_2018041700_fhr00_control +-rw-r--r-- 1 Hyun-Chul.Lee climate 632302296 Sep 8 20:44 bfg_2018041706_fhr00_control +-rw-r--r-- 1 Hyun-Chul.Lee climate 632302296 Sep 8 20:54 bfg_2018041712_fhr00_control +-rw-r--r-- 1 Hyun-Chul.Lee climate 632302296 Sep 8 21:05 bfg_2018041718_fhr00_control + Error: + NetCDF: Variable not found + + +2019070806 +[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] +ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2019070806.tar +HTAR: HTAR FAILED +###WARNING htar returned non-zero exit status. + 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2019070806.tar 2019070806/bfg_2019070806_fhr00_control + + + + + diff --git a/utilities/conv_gefs2datm/comp_f77_code.csh b/utilities/conv_gefs2datm/comp_f77_code.csh new file mode 100644 index 0000000..1aa8f31 --- /dev/null +++ b/utilities/conv_gefs2datm/comp_f77_code.csh @@ -0,0 +1,6 @@ +#! /bin/csh + +set name = $1 +#--- in hera +ifort -132 ${name}.f -o ${name} -L/apps/netcdf/4.6.1/intel/16.1.150/lib -lnetcdff + diff --git a/utilities/conv_gefs2datm/conv_gefs2datm.f b/utilities/conv_gefs2datm/conv_gefs2datm.f new file mode 100644 index 0000000..a9c3553 --- /dev/null +++ b/utilities/conv_gefs2datm/conv_gefs2datm.f @@ -0,0 +1,888 @@ +C---- convert surface fluxes data from GEFS to DATM +C---- hyun-chul.lee@noaa.gov +C---- 12/06/19 change input from prate_avesfc to prateb_avesfc : H-C +C---- ::: prate_avesfc is daily accumulated precipitation in [mm] +C---- ::: prateb_avesfc is precipitation rate in [Kg/m^3/s] + program conv_gefs2datm + implicit none + include 'netcdf.inc' + +C C This is the name of the data file we will read. + character*(*) FILE_NAME, FOUT_NAME + character STR*200,SATT*200 + parameter (FILE_NAME = "./gefs_input.nc") + parameter (FOUT_NAME = "./gefs_output.nc") + +C C We are reading + integer NX,NY,NT,nerr,NDIM,i,j,k,NSTR + parameter (NX = 1536, NY = 768, NT = 1, NDIM = 3) + real*4 :: latitude(NY), longitude(NX) + integer time(NT),NF_64BIT_OFFSET,vi_indx + real coef,Tcel,Hn,term1,term2,g0,Rsp + real, dimension(NX,NY,NT) :: + & dlwrf_avesfc, + & ulwrf_avesfc, + & dswrf_avesfc, + & vbdsf_avesfc, + & vddsf_avesfc, + & nbdsf_avesfc, + & nddsf_avesfc, + & uflx_avesfc, + & vflx_avesfc, + & shtfl_avesfc, + & lhtfl_avesfc, + & prateb_avesfc, + & cpofpsfc, + & ugrd10m, + & vgrd10m, + & pressfc, + & tmp_hyblev1sfc, + & spfh_hyblev1sfc, + & ugrd_hyblev1sfc, + & vgrd_hyblev1sfc, + & hgt_hyblev1sfc, + & icecsfc_in, + & landsfc, + & tmp2m, + & spfh2m, + & VAR1,VAR2 + +C C This will be the netCDF ID for the file and data variable. + integer :: ncid, + & vi_longitude, + & vi_latitude, + & vi_time_in, + & vi_dlwrf_avesfc, + & vi_ulwrf_avesfc, + & vi_dswrf_avesfc, + & vi_vbdsf_avesfc, + & vi_vddsf_avesfc, + & vi_nbdsf_avesfc, + & vi_nddsf_avesfc, + & vi_uflx_avesfc, + & vi_vflx_avesfc, + & vi_shtfl_avesfc, + & vi_lhtfl_avesfc, + & vi_prateb_avesfc, + & vi_cpofpsfc, + & vi_ugrd10m, + & vi_vgrd10m, + & vi_pressfc, + & vi_tmp_hyblev1sfc, + & vi_spfh_hyblev1sfc, + & vi_ugrd_hyblev1sfc, + & vi_vgrd_hyblev1sfc, + & vi_hgt_hyblev1sfc, + & vi_icecsfc_in, + & vi_landsfc, + & vi_tmp2m, + & vi_spfh2m + + +C ! for output + + real*4 :: lat(NY), lon(NX) +c integer :: time_in(NT) + real*8 :: time_in(NT),time_out(NT) +c real*8 :: timed(NT) + real, dimension(NX,NY,NT) :: + & DLWRF, + & ULWRF, + & DSWRF, + & vbdsf_ave, + & vddsf_ave, + & nbdsf_ave, + & nddsf_ave, + & dusfc, + & dvsfc, + & shtfl_ave, + & lhtfl_ave, + & totprcp_ave, + & u10m, + & v10m, + & hgt_hyblev1, + & psurf, + & tmp_hyblev1, + & spfh_hyblev1, + & ugrd_hyblev1, + & vgrd_hyblev1, + & q2m, + & t2m, + & slmsksfc, + & pres_hyblev1, + & precp, + & fprecp, + & icecsfc + + integer :: mcid,outdim(NDIM), outdimx,outdimy,outdimt, + & vi_t_dim, + & vi_y_dim, + & vi_x_dim, + & vi_lon, + & vi_lat, + & vi_time, + & vi_DLWRF, + & vi_ULWRF, + & vi_DSWRF, + & vi_vbdsf_ave, + & vi_vddsf_ave, + & vi_nbdsf_ave, + & vi_nddsf_ave, + & vi_dusfc, + & vi_dvsfc, + & vi_shtfl_ave, + & vi_lhtfl_ave, + & vi_totprcp_ave, + & vi_u10m, + & vi_v10m, + & vi_hgt_hyblev1, + & vi_psurf, + & vi_tmp_hyblev1, + & vi_spfh_hyblev1, + & vi_ugrd_hyblev1, + & vi_vgrd_hyblev1, + & vi_q2m, + & vi_t2m, + & vi_slmsksfc, + & vi_pres_hyblev1, + & vi_precp, + & vi_fprecp, + & vi_icecsfc + +C---------------------------------- + nerr = nf_open(FILE_NAME, NF_NOWRITE, ncid) + if (nerr .ne. nf_noerr) call handle_err(nerr) + +C ! Get the varid of the data variable, based on its name. + + nerr = nf_inq_varid(ncid , "lat" ,vi_latitude ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "lon" ,vi_longitude ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "time" ,vi_time_in ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + nerr = nf_inq_varid(ncid , "dlwrf_avesfc" ,vi_dlwrf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "ulwrf_avesfc" ,vi_ulwrf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "dswrf_avesfc" ,vi_dswrf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "vbdsf_avesfc" ,vi_vbdsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "vddsf_avesfc" ,vi_vddsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "nbdsf_avesfc" ,vi_nbdsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "nddsf_avesfc" ,vi_nddsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "uflx_avesfc" ,vi_uflx_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "vflx_avesfc" ,vi_vflx_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "shtfl_avesfc" ,vi_shtfl_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "lhtfl_avesfc" ,vi_lhtfl_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "prateb_avesfc" ,vi_prateb_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "cpofpsfc" ,vi_cpofpsfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "ugrd10m" ,vi_ugrd10m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "vgrd10m" ,vi_vgrd10m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "pressfc" ,vi_pressfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "tmp_hyblev1sfc" ,vi_tmp_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "spfh_hyblev1sfc" ,vi_spfh_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "ugrd_hyblev1sfc" ,vi_ugrd_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "vgrd_hyblev1sfc" ,vi_vgrd_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "hgt_hyblev1sfc" ,vi_hgt_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "icecsfc" ,vi_icecsfc_in ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "landsfc" ,vi_landsfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "tmp2m" ,vi_tmp2m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_inq_varid(ncid , "spfh2m" ,vi_spfh2m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + + +C ! Read the data. +C! call check( nf90_get_var(ncid, varid, data_in) ) + + + nerr = nf_get_var_real(ncid , vi_latitude, latitude ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_longitude, longitude ) + if (nerr .ne. nf_noerr) call handle_err(nerr) +c nerr = nf_get_var_int(ncid , vi_time_in, time_in ) + nerr = nf_get_var_double(ncid , vi_time_in, time_in ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_dlwrf_avesfc, dlwrf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_ulwrf_avesfc, ulwrf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_dswrf_avesfc, dswrf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_vbdsf_avesfc, vbdsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_vddsf_avesfc, vddsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_nbdsf_avesfc, nbdsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_nddsf_avesfc, nddsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_uflx_avesfc, uflx_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_vflx_avesfc, vflx_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_shtfl_avesfc, shtfl_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_lhtfl_avesfc, lhtfl_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_prateb_avesfc, prateb_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_cpofpsfc, cpofpsfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_ugrd10m, ugrd10m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_vgrd10m, vgrd10m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_pressfc, pressfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_tmp_hyblev1sfc, tmp_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_spfh_hyblev1sfc, spfh_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_ugrd_hyblev1sfc, ugrd_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_vgrd_hyblev1sfc, vgrd_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_hgt_hyblev1sfc, hgt_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_icecsfc_in, icecsfc_in ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_landsfc, landsfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_tmp2m, tmp2m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_get_var_real(ncid , vi_spfh2m, spfh2m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + +C------------------------------ + nerr = nf_close(ncid) + if (nerr .ne. nf_noerr) call handle_err(nerr) +C------------------------------ + + print *, "TIME= ", time_in +c print *, "Longitude= ", longitude + print *,"*** SUCCESS reading file ", FILE_NAME, "! " + + nerr = nf_create(FOUT_NAME, NF_CLOBBER, mcid) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + time_out = time_in + 12.0 + + print *, "TIME_OUT= ", time_out + +C: set dim + nerr = nf_def_dim(mcid, "lon", NX, vi_x_dim) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_dim(mcid, "lat", NY, vi_y_dim) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_dim(mcid, "time",NT,vi_t_dim) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + outdim(3) = vi_t_dim + outdim(2) = vi_y_dim + outdim(1) = vi_x_dim + outdimt = vi_t_dim + outdimy = vi_y_dim + outdimx = vi_x_dim + +C: define variables + + nerr = nf_def_var(mcid , "lon" ,NF_REAL, 1, vi_x_dim, vi_lon ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "lat" ,NF_REAL, 1, vi_y_dim, vi_lat ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "time" ,NF_DOUBLE, 1, vi_t_dim, vi_time ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "DLWRF" ,NF_REAL, NDIM, outdim, vi_DLWRF ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "ULWRF" ,NF_REAL, NDIM, outdim, vi_ULWRF ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "DSWRF" ,NF_REAL, NDIM, outdim, vi_DSWRF ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "vbdsf_ave" ,NF_REAL, NDIM, outdim, vi_vbdsf_ave ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "vddsf_ave" ,NF_REAL, NDIM, outdim, vi_vddsf_ave ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "nbdsf_ave" ,NF_REAL, NDIM, outdim, vi_nbdsf_ave ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "nddsf_ave" ,NF_REAL, NDIM, outdim, vi_nddsf_ave ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "dusfc" ,NF_REAL, NDIM, outdim, vi_dusfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "dvsfc" ,NF_REAL, NDIM, outdim, vi_dvsfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "shtfl_ave" ,NF_REAL, NDIM, outdim, vi_shtfl_ave ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "lhtfl_ave" ,NF_REAL, NDIM, outdim, vi_lhtfl_ave ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "totprcp_ave" ,NF_REAL, NDIM, outdim, vi_totprcp_ave ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "u10m" ,NF_REAL, NDIM, outdim, vi_u10m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "v10m" ,NF_REAL, NDIM, outdim, vi_v10m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "hgt_hyblev1" ,NF_REAL, NDIM, outdim, vi_hgt_hyblev1 ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "psurf" ,NF_REAL, NDIM, outdim, vi_psurf ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "tmp_hyblev1" ,NF_REAL, NDIM, outdim, vi_tmp_hyblev1 ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "spfh_hyblev1" ,NF_REAL, NDIM, outdim, vi_spfh_hyblev1 ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "ugrd_hyblev1" ,NF_REAL, NDIM, outdim, vi_ugrd_hyblev1 ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "vgrd_hyblev1" ,NF_REAL, NDIM, outdim, vi_vgrd_hyblev1 ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "q2m" ,NF_REAL, NDIM, outdim, vi_q2m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "t2m" ,NF_REAL, NDIM, outdim, vi_t2m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "slmsksfc" ,NF_REAL, NDIM, outdim, vi_slmsksfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "pres_hyblev1" ,NF_REAL, NDIM, outdim, vi_pres_hyblev1 ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "precp" ,NF_REAL, NDIM, outdim, vi_precp ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "fprecp" ,NF_REAL, NDIM, outdim, vi_fprecp ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_def_var(mcid , "icecsfc" ,NF_REAL, NDIM, outdim, vi_icecsfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + +C: put attribute + vi_indx = vi_lon + STR = "degrees_east" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_lon + STR = "Longitude" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_lon + nerr = nf_put_att_text(mcid, vi_indx, "modulo",1," ") + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_lat + STR = "degrees_north" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_lat + STR = "Latitude" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_time +C STR = "seconds since 1970-01-01 00.00.00" + STR = "hours since 1-1-1 00:00:00" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_time + STR = "gregorian" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"calendar",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_time + STR = "T" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"axis",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_DLWRF + STR = "W/m**2" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_DLWRF + STR = "surface downward longwave flux" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_ULWRF + STR = "W/m**2" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_ULWRF + STR = "surface upward longwave flux" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_DSWRF + STR = "W/m**2" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_DSWRF + STR = "averaged surface downward shortwave flux" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_vbdsf_ave + STR = "W/m**2" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_vbdsf_ave + STR = "Visible Beam Downward Solar Flux" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_vddsf_ave + STR = "W/m**2" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_vddsf_ave + STR = "Visible Diffuse Downward Solar Flux" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_nbdsf_ave + STR = "W/m**2" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_nbdsf_ave + STR = "Near IR Beam Downward Solar Flux" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_nddsf_ave + STR = "W/m**2" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_nddsf_ave + STR = "Near IR Diffuse Downward Solar Flux" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_dusfc + STR = "N/m**2" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_dusfc + STR = "surface zonal momentum flux" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_dvsfc + STR = "N/m**2" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_dvsfc + STR = "surface meridional momentum flux" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_shtfl_ave + STR = "w/m**2" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_shtfl_ave + STR = "surface sensible heat flux" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_lhtfl_ave + STR = "w/m**2" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_lhtfl_ave + STR = "surface latent heat flux" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_totprcp_ave + STR = "kg/m**2/s" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_totprcp_ave + STR = "surface precipitationrate" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_u10m + STR = "m/s" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_u10m + STR = "10 meter u wind" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_v10m + STR = "m/s" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_v10m + STR = "10 meter v wind" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_hgt_hyblev1 + STR = "m" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_hgt_hyblev1 + STR = "layer 1 height" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_psurf + STR = "Pa" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_psurf + STR = "surface pressure" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_tmp_hyblev1 + STR = "K" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_tmp_hyblev1 + STR = "layer 1 temperature" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_spfh_hyblev1 + STR = "kg/kg" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_spfh_hyblev1 + STR = "layer 1 specific humidity" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_ugrd_hyblev1 + STR = "m/s" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_ugrd_hyblev1 + STR = "layer 1 zonal wind" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_vgrd_hyblev1 + STR = "m/s" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_vgrd_hyblev1 + STR = "layer 1 meridional wind" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_q2m + STR = "kg/kg" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_q2m + STR = "2m specific humidity" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_t2m + STR = "K" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_t2m + STR = "2m temperature" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_slmsksfc + STR = "numerical" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_slmsksfc + STR = "sea-land-ice mask (0-sea, 1-land, 2-ice)" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_pres_hyblev1 + STR = "Pa" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_pres_hyblev1 + STR = "layer 1 pressure" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_precp + STR = "kg/m**2/s" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_precp + STR = "surface rain precipitation rate" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_fprecp + STR = "kg/m**2/s" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_fprecp + STR = "surface snow precipitation rate" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_icecsfc + STR = "numerical" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"units",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + vi_indx = vi_icecsfc + STR = "sea-ice coverage/fraction, proportion" + NSTR = len_trim(STR) + nerr = nf_put_att_text(mcid, vi_indx,"long_name",NSTR,trim(STR)) + if (nerr .ne. nf_noerr) call handle_err(nerr) + +C---------------- + nerr = nf_enddef(mcid) + if (nerr .ne. nf_noerr) call handle_err(nerr) +C---------------- + + +C: put vaiables + + nerr = nf_put_var_real(mcid , vi_lon, longitude ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_lat, latitude ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_double(mcid , vi_time, time_out ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_DLWRF, dlwrf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_ULWRF, ulwrf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_DSWRF, dswrf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_vbdsf_ave, vbdsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_vddsf_ave, vddsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_nbdsf_ave, nbdsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_nddsf_ave, nddsf_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_dusfc, uflx_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_dvsfc, vflx_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_shtfl_ave, shtfl_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_lhtfl_ave, lhtfl_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_totprcp_ave, prateb_avesfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_u10m, ugrd10m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_v10m, vgrd10m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_psurf, pressfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_tmp_hyblev1, tmp_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_spfh_hyblev1, spfh_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_ugrd_hyblev1, ugrd_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_vgrd_hyblev1, vgrd_hyblev1sfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_q2m, spfh2m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_t2m, tmp2m ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_slmsksfc, landsfc ) + if (nerr .ne. nf_noerr) call handle_err(nerr) +C-------------------------------------------------- + coef = 1.0 + g0 = 9.80665 +C--------- J/Kg/K + Rsp = 287.058 + do k=1,NT + do j=1,NY + do i=1,NX +C TCel = tmp2m(i,j,k)-273.15 +C if (TCel >= 0.0) then +C coef = 1.0 +C else if (TCel < -15.0) then +C coef = 0.0 +C else +C coef = (TCel + 15.0)/15.0 +C endif +C precp(i,j,k) = coef * prateb_avesfc(i,j,k) +C fprecp(i,j,k) = (1.0 - coef) * prateb_avesfc(i,j,k) +C +C cpofpsfc(i,j,k) is the ratio of ice precipitation/total precipitation + coef = cpofpsfc(i,j,k) + if (coef < 0.0) then + coef = 0.0 + else if (coef > 1.0) then + coef = 1.0 + endif + fprecp(i,j,k) = coef * prateb_avesfc(i,j,k) + precp(i,j,k) = (1.0 - coef) * prateb_avesfc(i,j,k) +C------------ + Hn = hgt_hyblev1sfc(i,j,k) + hgt_hyblev1(i,j,k) = Hn +C------------ + term1 = -1.0*g0*Hn/(Rsp*tmp_hyblev1sfc(i,j,k)) + term2 = exp(term1) + pres_hyblev1(i,j,k)=pressfc(i,j,k)*term2 + enddo + enddo + enddo +C-------------------------------------------------- + nerr = nf_put_var_real(mcid , vi_hgt_hyblev1, hgt_hyblev1 ) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_pres_hyblev1, pres_hyblev1) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_precp, precp) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_fprecp, fprecp) + if (nerr .ne. nf_noerr) call handle_err(nerr) + nerr = nf_put_var_real(mcid , vi_icecsfc, icecsfc_in ) + if (nerr .ne. nf_noerr) call handle_err(nerr) +C: close + nerr = nf_close(mcid) + if (nerr .ne. nf_noerr) call handle_err(nerr) + + end + +C======================================================= + subroutine handle_err(errcode) + implicit none + include 'netcdf.inc' + integer errcode + + print *, 'Error: ', nf_strerror(errcode) + stop 2 + end + + +C======================================================= + subroutine convolt(NX,NY,NT,invar,outvar,conmuti) + implicit none + integer I,J,K,NX,NY,NT + real conmuti + real :: invar(NT,NY,NX), outvar(NX,NY,NT) +C: set output variables + + do k=1,NT + do j=1,NY + do i=1,NX + outvar(i,j,k) = invar(k,j,i) * conmuti + enddo + enddo + enddo + return + end + + diff --git a/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh b/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh new file mode 100644 index 0000000..3e14a6b --- /dev/null +++ b/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh @@ -0,0 +1,52 @@ +#! /bin/csh -x +#-- hyun-chul.lee@noaa.gov + +module load NetCDF/4.5.0 + +set hdir = /scratch2/NCEPDEV/climate/Hyun-Chul.Lee/GEFS +set sdir = ${hdir}/Reanal +set tdir = ${hdir}/GEFS2DATM +set wdir = ${hdir}/Work +set grb2 = /home/Hyun-Chul.Lee/Tools/wgrib2 + +#set symd = 20180101 +#set eymd = 20180131 + +#set symd = 20130701 +#set eymd = 20130731 + +#set symd = 20171001 +#set eymd = 20171031 + +#set symd = 20000102 +#set eymd = 20031231 + +#set symd = 20040102 +#set eymd = 20111231 + +#set symd = 20120101 +set symd = 20130529 +set eymd = 20190831 + + +cd $wdir + +/bin/cp ${hdir}/conv_gefs2datm conv_gefs2datm + +set ymd = $symd +while ($ymd <= $eymd) + foreach hh (00 06 12 18) + set ymdh = ${ymd}${hh} + echo ${ymdh} + /bin/cp ${sdir}/bfg_${ymdh}_fhr00_control.nc gefs_input.nc + + ./conv_gefs2datm +ncdump gefs_output.nc | sed -e "5s#^.time = 1 ;#time = UNLIMITED ; // (1 currently)#" | ncgen -o gefs_output2.nc + mv gefs_output2.nc gefs_output.nc + /bin/mv gefs_output.nc ${tdir}/gefs.${ymdh}.nc + + rm gefs_input.* + end + set ymd = `date -d "${ymd} 1 day" +%Y%m%d` +end + From b26820b4202cb46648975158d6ed9a971d2388b7 Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Thu, 12 Mar 2020 15:07:38 -0400 Subject: [PATCH 03/16] Update 0readme --- utilities/conv_gefs2datm/0readme | 11 +++++++++++ 1 file changed, 11 insertions(+) diff --git a/utilities/conv_gefs2datm/0readme b/utilities/conv_gefs2datm/0readme index cef4e24..b9a90f4 100644 --- a/utilities/conv_gefs2datm/0readme +++ b/utilities/conv_gefs2datm/0readme @@ -1 +1,12 @@ conv_gefs2datm + + +Input data (GEFS) sources + +in HPSS : /3year/NCEPDEV/GEFSRR/gefsrr_replay_${YYYYstream}stream/${YYYYMMDDHH}.tar +in WCOSS : (venus)/gpfs/dell2/emc/modeling/noscrub/Hyun-Chul.Lee/GEFS/Reanal + +Output data (DATM) archive +in HPSS : /5year/NCEPDEV/marineda/DATM_INPUT/GEFS/${YYYYMM}/gefs.${YYYYMMDDHH}.nc +in Hera : /scratch2/NCEPDEV/marineda/DATM_INPUT/GEFS/${YYYYMM}/gefs.${YYYYMMDDHH}.nc + From df464b99b573b6eda23312480ba46b0896c5040f Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Thu, 12 Mar 2020 15:31:50 -0400 Subject: [PATCH 04/16] Update 0readme --- utilities/conv_gefs2datm/0readme | 18 ++++++++++++++++++ 1 file changed, 18 insertions(+) diff --git a/utilities/conv_gefs2datm/0readme b/utilities/conv_gefs2datm/0readme index b9a90f4..2f11dfe 100644 --- a/utilities/conv_gefs2datm/0readme +++ b/utilities/conv_gefs2datm/0readme @@ -10,3 +10,21 @@ Output data (DATM) archive in HPSS : /5year/NCEPDEV/marineda/DATM_INPUT/GEFS/${YYYYMM}/gefs.${YYYYMMDDHH}.nc in Hera : /scratch2/NCEPDEV/marineda/DATM_INPUT/GEFS/${YYYYMM}/gefs.${YYYYMMDDHH}.nc + +To conver from GEFS to DATM (netCDF), + +1) compile code +csh ./comp_f77_code.csh conv_gefs2datm + +2) modify the run script of conv_gefs2datm.fort.csh + +set sdir = !<-- source of GEFS +set tdir = !<-- output dir +set wdir = !<-- work dir +etc + +set symd = 20130529 !<-- start date +set eymd = 20190831 !<-- end date + +3) run the script +csh ./conv_gefs2datm.fort.csh From c0bf1801c19a7555463f4a86771101b005cccf59 Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Thu, 12 Mar 2020 15:45:51 -0400 Subject: [PATCH 05/16] Update 0readme --- utilities/conv_gefs2datm/0readme | 11 ++++++++--- 1 file changed, 8 insertions(+), 3 deletions(-) diff --git a/utilities/conv_gefs2datm/0readme b/utilities/conv_gefs2datm/0readme index 2f11dfe..96af661 100644 --- a/utilities/conv_gefs2datm/0readme +++ b/utilities/conv_gefs2datm/0readme @@ -1,17 +1,17 @@ conv_gefs2datm -Input data (GEFS) sources +- Input data (GEFS) sources in HPSS : /3year/NCEPDEV/GEFSRR/gefsrr_replay_${YYYYstream}stream/${YYYYMMDDHH}.tar in WCOSS : (venus)/gpfs/dell2/emc/modeling/noscrub/Hyun-Chul.Lee/GEFS/Reanal -Output data (DATM) archive +- Output data (DATM) archive in HPSS : /5year/NCEPDEV/marineda/DATM_INPUT/GEFS/${YYYYMM}/gefs.${YYYYMMDDHH}.nc in Hera : /scratch2/NCEPDEV/marineda/DATM_INPUT/GEFS/${YYYYMM}/gefs.${YYYYMMDDHH}.nc -To conver from GEFS to DATM (netCDF), +- To conver from GEFS to DATM (netCDF), 1) compile code csh ./comp_f77_code.csh conv_gefs2datm @@ -28,3 +28,8 @@ set eymd = 20190831 !<-- end date 3) run the script csh ./conv_gefs2datm.fort.csh + +- More comprensive information can be found + +https://docs.google.com/spreadsheets/d/1M32BvPgXuPqewQWO5cwsECLdcTj9o1x2Rq7GDvSpltk/edit#gid=0 + From 26061f33a47173d55775c9748ed1eadb32d35fb8 Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Fri, 13 Mar 2020 10:07:12 -0400 Subject: [PATCH 06/16] Update conv_gefs2datm.fort.csh --- .../conv_gefs2datm/conv_gefs2datm.fort.csh | 29 +++++-------------- 1 file changed, 7 insertions(+), 22 deletions(-) diff --git a/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh b/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh index 3e14a6b..7b1e55f 100644 --- a/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh +++ b/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh @@ -1,31 +1,16 @@ #! /bin/csh -x -#-- hyun-chul.lee@noaa.gov -module load NetCDF/4.5.0 +module load netcdf/4.6.1 -set hdir = /scratch2/NCEPDEV/climate/Hyun-Chul.Lee/GEFS -set sdir = ${hdir}/Reanal -set tdir = ${hdir}/GEFS2DATM -set wdir = ${hdir}/Work +set hdir = #<-- set home dir +set sdir = #<-- set input source dir +set tdir = #<-- set dir for output +set wdir = #<-- set work dir set grb2 = /home/Hyun-Chul.Lee/Tools/wgrib2 -#set symd = 20180101 -#set eymd = 20180131 - -#set symd = 20130701 -#set eymd = 20130731 - -#set symd = 20171001 -#set eymd = 20171031 - -#set symd = 20000102 -#set eymd = 20031231 - -#set symd = 20040102 -#set eymd = 20111231 - -#set symd = 20120101 +#<-- set intial date set symd = 20130529 +#<-- set end date set eymd = 20190831 From 26ff746cb298b038a15e18ab1263908c72b4b39e Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Fri, 13 Mar 2020 11:00:05 -0400 Subject: [PATCH 07/16] Update conv_gefs2datm.fort.csh --- utilities/conv_gefs2datm/conv_gefs2datm.fort.csh | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh b/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh index 7b1e55f..58a1e80 100644 --- a/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh +++ b/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh @@ -1,19 +1,18 @@ #! /bin/csh -x module load netcdf/4.6.1 +module load wgrib2/2.0.8 set hdir = #<-- set home dir set sdir = #<-- set input source dir set tdir = #<-- set dir for output set wdir = #<-- set work dir -set grb2 = /home/Hyun-Chul.Lee/Tools/wgrib2 #<-- set intial date set symd = 20130529 #<-- set end date set eymd = 20190831 - cd $wdir /bin/cp ${hdir}/conv_gefs2datm conv_gefs2datm From a2a8aa36c2f3a6bad1b8dd7b93cf766a3d6f60d3 Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Fri, 13 Mar 2020 11:33:00 -0400 Subject: [PATCH 08/16] Update conv_gefs2datm.fort.csh --- utilities/conv_gefs2datm/conv_gefs2datm.fort.csh | 1 - 1 file changed, 1 deletion(-) diff --git a/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh b/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh index 58a1e80..32c6b6a 100644 --- a/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh +++ b/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh @@ -1,7 +1,6 @@ #! /bin/csh -x module load netcdf/4.6.1 -module load wgrib2/2.0.8 set hdir = #<-- set home dir set sdir = #<-- set input source dir From 3051ebf70ff49c91b8cc76594f9f96166a659d17 Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Fri, 13 Mar 2020 11:50:48 -0400 Subject: [PATCH 09/16] Delete 0Read_me_missing_GEFS --- .../conv_gefs2datm/0Read_me_missing_GEFS | 233 ------------------ 1 file changed, 233 deletions(-) delete mode 100644 utilities/conv_gefs2datm/0Read_me_missing_GEFS diff --git a/utilities/conv_gefs2datm/0Read_me_missing_GEFS b/utilities/conv_gefs2datm/0Read_me_missing_GEFS deleted file mode 100644 index 60e47ef..0000000 --- a/utilities/conv_gefs2datm/0Read_me_missing_GEFS +++ /dev/null @@ -1,233 +0,0 @@ -#---- Followings are missing or failed in opening GEFS NEMSio data from HPSS archive -#---- hyun-chul.lee@noaa.gov - -Followings were all recovered from errors, by gary.bates@noaa.gov - - -2000103100 --> recovered (01/22) -2000103118 --> recovered (01/23) -2001121118 --> recovered (01/29) -2004062712 --> recovered (01/29) -2004062718 --> recovered (01/29) -2004100500 --> recovered (01/29) -2004100518 --> recovered (01/29) -2005080512 --> recovered (01/29) -2005082000 --> recovered (01/29) -2005082200 --> recovered (01/29) -2009080612 --> recovered (02/04) -2009080618 --> recovered (02/04) -2009080700 --> recovered (02/04) -2010121512 --> recovered (02/04) -2012120212 --> recovered (01/23) -2013051800 --> recovered (02/04) -2015022712 --> recovered (02/04) -2015083018 --> recovered (02/04) -2016043000 --> recovered (02/04) -2017031512 --> recovered (02/04) -2017031618 --> recovered (02/04) -2017031812 --> recovered (02/04) -2018041700 --> recovered (01/22) -2019070806 --> recovered (02/04) - - -#---- massages in error - - -2000103100 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2000103100.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2000103100.tar 2000103100/bfg_2000103100_fhr00_control - - -2000103118 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2000103118.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2000103118.tar 2000103118/bfg_2000103118_fhr00_control - - -2001121118 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2001121118.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_1999stream/2001121118.tar 2001121118/bfg_2001121118_fhr00_control - - -2004062712 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004062712.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004062712.tar 2004062712/bfg_2004062712_fhr00_control - - -2004062718 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004062718.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004062718.tar 2004062718/bfg_2004062718_fhr00_control - - -2004100500 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004100500.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004100500.tar 2004100500/bfg_2004100500_fhr00_control - - -2004100518 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004100518.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2004100518.tar 2004100518/bfg_2004100518_fhr00_control - - -2005080512 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005080512.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005080512.tar 2005080512/bfg_2005080512_fhr00_control - - -2005082000 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005082000.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005082000.tar 2005082000/bfg_2005082000_fhr00_control - - -2005082200 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005082200.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2003stream/2005082200.tar 2005082200/bfg_2005082200_fhr00_control - - -2009080612 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080612.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080612.tar 2009080612/bfg_2009080612_fhr00_control - - -2009080618 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080618.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080618.tar 2009080618/bfg_2009080618_fhr00_control - - -2009080700 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080700.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2009080700.tar 2009080700/bfg_2009080700_fhr00_control - - -2010121512 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2010121512.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2007stream/2010121512.tar 2010121512/bfg_2010121512_fhr00_control - - -2012120212 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2012120212.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2012120212.tar 2012120212/bfg_2012120212_fhr00_control - - -2013051800 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2013051800.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2013051800.tar 2013051800/bfg_2013051800_fhr00_control - - -2015022712 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] Archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2015022712.tar size = 23691833856 -ERROR: is inconsistent with size of 23691907584 from Index File /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2015022712.tar.idx -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2015022712.tar 2015022712/bfg_2015022712_fhr00_control - - -2015083018 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2015083018.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2011stream/2015083018.tar 2015083018/bfg_2015083018_fhr00_control - - -2016043000 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2016043000.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2016043000.tar 2016043000/bfg_2016043000_fhr00_control - - -2017031512 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031512.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031512.tar 2017031512/bfg_2017031512_fhr00_control - - -2017031618 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -Opening archive file -ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031618.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031618.tar 2017031618/bfg_2017031618_fhr00_control - - -2017031812 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] error -5 opening HPSS archive file /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031812.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2017031812.tar 2017031812/bfg_2017031812_fhr00_control - - -2018041700 --rw-r--r-- 1 Hyun-Chul.Lee climate 262144 Dec 26 11:41 bfg_2018041700_fhr00_control --rw-r--r-- 1 Hyun-Chul.Lee climate 632302296 Sep 8 20:44 bfg_2018041706_fhr00_control --rw-r--r-- 1 Hyun-Chul.Lee climate 632302296 Sep 8 20:54 bfg_2018041712_fhr00_control --rw-r--r-- 1 Hyun-Chul.Lee climate 632302296 Sep 8 21:05 bfg_2018041718_fhr00_control - Error: - NetCDF: Variable not found - - -2019070806 -[connecting to hpsscore1.fairmont.rdhpcs.noaa.gov/1217] -ERROR: [FATAL] no such HPSS archive file: /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2019070806.tar -HTAR: HTAR FAILED -###WARNING htar returned non-zero exit status. - 72 = /apps/hpss/bin/htar -T 4 -xvf /3year/NCEPDEV/GEFSRR/gefsrr_replay_2015stream/2019070806.tar 2019070806/bfg_2019070806_fhr00_control - - - - - From 2afe15d0448b0d58c300ca05124b95fee3d3108a Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Mon, 16 Mar 2020 10:15:10 -0400 Subject: [PATCH 10/16] Update comp_f77_code.csh --- utilities/conv_gefs2datm/comp_f77_code.csh | 2 ++ 1 file changed, 2 insertions(+) diff --git a/utilities/conv_gefs2datm/comp_f77_code.csh b/utilities/conv_gefs2datm/comp_f77_code.csh index 1aa8f31..6c3bc7e 100644 --- a/utilities/conv_gefs2datm/comp_f77_code.csh +++ b/utilities/conv_gefs2datm/comp_f77_code.csh @@ -1,4 +1,6 @@ #! /bin/csh +module load intel/19.0.4.243 +module load impi/2019.0.4 set name = $1 #--- in hera From b361fefd9173124b9fa53d9f88daf14ee49c1712 Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Mon, 16 Mar 2020 11:38:54 -0400 Subject: [PATCH 11/16] Update comp_f77_code.csh --- utilities/conv_gefs2datm/comp_f77_code.csh | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/utilities/conv_gefs2datm/comp_f77_code.csh b/utilities/conv_gefs2datm/comp_f77_code.csh index 6c3bc7e..312bb39 100644 --- a/utilities/conv_gefs2datm/comp_f77_code.csh +++ b/utilities/conv_gefs2datm/comp_f77_code.csh @@ -4,5 +4,5 @@ module load impi/2019.0.4 set name = $1 #--- in hera -ifort -132 ${name}.f -o ${name} -L/apps/netcdf/4.6.1/intel/16.1.150/lib -lnetcdff +ifort $FFLAGS ${name}.f -o ${name} -I$NETCDF/include -L$NETCDF/lib -lnetcdff From dd2ab4ab0ffcb21e52d13c77856a18866a55ffc5 Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Mon, 16 Mar 2020 11:44:27 -0400 Subject: [PATCH 12/16] Update comp_f77_code.csh --- utilities/conv_gefs2datm/comp_f77_code.csh | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/utilities/conv_gefs2datm/comp_f77_code.csh b/utilities/conv_gefs2datm/comp_f77_code.csh index 312bb39..786b85f 100644 --- a/utilities/conv_gefs2datm/comp_f77_code.csh +++ b/utilities/conv_gefs2datm/comp_f77_code.csh @@ -1,5 +1,5 @@ #! /bin/csh -module load intel/19.0.4.243 +module load intel/18.0.5.274 module load impi/2019.0.4 set name = $1 From 2cb7fc193c07f3e8b41b138c7013c63df96fdc14 Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Mon, 16 Mar 2020 12:10:02 -0400 Subject: [PATCH 13/16] Update conv_gefs2datm.f --- utilities/conv_gefs2datm/conv_gefs2datm.f | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/utilities/conv_gefs2datm/conv_gefs2datm.f b/utilities/conv_gefs2datm/conv_gefs2datm.f index a9c3553..2adf861 100644 --- a/utilities/conv_gefs2datm/conv_gefs2datm.f +++ b/utilities/conv_gefs2datm/conv_gefs2datm.f @@ -17,7 +17,8 @@ program conv_gefs2datm integer NX,NY,NT,nerr,NDIM,i,j,k,NSTR parameter (NX = 1536, NY = 768, NT = 1, NDIM = 3) real*4 :: latitude(NY), longitude(NX) - integer time(NT),NF_64BIT_OFFSET,vi_indx +C integer time(NT),NF_64BIT_OFFSET,vi_indx + integer time(NT),vi_indx real coef,Tcel,Hn,term1,term2,g0,Rsp real, dimension(NX,NY,NT) :: & dlwrf_avesfc, From 3089bd4f173c7da8a5842476407e740ed1ccc360 Mon Sep 17 00:00:00 2001 From: Hyun-Chul Lee -NOAA Date: Mon, 16 Mar 2020 12:12:26 -0400 Subject: [PATCH 14/16] Update comp_f77_code.csh --- utilities/conv_gefs2datm/comp_f77_code.csh | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/utilities/conv_gefs2datm/comp_f77_code.csh b/utilities/conv_gefs2datm/comp_f77_code.csh index 786b85f..3d443e9 100644 --- a/utilities/conv_gefs2datm/comp_f77_code.csh +++ b/utilities/conv_gefs2datm/comp_f77_code.csh @@ -1,8 +1,8 @@ #! /bin/csh module load intel/18.0.5.274 module load impi/2019.0.4 - +module load netcdf/4.6.1 set name = $1 #--- in hera -ifort $FFLAGS ${name}.f -o ${name} -I$NETCDF/include -L$NETCDF/lib -lnetcdff +ifort -132 $FFLAGS ${name}.f -o ${name} -I$NETCDF/include -L$NETCDF/lib -lnetcdff From 757dd7bbecec8f52bde6dde15511ae636bbf3288 Mon Sep 17 00:00:00 2001 From: flampouris Date: Mon, 16 Mar 2020 13:18:15 -0400 Subject: [PATCH 15/16] deleted 0readme replaced README.md --- utilities/conv_gefs2datm/0readme | 35 -------------------------------- 1 file changed, 35 deletions(-) delete mode 100644 utilities/conv_gefs2datm/0readme diff --git a/utilities/conv_gefs2datm/0readme b/utilities/conv_gefs2datm/0readme deleted file mode 100644 index 96af661..0000000 --- a/utilities/conv_gefs2datm/0readme +++ /dev/null @@ -1,35 +0,0 @@ -conv_gefs2datm - - -- Input data (GEFS) sources - -in HPSS : /3year/NCEPDEV/GEFSRR/gefsrr_replay_${YYYYstream}stream/${YYYYMMDDHH}.tar -in WCOSS : (venus)/gpfs/dell2/emc/modeling/noscrub/Hyun-Chul.Lee/GEFS/Reanal - -- Output data (DATM) archive -in HPSS : /5year/NCEPDEV/marineda/DATM_INPUT/GEFS/${YYYYMM}/gefs.${YYYYMMDDHH}.nc -in Hera : /scratch2/NCEPDEV/marineda/DATM_INPUT/GEFS/${YYYYMM}/gefs.${YYYYMMDDHH}.nc - - -- To conver from GEFS to DATM (netCDF), - -1) compile code -csh ./comp_f77_code.csh conv_gefs2datm - -2) modify the run script of conv_gefs2datm.fort.csh - -set sdir = !<-- source of GEFS -set tdir = !<-- output dir -set wdir = !<-- work dir -etc - -set symd = 20130529 !<-- start date -set eymd = 20190831 !<-- end date - -3) run the script -csh ./conv_gefs2datm.fort.csh - -- More comprensive information can be found - -https://docs.google.com/spreadsheets/d/1M32BvPgXuPqewQWO5cwsECLdcTj9o1x2Rq7GDvSpltk/edit#gid=0 - From 2ee0848015822e69c5d329822a84438ab0c71d63 Mon Sep 17 00:00:00 2001 From: flampouris <8270831+flampouris@users.noreply.github.com> Date: Mon, 16 Mar 2020 13:31:03 -0400 Subject: [PATCH 16/16] Updates for the PR7 to go through Necessary updates --- utilities/conv_gefs2datm/README.md | 32 +++++++++++++++++++ utilities/conv_gefs2datm/comp_f77_code.csh | 11 ++++--- utilities/conv_gefs2datm/conv_gefs2datm.f | 1 - .../conv_gefs2datm/conv_gefs2datm.fort.csh | 12 +++---- 4 files changed, 44 insertions(+), 12 deletions(-) create mode 100644 utilities/conv_gefs2datm/README.md diff --git a/utilities/conv_gefs2datm/README.md b/utilities/conv_gefs2datm/README.md new file mode 100644 index 0000000..3d99c4f --- /dev/null +++ b/utilities/conv_gefs2datm/README.md @@ -0,0 +1,32 @@ +# How to use conv_gefs2datm + + +## Input data (GEFS) sources + +At HPSS : `/3year/NCEPDEV/GEFSRR/gefsrr_replay_${YYYYstream}stream/${YYYYMMDDHH}.tar` +At WCOSS : `(venus)/gpfs/dell2/emc/modeling/noscrub/Hyun-Chul.Lee/GEFS/Reanal` + +## Output data (DATM) archive +At HPSS : `/5year/NCEPDEV/marineda/DATM_INPUT/GEFS/${YYYYMM}/gefs.${YYYYMMDDHH}.nc` +At Hera : `/scratch2/NCEPDEV/marineda/godas_input/DATM_INPUT/GEFS/${YYYYMM}/gefs.${YYYYMMDDHH}.nc` + +## To conver from GEFS to DATM (netCDF), + +1) Compile the code +`csh ./comp_f77_code.csh conv_gefs2datm` + +2) Modify the run script of conv_gefs2datm.fort.csh +`set sdir = !<-- source of GEFS +set tdir = !<-- output dir +set wdir = !<-- work dir ` + +Define the starting and ending date at YYYYMMDD format + `set symd = 20130529 !<-- start date` + `set eymd = 20190831 !<-- end date` + +3) Run the script +`csh ./conv_gefs2datm.fort.csh` + +## More information can be found + +https://docs.google.com/spreadsheets/d/1M32BvPgXuPqewQWO5cwsECLdcTj9o1x2Rq7GDvSpltk/edit#gid=0 diff --git a/utilities/conv_gefs2datm/comp_f77_code.csh b/utilities/conv_gefs2datm/comp_f77_code.csh index 3d443e9..f97fce9 100644 --- a/utilities/conv_gefs2datm/comp_f77_code.csh +++ b/utilities/conv_gefs2datm/comp_f77_code.csh @@ -1,8 +1,9 @@ #! /bin/csh -module load intel/18.0.5.274 -module load impi/2019.0.4 -module load netcdf/4.6.1 +module load intel +module load netcdf + set name = $1 -#--- in hera -ifort -132 $FFLAGS ${name}.f -o ${name} -I$NETCDF/include -L$NETCDF/lib -lnetcdff +#--- in hera +set FFLAGS='-extend-source 132' +ifort $FFLAGS ${name}.f -o ${name} -I$NETCDF/include -L$NETCDF/lib -lnetcdf \ No newline at end of file diff --git a/utilities/conv_gefs2datm/conv_gefs2datm.f b/utilities/conv_gefs2datm/conv_gefs2datm.f index 2adf861..1ddeb7f 100644 --- a/utilities/conv_gefs2datm/conv_gefs2datm.f +++ b/utilities/conv_gefs2datm/conv_gefs2datm.f @@ -17,7 +17,6 @@ program conv_gefs2datm integer NX,NY,NT,nerr,NDIM,i,j,k,NSTR parameter (NX = 1536, NY = 768, NT = 1, NDIM = 3) real*4 :: latitude(NY), longitude(NX) -C integer time(NT),NF_64BIT_OFFSET,vi_indx integer time(NT),vi_indx real coef,Tcel,Hn,term1,term2,g0,Rsp real, dimension(NX,NY,NT) :: diff --git a/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh b/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh index 32c6b6a..8d01639 100644 --- a/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh +++ b/utilities/conv_gefs2datm/conv_gefs2datm.fort.csh @@ -1,16 +1,17 @@ #! /bin/csh -x -module load netcdf/4.6.1 +module load intel +module load netcdf set hdir = #<-- set home dir set sdir = #<-- set input source dir set tdir = #<-- set dir for output set wdir = #<-- set work dir -#<-- set intial date +#<-- set intial date YYYYMMDD set symd = 20130529 -#<-- set end date -set eymd = 20190831 +#<-- set end date YYYYMMDD +set eymd = 20190831 cd $wdir @@ -31,5 +32,4 @@ ncdump gefs_output.nc | sed -e "5s#^.time = 1 ;#time = UNLIMITED ; // (1 current rm gefs_input.* end set ymd = `date -d "${ymd} 1 day" +%Y%m%d` -end - +end \ No newline at end of file