/[MITgcm]/MITgcm_contrib/dgoldberg/depth_control_no_nsa/README
ViewVC logotype

Contents of /MITgcm_contrib/dgoldberg/depth_control_no_nsa/README

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


Revision 1.3 - (show annotations) (download)
Fri Dec 8 14:38:55 2017 UTC (7 years, 7 months ago) by dgoldberg
Branch: MAIN
Changes since 1.2: +5 -3 lines
update README re: grdchk results

1 NOTES ON SETUP
2
3 1. cg2d_no_nsa.F should be exactly like cg2d.F -- but EXCH_S3D needed to be replaced by EXCH_XY_RL o/w the forward model crashed
4
5 2. cg2d.flow sets adjoint of cg2d_no_nsa to adcg2d_no_nsa -- this is the only reason cg2d_no_nsa.F even needs to exist
6
7 3. solve_for_pressure.F modified to call the above, only when ALLOW_DIFFERENTIATE_CG2D_MATRIX is defined and ALLOW_CG2D_NSA is undefined (in ECCO_CPPOPTIONS).
8
9 4. adcg2d_no_nsa finds adjoint of cg2d_b_ad as well as ac2d_ad etc. Needs common block from ad_taf_output.f to affect the latter. harcoded.
10
11 5. if line 312 of solve_for_pressure.F not included, mitgcmuv_ad will crash in forward pass (ac2d etc takes on different values -- don't know why)
12
13 6. line 212 of obcs_fields_load.F commented; if not, there are extensive recomputation errors in taf_ad.log i cannot fix.
14
15 7. shelfice_ad_diff.list in "main_branch_files" were modifief but do not seem to be required.
16
17 8. output below for depth control. with temp0 or salt0 control, grdchk is accurate with NONLIN_FRSURF and nonlinearfreesurface=4.
18
19 OUTPUT WITH NONLIN_FRSURF OFF:
20
21 grad-res -------------------------------
22 grad-res proc # i j k bi bj iobc fc ref fc + eps fc - eps
23 grad-res proc # i j k bi bj iobc adj grad fd grad 1 - fd/adj
24 grad-res -------------------------------
25 grad-res 0 1 4 13 1 1 1 1 -9.00753618943E+01 -9.00753618941E+01 -9.00753618945E+01
26 grad-res 0 1 1 10 0 1 1 1 2.07265124500E-08 2.07258210594E-08 3.33577897197E-05
27 grad-res -------------------------------
28 grad-res 0 2 1 14 1 1 1 1 -9.00753618943E+01 -9.00753614330E+01 -9.00753623556E+01
29 grad-res 0 2 2 11 0 1 1 1 4.61308022316E-05 4.61298220955E-05 2.12468893396E-05
30 grad-res -------------------------------
31 grad-res 0 3 2 14 1 1 1 1 -9.00753618943E+01 -9.00753618943E+01 -9.00753618942E+01
32 grad-res 0 3 3 12 0 1 1 1 -4.08823433387E-09 -4.08491018788E-09 8.13100646840E-04
33 grad-res -------------------------------
34 grad-res 0 4 3 14 1 1 1 1 -9.00753618943E+01 -9.00753618940E+01 -9.00753618946E+01
35 grad-res 0 4 4 13 0 1 1 1 3.04480641401E-08 3.04495983983E-08 -5.03893496338E-05
36 grad-res -------------------------------
37 grad-res 0 5 4 14 1 1 1 1 -9.00753618943E+01 -9.00753618940E+01 -9.00753618945E+01
38 grad-res 0 5 5 14 0 1 1 1 2.48672264310E-08 2.48682852089E-08 -4.25772386625E-05
39 grad-res -------------------------------
40 grad-res 0 6 5 14 1 1 1 1 -9.00753618943E+01 -9.00753618918E+01 -9.00753618967E+01
41 grad-res 0 6 6 15 0 1 1 1 2.44778868273E-07 2.44781261927E-07 -9.77883928943E-06
42
43 OUTPUT WITH NONLIN_FRSURF ON:
44
45 grad-res -------------------------------
46 grad-res proc # i j k bi bj iobc fc ref fc + eps fc - eps
47 grad-res proc # i j k bi bj iobc adj grad fd grad 1 - fd/adj
48 grad-res -------------------------------
49 grad-res 0 1 4 13 1 1 1 1 -9.00753618943E+01 -9.00753618941E+01 -9.00753618945E+01
50 grad-res 0 1 1 10 0 1 1 1 4.88519940983E-08 2.07258210594E-08 5.75742578334E-01
51 grad-res -------------------------------
52 grad-res 0 2 1 14 1 1 1 1 -9.00753618943E+01 -9.00753614330E+01 -9.00753623556E+01
53 grad-res 0 2 2 11 0 1 1 1 4.61308022316E-05 4.61298220955E-05 2.12468893396E-05
54 grad-res -------------------------------
55 grad-res 0 3 2 14 1 1 1 1 -9.00753618943E+01 -9.00753618943E+01 -9.00753618942E+01
56 grad-res 0 3 3 12 0 1 1 1 -1.02457536117E-08 -4.08491018788E-09 6.01307005545E-01
57 grad-res -------------------------------
58 grad-res 0 4 3 14 1 1 1 1 -9.00753618943E+01 -9.00753618940E+01 -9.00753618946E+01
59 grad-res 0 4 4 13 0 1 1 1 1.13033514034E-07 3.04495983983E-08 7.30614423001E-01
60 grad-res -------------------------------
61 grad-res 0 5 4 14 1 1 1 1 -9.00753618943E+01 -9.00753618940E+01 -9.00753618945E+01
62 grad-res 0 5 5 14 0 1 1 1 6.65438158395E-08 2.48682852089E-08 6.26287057706E-01
63 grad-res -------------------------------
64 grad-res 0 6 5 14 1 1 1 1 -9.00753618943E+01 -9.00753618918E+01 -9.00753618967E+01
65 grad-res 0 6 6 15 0 1 1 1 2.52987938174E-07 2.44781261927E-07 3.24390020583E-02
66
67

  ViewVC Help
Powered by ViewVC 1.1.22