;******************************************* ;******************************************* ; load ncl and user contributed libraries 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" load "$NCARG_ROOT/lib/ncarg/nclscripts/niwa/mapareas.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/niwa/maplabels.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/niwa/EcoConnectUtils.ncl" ;******************************************* begin ;******************************************* in = addfile("/data/riwrh/yangy/DRY_SOIL/DATA/2004010200/met00_2004010200_utc_nzlam_12.nc", "r") ter = in->orog_model ; terrain in0 = addfile("/data/riwrh/yangy/DRY_SOIL/DATA/2004010200/met02_2004010200_utc_nzlam_12.nc", "r") tlat = in0->latitude ; 2D true latitude tlon = in0->longitude ; 2D true longitude ter@lat2d = tlat ter@lon2d = tlon ;************************************************ ; plotting parameters ;************************************************ res2 = True ; plot resource (clouds) EcoConnectMapRes(res2) ; define a set of standard map resources res2@tiMainOn = False res2@mpGeophysicalLineColor = "black" ; map outline colour res2@cnFillOn = True res2@cnLinesOn = True res2@mpLabelFontColor = "burlywood" res2@mpGridLineColor = "burlywood" ; grid line colour res2@gsnSpreadColors = True ; we will define explicit fill colours res2@gsnSpreadColorStart = 70 ; start colour res2@gsnSpreadColorEnd = 190 ; end colour res2@cnFillMode = "AreaFill" res2@lbLabelStride = 1 ; colour bar label stride res2@lbTitleString = "Terrain height (m)" res2@cnSmoothingOn = False res2@cnSmoothingDistanceF = 0.03 res2@mpFillDrawOrder = "PreDraw" res2@lbTitleString = "Terrain height (m)" res2@mpGridAndLimbOn = False res2@cnLineColor = "black" ; white, dark green (eco_colourmap_1) res2@cnLevelSelectionMode = "ExplicitLevels" ; set explicit contour levels res2@cnLevels = (/ 0.1,100.,200.,300.,400.,500.,600./) res2@cnLineThicknesses = (/ 1,1,1,1,1,1,1/) res2@gsnContourLineThicknessesScale = 1.5 res2@cnLinesOn = True ; turn off the contour lines res2@cnLineLabelsOn = True res2@cnLineLabelInterval = 1 res2@cnInfoLabelOn = False res2@mpDataBaseVersion = "MediumRes" ; Will be fine for this area & reduces file size ;************************************************ ; e-nzlam12-lc area ;************************************************ wks = gsn_open_wks ("eps","Plot1") ; open workstation gsn_define_colormap(wks,"eco_colourmap_1") ; grey scale for clouds res2@vpHeightF = 0.705 ; view port width res2@mpProjection = "LambertConformal" res2@mpLimitMode = "Corners" ; choose range of map res2@mpLambertParallel1F = -43.012 ; map parameter res2@mpLambertParallel2F = -38.024 ; map parameter res2@mpLambertMeridianF = 172.911 ; map parameter res2@mpLeftCornerLatF = -39.8000 ; BLC -48.0000 165.0000 res2@mpLeftCornerLonF = 173.5 res2@mpRightCornerLatF = -38.8 ; TRC -33.0326 180.0035 res2@mpRightCornerLonF = 174.5 res2@mpGridLonSpacingF = 2.5 ; longitude grid line spacing res2@mpGridLatSpacingF = 2.5 ; res2@mpDataBaseVersion = "MediumRes" ; Will be fine for this area & reduces file size res2@gsnLeftString = "Plot1 " res2@mpOutlineBoundarySets ="NoBoundaries" plot = gsn_csm_contour_map(wks,ter,res2) ; generate the base map & overlay draw(plot) ; finally draw plot and advance frame maximize_output(wks,True) delete(plot) ; free resources delete(wks) ;************************************************ ; e-nz1-lc area ;************************************************ end