没有合适的资源？快使用搜索试试~ 我知道了~

首页idl在气象学中的实用应用：工具与实例解析

idl在气象学中的实用应用：工具与实例解析

需积分: 10 0 下载量 123 浏览量更新于2024-07-18 收藏 31.92MB PDF 举报

身份认证购VIP最低享 7 折!

30元优惠券

试读

135页

"Using IDL in Meteorology 5: An In-Depth Guide for Practitioners" 本指南深入探讨了如何在气象学领域中有效利用Interactive Data Language (IDL)。IDL,作为一种强大的数值计算和数据处理工具，已经成为气象学家不可或缺的辅助工具。本篇内容着重于提供实际应用案例和教程，旨在帮助读者理解并掌握如何通过IDL进行气象数据的分析、模型模拟以及可视化。首先，作者Andy Heaps，来自雷丁大学气象学部门，对NCAS（国家中心大气研究）表示感谢，他们的支持使得这本指南得以改进。众多气象专家如Chris Bell, Paul Berrisford, Charles Boschetto等人的专业知识和建议在本书中发挥了关键作用。他们不仅提供了数据和示例图，还积极参与讨论，解答了棘手的问题，从而丰富了内容。书中涵盖了核心的IDL脚本和程序，包括但不限于数据处理、天气模型的运行、气候数据分析、以及气象图像的生成，如Roger Brugge的ephigram绘制代码。这些实用的工具帮助读者在气象研究中提高效率，处理复杂的数据集，比如气温、风速、降水量等多变量数据。此外，指南还提到了pvwave-idl新sgroup和ITT-VIS支持的重要性，它们为用户提供了丰富的讨论平台和技术支援，解决了在实际工作中遇到的难题。一个生动的例子是介绍了一张照片，展示了在Tempelfjorden, Spitsbergen的雪地摩托队前往Noorderlicht schooner的情景，这与IDL在实时天气监测中的应用相呼应。 "Using IDL in Meteorology 5"是一本实用的教程，它详细介绍了如何运用IDL技术解决气象学中的各种问题，无论是基本的数据预处理，还是高级的气候模型分析，都能为气象工作者提供宝贵的指导和参考。通过阅读此书，读者将能提升自己的气象数据分析技能，进而推动气象预测的精确性和时效性。

资源详情

资源推荐

PRO ex12

d=NCREAD('orca2.nc')

PSOPEN, YPLOTS=2, /PORTRAIT

CS, SCALE=1

MAP

LEVS, MIN=-2, MAX=30, STEP=2

CON, FIELD=d.sst, X=d.longitude, Y=d.latitude, CB_TITLE='Temperature (Celsius)'

POS, YPOS=2

MAP, LAND=12, OCEAN=8, /DRAW

GPLOT, X=d.longitude, Y=d.latitude, SYM=1, /NOLINES, SIZE=10

AXES

PSCLOSE

END

The top plot shows the Orca2 grid point locations with two northern poles and one southern pole. The bottom

plot shows the sea surface temperature on this grid being contoured and plotted in the usual manner. Further

discussion of contouring irregular grids can be found in the guide routines discussion section after these

example plots.

PRO ex13

PSOPEN, SPACE2=800

LEVS, MIN=-32, MAX=32, STEP=4

CS, SCALE=1

GSET, XMIN=-90, XMAX=90, YMIN=1000, YMAX=0

CON, F=SF('gdata.nc', 'u', MEAN='lon'), TITLE='Jan 1987 - Zonal Mean Zonal Wind', $

NEGATIVE_STYLE=2, ZERO_THICK=200,CB_TITLE='ms!E-1!N'

AXES, XSTEP=30, XTITLE='Latitude', YSTEP=-250, YTITLE='Pressure (mb)'

PSCLOSE

END

The plotting region was set using GSET before contouring the data with negative contours dashed and a

thicker zero contour line. The data was averaged over the longitudes with SF('gdata.nc', 'u', MEAN='lon')

to get the zonal mean zonal wind.

The plot to colour bar spacing was increased to 800 pixels with PSOPEN, SPACE2=800 so that the colour

bar didn’t overwrite the x-axis title of ‘Latitude’.

A separate call was made to AXES to label them with the specified axis labels. Note that a YSTEP= -250

was used as this is a pressure axis with YMIN=1000 and YMAX=0.

To zonally average the zonal wind we could also have used

d=NCREAD('gdata.nc')

TOTAL(d.u,1)/N_ELEMENTS(d.lon)

Thhis sums the array over the first dimension of longitude and then divides by the number of longitudes to get

the zonal mean zonal wind.

PRO ex14

d=NCREAD('gdata.nc')

zonalu=SF('gdata.nc', 'u', MEAN='lon')

PSOPEN

CS, COLS=[464]

LEVS, MANUAL=[-200, 0], /EXACT

GSET, XMIN=-90, XMAX=90, YMIN=1000, YMAX=0.3, /YLOG

CON, FIELD=zonalu, TITLE='Jan 1987 - Zonal Mean Zonal Wind', $

CB_TITLE='ms!E-1!N', /NOLINES, /NOCOLBAR

LEVS, MIN=-100, MAX=35, STEP=5

CON, FIELD=zonalu, NEGATIVE_STYLE=2, ZERO_THICK=200, /NOFILL

AXES, XSTEP=30, XTITLE='Latitude', YVALS=['1000', '300', '100', '30', '10', '3', '1', '0.3'], $

YTITLE='Pressure (mb)', /NORIGHT

ylabels=[0, 10, 20, 30, 40, 50, 56]

yvals=PCON(ylabels, /TO_MB)

AXES, /ONLYRIGHT, YVALS=yvals, YLABELS=ylabels, YTITLE='Height (km)'

PSCLOSE

END

Here we made a logarithmic y-axis by adding the /YLOG keyword to GSET.

In this plot we used a single shade of grey for the negative values. A second call to CON was made to overlay

the contour levels. Shading with grey is a common way of avoiding paying journal colour page charges.

The right y-axis was labelled by height rather than pressure by using the PCON function to convert pressure

to height. Be careful to not do the same for troposphere plots as the pressure height relation is only reliable

for the stratosphere. The /NORIGHT and /ONLYRIGHT keywords to the AXES routine were used to plot

the right y-axis in the correct way.

PRO ex15

d=NCREAD('stipple.nc')

PSOPEN

CS, SCALE=11, NCOLS=14

MAP, /ROBINSON

LEVS, MIN=-30, MAX=30, STEP=5

CON, FIELD=d.climatology, X=d.longitude, Y=d.latitude, CB_TITLE='Temperature [Celcius]'

PFILL, FIELD=d.climatology-d.jan1963, X=d.longitude, Y=d.latitude, MIN=1, MAX=20

PSCLOSE

END

In example 15 the PFILL procedure is used to show and anomaly field on top of a contour plot. The contour

plot is that of the 1957-2002 January climatology. The stippled field is that where the January 1963

temperature is one degree or less below this climatology.

Stippling is commonly used to indicate anomalies and deviations.

[PFILL has replaced the STIPPLE routine in the guide 4 routine library.]

PRO ex16

PSOPEN, SPACE2=800

f=SF('gdata.nc', 'u', MEAN='lon')

LEVS, MIN=-32, MAX=32, STEP=4

CS, SCALE=1

GSET, XMIN=-90, XMAX=90, YMIN=1000, YMAX=0

CON, F=f, TITLE='Jan 1987 - Zonal Mean Zonal Wind', $

NEGATIVE_STYLE=2, ZERO_THICK=200,CB_TITLE='ms!E-1!N', /NOFILL

PFILL, FIELD=f.data, X=f.lat, Y=f.p, MIN=28, MAX=100, STYLE=6

AXES, XSTEP=30, XTITLE='Latitude', YSTEP=-250, YTITLE='Pressure (mb)'

PSCLOSE

END

In example 16 a cross hatch pattern was used to indicate the areas where the zonal wind was above 28m/s.

剩余134页未读，继续阅读

#!/bin/bash test "$BASH_SOURCE" = "" && echo "This script can be sourced only from bash" && return SCRIPT_SOURCE=$BASH_SOURCE test "$SCRIPT_SOURCE" = "$0" && echo "Script is being run, should be sourced" && exit 1 SCRIPT_DIR=$(dirname "${SCRIPT_SOURCE}") ZEEKR_BASE=$(cd "${SCRIPT_DIR}"; pwd) ZEEKR_TARGET=$ZEEKR_BASE/target/ ZEEKR_PROTO_IDL_PATH=$ZEEKR_TARGET/usr/idl ZEEKR_BIN_PATH=$ZEEKR_TARGET/aarch64/usr/bin ZEEKR_LIB_PATH=$ZEEKR_TARGET/aarch64/usr/lib ZEEKR_INCLUDE_PATH=$ZEEKR_TARGET/usr/include ZEEKR_CMAKE_INSTALL_PATH=$ZEEKR_TARGET/aarch64/usr/cmake ZEEKR_PROTOC_EXEC=$ZEEKR_BASE/host/protoc/protoc ZEEKR_PREBUILD=$ZEEKR_TARGET/aarch64/prebuild MAKEFLAGS=-I$ZEEKR_TARGET/usr/include PATH=$ZEEKR_TARGET/aarch64/usr/bin:$ZEEKR_TARGET/aarch64/usr/lib:$PATH export ZEEKR_TARGET ZEEKR_PROTO_IDL_PATH ZEEKR_BIN_PATH ZEEKR_LIB_PATH ZEEKR_INCLUDE_PATH ZEEKR_CMAKE_INSTALL_PATH ZEEKR_PREBUILD ZEEKR_PROTOC_EXEC MAKEFLAGS PATH echo ZEEKR_PROTO_IDL_PATH=$ZEEKR_PROTO_IDL_PATH echo ZEEKR_BIN_PATH=$ZEEKR_BIN_PATH echo ZEEKR_LIB_PATH=$ZEEKR_LIB_PATH echo ZEEKR_INCLUDE_PATH=$ZEEKR_INCLUDE_PATH echo ZEEKR_CMAKE_INSTALL_PATH=$ZEEKR_CMAKE_INSTALL_PATH echo ZEEKR_PROTOC_EXEC=$ZEEKR_PROTOC_EXEC echo ZEEKR_PREBUILD=$ZEEKR_PREBUILD #Do not edit past this line

- `ZEEKR_PROTO_IDL_PATH`：IDL文件的路径。 - `ZEEKR_BIN_PATH`：二进制文件的路径。 - `ZEEKR_LIB_PATH`：库文件的路径。 - `ZEEKR_INCLUDE_PATH`：头文件的路径。 - `ZEEKR_CMAKE_INSTALL_PATH`：CMake安装路径。...

class RwsQmiClientFactory { public: static RwsQmiClientFactory &getInstance(); ~RwsQmiClientFactory(); template<typename T> std::shared_ptr<QmiClient> getRwsQmiClient(int phoneId, IdlServiceObjectType type) { LOG(INFO, FUNCTION); std::lock_guard<std::mutex> lock(qmiClientFactoryMutex_); std::shared_ptr<QmiClient> rwsQmiClient = getQmiClient2Map(phoneId, type); if (rwsQmiClient != nullptr){ LOG(INFO, "the qmiClient has exist in map: ", phoneId); return rwsQmiClient; } auto iter = s_service_idl_map.find(type); auto verIter = s_service_version_map.find(type); if ((iter == s_service_idl_map.end()) || (verIter == s_service_version_map.end())) { LOG(INFO, "the qmiClient type not exist in map: ", phoneId); return nullptr; } qmi_idl_service_object_type idlServiceObject = iter->second(verIter->second.idl_maj_version, verIter->second.idl_min_version, verIter->second.library_version); auto rwsClient = std::make_shared<T>(phoneId); telux::common::Status status = rwsClient->init(idlServiceObject); if (status != telux::common::Status::SUCCESS){ LOG(ERROR, "rwsQmiClient init failed.: ", phoneId); return nullptr; } rwsQmiClient = std::dynamic_pointer_cast<QmiClient>(rwsClient); addQmiClient2Map(phoneId, type, rwsQmiClient); return rwsQmiClient; }} 解释下这段代码，并说明下这段代码中的函数getRwsQmiClient在什么时候会被编译到和执行

函数的作用是检查是否已经存在具有给定电话ID和IDL服务类型的QmiClient实例。如果实例已经存在，则返回该实例的共享指针。否则，使用给定的电话ID和IDL服务类型创建一个新的RwsQmiClient实例，并将其添加到映射中。...

解释以下代码pro envi_rotate_define_buttons , btninfo compile_opt idl2 ENVI_DEFINE_MENU_BUTTON,btnInfo,$ Value = 'rotate',$ UVALUE = ' ',$ Ref_Value = 'Window',$ Position = 'after',$ Event_Pro = 'envi_revolve' end Pro envi_rotate Compile_opt idl2 ENVI,/restore_base_save_files ENVI_BATCH_INIT ENVI_OPEN_FILE, fileIn, R_FID=fid ENVI_FILE_QUERY, fid,dims=dims,ns=ns,nl=nl,nb=nb pos=indgen(nb) base = widget_auto_base(title='setting') sb=widget_base(base, /row, /frame) wp=widget_param(sb, prompt='ROT_TYPE',dt=3,uvalue='ROT_TYPE',/auto) sb=widget_base(base, /row, /frame) wf=widget_outfm(sb, uvalue='outf', PROMPT='Enter output File',/auto) result = auto_wid_mng(base) if (result.accept eq 0) then return ROT_TYPE=(result.ROT_TYPE) outfile=result.outf.name envi_doit,'ROTATE_DOIT',fid=fid,dims=dims,pos=pos,$ r_fid=r_fid,out_name=outfile,$ ROT_TYPE=ROT_TYPE end

这段代码是一段IDL程序，用于实现ENVI软件中的图像旋转功能。其中，pro envi_rotate_define_buttons定义了旋转功能的按钮信息；pro envi_rotate定义了旋转功能的具体实现过程，包括打开文件、获取文件信息、设置旋转...

粉丝: 0
资源: 7

上传资源快速赚钱

我的内容管理展开

我的资源快来上传第一个资源

我的收益

登录查看自己的收益

我的积分登录查看自己的积分

我的C币登录后查看C币余额

我的收藏

我的下载

下载帮助

idl在气象学中的实用应用：工具与实例解析

IDL与C#的数据传递方法研究COM_IDL_CONNECT

nn.zip_IDL NN_IDL批量_IDL系统开发_idl系统_校正

idl_function.zip_IDL数学函数_IDl_idl convol_idl function_idl中functio

TAO_IDL 编译器

怎么使用file:///usr/local/bin/opendds_idl

高通，QMI消息机制中 qmi_idl_service_message_table_entry 怎么新增一个项

详细解释下c++ 中的这段代码的意思 ： std::function<qmi_idl_service_object_type (int32_t, int32_t, int32_t)>

IDL_IDLBridge

怎么在python中读取idl文件

输出helloworld的CORBA程序设计思路

为什么我的IDL显示READ_JPEG未定义，如何解决

ACF attributes in the IDL file need the /app_config switch : [implicit_handle]

使用IDL读取jpeg格式图片文件

可以直接使用FLIP么 为什么我的idl显示其未定义

error MIDL2150 : ACF attributes in the IDL file need the /app_config switch : [implicit_handle]

在IDL编写器中使用ENVI_SETUP_HEAD函数

举个python 读取idl文件的例子

最新资源

详细解释下c++ 中的这段代码的意思： std::function<qmi_idl_service_object_type (int32_t, int32_t, int32_t)>

可以直接使用FLIP么为什么我的idl显示其未定义