* Add parallel capability for ce0l.
* Small bug in grid_noro fixed (smoothed topography was used instead of unsmoothed one for geopotential computation at north pole).
* Removed average of mass at poles in etat0dyn_netcdf after start_init_dyn => different results in the zoomed grid case.
* ok_etat0=n and ok_limit=y combination now works fine (if no initial state is needed, but only limit.nc file). This required:
- to move grid_noro0 and start_init_noro0 subroutines from etat0dyn_netcdf.F90 to limit_netcdf.F90
- to create init_ssrf_m.F90 file, so that sub-surfaces can be initialized from limit_netcdf.F90 without any etat0*_netcdf routines call).
* Simplified somehow the corresponding code, in particular:
1) removed obsolete flags "oldice".
2) removed flag "ibar": barycentric interpolation is used everywhere (except in start_init_subsurf, still calling grille_m - to be changed soon).
3) removed useless CPP_PHY precompilation directives, considering the possibility to run ce0l without physics is useless (ce0l is dedicated to Earth physics).
1412 lines of code changed in 40 files:
Fixed regression from revision 2333. Add argument rot in call to
physiq from lmdz1d.
7 lines of code changed in 1 file:
Bug fix. maase, pbaru and pbarv cannot be intent(in) in caldyn0
because they are modified in caldyn0 (masse is computed by massdair
and pbaru, pbarv are computed by flumass). Also masse is already
computed in etat0dyn_netcdf, and the difference with the computation
in caldyn0 is the averaging at the poles. If we do not want to lose
this averaging, we should not overwrite masse with an output argument
from caldyn0. So the simplest thing to do is to downgrade masse from
argument to local variable of caldyn0. We also simply downgrade pbaru,
pbarv from arguments to local variables of caldyn0.
6 lines of code changed in 2 files:
New parameterization of gravity wave drag due to front/jet systems, by
A. de la Camara and F. Lott. The new Camara-Lott parameterization
replaces the Hines parameterization so it is activated if not ok_hines
and ok_gwd_rando.
Also changed distribution of phase speeds in FLOTT_GWD_rando, from
uniform to Gaussian. Bug fix in sugwd_strato. Bug fix in the arguments
of the call to add_phys_tend for methane oxydation.
For the new Camara-Lott parameterization, we need to compute relative
vorticity in calfis and pass it as a new argument "rot" to
physiq. Interpolation of relative vorticity to the physics grid is not
optimal for now: it is not weighted by cell areas.
Alvaro de la Camara, Fran\c{}cois Lott
802 lines of code changed in 15 files:
Corrections pour recuperer les cas 1D bomex et cindy.
215 lines of code changed in 6 files:
Fixed regression from revision 2315: comvert.h was replaced by
vertical_layers_mod in test_disvert, but variables ap, bp, preff of
vertical_layers_mod were not defined. So, in main program ce0l, moved
call to test_disvert after call to Init_Phys_lmdz, and inserted in
between them calls to infotrac_init and iniphysiq (required). Had then
to remove the call to infotrac_init in etat0dyn_netcdf. In main
program ce0l, had to remove the call to InitComgeomphy since this is
done in iniphysiq.
In main program ce0l: no need to use indice_sol_mod; removed
preprocessor tests on CPP_PHYS in ce0l.
133 lines of code changed in 3 files:
That one as well grrrrrr
27 lines of code changed in 2 files:
Needed here finally
17 lines of code changed in 2 files:
Special option for Dice 1D case with a new value of iflag_radia=2.
It allows to force LWUP computation via ground temperature
MPL
117 lines of code changed in 6 files:
Changes in output levels of some variables and
Fixes of units for Ale and Alp variables.
Cleaning of variable declarations in cv3_yield.
70 lines of code changed in 3 files:
Further code reorganization: adding "phy_common", a directory which should contain routines common (wrt structural nature of the underlying code/grid) to all LMDZ-related physics packages.
EM
5461 lines of code changed in 29 files:
Enabling aerosols for the case RRTM + NSW=2
In this case we use the old aerosol properties
32 lines of code changed in 2 files:
Tidying up the old aerosol optical properties
especially use naero_tot dimension wherever needed
78 lines of code changed in 3 files:
add seasalt accumulation mode in mass_solu_aero calcul
it will make some difference in ok_aie case
(Modification by Yves B., and Olivier B.)
4 lines of code changed in 2 files:
Bug fix in alboc_cd (more than 11 years old). rmu0 must not be
modified in alboc_cd. The corresponding actual argument in subroutine
surf_ocean is an intent(in) argument of surf_ocean. It probably worked
before because the compiler made a copy of the argument but I would
say the behavior was in principle compiler-dependent and
optimization-level-dependent.
148 lines of code changed in 2 files:
Create the "obsolete" directory where old and unused stuff should go. And move some obsolete routines there.
Minor correction in "print_control_mod": use "opened" argument to inquire instead of "exist".
EM
907 lines of code changed in 15 files:
Physics/dynamics separation: make an infotrac_phy module, which should be used from within the physics, and is initialized from infotrac (dynamics) via iniphysiq.
EM
300 lines of code changed in 36 files: