New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
ticket/0863 (diff) – NEMO

Changes between Version 19 and Version 20 of ticket/0863


Ignore:
Timestamp:
2011-10-21T16:30:57+02:00 (12 years ago)
Author:
mlelod
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • ticket/0863

    v19 v20  
    124124=== Time splitting external mode === 
    125125=== Coupling between ocean and sea ice === 
     126In 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 
    126128=== 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  diffusion routines have thus been modified to use these new variables.[[BR]]This could be done in many other routines. 
     129The 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. 
    128130 
    129131=== Isopycnal diffusion === 
     
    131133 
    132134=== Thickness weighted outputs === 
    133 Still to do. 
     135In 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: 
    134136 
     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]] 
    135149---- 
    136150== Testing ==