一月分のデータだけが含まれたファイルを使ってのシミュレーションができない…

まとめ：NEMOで一ヶ月分（例えば一月分）のoffline dynamicsを使ってオフラインでPISCESをシミュレートしたいんだがうまくいかなかった。データが一ヶ月分しか含まれていないということの認識はできているようだがsowaflcdやsoicecovなどの2次元空間のデータが正しく読み込まれていないようだ。とりあえず本日のところは未解決。

日平均のデータを使って一ヶ月ずつシミュレーションをするように言われたのでそのメモ。

time resolution = 1hour (3600sec)

24 hours in 1 day

24 time steps per day

24*31 = 744 time sptes

24*30 = 720 time steps

24*28 = 672 time steps (February)

1994年のラブラドール海をシミュレーションした物理モデルデータなので大気中二酸化炭素濃度は360ppmに設定。

オフラインデータの情報もnamelist内で編集。
   sn_tem  = 'Lab_T_94_01' ,   24           , 'votemper' ,  .true.    , .true. ,   'yearly'  , ''       , ''

そしてシミュレーションスタート！っとすぐにストップ！

ocean.outputを見てみると原因がよくわからんエラーが。。。

===>>> : E R R O R
         ===========

           iom_get_123d, file: ./Lab_T_94_01.nc, var: votemper
 start and count too big regarding to the size of the data,
 (istart(4) + icnt(4) - 1) =   365                                                                                                                            
オフラインのデータは一ヶ月分あるからそれを使って一ヶ月間のシミュレーションはできるはずなんだがな〜。

少し原因を調べてみるとわかったことがひとつ。
namelistでオフラインデータを定義するときの'yearly'という情報があることに気づいた。どうやらこれはオフラインデータは一年間分のデータです、ということを定義しているみたいだ。だからこれを'monthly'というオプションに変えてシミュレート再スタート。

すると今度のエラーはこちら：
===>>> : W A R N I N G
         ===============                                                        
 previous year/month/week/day file: ./Lab_T_94_01_m12 not present -> back to current year/month/week/day

 ===>>> : E R R O R
         ===========
 
                     iom_open ~~~                                                File ./Lab_T_94_01_m01.nc* not found
fld_init : time-interpolation for votemper read previous record =    1 at time =   -0.50 days

Lab_T_94_01.ncというファイル名を定義したのに語尾が_m01.ncになったファイルを勝手に読み込もうとしているみたい。多分月単位のデータはm01を一月、m02を二月、というふうに定義するように設定されているみたいだ。

だったらファイル名をLab_T_94_m01.ncに変更して読み込めるようにしようじゃないかとやってみたら次のエラーはこちら…

===>>> : E R R O R

         ===========

           iom_get_123d, file: ./Lab_T_94_m01.nc, var: somixhgt

 The file has only 2 spatial dimension

 it is impossible to read a 3D array from this file...

 ===>>> : E R R O R

         ===========

           iom_get_123d, file: ./Lab_T_94_m01.nc, var: somixhgt

 start and count too big regarding to the size of the data,

 (istart(3) + icnt(3) - 1) =    46                                                                                                                                                                                                                

  is larger than idimsz(3) =    31                                                                                                                                                                                                                

fld_init : time-interpolation for somixhgt read previous record =    1 at time =   -0.50 days

最初の次元に関するエラー文は理解できない（ファイルの空間次元数は２だが３は必要ないのでこのエラーはおかしい）。次のエラー文の４６というのはZ-coordinateのことをさしてるのでまたどこかでデータの読み込みが正しくされてないみたい。

明日教授に聞いてみよう。