| 1 |
dgoldberg |
1.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 |
dgoldberg |
1.2 |
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 |
dgoldberg |
1.3 |
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 |
dgoldberg |
1.2 |
|
| 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 |
dgoldberg |
1.3 |
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 |
dgoldberg |
1.2 |
|
| 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 |
dgoldberg |
1.3 |
|
Parent Directory
| 
Revision Log
| 
Revision Graph