Re: POP remapping

From: David Bailey <dbailey_at_nyahnyahspammersnyahnyah>
Date: Thu Jul 12 2012 - 11:18:27 MDT

Alex,

I have added special routines PopLatLon2 and PopLatLonV2 that use
different arguments. With these you can pass the angle array as an
argument for the remapping.

Dave

On 07/12/2012 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

-- 
--
David A Bailey                             email: dbailey_at_ucar.edu
National Center for Atmospheric Research   phone: 303-497-1737
PO Box 3000                                fax  : 303-497-1700
Boulder, CO 80307-3000
              http://www.cgd.ucar.edu/staff/dbailey

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jul 12 11:18:38 2012

This archive was generated by hypermail 2.1.8 : Wed Jul 18 2012 - 14:33:00 MDT