Changes between Version 9 and Version 10 of 2009WP/2009Stream3/LeapFrogC
- Timestamp:
- 2009-06-22T08:13:10+02:00 (15 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
2009WP/2009Stream3/LeapFrogC
v9 v10 39 39 40 40 ---- 41 == 0. Flow Chart==41 === 0. Flow Chart === 42 42 The flow chart of ''step'' used must be the one introduce in the trunk after the v3.1, to restore the key_vvl option. Its main characteristic is that the ''ssh'' is time-stepped first, followed immediately by the computation of the ''now'' vertical velocity. [[BR]] 43 43 … … 46 46 47 47 ---- 48 == 1. Forcing evaluation==48 === 1. Forcing evaluation === 49 49 50 50 In order to evaluate the forcing at kt+1/2, all the input data in the surface module should be interpolated in time at kt+1/2. This have been introduced in fldread.F90 module since release 3.0. There is no required additional action.[[BR]] … … 52 52 53 53 ---- 54 == 2. Forcing : time discretisation==54 === 2. Forcing : time discretisation === 55 55 Each time a forcing term appear in the 3 time evolution equation solved in NEMO/OPA (''ssh'', momentum, tracer), it is the mean of kt+1/2 and kt+1/2 which is used. Therefore, the surface module (including the sea-ice!) must now provide both kt+1/2 '''and''' kt-1/2 forcing fields, and all the module involving a forcing term have to be modified.[[BR]] 56 56 57 === • surface module (SBC)===58 59 ==== - sbc_oce.F90 ====60 61 add before field (denoted by a '_b' sufix) for all ocean forcing (utau_b, vtau_b, qns_b, qsr_b, emp_b, emps_b)57 ==== • '''surface module''' (SBC) ==== 58 59 ===== - sbc_oce.F90 ==== 60 61 add before field (denoted by a '_b' sufix) for all ocean forcing (utau_b, vtau_b, qns_b, qsr_b, emp_b, emps_b) 62 62 63 63 the following lines … … 176 176 [[BR]] 177 177 178 === • tracer forcing term (TRA)===178 ==== • '''tracer forcing term''' (TRA) ==== 179 179 180 180 ==== - trasbc.F90 ==== … … 327 327 328 328 [[BR]] 329 === • momentum forcing term (DYN)===329 ==== • '''momentum forcing term''' (DYN) ==== 330 330 331 331 ==== - dynzdf_exp.F90 ====