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

Timestamp:

2013-04-27T17:42:53+02:00 (11 years ago)

Author:

gm

Comment:

--

ticket/0927_Energy_diag_v3.5

@@ -12 +12 @@
 '''     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]]
 
-  (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
+  (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
 
 '''     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]]
@@ -25 +25 @@
 see 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]"]
 
-'''II. First revision of the eos '''(see''' '''Changeset [https://forge.ipsl.jussieu.fr/nemo/changeset/3893 3893] )
+'''II. First revision of the eos '''(see''' '''Changeset [https://forge.ipsl.jussieu.fr/nemo/changeset/3893 3893] )
 
 • introduce Vallis eos as the alternative to UNESCO type eos,
@@ -39 +39 @@
   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
 
-[[BR]] '''Pending issues''' :
+'''III. Second revision of the eos '''(see''' '''Changeset )
+
+• simplification in eosbn2 :
+
+ * 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)
+ * N!^2 computed using vertical interpolation of rab (alpha & beta) at w-point. 
+   this has been shown to have negligible effect on the precision at which N!^2 is computed
+ * suppression of the calculation of rrau in lk_zdfddm case. ==> zdfddm.F90 : the calculation of rrau is now local to zdfddm using zab_n
+
+• use of rab_b or rab_n when possible ( zdfkpp, zdfddm, trabbl)
+
+  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.
+
+• Freezing point: change the name of tfreez into eos_fzp, introduce the depth as a OPTIONAL argument, and allow 2D and 3D calculations.
+
+This will be usefull for under iceshelf seas and can already be used in asminc.F90
+
+----
+'''Verification to be shown:'''
+
+- N!^2 negligible effect of use of interpolation of alpha & beta at T-point instead re-calculation of alpha & beta at w-point ?
+
+- 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
+
+----
+• '''Pending issues''' :
 
  * diagnostic of trends missing in the following modules :                             '''==>>  To be done !!!!''' '''
@@ -59 +84 @@
  * 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
 
-----
 '''Fabien addition (copy from the v3.4 page):'''