wgrib2 Z__C_RJTD_20191226000000_MSM_GPV_Rjp_Lsurf_FH00-15_grib2.bin | grep “TMP” | wgrib2 -i -no_header Z__C_RJTD_20191226000000_MSM_GPV_Rjp_Lsurf_FH00-15_grib2.bin -no_append -ieee output.bin
wgrib2 Z__C_RJTD_20171205000000_MSM_GPV_Rjp_Lsurf_FH00-15_grib2.bin | grep “APCP” | wgrib2 -i -no_header Z__C_RJTD_20171205000000_MSM_GPV_Rjp_Lsurf_FH00-15_grib2.bin -no_append -ieee outputAPCP.bin
wgrib2 Z__C_RJTD_20171205000000_MSM_GPV_Rjp_Lsurf_FH00-15_grib2.bin | grep “PRMSL” | wgrib2 -i -no_header Z__C_RJTD_20171205000000_MSM_GPV_Rjp_Lsurf_FH00-15_grib2.bin -no_append -ieee outputPRMSL.bin
$ vi grads.ctl
dset output.bin
options big_endian
undef 9.999E+20
title MSM-GPV temperature
xdef 481 linear 120.0 0.0625
ydef 505 linear 22.4 0.05
zdef 1 levels 1000.0
tdef 16 linear 00Z26DEC2019 1hr
vars 1
tmp 1 0 temperature
endvars
open grads.ctl
set mpdset hires
set t 1
set gxout grfill
set gxout contour
open HECB88.ctl
set lat 25 50
set lon 120 150
enable print tmo
set gxout shaded
set clevs 0.5 5 10 20 30 50 80 120 200
set ccols 15 5 4 3 7 8 6 2 9
set csmooth on
d apcpsfc.1
ga-> set clevs 0 5 10 15 20
ga-> set ccols 1 2 3 4 5
outputAPCP.bin
outputPRMSL.bin
outputTMP.bin
gradsPRMSL.ctl
gradsAPCP.ctl
gradsTMP.ctl
💚ある値以下（以上）をマスクして表示
IT memo/GrADS memo． †
d maskout(var,var-10) # 10以上の値を表示
💚ファイル番号を指定して描画する
デフォルトのファイル番号を変更する
GrADSのTips
ga-> d U.1
ga-> d V.2
ga-> set dfile 1 // デフォルトのファイルを1に指定
ga-> d u // Uを描画
ga-> set dfile 2 // デフォルトのファイルを2に指定
ga-> d v
| 記号 | 物理量 |
|---|---|
| PRMSL | 海面更生気圧 |
| PRES | 地上気圧 |
| UGRD | 東西風速 |
| VGRD | 南北風速 |
| TMP | 気温 |
| RH | 相対湿度 |
| LCDC | 下層雲量 |
| MCDC | 中層雲量 |
| HCDC | 上層雲量 |
| TCDC | 全雲量 |
| APCP | 降水量 |
| DSWRF | 日射量 |
open gradsPRMSL.ctl
set mpdset hires
set t 1
set gxout contour
d prmsl
open gradsAPCP.ctl
set gxout grfill
d maskout(prmsl.2,prmsl.2-0.01)
d maskout(apcp.2,apcp.2-0.01)
/home/office70/tmp/gradsdods -b
printim /home/office70/www/office70/images/test.gif
open /home/office70/scripts/grads.ctl
set t 1
set display color white
set gxout grfill
d A
printim /home/office70/www/office70/images/test.gif
MAC版
‘open gradsPRMSL.ctl’
‘open gradsAPCP.ctl’
‘set mpdset hires’
‘set t 1’
#——————–
‘set gxout contour’
# 0.01 倍して hPa に換算
‘prmslh=0.01*prmsl.1’
‘d prmslh’
#——————–
‘set gxout shaded’
‘set clevs 0.05 0.5 5 10 20 30 50 80 120 200’
‘set ccols 15 5 4 3 7 8 6 2 9’
‘set csmooth on’
‘d maskout(apcp.2,apcp.2-0.05)’
‘gxprint out.png’
サーバー版
/home/office70/scripts
以下MAC版を、２か所変更！！
‘open /home/office70/scripts/grads.ctl’
‘printim /home/office70/www/office70/images/test.gif’
２０２２ 令和４年 ７月９日（土曜日）
Shared object “libgfortran.so.3” not found, required by “wgrib2”
/home/office70/scripts/grib2tobin_PRMSL.sh: line 2: grep: command not found
Shared object “libgfortran.so.3” not found, required by “wgrib2”
find 検索先 -name ファイル名
% find ~/ -name libgfortran.so.3
/home/office70/local/lib/libgfortran.so.3
which curl
/usr/local/bin/curl
which gfortran
/home/office70/local/bin/gfortran
/home/office70/local/bin:/home/office70/local/lib/:/sbin:/bin:/usr/sbin:/usr/bin:/usr/local/sbin:/usr/local/bin:/home/office70/bin
% chsh -s /usr/local/bin/bash
Password:
chsh: Permission denied
/bin/csh
setenv EDITOR vi
setenv PAGER more
setenv BLOCKSIZE K
setenv PKG_DBDIR ~/db/pkg
setenv C_INCLUDE_PATH $HOME/local/include
setenv LIBRARY_PATH $HOME/local/lib
setenv LD_LIBRARY_PATH $HOME/local/lib
setenv CC /home/office70/local/bin/gcc
setenv FC /home/office70/local/bin/gfortran
setenv GADDIR /home/office70/tmp/grads-2.0.2/data
setenv GASCRP /home/office70/tmp/grads-2.0.2/scripts
/bin/csh /home/office70/scripts/grib2tobin_PRMSL.sh
PATH=/home/office70/local/bin:/home/office70/local/lib/:/sbin:/bin:/usr/sbin:/usr/bin:/usr/local/sbin:/usr/local/bin:/home/office70/bin
/home/office70/local/grib2/wgrib2/wgrib2 /home/office70/downloads/downloaddata_Lsurf.bin | grep “APCP” | /home/office70/local/grib2/wgrib2/wgrib2 -i -no_header /home/office70/downloads/downloaddata_Lsurf.bin -no_append -ieee /home/office70/downloads/outputAPCP.bin