Changes between Version 2 and Version 3 of Working Groups/TAM/Reference Manual/directmodel
 Timestamp:
 20100219T11:30:54+01:00 (11 years ago)
Legend:
 Unmodified
 Added
 Removed
 Modified

Working Groups/TAM/Reference Manual/directmodel
v2 v3 14 14 15 15 == 1  NonLinear Trajectory == 16 To evaluate the tangent and/or the adjoint at a time t, some actual parameters of the direct model are needed at the corresponding time t. To fulfill this purpose a subroutine, '' asm_trj_wri'' is added to write the trajectory. It is controlled with two logical parameters: ''lk_asminc'' and ''ln_trjwri''. The trajectory is saved with a certain saving frequency and only both logical parameter areset to TRUE.16 To evaluate the tangent and/or the adjoint at a time t, some actual parameters of the direct model are needed at the corresponding time t. To fulfill this purpose a subroutine, ''tamtrj_wri'' is added to write the trajectory. It is controlled with one logical parameter: ''ln_trjwri''. The trajectory is saved upon a given frequency "nittrjfrq" only if the logical parameter is set to TRUE. 17 17 18 === Subroutine: asm_trj_wri ===19 A new subdirectory is created into /OPA/ named ''' ASM'''. The file is '''asmtrj.F90'''.20 The purpose of this routine is to write to file the model state trajectory for use with 4DVar. It is called in ''opa_model'' routine to save the initial state needed for the TAM and at the end of the ''stp'' routine.18 === Subroutine: tam_trj_wri === 19 A new subdirectory is created into /OPA/ named '''TAM'''. The file is '''tamtrj.F90'''. 20 The purpose of this routine is to write to file the model state trajectory for use with 4DVar. It is called in ''opa_model'' routine to initialize tamtrj routine TAM and at in the ''stp'' routine. 21 21 22 22 Currently the following parameters are saved: … … 45 45 hmlp  key_tradmp  46 46 47 ''note: is it worth to add why (and where) there will be a need of these variables later in TAM?'' 47 ==== Specific management of the initial state ==== 48 For the initial state of the model, we call the tam_trj_wri in 'step' routine just before the tracer block. We,then, have access to some parameters (tke and sbc for instance) as they are not yet computed in opa (outside 'step'). 48 49 49 ==== Specific management of the initial state ====50 For the initial state of the model, some parameters need to be computed in ''asmtrj'' as they are not yet computed in opa.51 52 ''note: we still have zero value in emp/emps parameters. It needs to be solved quickly.''53 50 54 51 ==== Linear Interpollation ==== 55 For efficiency matters, the whole trajectory is of course not saved. Only a predefined steptime frequency is saved. The frequency ''nittrjfr eq'' is controlled by namelist. The intermediate steptime will be estimated in TAM by linear interpollation.52 For efficiency matters, the whole trajectory is of course not saved. Only a predefined steptime frequency is saved. The frequency ''nittrjfrq'' is controlled by namelist. The intermediate steptime will be estimated in TAM by linear interpollation. 56 53 57 54 ==== Namelist ==== 58 Parameters related to the model state saving is controlled by the namelist '''nam _asminc''' with:55 Parameters related to the model state saving is controlled by the namelist '''namtam''' with: 59 56 60 57 * ln_trjwri: Logical switch for writing out state trajectory 61 * nittrjfrq: frequency of the trajectory used in 4DVAR58 * nittrjfrq: saving frequency of the trajectory 62 59 63 60 == 2  Routine modification == 64 61 65 62 * '''opa.F90''': 66 * include output trajectory fields (USE asmtrj and CALL asm_trj_wri)63 * include initilialization of trajectory saving handling (USE tamtrj) 67 64 * '''step.F90''': 68 * include output trajectory fields (CALL asm_trj_wri)69 * add temporary variables to save ta/sa (used as workspace in standard)65 * include calls for trajectory saving (CALL tam_trj_wri) 66 * add temporary variables (zta, zsa) to save ta/sa (used as workspace in dynamic block) 70 67 71 68 == 3  Declaration Changes ==