-
Notifications
You must be signed in to change notification settings - Fork 0
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Lorenz96でのデータ同化実装とv1のリファクタリング #33
Conversation
@@ -311,7 +310,7 @@ program lorenz63 | |||
write(6,*) ' TRUTH = ', x_true(it), y_true(it), z_true(it) | |||
write(6,*) ' PREDICT = ', x_sim(it), y_sim(it), z_sim(it) | |||
write(6,*) ' OBSERVE = ', x_obs(it/obs_interval), y_obs(it/obs_interval), z_obs(it/obs_interval) | |||
write(6,*) ' ANALYSIS = ', x_anl(it), y_anl(it), z_anl(it) | |||
write(6,*) ' ANALYSIS (BEFORE) = ', x_anl(it), y_anl(it), z_anl(it) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
logで確認した方が早いので、効率性を意識した出力形式にした。
@@ -162,7 +162,7 @@ def lorenz_rmse_draw(self, rmse_sim:list, rmse_obs:list, rmse_da:list, | |||
ax1.set_ylabel('RMSE') | |||
ax1.set_xlim(0, 1) | |||
ax1.set_ylim(0, 10) | |||
ax1.set_title('Lorenz(1963) RMSE infla. ***', loc='left') | |||
ax1.set_title('Lorenz(1963) RMSE infla. 0.0d0', loc='left') |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
lorenz63では、共分散膨張スキームを活用したくてもfiler divergenceは生じなかった。
@@ -21,16 +23,21 @@ oneday=0.2d0 | |||
spinup_period=365 | |||
normal_period=40 | |||
|
|||
da_method='KF' | |||
da_method='EnKF' |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
EnKFの導入(sh_ver)
#24 #31 と関連。
新たに実装した点(重要度が高い順に)
※行列計算には、
xt_vec(nx,1), yt_vec(ny,1)
のような一時的な行列フォーマットを使用した方が良いと分かった。