NCL提取WRFout的经纬度信息

杨雨轩

我想用NCL提取wrfout中变量垂直速度对应的经纬度，将其存为二维数组，并保存到TXT中，不知道怎么写，请各位知道

1649518749

同问。就是怎么提取资料。有时候要做月平均，年平均。

weiguang1233

表示密切关注回复

杨雨轩

在同学的帮助下，提取出来了，但并不是用垂直速度，因为垂直速度维度中并没有经纬度信息（它是格点），而是采用XLAT和XLONG，并将经纬度分别以二维数组的形式存入两个txt文件，NCL脚本如下：
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/wrf/WRFUserARW.ncl"

begin

  a=addfile("/yhhome/home/yangyuxuan/trunk/PODEn4DVar/data/rdsample/wrfout331/wrfout_d01_2010-07-08_00:00:00.nc","r")
  xlat=a->XLAT(0,:,:)
  xlong=a->XLONG(0,:,:)
  asciiwrite("XLAT_M.txt",xlat)
  asciiwrite("XLONG_M.txt",xlong)
  
end

yeah...

杨雨轩 发表于 2016-9-16 21:14
在同学的帮助下，提取出来了，但并不是用垂直速度，因为垂直速度维度中并没有经纬度信息（它是格点），而是 ...

您好 我想请教一下
我们读取一般的文件的时候  经度为一维  纬度为一维
为什么在读取wrfout的时候 要用xlat=a->XLAT(0,:,:)
  xlong=a->XLONG(0,:,:)来读呢？我不理解这个  

杨雨轩

yeah... 发表于 2017-7-3 21:06
您好 我想请教一下
我们读取一般的文件的时候  经度为一维  纬度为一维
为什么在读取wrfout的时候 要用 ...

wrfout数据集里的经纬度信息是三维的，（time,lat,lon）,这是由它的投影方式决定的，但虽然它是三维的，但每个时刻的经纬度实质上是一样的，所以经纬度其实是二维的，在NCL里的好多函数所需的经纬度是一维的，一般是在二维的基础上进行处理。

Thaliajin

{:eb513:}

mona

您好，我想请问，若提取经纬度范围的数据，xlat=a->XLAT(0,:,:)，该怎么写，时间就一个时刻的，所以时间是不变的。

mona

杨雨轩 发表于 2016-9-16 21:14
在同学的帮助下，提取出来了，但并不是用垂直速度，因为垂直速度维度中并没有经纬度信息（它是格点），而是 ...

xlat=a->XLAT(0,:,:)，如何设定经纬度范围，比如我想提取经度为72°-136°，纬度为2°-54°，那我该怎么写？求助

mona

杨雨轩 发表于 2016-9-16 21:14
在同学的帮助下，提取出来了，但并不是用垂直速度，因为垂直速度维度中并没有经纬度信息（它是格点），而是 ...

xlat=a->XLAT(0,:,:)，如何设定经纬度范围，比如我想提取经度为72°-136°，纬度为2°-54°，那我该怎么写？求助