[2018-04-03_17:19]
[/work2/am/2017.KYUSHU-HOKUBU.HEAVY.RAIN/FNL0.25.FLUX.CONVERT]
[am@localhost]
$ fnl0.25.flux.convert.sh
Grid Analysis and Display System (GrADS) Version 2.1.a3
Copyright (C) 1988-2015 by the Institute for Global Environment and Society (IGES)
GrADS comes with ABSOLUTELY NO WARRANTY
See file COPYRIGHT for more information
Config: v2.1.a3 little-endian readline grib2 netcdf hdf4-sds hdf5 opendap-grids,stn geotiff shapefile cairo
Issue 'q config' command for more detailed configuration information
GX Package Initialization: Size = 11 8.5
Running in Batch mode
Default file number is: 1
X is varying Lon = 115 to 145 X = 461 to 581
Y is varying Lat = 20 to 45 Y = 441 to 541
Z is fixed Lev = 1000 Z = 1
T is varying Time = 00Z30MAY2017 to 18Z13JUL2017 T = 1 to 180
E is fixed Ens = 1 E = 1
COMPUTING CLOUD WATER PATH (CWP)...
OUTPUT: NCEP_fnl0.25_CWP.nc
COMPUTING SHORT-WAVE RADIATION (SWR) ...
OUTPUT: NCEP_fnl0.25_SWR.nc
COMPUTING LONG-WAVE RADIATION (LWR) ...
OUTPUT: NCEP_fnl0.25_LWR.nc
COMPUTING SENSIBLE HEAT FLUX (SHF) ...
OUTPUT: NCEP_fnl0.25_SHF.nc
COMPUTING LATENT HEAT FLUX (LHF) ...
OUTPUT: NCEP_fnl0.25_LHF.nc
COMPUTING NET HEAT FLUX (NET) ...
OUTPUT: NCEP_fnl0.25_NET.nc
Default file number is: 1
X is varying Lon = 115 to 145 X = 461 to 581
Y is varying Lat = 20 to 45 Y = 441 to 541
Z is fixed Lev = 1000 Z = 1
T is varying Time = 00Z30MAY2017 to 18Z13JUL2017 T = 1 to 180
E is fixed Ens = 1 E = 1
GX Package Terminated
[2018-04-03_17:21]
[/work2/am/2017.KYUSHU-HOKUBU.HEAVY.RAIN/FNL0.25.FLUX.CONVERT]
[am@localhost]
$ ncdump -h NCEP_fnl0.25_CWP.nc
netcdf NCEP_fnl0.25_CWP {
dimensions:
lon = 121 ;
lat = 101 ;
time = 180 ;
variables:
double lon(lon) ;
lon:units = "degrees_east" ;
lon:long_name = "Longitude" ;
double lat(lat) ;
lat:units = "degrees_north" ;
lat:long_name = "Latitude" ;
double time(time) ;
time:long_name = "Time" ;
time:units = "minutes since 2017-05-30 00:00" ;
double cwp(time, lat, lon) ;
cwp:_FillValue = -999000000. ;
}
[2018-04-03_17:24]
[/work2/am/2017.KYUSHU-HOKUBU.HEAVY.RAIN/FNL0.25.FLUX.CONVERT]
[am@localhost]
$ ncdump -h NCEP_fnl0.25_SWR.nc
netcdf NCEP_fnl0.25_SWR {
dimensions:
lon = 121 ;
lat = 101 ;
time = 180 ;
variables:
double lon(lon) ;
lon:units = "degrees_east" ;
lon:long_name = "Longitude" ;
double lat(lat) ;
lat:units = "degrees_north" ;
lat:long_name = "Latitude" ;
double time(time) ;
time:long_name = "Time" ;
time:units = "minutes since 2017-05-30 00:00" ;
double swr(time, lat, lon) ;
swr:_FillValue = -999000000. ;
}
[2018-04-03_17:24]
[/work2/am/2017.KYUSHU-HOKUBU.HEAVY.RAIN/FNL0.25.FLUX.CONVERT]
[am@localhost]
$ ncdump -h NCEP_fnl0.25_LWR.nc
netcdf NCEP_fnl0.25_LWR {
dimensions:
lon = 121 ;
lat = 101 ;
time = 180 ;
variables:
double lon(lon) ;
lon:units = "degrees_east" ;
lon:long_name = "Longitude" ;
double lat(lat) ;
lat:units = "degrees_north" ;
lat:long_name = "Latitude" ;
double time(time) ;
time:long_name = "Time" ;
time:units = "minutes since 2017-05-30 00:00" ;
double lwr(time, lat, lon) ;
lwr:_FillValue = -999000000. ;
}
[2018-04-03_17:24]
[/work2/am/2017.KYUSHU-HOKUBU.HEAVY.RAIN/FNL0.25.FLUX.CONVERT]
[am@localhost]
$ ncdump -h NCEP_fnl0.25_SHF.nc
netcdf NCEP_fnl0.25_SHF {
dimensions:
lon = 121 ;
lat = 101 ;
time = 180 ;
variables:
double lon(lon) ;
lon:units = "degrees_east" ;
lon:long_name = "Longitude" ;
double lat(lat) ;
lat:units = "degrees_north" ;
lat:long_name = "Latitude" ;
double time(time) ;
time:long_name = "Time" ;
time:units = "minutes since 2017-05-30 00:00" ;
double shf(time, lat, lon) ;
shf:_FillValue = -999000000. ;
}
[2018-04-03_17:24]
[/work2/am/2017.KYUSHU-HOKUBU.HEAVY.RAIN/FNL0.25.FLUX.CONVERT]
[am@localhost]
$ ncdump -h NCEP_fnl0.25_LHF.nc
netcdf NCEP_fnl0.25_LHF {
dimensions:
lon = 121 ;
lat = 101 ;
time = 180 ;
variables:
double lon(lon) ;
lon:units = "degrees_east" ;
lon:long_name = "Longitude" ;
double lat(lat) ;
lat:units = "degrees_north" ;
lat:long_name = "Latitude" ;
double time(time) ;
time:long_name = "Time" ;
time:units = "minutes since 2017-05-30 00:00" ;
double lhf(time, lat, lon) ;
lhf:_FillValue = -999000000. ;
}
[2018-04-03_17:24]
[/work2/am/2017.KYUSHU-HOKUBU.HEAVY.RAIN/FNL0.25.FLUX.CONVERT]
[am@localhost]
$ ncdump -h NCEP_fnl0.25_NET.nc
netcdf NCEP_fnl0.25_NET {
dimensions:
lon = 121 ;
lat = 101 ;
time = 180 ;
variables:
double lon(lon) ;
lon:units = "degrees_east" ;
lon:long_name = "Longitude" ;
double lat(lat) ;
lat:units = "degrees_north" ;
lat:long_name = "Latitude" ;
double time(time) ;
time:long_name = "Time" ;
time:units = "minutes since 2017-05-30 00:00" ;
double net(time, lat, lon) ;
net:_FillValue = -999000000. ;
}
[2018-04-03_17:24]
[/work2/am/2017.KYUSHU-HOKUBU.HEAVY.RAIN/FNL0.25.FLUX.CONVERT]
[am@localhost]
$ srcdump.sh fnl0.25.flux.convert.sh
#------------------------------
# List of the following files:
#------------------------------
fnl0.25.flux.convert.sh
#------------------------------
# Machine info
#------------------------------
localhost
/work2/am/2017.KYUSHU-HOKUBU.HEAVY.RAIN/FNL0.25.FLUX.CONVERT
Tue Apr 3 17:24:52 JST 2018
#======================
# fnl0.25.flux.convert.sh
#======================
#!/bin/bash
gs=$(basename $0 .sh).gs
figdir=Fig_$(basename $0 .sh)
mkdir -vp $figdir
dataset_name='NCEP_fnl0.25'
ctl=fnl0.25.ctl
export LANG=C
host=$(hostname)
cwd=$(pwd)
timestamp=$(date)
lonw=115
lone=145
lats=20
latn=45
datetimes=00Z30MAY2017
datetimee=18Z13JUL2017
cat <<EOF>$gs
'open ${ctl}'
'set time ${datetimes} ${datetimee}'
'set lon ${lonw} ${lone}'
'set lat ${lats} ${latn}'
'set z 1'
'q dims'
say result
title=CWP
say 'COMPUTING CLOUD WATER PATH ('title')...'
'define cwp=vint(lev(z=1),CLWMRprs,100)'
fname='${dataset_name}_'title'.nc'
say 'OUTPUT: 'fname
'set sdfwrite 'fname
'sdfwrite cwp'
title=SWR
say 'COMPUTING SHORT-WAVE RADIATION ('title') ...'
'define swr=DSWRFsfc-USWRFsfc'
fname='${dataset_name}_'title'.nc'
say 'OUTPUT: 'fname
'set sdfwrite 'fname
'sdfwrite swr'
title=LWR
say 'COMPUTING LONG-WAVE RADIATION ('title') ...'
'define lwr=DLWRFsfc-ULWRFsfc'
fname='${dataset_name}_'title'.nc'
say 'OUTPUT: 'fname
'set sdfwrite 'fname
'sdfwrite lwr'
title=SHF
say 'COMPUTING SENSIBLE HEAT FLUX ('title') ...'
'define shf=-SHTFLsfc'
fname='${dataset_name}_'title'.nc'
say 'OUTPUT: 'fname
'set sdfwrite 'fname
'sdfwrite shf'
title=LHF
say 'COMPUTING LATENT HEAT FLUX ('title') ...'
'define lhf=-LHTFLsfc'
fname='${dataset_name}_'title'.nc'
say 'OUTPUT: 'fname
'set sdfwrite 'fname
'sdfwrite lhf'
title=NET
say 'COMPUTING NET HEAT FLUX ('title') ...'
'define net=DSWRFsfc-USWRFsfc+DLWRFsfc-ULWRFsfc-SHTFLsfc-LHTFLsfc'
fname='${dataset_name}_'title'.nc'
say 'OUTPUT: 'fname
'set sdfwrite 'fname
'sdfwrite net'
'q dims'
say result
quit
EOF
grads -bcl "${gs}"
exit 0
#----------------------
# End of fnl0.25.flux.convert.sh
#----------------------