Changes between Version 3 and Version 4 of ticket/0927_Energy_diag_v3.5


Ignore:
Timestamp:
2013-04-27T17:42:53+02:00 (8 years ago)
Author:
gm
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • ticket/0927_Energy_diag_v3.5

    v3 v4  
    1212'''     Motivation: '''    output 3D trends of tracers, momentum, kinetic energy and potential energy.[[BR]] '''     Status :''' the extraction of trends terms exists, but not the 3D output of the trends   [[BR]] '''     Main tasks :'''         [[BR]] 
    1313 
    14   (1) implement the 3D output of tracers and momentum trends using iom_put [[BR]]                                               (2) compute and output the 3D trends of PE and KE. In particular introduce a consistent way of computing PE[[BR]]    (3) add Vallis equation of state in replacement of linear in T and in T-S eos (nn_eos=2)[[BR]]    (4) add the exact Jacket and McDougall     eos (nn_eos=-1) keep for backward compatibility the nn_eos=0 old case[[BR]]    (5) add TEOS10 equation of state[[BR]]    (6) validation + documentation 
     14  (1) implement the 3D output of tracers and momentum trends using iom_put [[BR]]                                                     (2) compute and output the 3D trends of PE and KE. In particular introduce a consistent way of computing PE[[BR]]          (3) add Vallis equation of state in replacement of linear in T and in T-S eos (nn_eos=2)[[BR]]          (4) add the exact Jacket and McDougall           eos (nn_eos=-1) keep for backward compatibility the nn_eos=0 old case[[BR]]          (5) add TEOS10 equation of state[[BR]]          (6) validation + documentation 
    1515 
    1616'''     Science Reviewer:''' NOCS guy?     [[BR]]     '''     System Reviewer:'''  NOCS guy?     [[BR]]    '''     Deadline:'''   automn 2013 [[BR]] '''     Priority:'''      high    [[BR]] '''     Depends on:''' - [[BR]] '''     Principal Investigator : '''      Gurvan Madec and Fabien Roquet (gurvan.madec@locean-ipsl.upmc.fr)  [[BR]] [[BR]] 
     
    2525see Changeset [https://forge.ipsl.jussieu.fr/nemo/changeset/3874 "[3874]"] ; [https://forge.ipsl.jussieu.fr/nemo/changeset/3876 "[3876]"] and [https://forge.ipsl.jussieu.fr/nemo/changeset/3878 "[3878]"] 
    2626 
    27 '''II. First revision of the eos '''(see''' '''Changeset [https://forge.ipsl.jussieu.fr/nemo/changeset/3893 3893] ) 
     27'''II. First revision of the eos '''(see''' '''Changeset [https://forge.ipsl.jussieu.fr/nemo/changeset/3893 3893] ) 
    2828 
    2929• introduce Vallis eos as the alternative to UNESCO type eos, 
     
    3939  NB:  rhd and rhop are kept for backward compatibility, but they are meant to become obsolete (especially rhop which could be easily removed, with minor changes in zdfkpp, zdfmxl and zdfric 
    4040 
    41 [[BR]] '''Pending issues''' : 
     41'''III. Second revision of the eos '''(see''' '''Changeset ) 
     42 
     43• simplification in eosbn2 : 
     44 
     45 * density anomaly computation in 3D and 2D cases regrouped in 1 routine (use of "fake" 3D/4D arrays in argument for 2D calculation, and systematic use of depth as an argument) 
     46 * N!^2 computed using vertical interpolation of rab (alpha & beta) at w-point.  
     47   this has been shown to have negligible effect on the precision at which N!^2 is computed 
     48 * suppression of the calculation of rrau in lk_zdfddm case. ==> zdfddm.F90 : the calculation of rrau is now local to zdfddm using zab_n 
     49 
     50• use of rab_b or rab_n when possible ( zdfkpp, zdfddm, trabbl) 
     51 
     52  trabbl case: use of interpolation of alpha and beta at T-point instead of re-calculation of alpha and beta with mean T, S and depth. 
     53 
     54• Freezing point: change the name of tfreez into eos_fzp, introduce the depth as a OPTIONAL argument, and allow 2D and 3D calculations. 
     55 
     56This will be usefull for under iceshelf seas and can already be used in asminc.F90 
     57 
     58---- 
     59'''Verification to be shown:''' 
     60 
     61- N!^2 negligible effect of use of interpolation of alpha & beta at T-point instead re-calculation of alpha & beta at w-point ? 
     62 
     63- bbl  negligible effect of use of interpolation of alpha & beta at T-point instead of re-calculation of alpha & beta with mean T, S and depth 
     64 
     65---- 
     66• '''Pending issues''' : 
    4267 
    4368 * diagnostic of trends missing in the following modules :                             '''==>>  To be done !!!!''' ''' 
     
    5984 * reshape trdtra so that T and S are treated separately in all trd routine (including mld diag...) create umask_i and vmask_i (2D) fields that mirror task_i field but for the velocity points 
    6085 
    61 ---- 
    6286'''Fabien addition (copy from the v3.4 page):''' 
    6387