Re: POP remapping

The error messsage says that NCL is looking for the file rot_gx1v6.nc

> ls: rot_*: No such file or directory
> (0)create_rotPopFileName: file not recognized/found: rot_gx1v6.nc

The original PopLatLonV was written 12 years ago. There were several
POP models. They either did not include the necessary rotation
angles for vector rotation or they had different variable names

   rotaName = (/"ANGLE", "angle", "UTAN", "utan" \
                ,"ROTA", "rota" , "ROTANG", "rotang" /)

---
Now, POP files have "ANGLE". I think you could create
your own 'rot_gx1v6.nc' via
Pick any "Gexp.cice.h.0001-03-???.nc" file
%> ncks -v ANGLE,ANGLET input_file.nc rot_gx1v6.nc
This needs to be done once.
On 7/12/12 10:49 AM, Alexander Semenov wrote:
> 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
> <http://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
> <http://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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk