/[MITgcm]/MITgcm/verification/carbon/DESCRIPTION
ViewVC logotype

Contents of /MITgcm/verification/carbon/DESCRIPTION

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (show annotations) (download)
Fri Jul 13 13:19:00 2001 UTC (22 years, 8 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint46b_post, checkpoint57t_post, checkpoint48f_post, checkpoint46k_post, checkpoint51k_post, checkpoint53f_post, checkpoint47j_post, checkpoint54a_pre, checkpoint55c_post, checkpoint53b_pre, checkpoint48d_pre, checkpoint51l_post, checkpoint51j_post, branch-exfmods-tag, checkpoint47e_post, checkpoint57m_post, checkpoint44h_pre, checkpoint47i_post, checkpoint52l_pre, checkpoint48i_post, checkpoint52e_pre, release1_p12, release1_p13, release1_p10, checkpoint57v_post, release1_p16, checkpoint57g_pre, checkpoint52j_post, release1_p15, release1_p11, checkpoint47f_post, checkpoint48d_post, checkpoint51o_pre, checkpoint57f_post, checkpoint46j_post, checkpoint47c_post, checkpoint50e_post, checkpoint52e_post, checkpoint50c_post, checkpoint57s_post, checkpoint51n_pre, checkpoint47d_post, checkpoint44e_post, checkpoint40pre3, checkpoint57j_post, checkpoint58b_post, checkpoint44f_pre, checkpoint47a_post, checkpoint46c_post, checkpoint57b_post, checkpoint46f_post, checkpoint46l_pre, checkpoint53c_post, checkpoint53d_post, checkpoint57f_pre, checkpoint46a_post, checkpoint48a_post, checkpoint55d_pre, checkpoint46n_post, checkpoint51f_pre, release1_beta1, checkpoint46d_pre, checkpoint57g_post, checkpoint48e_post, checkpoint46e_post, checkpoint57a_post, checkpoint48h_post, checkpoint55j_post, checkpoint56b_post, checkpoint50c_pre, release1-branch_tutorials, checkpoint57h_pre, checkpoint57y_post, release1_p14, checkpoint44g_post, branchpoint-genmake2, checkpoint46h_pre, checkpoint44h_post, checkpoint46l_post, chkpt44c_post, checkpoint57x_post, checkpoint52j_pre, checkpoint54a_post, checkpoint46e_pre, checkpoint43a-release1mods, branch-netcdf, checkpoint50d_pre, checkpoint55h_post, checkpoint45d_post, checkpoint51r_post, checkpoint52b_pre, checkpoint52n_post, checkpoint54b_post, checkpoint46j_pre, checkpoint45b_post, checkpoint46b_pre, checkpoint51i_post, checkpoint57e_post, checkpoint54d_post, checkpoint47h_post, checkpoint48c_post, chkpt44a_pre, release1-branch-end, checkpoint56c_post, checkpoint54e_post, release1_final_v1, checkpoint55b_post, checkpoint51e_post, checkpoint51b_post, checkpoint46, checkpoint57h_post, checkpoint51l_pre, checkpoint52m_post, checkpoint51c_post, checkpoint57y_pre, checkpoint55, checkpoint53a_post, checkpoint41, checkpoint55a_post, checkpoint57c_pre, checkpoint44, checkpoint45, checkpoint48, checkpoint49, checkpoint44f_post, checkpoint47b_post, checkpoint40pre2, checkpoint40pre5, checkpoint53b_post, checkpoint57o_post, checkpoint40pre6, checkpoint55g_post, checkpoint57r_post, checkpoint51o_post, checkpoint40pre8, checkpoint48g_post, release1_p17, checkpoint57k_post, checkpoint57d_post, release1_b1, checkpoint44b_post, checkpoint55f_post, checkpoint57i_post, checkpoint51q_post, checkpoint52l_post, checkpoint52k_post, chkpt44d_post, checkpoint58, checkpoint42, checkpoint57a_pre, checkpoint54, checkpoint57, checkpoint56, checkpoint51, checkpoint50, checkpoint53, checkpoint52, release1_p8, release1_p9, checkpoint50d_post, checkpoint52d_post, checkpoint46m_post, checkpoint46g_pre, release1_p2, release1_p3, release1_p4, checkpoint51b_pre, release1_p6, checkpoint52a_post, checkpoint57h_done, checkpoint47g_post, checkpoint52b_post, checkpoint53g_post, chkpt44a_post, checkpoint52f_post, checkpoint44b_pre, checkpoint57n_post, checkpoint52c_post, release1_p1, checkpoint58d_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint40pre4, checkpoint51h_pre, checkpoint58a_post, checkpoint46a_pre, checkpoint50g_post, checkpoint45c_post, checkpoint57q_post, release1_p5, checkpoint44e_pre, checkpoint51g_post, ecco_c52_e35, chkpt44c_pre, checkpoint40pre9, release1_p7, checkpoint57z_post, checkpoint54f_post, checkpoint51f_post, checkpoint46d_post, checkpoint48b_post, checkpoint50b_post, eckpoint57e_pre, checkpoint46g_post, checkpoint58c_post, release1_p13_pre, release1_p12_pre, checkpoint45a_post, checkpoint57c_post, checkpoint50f_post, checkpoint50a_post, checkpoint46c_pre, checkpoint50f_pre, checkpoint52a_pre, checkpoint43, checkpoint47d_pre, checkpoint51d_post, checkpoint40, checkpoint48c_pre, release1-branch_branchpoint, checkpoint46i_post, checkpoint51m_post, checkpoint51t_post, checkpoint53d_pre, checkpoint47, checkpoint55e_post, checkpoint40pre7, checkpoint54c_post, checkpoint50h_post, checkpoint52i_post, checkpoint51a_post, checkpoint46h_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51p_post, checkpoint51n_post, release1_chkpt44d_post, checkpoint55i_post, checkpoint51i_pre, checkpoint57l_post, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, checkpoint50b_pre, checkpoint56a_post, checkpoint51s_post, checkpoint55d_post
Branch point for: netcdf-sm0, branch-genmake2, release1_coupled, branch-nonh, tg2-branch, release1_final, checkpoint51n_branch, release1-branch, release1, release1_50yr, branch-exfmods-curt
Adding carbon sequestration verification experiment
(forward and adjoint integration over 3 timesteps with gradient checks).

1 Setup:
2 -----
3
4 The experiment is in verification/carbon/.
5 The code/ directory has the extra and modified code to
6 carry a passive tracer (variables: tr1, gTr1, gTr1Nm1).
7 Additional mod's are in pkg/ctrl/
8
9 The tracer fields is initialised at the surface (ini_tr1.F)
10 and then carried around by the flow field.
11 The cost funciton is cost_tracer.F
12 It basically accumulates all tracers reaching the surface,
13
14 The controls so far are
15 initial temp./salinity and surface fluxes,
16 also initial tracer (tr1) distribution
17 (see adjoint/makefile).
18
19 I added the common block /adtr1/ to addummy_in_stepping
20 to be able to monitor adtr1 in addition to adtheta, adsalt, ...
21
22 The model picks up from a 400 year spinup.
23 In a test phase it performs a 1 year forward/adjoint run
24 (timestep is 1 day).
25 dumpFreq is set to 30 days, yielding adjoint fields of
26 theta, salt, uvel, vvel, tr1
27 gtnm1,gsnm1, gunm1, gvnm1, gtr1nm1
28
29 If you plot the adjoint fields, e.g. at initial timestep index 144000
30 you will see that:
31 ADJtheta.0000144000.001.001.data
32 ADJsalt
33 ADJuvel
34 ADJvvel
35 have a checkerboard structure, whereas
36 ADJtr1.0000144000.001.001.data
37 seems to look ok.
38
39 I had the checkerboard structure also for a simpler
40 adjoint test case (see verification/vero/).
41 This is basically a box setup,
42 cost function is final temp. at one grid point.
43 The checkerboard structure seems to disappear, if you
44 set in "data"
45 tAlpha = 0.
46 and take a linear eqn. of state
47 eosType='LINEAR'

  ViewVC Help
Powered by ViewVC 1.1.22