如何获得CALIPSO数据,

发布网友 发布时间：2022-05-30 02:09

我来回答

共2个回答

热心网友 时间：2023-09-16 06:00

代码如下：

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/csm/skewt_func.ncl" 
load 
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl" 


begin
dir=("./DATA/")
fi=addfile(dir+"./CAL_LID_L2_VFM-ValStage1-V3-01.2007-02-06T16-29-28ZN.hdf","r")

lat 
= fi->Latitude(:,:)
lon = fi->Longitude(:,:)
time = 
fi->Profile_UTC_Time(:,:)
aerosol = 
fi->Feature_Classification_Flags(:,:)
aerosol@units 
="categroy"

wks  = gsn_open_wks("png","gsn_contour")    ; open a ps 
file
gsn_define_colormap(wks,"rainbow")

res                       = 
True
  res@gsnMaximize         = True
  res@gsnLeftString       = 
"Vertical Feature Mask"
  res@tiYAxisString          = " Altitude 
(km)"
  res@tiYAxisFontHeightF  = 0.02
  res@cnLinesOn           = 
False
  res@cnFillOn            = True
  res@cnFillMode          = 
"RasterFill"    ; Much faster!
  res@gsnAddCyclic        = 
False
  res@cnLevelSelectionMode = "ExplicitLevels"
  res@cnLevels        
     = (/0,1,2,3,4,5,6/)
  res@gsnSpreadColors     = 
True
  res@gsnSpreadColorStart = 1
  res@gsnSpreadColorEnd   = -3  ; don't 
include gray or white at end
;---Tickmark resources
  latstep              
    = 10
  res@tmXBMode             = "Explicit"
  res@tmXBValues        
   = lat(::latstep)
  res@tmXBLabels           = 
sprintf("%5.2f",lat(::latstep)) + "~C~" + \
                          
   sprintf("%5.2f",lon(::latstep))
;---Tickmark 
resources
  res@lbOrientation      = "Vertical"
  res@lbLabelAlignment   = 
"BoxCenters"     ; Default is "InteriorEdges"
  res@lbLabelAutoStride  = 
True
  res@lbLabelFontHeightF = 0.015
  res@lbLabelStrings     = (/"Not 
Determined","Clean Marine","Dust","Polluted Continental",\
                  
           "Clean Continental","Polluted Dust","Smoke","Other"/)
  plot = 
gsn_csm_contour(wks,aerosol(:,:),res)    ; Draw a contour plot.

end

热心网友 时间：2023-09-16 06:00

获得CALIPSO数据代码只有唯一的一组代码，如同长方形的面积公式就是长乘宽一样，代码如下：

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/csm/skewt_func.ncl"
load
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"

begin
dir=("./DATA/")
fi=addfile(dir+"./CAL_LID_L2_VFM-ValStage1-V3-01.2007-02-06T16-29-28ZN.hdf","r")

lat
= fi->Latitude(:,:)
lon = fi->Longitude(:,:)
time =
fi->Profile_UTC_Time(:,:)
aerosol =
fi->Feature_Classification_Flags(:,:)
aerosol@units
="categroy"

wks = gsn_open_wks("png","gsn_contour") ; open a ps
file
gsn_define_colormap(wks,"rainbow")

res =
True
res@gsnMaximize = True
res@gsnLeftString =
"Vertical Feature Mask"
res@tiYAxisString = " Altitude
(km)"
res@tiYAxisFontHeightF = 0.02
res@cnLinesOn =
False
res@cnFillOn = True
res@cnFillMode =
"RasterFill" ; Much faster!
res@gsnAddCyclic =
False
res@cnLevelSelectionMode = "ExplicitLevels"
res@cnLevels
= (/0,1,2,3,4,5,6/)
res@gsnSpreadColors =
True
res@gsnSpreadColorStart = 1
res@gsnSpreadColorEnd = -3 ; don't
include gray or white at end
;---Tickmark resources
latstep
= 10
res@tmXBMode = "Explicit"
res@tmXBValues
= lat(::latstep)
res@tmXBLabels =
sprintf("%5.2f",lat(::latstep)) + "~C~" + \

sprintf("%5.2f",lon(::latstep))
;---Tickmark
resources
res@lbOrientation = "Vertical"
res@lbLabelAlignment =
"BoxCenters" ; Default is "InteriorEdges"
res@lbLabelAutoStride =
True
res@lbLabelFontHeightF = 0.015
res@lbLabelStrings = (/"Not
Determined","Clean Marine","Dust","Polluted Continental",\

"Clean Continental","Polluted Dust","Smoke","Other"/)
plot =
gsn_csm_contour(wks,aerosol(:,:),res) ; Draw a contour plot.

end

数据就是数值，也就是我们通过观察、实验或计算得出的结果。数据有很多种，最简单的就是数字。数据也可以是文字、图像、声音等。数据可以用于科学研究、设计、查证等。