[2015年 8月 13日 木曜日 15:56:30 JST]
[~/TEACHING/Outreach/141002.Nagasaki.Minami/Velocity/ECCO]
[am@aofd165]
流速ベクトルの大きさ(m/s)の水平分布
流速ベクトルの大きさ(m/s)の水平分布
::::::::::::::
get.ecco.uv.sh
::::::::::::::
#!/bin/sh
wget -nc ftp://ecco2.jpl.nasa.gov/data1/cube/cube92/lat_lon/quart_90S_90N/readme.txt
urlu=ftp://ecco2.jpl.nasa.gov/data1/cube/cube92/lat_lon/quart_90S_90N/UVEL_monthly.nc/
urlv=ftp://ecco2.jpl.nasa.gov/data1/cube/cube92/lat_lon/quart_90S_90N/VVEL_monthly.nc/
#wget $urlu/UVEL.1440x720x50.201212.nc
wget $urlv/VVEL.1440x720x50.201212.nc
::::::::::::::
pl.uv.ncl
::::::::::::::
;
; Plot Current Speed
;
;
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
fu=addfile("UVEL.1440x720x50.201212.nc","r")
fv=addfile("VVEL.1440x720x50.201212.nc","r")
u = fu->UVEL
v = fv->VVEL
printVarSummary(u)
wks = gsn_open_wks("ps","uv_201212")
gsn_define_colormap(wks,"BlueDarkRed18")
umag=u(0,0,:,:)
umag=sqrt(u(0,0,:,:)^2+v(0,0,:,:)^2)
umag@_FillValue=-9.9999998e+22
do j=0,720-1
do i=0,1440-1
if(umag(j,i) .eq. 0)then
umag(j,i)=umag@_FillValue
end if
end do
end do
;print(umag)
;gsn_define_colormap(wks,"MPL_gnuplot2")
res = True
res@mpFillOn = False
res@gsnFrame = False
res@vpWidthF = 0.5*1.9
res@vpHeightF = 0.32*1.9
res@vpXF = .1*0.1
; res@vpYF = .75
res@mpProjection = "Robinson" ; choose projection
res@mpFillOn = False
res@mpGridAndLimbOn = True ; turn on lat/lon lines
res@mpPerimOn = False
res@mpGridLatSpacingF = 90 ; change latitude line spacing
res@mpGridLonSpacingF = 180. ; change longitude line spacing
res@mpGridLineColor = "transparent" ; trick ncl into drawing
; perimeter
res@mpCenterLonF = 200.
res@cnFillOn = True ; color plot desired
res@cnLinesOn = False ; turn off/on contour lines
res@cnLineLabelsOn = False ; turn off contour lines
res@cnLineColor = "Background"
res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
res@cnMinLevelValF = 0.1 ; set min contour level
res@cnMaxLevelValF = 1.3 ; set max contour level
res@cnLevelSpacingF = 0.2 ; set contour spacing
res@lbLabelFontHeightF = 0.02
res@pmLabelBarHeightF = 0.07
res@pmLabelBarWidthF = 0.6
res@pmLabelBarOrthogonalPosF = 0.025
res@gsnLeftString = ""
res@gsnCenterString = ""
res@gsnRightString = ""
res@mpDataBaseVersion = "Ncarg4_1" ; use GMT coastline
plot = gsn_csm_contour_map(wks,umag,res) ; Draw a contour plot.
end
::::::::::::::
pl.uv.pacific.ncl
::::::::::::::
;
; Plot Current Speed
;
;
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
fu=addfile("UVEL.1440x720x50.201212.nc","r")
fv=addfile("VVEL.1440x720x50.201212.nc","r")
u = fu->UVEL
v = fv->VVEL
printVarSummary(u)
wks = gsn_open_wks("ps","uv_201212.Pacific")
gsn_define_colormap(wks,"BlueDarkRed18")
umag=u(0,0,:,:)
umag=sqrt(u(0,0,:,:)^2+v(0,0,:,:)^2)
umag@_FillValue=-9.9999998e+22
do j=0,720-1
do i=0,1440-1
if(umag(j,i) .eq. 0)then
umag(j,i)=umag@_FillValue
end if
end do
end do
;print(umag)
;gsn_define_colormap(wks,"MPL_gnuplot2")
res = True
res@mpFillOn = False
res@gsnFrame = False
res@vpWidthF = 0.5*1.9
res@vpHeightF = 0.32*1.9
res@vpXF = .1*0.1
; res@vpYF = .75
; res@mpProjection = "Robinson" ; choose projection
res@mpMaxLatF = 50 ; choose subregion
res@mpMinLatF = 15
res@mpMaxLonF = 155
res@mpMinLonF = 120
res@mpFillOn = False
res@mpGridAndLimbOn = True ; turn on lat/lon lines
res@mpPerimOn = False
res@mpGridLatSpacingF = 90 ; change latitude line spacing
res@mpGridLonSpacingF = 180. ; change longitude line spacing
res@mpGridLineColor = "transparent" ; trick ncl into drawing
; perimeter
res@mpCenterLonF = 200.
res@cnFillOn = True ; color plot desired
res@cnLinesOn = False ; turn off/on contour lines
res@cnLineLabelsOn = False ; turn off contour lines
res@cnLineColor = "Background"
res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
res@cnMinLevelValF = 0.1 ; set min contour level
res@cnMaxLevelValF = 1.3 ; set max contour level
res@cnLevelSpacingF = 0.2 ; set contour spacing
res@lbLabelFontHeightF = 0.02
res@pmLabelBarHeightF = 0.07
res@pmLabelBarWidthF = 0.6
res@pmLabelBarOrthogonalPosF = 0.2 ; 25
res@gsnLeftString = ""
res@gsnCenterString = ""
res@gsnRightString = ""
res@mpDataBaseVersion = "Ncarg4_1" ; use GMT coastline
plot = gsn_csm_contour_map(wks,umag,res) ; Draw a contour plot.
end