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	ifenty	1.1.2.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.