Hi.
I have a question about regridding of POP output to lat-lon grid. I
downloaded the POP weight file and I have it in the same directory where I
have the script. However, when I run it it looks like ncl can't find this
file or something....
Here's the error message:
pacman3 1162% ncl vector4.ncl
Copyright (C) 1995-2011 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.0.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
(0)
ls: map_*: No such file or directory
(0) create_rmpPopFileName: ./map_gx1v6_to_1x1d_bilin_da_100716.nc
ls: rot_*: No such file or directory
(0) create_rotPopFileName: file not recognized/found: rot_gx1v6.nc
(0) create_rotPopFileName: dir searched: ./
Here's my script:
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/popRemap.ncl"
;************************************************
begin
do i=15,31
in =
addfile("/import/wrkdir1/asemenov/archive/Gexp/ice/hist/Gexp.cice.h.0001-03-"+i+".nc","r")
u = in->uvel
v = in->vvel
remap = PopLatLonV(u,v,"gx1v6","1x1d","bilin","da","100716")
u_remap = remap(0,:,:,:)
v_remap = remap(1,:,:,:)
wks = gsn_open_wks("x11","pop2lat") ; open ps file
gsn_define_colormap(wks,"WhBlGrYeRe") ; choose colormap
i = NhlNewColor(wks,0.7,0.7,0.7) ; add gray to colormap
res = True
res@gsnPolar = "NH" ; choose hemisphere
res@mpLandFillColor = "grey" ; color of land
res@mpMinLatF = 50 ; adjust map limit
res@vcRefMagnitudeF = 5.0 ; vec magnitude
res@vcRefLengthF = 0.050 ; size of ref vector
res@vcMonoLineArrowColor = False
res@vcMinDistanceF = 0.007 ; thins arrows near pole
res@vcGlyphStyle = "CurlyVector" ; turn on curley vectors
res@gsnSpreadColors = True ; use full colormap
res@gsnSpreadColorStart = 2
res@gsnSpreadColorEnd = 97
res@tiMainString = "gx1v2 to 1x1d"
plot = gsn_csm_vector_map_polar(wks,u_remap(0,:,:), \
v_remap(0,:,:),res)
end do
end
-- regards ******************************************************* *Alex Semenov* * * PhD Student - Research Assistant International Arctic Research Center Department of atmospheric sciences University of Alaska Fairbanks 930 Koyukuk dr, 408c3 Fairbanks, AK, USA, 99775 work phone +19074742672
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jul 12 10:49:32 2012
This archive was generated by hypermail 2.1.8 : Wed Jul 18 2012 - 14:33:00 MDT