Changes between Version 3 and Version 4 of ticket/1296
- Timestamp:
- 2014-05-15T17:26:41+02:00 (10 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
ticket/1296
v3 v4 19 19 1) Damping climatologies can be different to climatological start files - use separate version of dia_tsd_read or implement new version?[[BR]]2) No grid dependence - controlled by mask file only[[BR]]3) Region entirely specified by mask file (both horizontal and vertical)[[BR]]4) Mask file can be used for varying coefficient (use values between 0 and 1) - i.e. 0 for no relxation or 1 for relaxation at the time scale specified in namelist[[BR]]5) Should work with VVL[[BR]]6) Make sure only one of tra_dmp or tra_rlx is used at one time. 20 20 21 '''Namelist requirements:''' 22 23 ''''''Namelist requirements:[[BR]]Although a new namelist section is needed this will eventually be a replacement for tra_dmp so still a simplification[[BR]]Needs:[[BR]]1) Mask file name[[BR]]2) Relaxation time scale (in seconds or timesteps?)[[BR]]3) Switch to damp in full column/only in mixed layer/only below mixed layer - this is current functionality that can't be reproduced using a maskfile.[[BR]]4) Mask variable name? 21 '''Namelist requirements:''[[BR]]'''''Although a new namelist section is needed this will eventually be a replacement for tra_dmp so still a simplification[[BR]]Needs:[[BR]]1) Mask file name[[BR]]2) Relaxation time scale (in seconds or timesteps?)[[BR]]3) Switch to damp in full column/only in mixed layer/only below mixed layer - this is current functionality that can't be reproduced using a maskfile.[[BR]]4) Mask variable name? 24 22 25 23 '''Subroutines:''' 26 24 27 - Copy of dia_tsd_read or modify dia_tsd_read to allow it to be used with different files (after initial conditions have been read)[[BR]]- tra_rlx_init: [[BR]] read the namelist, read the mask file, calculate restoration coefficient[[BR]]- tra_rlx:[[BR]] Calculate the trend at each point = rlx_mask * resto * (zts_dta - tsb) (where zts_dta is the climatological value, resto is the restoration coefficient) and add to total trend[[BR]]The climatology that we are relaxing to should be on the original model grid. Given the long time scales of the climatology compared to variations in SSH it seems reasonable not worry about the slight changes in vertical grid spacing.[[BR]]- Link to trends diagnostics to output the trends?25 * Copy of dia_tsd_read or modify dia_tsd_read to allow it to be used with different files (after initial conditions have been read)[[BR]]- tra_rlx_init: [[BR]] read the namelist, read the mask file, calculate restoration coefficient[[BR]]- tra_rlx:[[BR]] Calculate the trend at each point = rlx_mask * resto * (zts_dta - tsb) (where zts_dta is the climatological value, resto is the restoration coefficient) and add to total trend[[BR]] The climatology that we are relaxing to should be on the original model grid. Given the long time scales of the climatology compared to variations in SSH it seems reasonable not worry about the slight changes in vertical grid spacing.[[BR]]- Link to trends diagnostics to output the trends? 28 26 29 27 ----