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

Annotation of /MITgcm/verification/carbon/DESCRIPTION

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

Revision 1.1 - (hide annotations) (download)
Fri Jul 13 13:19:00 2001 UTC (22 years, 10 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 heimbach 1.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