問題タブ [netcdf4]

を使用して netcdf からデータをプロットしようとしていますBasemapが、緯度のインデックスが反転しているため、マップが上下逆になっていると思います。これを修正するにはどうすればよいですか? ありがとう！

netCDF4 ライブラリでデフォルトのチャンク設定を確認できますが、チャンク サイズを変更する方法がわかりません。

設定を手伝ってくれる人はいますか？

LandMask_NaNnetCDF ファイルの値を、ゼロより大きい場合に呼び出される NaN に変換したいと考えています。LandMaskただし、と numpy が NaN に変換するものとの間に型の不一致があるようです。以下のコードと情報を大歓迎します。

私が得るエラーは

しようとすると

また

NaN に変換しようとする前に、

取得する LandMask のタイプを確認するとき

<type 'netCDF4._netCDF4.Variable'>

変数の型を調べる方法がわかりません。

更新: nc 変数の詳細

次のコードを使用して、Fortran で既存の netCDF ファイルに新しい次元を追加しました -

ここで、このディメンションをゼロの値で埋められるようにしたいと考えています。このセットのカーディナリティは、私の場合のジオポテンシャルの高さの変数のカーディナリティと同じです。さらに、時間 (無制限)、緯度、経度、レベルの 3 つの次元があります。

fortran で netCDF API を調べたのですが、どの API を呼び出すべきかわかりません。次の API を使用すると、

ジオポテンシャルの高さの値を 0.0 で上書きしてしまいます (これは私の netCDF ファイルの唯一の変数です)。

これを行うにはどうすればよいですか？

コードは、この答えは機能します新しい次元を追加しますが、それは私がやりたいことをしません。

私の問題をさらに説明するために、次の定義があります。時間（無制限）、緯度、経度、レベルとして 4 つの次元と 1 つの可変ジオポテンシャルの高さです。私がやりたいことは、ディメンション dimz_Z を既存の可変ジオポテンシャルの高さに追加することです。

現在、上記のコードの答えは、netCDF ファイルが ncdump でどのように見えるかです。

しかし、ncks の出力から以下に示すように、dimz_Z 次元が欠落しています。

定義には表示されますが、スタンドアロンのディメンションとして表示され、hgt のディメンションの一部ではありません。

これはファイルの出力です

私が欲しいのは次のとおりです

私が友人から言われたのは、この次元を既存の変数に追加することはできないということです。代わりに、Fortran では、変数 hgt を 5 次元変数としてもう一度定義し、値をもう一度入力する必要があります。それは正しいですか、それとも代替手段はありますか? Fortran 77 netcdf API を使用しています。