 |
 |
 |
||||||||||
 |
|
|||||||||||
 |
 |
 |
||||||||||
Next: 3.19 Sensitivity of Air-Sea Up: 3.18 Global Ocean State Previous: 3.18.4 Compiling Contents 3.18.5 Running the estimationCopy the mitgcmuv_ad executable to input_ad and optim.x to the subdirectory input_ad/OPTIM. Move into input_ad/. The first iteration has be done "by hand". Check that the iteration number is set to 0 in data.optim and run the Mitgcm % ./mitgcmuv_ad
The output files adxx_hfluxm.0000000000.* and xx_hfluxm.0000000000.* contain
the sensitivity of the cost function to
At the first iteration, move these two files ecco_cost and ecco_ctrl to OPTIM/, open data.optim and check the iteration number is set to 0. The target cost function fmin needs to be specified: a rule of thumb suggest it should be set to about 0.95-0.90 times the value of the cost function at the first iteration. This value is only used at the first iteration and does not need to be updated afterwards. However, it implicitly specifies the "pace" at which the cost function is going down (if you are lucky and it goes down). More in the ECCO section maybe ? Once this is done, run the line-search algorithm: % ./optim.xwhich computes the updated perturbations for iteration 1 ecco_ctrl_1. The following iterations can be executed automatically using the shell script cycsh found in input_ad/. This script will take care of changing the iteration numbers in the data.optim, launch the adjoint model, clean and store the ouputs, move the ecco_cost* and ecco_ctrl* files, and run the line-search algotrithm. Edit cycsh to specify the prefix of the directories used to store the ouputs and the maximum number of iteration.
Next: 3.19 Sensitivity of Air-Sea Up: 3.18 Global Ocean State Previous: 3.18.4 Compiling Contents mitgcm-support@mitgcm.org
|
|||||
and the perturbations
to
and an input of the line-search. The
latter returns an updated perturbation in ecco_ctrl*n+1 to be used as an input of
the adjoint model at iteration n+1.
