POP remapping

From: Alexander Semenov <asemenov_at_nyahnyahspammersnyahnyah>
Date: Thu Jul 12 2012 - 10:49:22 MDT

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