1 |
SHELL = /bin/tcsh |
2 |
RM = rm -f |
3 |
LOC = $(PWD) |
4 |
TAMC = tamc |
5 |
TAF = ~fastopt/bin/taf |
6 |
DEST = . |
7 |
SOURCE_CODE = ../bin |
8 |
ADJOINT_SCRIPT = ../adjoint |
9 |
BASE = .. |
10 |
TAMCFLAGS = -reverse -i4 -r4 |
11 |
TAFFLAGS = -reverse -i4 -r4 |
12 |
PKG = $(BASE)/pkg |
13 |
|
14 |
DIFF_FLAGS = -toplevel the_main_loop \ |
15 |
-input ' xx_theta_dummy \ |
16 |
xx_salt_dummy \ |
17 |
xx_tr1_dummy \ |
18 |
xx_hflux_dummy \ |
19 |
xx_sflux_dummy \ |
20 |
xx_tauu_dummy \ |
21 |
xx_tauv_dummy \ |
22 |
xx_sss_dummy \ |
23 |
xx_sst_dummy \ |
24 |
xx_diffkr_dummy \ |
25 |
xx_kapgm_dummy' \ |
26 |
-output 'fc' |
27 |
|
28 |
TAMCFLAG = $(TAMCFLAGS) $(DIFF_FLAGS) -admark ad |
29 |
|
30 |
TAFFLAG = $(TAFFLAGS) $(DIFF_FLAGS) -admark ad |
31 |
|
32 |
|
33 |
# Source files for the MITgcmUV |
34 |
SRC_MODEL = \ |
35 |
$(PKG)/autodiff/read_write.flow \ |
36 |
$(PKG)/autodiff/read_write_fld.flow \ |
37 |
$(PKG)/autodiff/read_write_rec.flow \ |
38 |
$(PKG)/autodiff/diags.flow \ |
39 |
$(PKG)/autodiff/cg2d.flow \ |
40 |
$(PKG)/autodiff/debug.flow \ |
41 |
$(PKG)/autodiff/dummy_in_stepping.flow \ |
42 |
$(PKG)/autodiff/external_fields_load.flow \ |
43 |
$(PKG)/autodiff/eesupp.flow \ |
44 |
$(PKG)/autodiff/active_file.flow \ |
45 |
$(PKG)/autodiff/write_state.flow \ |
46 |
$(PKG)/autodiff/ecco_check_exp.flow \ |
47 |
$(PKG)/autodiff/exch_ad.flow \ |
48 |
$(PKG)/autodiff/exch_z_ad.flow \ |
49 |
$(PKG)/autodiff/exch_uv_ad.flow \ |
50 |
$(PKG)/autodiff/global_sum_ad.flow \ |
51 |
$(PKG)/autodiff/global_max_ad.flow \ |
52 |
$(PKG)/autodiff/initialisation.flow \ |
53 |
$(PKG)/autodiff/mdsio.flow \ |
54 |
$(PKG)/autodiff/monitor.flow \ |
55 |
$(PKG)/autodiff/checkpoint.flow \ |
56 |
$(PKG)/autodiff/print_message.flow \ |
57 |
$(PKG)/autodiff/open_copy_data_file.flow\ |
58 |
$(PKG)/autodiff/timeave.flow \ |
59 |
\ |
60 |
$(SOURCE_CODE)/calc_buoyancy.f \ |
61 |
$(SOURCE_CODE)/calc_common_factors.f \ |
62 |
$(SOURCE_CODE)/calc_diffusivity.f \ |
63 |
$(SOURCE_CODE)/calc_div_ghat.f \ |
64 |
$(SOURCE_CODE)/calc_grad_phi_surf.f \ |
65 |
$(SOURCE_CODE)/calc_gs.f \ |
66 |
$(SOURCE_CODE)/calc_gt.f \ |
67 |
$(SOURCE_CODE)/calc_gtr1.f \ |
68 |
$(SOURCE_CODE)/calc_ivdc.f \ |
69 |
$(SOURCE_CODE)/calc_mom_rhs.f \ |
70 |
$(SOURCE_CODE)/calc_phi_hyd.f \ |
71 |
$(SOURCE_CODE)/calc_viscosity.f \ |
72 |
$(SOURCE_CODE)/comm_stats.f \ |
73 |
$(SOURCE_CODE)/convect.f \ |
74 |
$(SOURCE_CODE)/convective_adjustment.f \ |
75 |
$(SOURCE_CODE)/convective_adjustment_ini.f \ |
76 |
$(SOURCE_CODE)/correction_step.f \ |
77 |
$(SOURCE_CODE)/cycle_tracer.f \ |
78 |
$(SOURCE_CODE)/different_multiple.f \ |
79 |
$(SOURCE_CODE)/do_fields_blocking_exchanges.f \ |
80 |
$(SOURCE_CODE)/do_the_model_io.f \ |
81 |
$(SOURCE_CODE)/dynamics.f \ |
82 |
$(SOURCE_CODE)/external_forcing.f \ |
83 |
$(SOURCE_CODE)/external_forcing_surf.f \ |
84 |
$(SOURCE_CODE)/find_alpha.f \ |
85 |
$(SOURCE_CODE)/find_rho.f \ |
86 |
$(SOURCE_CODE)/forward_step.f \ |
87 |
$(SOURCE_CODE)/freeze.f \ |
88 |
$(SOURCE_CODE)/get_offline_fields.f \ |
89 |
$(SOURCE_CODE)/gad_biharm_x.f \ |
90 |
$(SOURCE_CODE)/gad_biharm_y.f \ |
91 |
$(SOURCE_CODE)/gad_c2_adv_r.f \ |
92 |
$(SOURCE_CODE)/gad_c2_adv_x.f \ |
93 |
$(SOURCE_CODE)/gad_c2_adv_y.f \ |
94 |
$(SOURCE_CODE)/gad_c4_adv_r.f \ |
95 |
$(SOURCE_CODE)/gad_c4_adv_x.f \ |
96 |
$(SOURCE_CODE)/gad_c4_adv_y.f \ |
97 |
$(SOURCE_CODE)/gad_calc_rhs.f \ |
98 |
$(SOURCE_CODE)/gad_del2.f \ |
99 |
$(SOURCE_CODE)/gad_diff_r.f \ |
100 |
$(SOURCE_CODE)/gad_diff_x.f \ |
101 |
$(SOURCE_CODE)/gad_diff_y.f \ |
102 |
$(SOURCE_CODE)/gad_fluxlimit_adv_r.f \ |
103 |
$(SOURCE_CODE)/gad_fluxlimit_adv_x.f \ |
104 |
$(SOURCE_CODE)/gad_fluxlimit_adv_y.f \ |
105 |
$(SOURCE_CODE)/gad_grad_x.f \ |
106 |
$(SOURCE_CODE)/gad_grad_y.f \ |
107 |
$(SOURCE_CODE)/gad_u3_adv_r.f \ |
108 |
$(SOURCE_CODE)/gad_u3_adv_x.f \ |
109 |
$(SOURCE_CODE)/gad_u3_adv_y.f \ |
110 |
$(SOURCE_CODE)/grad_sigma.f \ |
111 |
$(SOURCE_CODE)/impldiff.f \ |
112 |
$(SOURCE_CODE)/ini_cartesian_grid.f \ |
113 |
$(SOURCE_CODE)/ini_cg2d.f \ |
114 |
$(SOURCE_CODE)/ini_communication_patterns.f \ |
115 |
$(SOURCE_CODE)/ini_cori.f \ |
116 |
$(SOURCE_CODE)/ini_curvilinear_grid.f \ |
117 |
$(SOURCE_CODE)/ini_depths.f \ |
118 |
$(SOURCE_CODE)/ini_fields.f \ |
119 |
$(SOURCE_CODE)/ini_forcing.f \ |
120 |
$(SOURCE_CODE)/ini_grid.f \ |
121 |
$(SOURCE_CODE)/ini_linear_phisurf.f \ |
122 |
$(SOURCE_CODE)/ini_masks_etc.f \ |
123 |
$(SOURCE_CODE)/ini_mixing.f \ |
124 |
$(SOURCE_CODE)/ini_p_ground.f \ |
125 |
$(SOURCE_CODE)/ini_pnh.f \ |
126 |
$(SOURCE_CODE)/ini_procs.f \ |
127 |
$(SOURCE_CODE)/ini_psurf.f \ |
128 |
$(SOURCE_CODE)/ini_salt.f \ |
129 |
$(SOURCE_CODE)/ini_spherical_polar_grid.f \ |
130 |
$(SOURCE_CODE)/ini_theta.f \ |
131 |
$(SOURCE_CODE)/ini_tr1.f \ |
132 |
$(SOURCE_CODE)/ini_dynvars.f \ |
133 |
$(SOURCE_CODE)/ini_vel.f \ |
134 |
$(SOURCE_CODE)/ini_vertical_grid.f \ |
135 |
$(SOURCE_CODE)/initialise_varia.f \ |
136 |
$(SOURCE_CODE)/integrate_for_w.f \ |
137 |
$(SOURCE_CODE)/modeldata_example.f \ |
138 |
$(SOURCE_CODE)/mom_calc_hfacz.f \ |
139 |
$(SOURCE_CODE)/mom_calc_ke.f \ |
140 |
$(SOURCE_CODE)/mom_cdscheme.f \ |
141 |
$(SOURCE_CODE)/mom_u_adv_uu.f \ |
142 |
$(SOURCE_CODE)/mom_u_adv_vu.f \ |
143 |
$(SOURCE_CODE)/mom_u_adv_wu.f \ |
144 |
$(SOURCE_CODE)/mom_u_bottomdrag.f \ |
145 |
$(SOURCE_CODE)/mom_u_del2u.f \ |
146 |
$(SOURCE_CODE)/mom_u_metric_nh.f \ |
147 |
$(SOURCE_CODE)/mom_u_metric_sphere.f \ |
148 |
$(SOURCE_CODE)/mom_u_rviscflux.f \ |
149 |
$(SOURCE_CODE)/mom_u_sidedrag.f \ |
150 |
$(SOURCE_CODE)/mom_u_xviscflux.f \ |
151 |
$(SOURCE_CODE)/mom_u_yviscflux.f \ |
152 |
$(SOURCE_CODE)/mom_v_adv_uv.f \ |
153 |
$(SOURCE_CODE)/mom_v_adv_vv.f \ |
154 |
$(SOURCE_CODE)/mom_v_adv_wv.f \ |
155 |
$(SOURCE_CODE)/mom_v_bottomdrag.f \ |
156 |
$(SOURCE_CODE)/mom_v_del2v.f \ |
157 |
$(SOURCE_CODE)/mom_v_metric_nh.f \ |
158 |
$(SOURCE_CODE)/mom_v_metric_sphere.f \ |
159 |
$(SOURCE_CODE)/mom_v_rviscflux.f \ |
160 |
$(SOURCE_CODE)/mom_v_sidedrag.f \ |
161 |
$(SOURCE_CODE)/mom_v_xviscflux.f \ |
162 |
$(SOURCE_CODE)/mom_v_yviscflux.f \ |
163 |
$(SOURCE_CODE)/packages_init_variables.f \ |
164 |
$(SOURCE_CODE)/plot_field.f \ |
165 |
$(SOURCE_CODE)/solve_for_pressure.f \ |
166 |
$(SOURCE_CODE)/state_summary.f \ |
167 |
$(SOURCE_CODE)/timestep.f \ |
168 |
$(SOURCE_CODE)/the_correction_step.f \ |
169 |
$(SOURCE_CODE)/the_main_loop.f \ |
170 |
$(SOURCE_CODE)/thermodynamics.f \ |
171 |
$(SOURCE_CODE)/timestep_tracer.f \ |
172 |
$(SOURCE_CODE)/swfrac.f \ |
173 |
\ |
174 |
$(SOURCE_CODE)/kpp_calc.f \ |
175 |
$(SOURCE_CODE)/kpp_calc_diff.f \ |
176 |
$(SOURCE_CODE)/kpp_calc_visc.f \ |
177 |
$(SOURCE_CODE)/kpp_do_diags.f \ |
178 |
$(SOURCE_CODE)/kpp_init.f \ |
179 |
$(SOURCE_CODE)/kpp_open_diags.f \ |
180 |
$(SOURCE_CODE)/kpp_routines.f \ |
181 |
$(SOURCE_CODE)/kpp_transport_s.f \ |
182 |
$(SOURCE_CODE)/kpp_transport_t.f \ |
183 |
\ |
184 |
$(SOURCE_CODE)/gmredi_calc_diff.f \ |
185 |
$(SOURCE_CODE)/gmredi_calc_tensor.f \ |
186 |
$(SOURCE_CODE)/gmredi_diags.f \ |
187 |
$(SOURCE_CODE)/gmredi_init.f \ |
188 |
$(SOURCE_CODE)/gmredi_rtransport.f \ |
189 |
$(SOURCE_CODE)/gmredi_slope_limit.f \ |
190 |
$(SOURCE_CODE)/gmredi_xtransport.f \ |
191 |
$(SOURCE_CODE)/gmredi_ytransport.f \ |
192 |
\ |
193 |
$(SOURCE_CODE)/ctrl_map_ini.f \ |
194 |
$(SOURCE_CODE)/ctrl_map_forcing.f \ |
195 |
\ |
196 |
$(SOURCE_CODE)/cost_final.f \ |
197 |
$(SOURCE_CODE)/cost_init.f \ |
198 |
$(SOURCE_CODE)/cost_tile.f \ |
199 |
$(SOURCE_CODE)/cost_test.f \ |
200 |
$(SOURCE_CODE)/cost_tracer.f \ |
201 |
\ |
202 |
|
203 |
|
204 |
#-------------------------------------------------------- |
205 |
# rules |
206 |
#-------------------------------------------------------- |
207 |
|
208 |
adcode: |
209 |
cd $(DEST) ; cat $(SRC_MODEL) >! tamc_code_ecco.f |
210 |
admodeltamc: |
211 |
cd $(DEST) ; $(TAMC) $(TAMCFLAG) tamc_code_ecco.f \ |
212 |
>! tamc_code_ecco_ad.prot |
213 |
|
214 |
admodeltaf: |
215 |
cd $(DEST) ; $(TAF) $(TAFFLAG) tamc_code_ecco.f |
216 |
|
217 |
adchange: |
218 |
${ADJOINT_SCRIPT}/adjoint_ecco_sed.com \ |
219 |
>! tamc_code_ecco_sed_ad.f ; \ |
220 |
cp $(PKG)/autodiff/adjoint_model.F \ |
221 |
$(PKG)/autodiff/adjoint_model.F_old ; \ |
222 |
cat tamc_code_ecco_sed_ad.f >> $(PKG)/autodiff/adjoint_model.F |
223 |
adrestore: |
224 |
cp $(PKG)/autodiff/adjoint_model.F_template \ |
225 |
$(PKG)/autodiff/adjoint_model.F |
226 |
|
227 |
adtamc: adcode admodeltamc |
228 |
|
229 |
adtaf: adcode admodeltaf |