/[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.3 by dgoldberg, Fri Dec 8 14:38:55 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 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    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22