| 1 |
Starting a configuration for OpenAD |
| 2 |
19-Aug-2005, heimbach@mit.edu, utke@mcs.anl.gov, cnh@mit.edu |
| 3 |
############################################################ |
| 4 |
|
| 5 |
This experiment is derived from global_ocean.90x40x15, |
| 6 |
but excludes packages gmredi, kpp. |
| 7 |
|
| 8 |
STEP 1: |
| 9 |
###### |
| 10 |
To bypass incomplete canonicalizer, convert COMMON blocks to MODULES. |
| 11 |
Successfully completed. |
| 12 |
The built process needed to be modified, and some routines |
| 13 |
needed changes. Most changes were commited to default routines, |
| 14 |
the remaining changes are kept in code/ for now. |
| 15 |
|
| 16 |
To build: |
| 17 |
-------- |
| 18 |
o short version: |
| 19 |
--------------- |
| 20 |
make makefile ; make depend ; make cb2m ; make makefile ; make small_f ; make allmods ; make |
| 21 |
|
| 22 |
o long version: |
| 23 |
-------------- |
| 24 |
# (chdir to built/; assume we are in built/) |
| 25 |
# |
| 26 |
# generate makefile using ifort and modified genmake2 in ../code/ |
| 27 |
../code/genmake2 -of ../../../tools/build_options/linux_ia32_ifort -mods ../code |
| 28 |
# |
| 29 |
# make dependencies |
| 30 |
make depend |
| 31 |
# |
| 32 |
# invoke script to convert COMMON block headers to MODULE headers |
| 33 |
# converts FILE.h to FILE_mod.h which uses newi module FILE_mod.F90 |
| 34 |
make cb2m |
| 35 |
# |
| 36 |
# re-generate makefile which takes into account newly created files |
| 37 |
# FILE_mod.h, FILE_mod.F90 |
| 38 |
make makefile |
| 39 |
# |
| 40 |
# make .f, .f90 |
| 41 |
make small_f |
| 42 |
# |
| 43 |
# first compile all module files .f90 |
| 44 |
make allmods |
| 45 |
# |
| 46 |
# compile everything else |
| 47 |
make |
| 48 |
|
| 49 |
To clean: |
| 50 |
-------- |
| 51 |
# Since soft links get overwritten, for now do: |
| 52 |
make CLEAN ; rm *.F *.F90 *.h |
| 53 |
|
| 54 |
to run: |
| 55 |
------ |
| 56 |
# paremeter files are in input/ |
| 57 |
ln -s ../input/* . |
| 58 |
# initial and forcing fields are elsewhere; |
| 59 |
# get these from verif. exp. global_ocean.90x40x15/input/*.bin |
| 60 |
ln -s ../../global_ocean.90x40x15/input/*.bin . |
| 61 |
./mitgcmuv >! output.txt |
| 62 |
|
| 63 |
STEP 2: |
| 64 |
###### |
| 65 |
Generate code for AD-related routines. |
| 66 |
|
| 67 |
Similar to step 1, but look in code_ad/ instead of code/ |
| 68 |
|
| 69 |
To build: |
| 70 |
-------- |
| 71 |
o short version: |
| 72 |
--------------- |
| 73 |
|
| 74 |
../code_ad/genmake2 -of ../../../tools/build_options/linux_ia32_ifort -adof ../../../tools/adjoint_options/adjoint_f95 -mods ../code_ad |
| 75 |
|
| 76 |
make adall |
| 77 |
|
| 78 |
STEP 3: |
| 79 |
###### |
| 80 |
Compare reference outputs. |
| 81 |
|
| 82 |
# uncompress outputs in results_... directories: |
| 83 |
gunzip ../results_shallow_*/output.txt.gz |
| 84 |
|
| 85 |
# extract gradients from TAF-generated adjoint: |
| 86 |
grep grad-res ../results_shallow_taf/output.txt |
| 87 |
|
| 88 |
# extract gradients from OpenAD-generated adjoint: |
| 89 |
grep ' salt/theta ' ../results_shallow_openad/output.txt |
| 90 |
|
| 91 |
The two positions tested via gradient checks (results_shallow_taf) are |
| 92 |
(i,j,k) = (43,2,1) and (44,2,1) |
| 93 |
Find corresponding locations in ' salt/theta ' output |
| 94 |
(which outputs whole theta, salt gradient) |
Parent Directory
| 
Revision Log
| 
Revision Graph