Changes between Version 19 and Version 20 of ticket/0863
- Timestamp:
- 2011-10-21T16:30:57+02:00 (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
ticket/0863
v19 v20 124 124 === Time splitting external mode === 125 125 === Coupling between ocean and sea ice === 126 In the vvl case `fse3t_m` is now a statement function for `e3t_m` whis is a 2D variable. In non-vvl cases `fse3t_m(:,:)` stands for fse3t_0(:,:,1).[[BR]]`fse3t_m` is treated in the vvl case as all other `*_m` variables in `sbcssm` module. 127 126 128 === Tracer lateral diffusion === 127 The 2D fields `e2_1u = e2u / e1u` and `e1_2v = e1v / e2v` are used in the computation of the thickness diffusion term. They were allready computed and saved locally in `traldf_lap` module under the names e1ur and e2vr. Because they are used in all other tracer diffusion routines (and now the thickness diffusion term in `dom_vvl_sf_nxt`), they are now declared in `dom_oce` and computed in `domain` just after all e12x and e12x_1 fields. The different tracer lateral 129 The 2D fields `e2_1u = e2u / e1u` and `e1_2v = e1v / e2v` are used in the computation of the thickness diffusion term. They were allready computed and saved locally in `traldf_lap` module under the names e1ur and e2vr. Because they are used in all other tracer diffusion routines (and now the thickness diffusion term in `dom_vvl_sf_nxt`), they are now declared in `dom_oce` and computed in `domain` just after all e12x and e12x_1 fields. The different tracer lateral diffusion routines have thus been modified to use these new variables.[[BR]]This could be done in many other routines. 128 130 129 131 === Isopycnal diffusion === … … 131 133 132 134 === Thickness weighted outputs === 133 Still to do. 135 In vvl case, the "toce" variable in the netcdf file does not contain temperature anymore but the heat content fse3t_n * tn.[[BR]]The same thing is done for salinity and horizontal velocities (only in flux form advection but I'm not sure if it has to be done in any vvl case).[[BR]]3 other variables are created in netcdf outputs: 134 136 137 * e3t_n is the vertical scale factor at t-points (_n suffix indicates that it is different from the background scale factor).[[BR]]e3u and e3v and e3w are not stored, they can be recomputed off-line. 138 * dept_n is the depth at t-points. Depth at u, v and w points should be computed off-line also. 139 * e3tdef is the squared deformation ratio of e3t given in %!^2.[[BR]]During the post treatment, the square root of this quantity gives the variance of the level thickness in %. 140 141 {{{ 142 z_e3t_def(:,:,:) = ( ( fse3t_n(:,:,:) - fse3t_0(:,:,:) ) / fse3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 143 CALL iom_put( "e3t_n" , fse3t_n (:,:,:) ) 144 CALL iom_put( "dept_n" , fsdept_n (:,:,:) ) 145 CALL iom_put( "e3tdef" , z_e3t_def(:,:,:) ) 146 147 }}} 148 [[BR]] 135 149 ---- 136 150 == Testing ==