* Fix for unadvected tracers (iadv==0)
* The key %isH2Ofamily, from the derived type "trac_type", is replaced with the more general
key %isInPhysics, which is TRUE for tracers both in "qx" and "tr_seri".
Currently, FALSE for tracers descending on H2O (isotopes and tagging tracers included). Could be set to FALSE
for interactive CO2 (type_trac=='inco') or ice supersaturated cloud content (tranfered to "rneb_seri")
* Parser routines taken from version 10 of https://svn.lmd.jussieu.fr/tracers-parser
41 lines of code changed in 20 files:
Some diagnostics are subject to the flag being activated.
23 lines of code changed in 1 file:
Fix in strings_mod ; parser routines taken from version 9 of https://svn.lmd.jussieu.fr/tracers-parser
6 lines of code changed in 2 files:
FIx for last version: function version of strReduce crashes with gfortran (but not for ifort) => back to subroutine version.
22 lines of code changed in 2 files:
Fixes mainly for isotopes (more to be done).
Fix (to be confirmed) in physiq to avoid attempting to send a non-transported (iadv==0) tracer to the physics.
79 lines of code changed in 4 files:
Added missing initialisation
LF
2 lines of code changed in 1 file:
Corrections to r4059/r4062 so that phylmdiso compiles
LF
1034 lines of code changed in 4 files:
* minor fixes (unused variables suppressed, comas after a WRITE() statement, etc.)
* parser routines taken from version 7 of https://svn.lmd.jussieu.fr/tracers-parser
* few changes in infotrac, and few fixes of (at least) the sequential version:
- uadv and vadv were deallocated twice (fix was lost by mistake just before last commit)
- in [( dum(im), im=1, nm)] implicit loops, ifort evaluates "dum(im)" even if nm==0,
resulting in a crash, "im" being unitialized.
66 lines of code changed in 6 files:
Except if the bench fails, last commit before switching to the new parser.
* parser routines taken from version 6 of https://svn.lmd.jussieu.fr/tracers-parser
* significant changes in "infotrac", now close to its target version.
Modified old code is active (hard-coded flag: lOldCode=.TRUE.).
New code will be activated in the next commit, then old code will be suppressed.
* probably problems with the isotopic version, to be fixed after tgcc maintenance...
973 lines of code changed in 8 files:
more keys for plane effects to avoid pb on outputs
65 lines of code changed in 6 files:
petit oubli
1 lines of code changed in 1 file:
add changes made to physiq_mod in phylmd
25 lines of code changed in 1 file:
Audran Borella's parametrisation for ice supersaturation
activated with flag_ice_sursat (FALSE by default)
1208 lines of code changed in 15 files:
Fix the call to vlspltgen_loc in advtrac_loc (one wrong additional calling argument).
7 lines of code changed in 1 file:
Temporarly removing the test on nqtottr.
Will be included later, when tracers(:)%iso_iName will be defined.
2 lines of code changed in 1 file:
Most of the changes are intended to help to eventually remove the constraints about the tracers assumptions, in particular water tracers.
- Remove index tables itr_indice and niadv, replaced by tracers(:)%isAdvected and tracers(:)%isH2OFamily.
Most of the loops are now from 1 to nqtot:
* DO iq=nqo+1,nqtot loops are replaced with:
DO iq=1,nqtot
IF(tracers(iq)%isH2Ofamily) CYCLE
* DO it=1,nbtr; iq=niadv(it+nqo)
and DO it=1,nqtottr; iq=itr_indice(it) loops are replaced with:
it = 0
DO iq = 1, nqtot
IF(.NOT.tracers(iq)%isAdvected .OR. tracers(iq)%isH2Ofamily) CYCLE
it = it+1
- Move some StratAer related code from infotrac to infotrac_phy
- Remove "nqperes" variable:
DO iq=1,nqpere loops are replaced with:
DO iq=1,nqtot
IF(tracers(iq)%parent/='air') CYCLE
- Cosmetic changes (justification, SELECT CASE instead of multiple IF...) mostly in advtrac* routines.
1454 lines of code changed in 34 files:
Check that the value of `iflag_top_bound` is valid
This check is useful because if `iflag_top_bound` is not between 0 and
2 then the program runs with an undefined value of lambda in procedure
`top_bound_loc`.
4 lines of code changed in 2 files:
Only compute pk if we are going to need it
In revsion revision 4042, wa added the computation of pk for
guide_plevs == 1. We can restrict a little better the cases where we
need pk.
10 lines of code changed in 1 file: