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

Diff of /MITgcm_contrib/dgoldberg/depth_control_no_nsa/README

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

revision 1.1 by dgoldberg, Thu Dec 7 23:21:12 2017 UTC revision 1.2 by dgoldberg, Thu Dec 7 23:30:59 2017 UTC
# Line 2  NOTES ON SETUP Line 2  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  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    
 2. solve_for_pressure.F modified to call the above, when  
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 (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" may be required.
16    
17    OUTPUT WITH NONLIN_FRSURF OFF:
18    
19    grad-res -------------------------------
20     grad-res  proc    #    i    j    k   bi   bj iobc       fc ref            fc + eps           fc - eps
21     grad-res  proc    #    i    j    k   bi   bj iobc      adj grad            fd grad          1 - fd/adj
22    grad-res -------------------------------
23     grad-res     0    1    4   13    1    1    1    1  -9.00753618943E+01 -9.00753618941E+01 -9.00753618945E+01
24     grad-res     0    1    1   10    0    1    1    1   2.07265124500E-08  2.07258210594E-08  3.33577897197E-05
25    grad-res -------------------------------
26     grad-res     0    2    1   14    1    1    1    1  -9.00753618943E+01 -9.00753614330E+01 -9.00753623556E+01
27     grad-res     0    2    2   11    0    1    1    1   4.61308022316E-05  4.61298220955E-05  2.12468893396E-05
28    grad-res -------------------------------
29     grad-res     0    3    2   14    1    1    1    1  -9.00753618943E+01 -9.00753618943E+01 -9.00753618942E+01
30     grad-res     0    3    3   12    0    1    1    1  -4.08823433387E-09 -4.08491018788E-09  8.13100646840E-04
31    grad-res -------------------------------
32     grad-res     0    4    3   14    1    1    1    1  -9.00753618943E+01 -9.00753618940E+01 -9.00753618946E+01
33     grad-res     0    4    4   13    0    1    1    1   3.04480641401E-08  3.04495983983E-08 -5.03893496338E-05
34    grad-res -------------------------------
35     grad-res     0    5    4   14    1    1    1    1  -9.00753618943E+01 -9.00753618940E+01 -9.00753618945E+01
36     grad-res     0    5    5   14    0    1    1    1   2.48672264310E-08  2.48682852089E-08 -4.25772386625E-05
37    grad-res -------------------------------
38     grad-res     0    6    5   14    1    1    1    1  -9.00753618943E+01 -9.00753618918E+01 -9.00753618967E+01
39     grad-res     0    6    6   15    0    1    1    1   2.44778868273E-07  2.44781261927E-07 -9.77883928943E-06
40    
41    OUTPUT WITH NONLIN_FRSURF ON:
42    
43    grad-res -------------------------------
44     grad-res  proc    #    i    j    k   bi   bj iobc       fc ref            fc + eps           fc - eps
45     grad-res  proc    #    i    j    k   bi   bj iobc      adj grad            fd grad          1 - fd/adj
46    grad-res -------------------------------
47     grad-res     0    1    4   13    1    1    1    1  -9.00753618943E+01 -9.00753618941E+01 -9.00753618945E+01
48     grad-res     0    1    1   10    0    1    1    1   4.88519940983E-08  2.07258210594E-08  5.75742578334E-01
49    grad-res -------------------------------
50     grad-res     0    2    1   14    1    1    1    1  -9.00753618943E+01 -9.00753614330E+01 -9.00753623556E+01
51     grad-res     0    2    2   11    0    1    1    1   4.61308022316E-05  4.61298220955E-05  2.12468893396E-05
52    grad-res -------------------------------
53     grad-res     0    3    2   14    1    1    1    1  -9.00753618943E+01 -9.00753618943E+01 -9.00753618942E+01
54     grad-res     0    3    3   12    0    1    1    1  -1.02457536117E-08 -4.08491018788E-09  6.01307005545E-01
55    grad-res -------------------------------
56     grad-res     0    4    3   14    1    1    1    1  -9.00753618943E+01 -9.00753618940E+01 -9.00753618946E+01
57     grad-res     0    4    4   13    0    1    1    1   1.13033514034E-07  3.04495983983E-08  7.30614423001E-01
58    grad-res -------------------------------
59     grad-res     0    5    4   14    1    1    1    1  -9.00753618943E+01 -9.00753618940E+01 -9.00753618945E+01
60     grad-res     0    5    5   14    0    1    1    1   6.65438158395E-08  2.48682852089E-08  6.26287057706E-01
61    grad-res -------------------------------
62     grad-res     0    6    5   14    1    1    1    1  -9.00753618943E+01 -9.00753618918E+01 -9.00753618967E+01
63     grad-res     0    6    6   15    0    1    1    1   2.52987938174E-07  2.44781261927E-07  3.24390020583E-02
64    
65    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22