Changes between Version 66 and Version 67 of ticket/0927_Energy_diag
- Timestamp:
- 2013-04-24T15:55:40+02:00 (11 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
ticket/0927_Energy_diag
v66 v67 7 7 '''Branch v3.4''' : [https://forge.ipsl.jussieu.fr/nemo/browser/branches/2012/dev_r3309_LOCEAN12_Ediag 2012/dev_r3309_LOCEAN12_Ediag] 8 8 9 '''Branch v3.5''' : [https://forge.ipsl.jussieu.fr/nemo/browser/branches/2013/dev_r3858_CNRS3_Ediag 2013/dev_r3858_CNRS3_Ediag] 9 '''The development has been transfer to v3.5''' 10 11 see [https://forge.ipsl.jussieu.fr/nemo/wiki/ticket/927_Energy_diag_v3.5 v3.5 Ediag wiki page] 12 13 and Branch v3.5''' : [https://forge.ipsl.jussieu.fr/nemo/browser/branches/2013/dev_r3858_CNRS3_Ediag 2013/dev_r3858_CNRS3_Ediag]''' 10 14 11 15 ---- … … 14 18 ''' 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]] 15 19 16 (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 [[BR]](3) validation + documentation [[BR]]17 18 ''' 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]]20 (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 [[BR]] (3) validation + documentation [[BR]] 21 22 ''' 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]] 19 23 20 24 ---- … … 26 30 see [https://forge.ipsl.jussieu.fr/nemo/log/branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC?rev=3316 revision 3316] 27 31 28 ''' trdmod_oce''' module and ''' namtrd namelist'''[[BR]]32 ''' trdmod_oce''' module and ''' namtrd namelist''' [[BR]] 29 33 30 34 logical flags added in namlist namtrd which now controls what is done with the trends. [[BR]] All the types of treatment of a given trend are available at the same time. The memory requirement will only increase due to the time averaged arrays defined in IOM. … … 42 46 Add these new logical in the namelist. '''==>>> CAUTION only in the ORCA2_LIM directory''' [[BR]] NB: here is the new name set in revision number 3318[[BR]] 43 47 44 ''' trdtra'''module[[BR]]48 ''' trdtra''' module[[BR]] 45 49 46 50 Add a systematic mask of the trend.[[BR]] Change the comments to better describe the purpose of this module. Its purpose is: [[BR]] 47 51 48 'TRA' case: to regroup T & S trends and send them to trd_mod, with, in case of advection, transform the incoming advective fluxes into advctive trend (U.grad[T])[[BR]] 'TRC' case: send trend to ted_mod_trc, with, in case of advection, transform the incoming advective fluxes into advective trend [[BR]] all cases : mask the trend (''' ===>>> PROBABLY add in the module a lbc_lnk so that the trend is defined everywhere''')49 50 ''' dynadv_cen2 and _ubs'''modules[[BR]]52 'TRA' case: to regroup T & S trends and send them to trd_mod, with, in case of advection, transform the incoming advective fluxes into advctive trend (U.grad[T])[[BR]] 'TRC' case: send trend to ted_mod_trc, with, in case of advection, transform the incoming advective fluxes into advective trend [[BR]] all cases : mask the trend (''' ===>>> PROBABLY add in the module a lbc_lnk so that the trend is defined everywhere''' ) 53 54 ''' dynadv_cen2 and _ubs''' modules[[BR]] 51 55 52 56 change jpdyn_trd_had into jpdyn_trd_keg. Now in flux form _keg corresponds to the horizontal advection trends and _rvo to the metric terms[[BR]] [[BR]] 53 57 54 ''' dynnxt'''module[[BR]]58 ''' dynnxt''' module[[BR]] 55 59 56 60 add the output using sum of the total dyn trend (except asselin time filter) ("utrd_tot", "vtrd_tot") and of the asselin time filter trend ("utrd_atf", "vtrd_atf") but with a shift by one time step[[BR]] [[BR]] 57 61 58 ''' dynvor, trdvor and trdmod_oce'''modules[[BR]]62 ''' dynvor, trdvor and trdmod_oce''' modules[[BR]] 59 63 60 64 suppress the call to trd_mod in the jpdyn_trd_dat case (computation of beta.V) add add the calculation of beta.V term in 'trdvor' in jpvor_pvo case. And obviously suppress jpdyn_trd_dat from trdmod_oce [[BR]] Also suppress the jpdyn_trd_had case horizontal advection for the dynamics is 'keg' + 'vor' ; in case of flux form, 'had' is put in 'keg' and the metric terms is put in 'vor' [[BR]] there is now only 10 trends on the dynamics instead of 12[[BR]] [[BR]] 61 65 62 ''' trdmod'''module[[BR]]66 ''' trdmod''' module[[BR]] 63 67 64 68 1- introduce the new logical namelist parameters[[BR]] 2- introduce new subroutines : '''trd_budget''' : computation of the domain averaged T,T^2^, PE, KE trends formerly computes in trd_mod routine)[[BR]] ''' trd_3Diom''': output of the 3D trends using IOM [[BR]] [[BR]] 65 69 66 ''' trdicp'''module[[BR]]70 ''' trdicp''' module[[BR]] 67 71 68 72 add in trd_twr routine the computation of the vertical diffusive trend on T & S in case of iso-neutral diffusion (ln_traldf_iso=T). These trends ("ttrd_zdfp", "strd_zdfp") name zdfp for "PURE" vertical diffusion trends are output so that by difference with "zdf" trends we can access to the vertical contribution of the iso-neutral operator [[BR]] [[BR]] 69 73 70 ''' iodef.xml'''file[[BR]]74 ''' iodef.xml''' file[[BR]] 71 75 72 76 add all the trends nick name : '''==>>> CAUTION only in the ORCA2_LIM directory''' [[BR]] [[BR]] … … 75 79 ''' Step II : simplification of the structure ''' 76 80 77 ''' trdicp and trdicp_oce'''modules[[BR]]81 ''' trdicp and trdicp_oce''' modules[[BR]] 78 82 79 83 see [https://forge.ipsl.jussieu.fr/nemo/log/branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC?rev=3317 revision 3317] … … 81 85 1- suppress trdicp_oce module (put required variables in trdmod_oce parameters)[[BR]] 2- use the jptra_trd_... instead of jpicp_... 3- move trd_budget subroutine from trdmod to trdicp.F90. Suppress trd_icp routine (i.e. trd_2d, trd_3d) as the work is now simply done in trd_budget.[[BR]] 4- rescan all dyn/tra trend indices. add jptra_trd_zdfp for "PURE Kz dissusive trend when ln_traldf_iso=T (see also changes in tranxt where a call to trd_tra is done just before the swap, so that PURE Kz trends can be diagnosed in tra_trd. See also xml file). 5- jptra_trd_nsr and jptra_trd_cdt are 3D trends as their incorporate both surface forcing AND runoff, the later being possibly spread in depth ('''==>>> probably to be changed''')[[BR]] 6- suppress the key_trddyn key trdtra from trdicp[[BR]] 82 86 83 ''' trdtra - trddyn - trdmod - trdvor''' modules : impact on ''' TRA''' , ''' DYN''' , ''' TRD''' and almost all ''' TOP_SRC'''modules[[BR]]87 ''' trdtra - trddyn - trdmod - trdvor''' modules : impact on ''' TRA''' , ''' DYN''' , ''' TRD''' and almost all ''' TOP_SRC''' modules[[BR]] 84 88 85 89 see [https://forge.ipsl.jussieu.fr/nemo/log/branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC?rev=3318 revision 3318] … … 96 100 ''' Step III : addition of 3D KE diagnostics ''' 97 101 98 ''' trdken'''module[[BR]]102 ''' trdken''' module[[BR]] 99 103 100 104 see [https://forge.ipsl.jussieu.fr/nemo/log/branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC?rev=3325 revision 3325] … … 111 115 ''' Step IV : addition of 3D PE diagnostics ''' 112 116 113 ''' trdpen'''module[[BR]]117 ''' trdpen''' module[[BR]] 114 118 115 119 see [https://forge.ipsl.jussieu.fr/nemo/log/branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC?rev=3326 revision 3326] and [https://forge.ipsl.jussieu.fr/nemo/log/branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC?rev=3327 revision 3327] 0- correct a bug in trdtra and trdken … … 125 129 NB: runoff put as a source term in the divergence computation since v3.3 (see divcur), implication for source term in PE, but also tracer trends ? '''Think about that!''''''''' 126 130 127 [[BR]] [[BR]] [[BR]] 128 ---- 129 '''Step IV : port to v3.5 beta''' 130 131 132 [[BR]]april 17 2013: 133 134 • create a new branch from v3.5 : [https://forge.ipsl.jussieu.fr/nemo/browser/branches/2013/dev_r3858_CNRS3_Ediag 2013/dev_r3858_CNRS3_Ediag][[BR]] [[BR]] • merge [https://forge.ipsl.jussieu.fr/nemo/browser/branches/2012/dev_r3309_LOCEAN12_Ediag 2012/dev_r3309_LOCEAN12_Ediag] with [https://forge.ipsl.jussieu.fr/nemo/browser/branches/2013/dev_r3858_CNRS3_Ediag 2013/dev_r3858_CNRS3_Ediag][[BR]] [[BR]][[BR]][[BR]][[BR]] 135 ---- 136 '''Step IV : port to v''' 137 138 [[BR]] '''Pending issues''' : 139 140 atmospheric pressure gradient trend not taken into account (see dynspg.F90 ''' To be done !!!!''' [[BR]] ''' 141 142 kpp non-local trend put in zdf trends !!! this will not work ! a additional trend term should be add ''' To be done !!!!''' [[BR]] ''' 131 ---- 132 '''Pending issues''' : 133 134 atmospheric pressure gradient trend not taken into account (see dynspg.F90 ''' To be done !!!!''' [[BR]] ''' 135 136 kpp non-local trend put in zdf trends !!! this will not work ! a additional trend term should be add ''' To be done !!!!''' [[BR]] ''' 143 137 144 138 problems to be solved: vvl case for tracer sad trends ; flux form case for had (keg) and zad momentum trends[[BR]] … … 148 142 create the momentum diag over the ML 149 143 150 reshape trdtra so that T and S are treated separately in all ted routine (including mld diag...) create umask_i and vmask_i (2D) fields that mirror task_i field but for the velocity points ''''' ''''''144 reshape trdtra so that T and S are treated separately in all ted routine (including mld diag...) create umask_i and vmask_i (2D) fields that mirror task_i field but for the velocity points ''''' ''''' ' 151 145 152 146 ----