verification/1D_ocean_ice_column/README_11K_TIME_STEP_SIMULATION.TXT

HOW TO RUN AND ANALYZE A 11,000 TIME STEP EXPERIMENT
JAN 23, 2013

MITgcm checkpoint 64c.  
seaice_growth.F version: 1.185 (2012/12/27)

1) Make the following changes to the input_ad/data* files 
for the 11,000 time step simulation.

a) data 

>  nTimeSteps= 11000,
>  monitorFreq=0.
>  adjMonitorFreq = 1.,

b) data.autodiff 

>   useKPPinAdMode = .FALSE.,

c) data.diagnostics

>   frequency(1)  = 86400,
>   fields(1,1)   = 'EXFuwind','EXFvwind','EXFlwdn', 'EXFswdn','EXFatemp','EXFaqh','EXFevap','EXFpreci','EXFempmr','EXFhl   ','EXFhs   ',
>   filename(1)   = 'EXF_avg_daily',

>   frequency(2)  = 86400,
>   fields(1:3,2)   = 'SIarea  ','SIheff  ','SIhsnow ',
>   filename(2)   = 'ICE_avg_daily',

d) data.pkg

>   useGrdchk = .FALSE.,
>   useMNC    = .TRUE,
>   useDiagnostics = .TRUE.,


2) Make the following changes to packages.conf in code_ad:

a) packages.conf 

> diagnostics
> mnc


3) After running, use the two matlab scripts (analyze_forward_run_general.m, analyze_adjoint_run_general.m) in ./matlab_scripts to plot the results.  Be sure that the helper script, 'tenPercentAboveBelowLeftRight.m' is in your matlab path. 
1	HOW TO RUN AND ANALYZE A 11,000 TIME STEP EXPERIMENT
2	JAN 23, 2013
3
4	MITgcm checkpoint 64c.
5	seaice_growth.F version: 1.185 (2012/12/27)
6
7	1) Make the following changes to the input_ad/data* files
8	for the 11,000 time step simulation.
9
10	a) data
11
12	> nTimeSteps= 11000,
13	> monitorFreq=0.
14	> adjMonitorFreq = 1.,
15
16	b) data.autodiff
17
18	> useKPPinAdMode = .FALSE.,
19
20	c) data.diagnostics
21
22	> frequency(1) = 86400,
23	> fields(1,1) = 'EXFuwind','EXFvwind','EXFlwdn', 'EXFswdn','EXFatemp','EXFaqh','EXFevap','EXFpreci','EXFempmr','EXFhl ','EXFhs ',
24	> filename(1) = 'EXF_avg_daily',
25
26	> frequency(2) = 86400,
27	> fields(1:3,2) = 'SIarea ','SIheff ','SIhsnow ',
28	> filename(2) = 'ICE_avg_daily',
29
30	d) data.pkg
31
32	> useGrdchk = .FALSE.,
33	> useMNC = .TRUE,
34	> useDiagnostics = .TRUE.,
35
36
37	2) Make the following changes to packages.conf in code_ad:
38
39	a) packages.conf
40
41	> diagnostics
42	> mnc
43
44
45	3) After running, use the two matlab scripts (analyze_forward_run_general.m, analyze_adjoint_run_general.m) in ./matlab_scripts to plot the results. Be sure that the helper script, 'tenPercentAboveBelowLeftRight.m' is in your matlab path.