[ncl-talk] Eady _growth_rate

[Sri Nandini]

​Hello

I am trying to plot the eady growth rate as a measure of baroclinic instability. I cannot due to: 

fatal:conform: The array to be conformed must have the same number of dimensions as indicated by the length of the last argument
I cannot conform the latitude to same dimensions as the th variable.

Below is my script, would be grateful for any suggestions.


load "lat_lon_masked_lambert.ncl"

; ==============================================================
; Open the file: Read only the user specified period first observations then model
; ==============================================================
   f= addfile("SLP_PHIS_Z3_PI.nc", "r") ;
   T41    = f->Z3(:,:,:,:)

   printVarSummary(T41)                                 ; ; [Time|12]x[ilev | 26] x [lat | 96] x [lon | 144]
   T41 at _FillValue = -9.96921e+36     

;==============================================================

   f1= addfile("totalwinds_PI.nc", "r") ;
   T411    = f1->U(:,:,:,:)

   printVarSummary(T411)                                 ; ; [Time|12]x[ilev | 26] x [lat | 96] x [lon | 144]
   T411 at _FillValue = -9.96921e+36     


;==============================================================
   f2= addfile("th_PI.nc", "r") ;
   T412    = f2->TH(:,:,:,:)

   printVarSummary(T412)                                 ; ; [Time|12]x[ilev | 27] x [lat | 96] x [lon | 144]
   T412 at _FillValue = -9.96921e+36     

;==============================================================
   
;    Read latitudes 
;    The 'eady_growth_rate' function requires that 'lat' and 'th' agree
;    Use 'conform' the propogate the lat values
;==============================================================

   xlat = f->lat                               ;  [lat | 96] 
   printVarSummary(xlat)

   XLAT = conform(T412, xlat, (/0,2,3/))     ; problem here
   printVarSummary(XLAT)

   egr = eady_growth_rate(aveX2, aveX1, aveX, XLAT, 0,  1)
   printVarSummary(egr)
   printMinMax(egr, 0)

;==============================================================
   wks = gsn_open_wks("pdf","Eady")              ; send graphics to PNG file

;---Set some basic plot options

   res               = True
   res at gsnMaximize   = True       ; maximize plot in frame
   res at gsnAddCyclic  = False

   res at cnFillOn      = True  
   res at cnLinesOn     = False
  ;res at cnFillMode    = "RasterFill"                 ; slow here
   res at cnFillMode    = "CellFill"                   ; faster


  minlat = 25.                          ; min lat to mask
  maxlat = 80.                          ; max lat to mask
  minlon = -10.                          ; min lon to mask
  maxlon =  110.                          ; max lon to mask
  res at mpProjection = "LambertConformal"            ; choose projection

;---masked plot
  res at gsnAddCyclic = True                ; regional plot

  res at mpMinLatF = minlat              ; min lat to mask
  res at mpMaxLatF = maxlat              ; max lat to mask
  res at mpMinLonF = minlon              ; min lon to mask
  res at mpMaxLonF = maxlon              ; max lon to mask

  res at gsnMaskLambertConformal = True                ; turn on lc masking

; specify a level or levels ... within boundary layer

   klStrt = 5
   klLast = 5
   nt     = 0

;--- Eady growth rate (1/day)

   egr        = egr*86400
   egr at units  = "1/day"

   res at cnFillPalette        = "precip2_17lev"
   res at cnLevelSelectionMode = "ManualLevels"     ; set manual contour levels
   res at cnMinLevelValF       =  0.5               ; set min contour level
   res at cnMaxLevelValF       =  4.0               ; set max contour level
   res at cnLevelSpacingF      =  0.25              ; set contour spacing

   do kl=klStrt,klLast
      res at gsnCenterString = "znu="+znu(nt,kl)
      res at gsnRightString  =  egr at units
      contour = gsn_csm_contour_map(wks,  egr(nt,kl,:,:),res)
   end do






-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180817/2954bc8b/attachment.html>