r2291-r2396
------------------------------------------------------------------------
r2396 | fairhead | 2015-11-18 15:12:46 +0100 (Wed, 18 Nov 2015) | 4 lines
Modifications from 2274 that were somehow lost
LF
------------------------------------------------------------------------
r2395 | emillour | 2015-11-18 13:25:35 +0100 (Wed, 18 Nov 2015) | 6 lines
Bug fix for aquaplanets: "rlat" and "rlon" were uninitialized when written to startphy.nc; phyredem
should write "longitude_deg" and "latitude_deg" to the file. Overall "rlat" and "rlon" should not b
e used and "latitude_deg" and "longitude_deg" should be used instead in the physics. To be further c
leaned up.
Some cleanup on that mater also made in phydev.
Note that this change will make bench test results different for longitudes and latitudes, because o
f roundoff effects.
EM
------------------------------------------------------------------------
r2394 | acozic | 2015-11-18 12:41:49 +0100 (Wed, 18 Nov 2015) | 7 lines
to fit with inca we need to increase dimension of naero_grp in aero_mod.
When lmdz is running without inca dim 10 to 12 will be empty and there is no change in results
these modifications only impact the old physics coupled with inca
------------------------------------------------------------------------
r2393 | jyg | 2015-11-18 12:25:20 +0100 (Wed, 18 Nov 2015) | 6 lines
Add various intializations of arrays in lmdz1d.F90
and in the convection scheme. Add output variables
for boundary layer splitting.
------------------------------------------------------------------------
r2392 | fhourdin | 2015-11-11 19:42:11 +0100 (Wed, 11 Nov 2015) | 4 lines
Retour a 1+1=2. Probleme dans thermcell_alim.
Bug fixing
------------------------------------------------------------------------
r2391 | emillour | 2015-11-10 09:51:52 +0100 (Tue, 10 Nov 2015) | 10 lines
Fix some minor anomalies spotted by the xlf compiler:
- infotrac : wrongly giving integer values to logicals
- surf_ocean_mod and yamada_c : should be a space between "stop" and message; but avoid using stop,
use abort_physic routine instead.
- readchlorophyll: using isnan() is not standard; compare the variable to itself instead (will retur
n .false. if NaN).
While at it, also added some missing "only" clauses when using modules.
EM
------------------------------------------------------------------------
r2390 | musat | 2015-11-09 16:35:21 +0100 (Mon, 09 Nov 2015) | 2 lines
Correction noms variables diagnostiques ptstar, pt0
------------------------------------------------------------------------
r2389 | musat | 2015-11-09 16:31:08 +0100 (Mon, 09 Nov 2015) | 2 lines
Ajout variables ptstar, pt0 utilisees dans le diagnostic de la slp
------------------------------------------------------------------------
r2388 | fhourdin | 2015-11-07 10:37:38 +0100 (Sat, 07 Nov 2015) | 3 lines
Ajout de la routine thermcell_alim calculant l'alimentation laterale.
------------------------------------------------------------------------
r2387 | fhourdin | 2015-11-07 10:27:40 +0100 (Sat, 07 Nov 2015) | 7 lines
Poursuite du travail sur les thermiques.
Extraction d'une routine calculant l'alimentation laterale à la base.
+ une correction de Jean-Yves sur des initialisations maquantes dans
thermcell_main (mais qui sont maintenant dans thermcell_alp).
------------------------------------------------------------------------
r2386 | musat | 2015-11-06 18:31:30 +0100 (Fri, 06 Nov 2015) | 2 lines
Ajout routines calcul slp
------------------------------------------------------------------------
r2385 | musat | 2015-11-06 18:30:34 +0100 (Fri, 06 Nov 2015) | 3 lines
Correction calcul sea level pressure (slp) cf. Arpege-IFS
Ajout initialisation cldfrarad
------------------------------------------------------------------------
r2384 | fhourdin | 2015-11-05 18:49:55 +0100 (Thu, 05 Nov 2015) | 4 lines
Séparation de la partie de calcul de ALP dans thermcell_main
Spliting thermcell_main.F90 in 2, isolating the ALP computation.
------------------------------------------------------------------------
r2383 | jbmadeleine | 2015-11-02 14:53:48 +0100 (Mon, 02 Nov 2015) | 1 line
Modified iniphysiq_mod.F90 that includes nqtot
------------------------------------------------------------------------
r2382 | jbmadeleine | 2015-11-02 14:52:09 +0100 (Mon, 02 Nov 2015) | 6 lines
>>> Updated phymar routines and made sure it compiles
- Changed X___AT in Mod_PHY_AT_kkl.f90 from REAL*16 (obsolete) to double precision
- The number of tracers nqtot is now given to MAR via the infotrac_phy routine
- Changed iim and jjm into ii and jj in dynlonlat_phylonlat/phymar/iniphysiq_mod.F90
- Added the declaration of the rotational in phymar/physiq.F90
------------------------------------------------------------------------
r2380 | musat | 2015-10-27 16:59:53 +0100 (Tue, 27 Oct 2015) | 4 lines
Bug fix for 1D cases concerning outputs on NMC levels when running with
instal1D.sh
IM
------------------------------------------------------------------------
r2379 | musat | 2015-10-23 14:21:06 +0200 (Fri, 23 Oct 2015) | 3 lines
2nd Bug fix for monthly outputs on NMC levels when running with install.sh
IM
------------------------------------------------------------------------
r2378 | musat | 2015-10-22 17:58:25 +0200 (Thu, 22 Oct 2015) | 3 lines
Bug fix for monthly outputs on NMC levels when running with install.sh
IM
------------------------------------------------------------------------
r2377 | musat | 2015-10-20 15:48:18 +0200 (Tue, 20 Oct 2015) | 7 lines
Correct zfull and zhalf CMIP6' diagnostics :
zhalf correspond to the interface levels (and paprs) and include
both the top of the model atmosphere and surface levels
zfull is the actual height above mean sea level and corresponds to the
geopotential height (and pplay levels)
ASima/AI/IM
------------------------------------------------------------------------
r2376 | jyg | 2015-10-20 11:57:39 +0200 (Tue, 20 Oct 2015) | 4 lines
Add initialization of ice fraction below cloud
base in cv3_undilute2 in cv3_routines.F90.
------------------------------------------------------------------------
r2375 | emillour | 2015-10-18 08:38:58 +0200 (Sun, 18 Oct 2015) | 5 lines
Fix in the computation of the date; the convention is that it corresponds to the time at the end of
the current dynamics or physics step (except when in backward Matsuno step where it remains unchange
d as it was correctly updated during the forward part of the step).
Note that the relationship between itau and date is a bit tricky as itau is incremented between Mats
uno forward and backward steps (and from a leapfrog step to the next) but unchanged from Matsuno bac
ward step to leapfrog step.
Because of change in date when calling physics, bench results will change with this revision.
EM
------------------------------------------------------------------------
r2374 | jyg | 2015-10-13 19:43:21 +0200 (Tue, 13 Oct 2015) | 8 lines
Creation of a new closure routine for Emanuel
convective scheme: cv3p2_closure.F90 (called when
iflag_clos=3) is a cleaned up and reordered
version of cv3p1_closure.F90. cv3p1_closure.F90
(called when iflag_clos=2) is kept for numerical
compatibility with earlier versions.
------------------------------------------------------------------------
r2373 | jyg | 2015-10-13 19:28:01 +0200 (Tue, 13 Oct 2015) | 5 lines
From MPL:
Correction of some missing initializations in
dyn1D and of two 1D case bugs.
------------------------------------------------------------------------
r2372 | acozic | 2015-10-13 14:41:47 +0200 (Tue, 13 Oct 2015) | 3 lines
Change call to inca initialisation to fit with new sections dynamique/physic
------------------------------------------------------------------------
r2371 | fairhead | 2015-10-13 12:14:39 +0200 (Tue, 13 Oct 2015) | 2 lines
Bug correction
------------------------------------------------------------------------
r2370 | fairhead | 2015-09-28 15:14:06 +0200 (Mon, 28 Sep 2015) | 2 lines
Rajout de la propriete Id pour un meilleur suivi des versions des physiq.def
------------------------------------------------------------------------
r2369 | oboucher | 2015-09-23 22:39:13 +0200 (Wed, 23 Sep 2015) | 3 lines
A correction on recent changes to rmu0 calculation
radpas was hard-coded in one place
------------------------------------------------------------------------
r2368 | oboucher | 2015-09-22 17:26:15 +0200 (Tue, 22 Sep 2015) | 3 lines
Removing radsol from output list of radlwsw_inca.F90
Anne Cozic will synchronize INCA routine accordingly
------------------------------------------------------------------------
r2367 | musat | 2015-09-22 11:02:39 +0200 (Tue, 22 Sep 2015) | 4 lines
Correction sur les sortiesi Cosp.
(Remise a zero des valeurs non definies)
------------------------------------------------------------------------
r2366 | oboucher | 2015-09-21 22:41:04 +0200 (Mon, 21 Sep 2015) | 4 lines
Second batch of changes for diagnosing SW radiative
diagnostics at the physics timestep resolution
radsol is now computed in physiq.F90
------------------------------------------------------------------------
r2365 | emillour | 2015-09-17 13:19:54 +0200 (Thu, 17 Sep 2015) | 3 lines
Arch files for Occigen; possibly to be improved later.
EM
------------------------------------------------------------------------
r2364 | jyg | 2015-09-17 13:00:25 +0200 (Thu, 17 Sep 2015) | 2 lines
Commenting out a spurious print in cv3_routines.F90
------------------------------------------------------------------------
r2363 | acozic | 2015-09-16 16:44:59 +0200 (Wed, 16 Sep 2015) | 3 lines
make some modification to allows compilation with inca model
------------------------------------------------------------------------
r2362 | acozic | 2015-09-16 14:50:44 +0200 (Wed, 16 Sep 2015) | 7 lines
several correction to use LMDZ coupled with INCA
1- add an "use" in gcm to define "klon_glo" variable
2- correct an error on a test in infotrac
3- add parents to inca's tracer in infotrac
4- change calendar name for "leap calendar" in wxios to fit with rev 2229 on gcm
------------------------------------------------------------------------
r2361 | dcugnet | 2015-09-07 17:50:29 +0200 (Mon, 07 Sep 2015) | 5 lines
In etat0dyn: removed few useless lines: "masque" is always known because etat0dyn is called after e
tat0phys.
In grid_atob: shape error in routine fine2coarse now fixed: "msk" argument and local variable mask m
ust have the dimensions of the output grid. Working unit of dist_sphe is no longer degree, but radia
n.
------------------------------------------------------------------------
r2360 | oboucher | 2015-09-05 12:06:39 +0200 (Sat, 05 Sep 2015) | 4 lines
Small structural changes to rmu0 calculations
Now the ocean albedo sees the high-resolution rmu0
if iflag_cycle_diurne=2
------------------------------------------------------------------------
r2359 | oboucher | 2015-09-04 18:48:31 +0200 (Fri, 04 Sep 2015) | 9 lines
Recoded diurnal cycle in physiq.F90
iflag_diurnal_cycle is the new flag
=0 for no diurnal cycle
=1 for old diurnal cycle
=2 for new diurnal cycle with updates
in between two calls to radiation according
to how rmu0 has changed over the previous
timestep of physiq
------------------------------------------------------------------------
r2358 | emillour | 2015-09-04 09:10:16 +0200 (Fri, 04 Sep 2015) | 4 lines
- Correction of xml files to follow up on the new variables added in rev 2333.
- Added option '-libphy' to makelmdz_fcm to generate a given physics package library (without compil
ing dynamics or a main program).
EM
------------------------------------------------------------------------
r2357 | lguez | 2015-08-31 18:25:19 +0200 (Mon, 31 Aug 2015) | 11 lines
The default value for gwd_rando_ruwmax is changed. This is ok because
the Lott GWD rando parameterization has not been much used yet.
New parameters GWD_FRONT_RUWMAX, GWD_FRONT_SAT, that can be chosen at
run-time, for frontal gravity waves.
New paramters sso_gkdrag, sso_grahil, sso_grcrit, sso_gfrcri,
sso_gkwake, sso_gklift that can be chosen at run-time, for orographic
gravity waves. The default values for those parameters are those that
were hard-coded before this revision.
------------------------------------------------------------------------
r2356 | jghattas | 2015-08-31 15:48:29 +0200 (Mon, 31 Aug 2015) | 3 lines
Corrected bug: sollw was written twice to the restart file and sollwdown was never written.
------------------------------------------------------------------------
r2355 | emillour | 2015-08-28 13:04:00 +0200 (Fri, 28 Aug 2015) | 3 lines
Bug fix on bug introduced in previous bug fix (klon_glo is not set before call to init_physics_distr
ibution, but has a correct value after).
EM
------------------------------------------------------------------------
r2354 | emillour | 2015-08-28 11:38:53 +0200 (Fri, 28 Aug 2015) | 3 lines
Corrections for the 1D model to fix initializations broken with previous revision.
EM+MPL
------------------------------------------------------------------------
r2353 | emillour | 2015-08-26 14:47:43 +0200 (Wed, 26 Aug 2015) | 3 lines
Correction for ce0l to compile and work in serial mode.
EM
------------------------------------------------------------------------
r2352 | emillour | 2015-08-26 08:08:40 +0200 (Wed, 26 Aug 2015) | 3 lines
Fix a circular dependency introduced in previous commit: occurs when routines in phy_common use modu
le print_control. Make them 'include "iniprint.h"' instead (not an ideal solution but will have to d
o for now), and therfore move iniprint.h to "misc" directory.
EM
------------------------------------------------------------------------
r2351 | emillour | 2015-08-25 17:14:59 +0200 (Tue, 25 Aug 2015) | 7 lines
More on physics/dynamics separation and cleanup:
- Set things up so that all physics-related initializations are done via iniphysiq.
- Created a "geometry_mod.F90" module in phy_common to store information on the loacl grid (i.e. rep
laces comgeomphy) and moreover give these variables more obvious names (e.g.: rlond => longitude, rl
atd => latitude, airephy => cell_area).
- removed obsolete comgeomphy.h and comgeomphy.F90
EM
------------------------------------------------------------------------
r2350 | emillour | 2015-08-25 13:40:19 +0200 (Tue, 25 Aug 2015) | 3 lines
Corrections for 1D case where nbp_lon=nbp_lat=1.
EM + MPL
------------------------------------------------------------------------
r2349 | emillour | 2015-08-25 07:51:24 +0200 (Tue, 25 Aug 2015) | 3 lines
Fix call to iniphysiq in ce0l.
EM
------------------------------------------------------------------------
r2348 | emillour | 2015-08-24 14:41:28 +0200 (Mon, 24 Aug 2015) | 3 lines
Added missing declarations required when CPP_VEGET preprocessing flag is defined (i.e. for LMDZOR).
EM
------------------------------------------------------------------------
r2347 | emillour | 2015-08-24 09:26:48 +0200 (Mon, 24 Aug 2015) | 4 lines
Make iniphysiq a module.
Fix call to iniphysiq in lmdz1d (missing arguments and arrays of wrong sizes).
EM
------------------------------------------------------------------------
r2346 | emillour | 2015-08-21 17:13:46 +0200 (Fri, 21 Aug 2015) | 6 lines
Physics/dynamics separation:
- remove all references to dimensions.h from physics. nbp_lon (==iim) , nbp_lat (==jjm+1) and nbp_le
v (==llm) from mod_grid_phy_lmdz should be used instead.
- added module regular_lonlat_mod in phy_common to store information about the global (lon-lat) grid
cell boundaries and centers.
EM
------------------------------------------------------------------------
r2345 | emillour | 2015-08-21 11:57:36 +0200 (Fri, 21 Aug 2015) | 6 lines
Physics/dynamics separation:
- move test_disvert_m to dynlonlat_phylonlat/phylmd since it is only used by ce0l and relies on dyna
mics.
- put "config_inca" in tracinca_mod so physics routines can get the info from there rather than from
control_mod.
- get rid of references to "control_mod" from within the physics.
EM
------------------------------------------------------------------------
r2344 | emillour | 2015-08-21 09:23:13 +0200 (Fri, 21 Aug 2015) | 3 lines
Physics/dynamics separation: get rid of all the 'include "temps.h"' in the physics; variables in mod
ule time_phylmdz_mod must be used instead. Also added JD_cur, JH_cur and JD_ref in module phys_cal_m
od, in preparation for having physics handle its calendar internally.
EM
------------------------------------------------------------------------
r2343 | emillour | 2015-08-20 12:02:53 +0200 (Thu, 20 Aug 2015) | 4 lines
Another step towards a clean separation between physics and dynamics: adapted read_pstoke.F90, read_
pstoke0.F90 initphysto.F90 and phystokenc.F90 (now module phystokenc_mod.F90) to not explicitely inc
lude/use file/modules from the dynamics.
In the process, added module "time_phylmdz_mod.F90" in the physics, which contains the information o
therwise found in "temps.h" (which is in the dynamics) and should be used instead.
EM
------------------------------------------------------------------------
r2342 | emillour | 2015-08-19 15:21:38 +0200 (Wed, 19 Aug 2015) | 7 lines
Follow-up for a cleaner separation between dynamics and physics:
- "write_field" is called from physics and dynamics but has no dependence on either so it should be
in "misc"
- "write_field_phy" is common to all physics, so it goes in "phy_common"
- "init_ssrf_m" and "limit_netcdf" are only used by ce0l, so these should be in "dynlonlat_phylonlat
/phylmd"
- "q_sat" is called from physics and dynamics but has no dependence on either so it should be in "misc"
EM
------------------------------------------------------------------------
r2341 | musat | 2015-08-19 14:56:45 +0200 (Wed, 19 Aug 2015) | 5 lines
Bug fix : Add !$OMP THREADPRIVATE for
zustr_gwd_hines, zvstr_gwd_hines
zustr_gwd_front, zvstr_gwd_front
IM
------------------------------------------------------------------------
r2340 | emillour | 2015-08-18 17:27:50 +0200 (Tue, 18 Aug 2015) | 3 lines
Bug fix: Adding the tendencies from wakes to dry adjustement tendencies should only be done if wakes
are computed.
EM + JYG
------------------------------------------------------------------------
r2339 | musat | 2015-08-18 17:08:09 +0200 (Tue, 18 Aug 2015) | 3 lines
Bug on l initialization for debug mode
IM
------------------------------------------------------------------------
r2338 | dcugnet | 2015-08-05 17:37:29 +0200 (Wed, 05 Aug 2015) | 1 line
------------------------------------------------------------------------
r2337 | dcugnet | 2015-08-03 13:06:55 +0200 (Mon, 03 Aug 2015) | 1 line
------------------------------------------------------------------------
r2336 | dcugnet | 2015-07-31 19:22:21 +0200 (Fri, 31 Jul 2015) | 12 lines
* Add parallel capability for ce0l.
* Small bug in grid_noro fixed (smoothed topography was used instead of unsmoothed one for geopote
ntial 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 li
mit.nc file). This required:
- to move grid_noro0 and start_init_noro0 subroutines from etat0dyn_netcdf.F90 to limit_netcdf.F
90
- 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_subsu
rf, still calling grille_m - to be changed soon).
3) removed useless CPP_PHY precompilation directives, considering the possibility to run ce0l wi
thout physics is useless (ce0l is dedicated to Earth physics).
------------------------------------------------------------------------
r2335 | lguez | 2015-07-28 16:01:58 +0200 (Tue, 28 Jul 2015) | 3 lines
Fixed regression from revision 2333. Add argument rot in call to
physiq from lmdz1d.
------------------------------------------------------------------------
r2334 | lguez | 2015-07-24 19:32:31 +0200 (Fri, 24 Jul 2015) | 10 lines
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.
------------------------------------------------------------------------
r2333 | lguez | 2015-07-22 16:14:51 +0200 (Wed, 22 Jul 2015) | 16 lines
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
------------------------------------------------------------------------
r2332 | fhourdin | 2015-07-21 17:40:58 +0200 (Tue, 21 Jul 2015) | 3 lines
Corrections pour recuperer les cas 1D bomex et cindy.
------------------------------------------------------------------------
r2331 | lguez | 2015-07-17 14:17:02 +0200 (Fri, 17 Jul 2015) | 12 lines
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.
------------------------------------------------------------------------
r2330 | fairhead | 2015-07-16 16:58:30 +0200 (Thu, 16 Jul 2015) | 2 lines
That one as well grrrrrr
------------------------------------------------------------------------
r2329 | fairhead | 2015-07-16 16:18:16 +0200 (Thu, 16 Jul 2015) | 2 lines
Needed here finally
------------------------------------------------------------------------
r2328 | fairhead | 2015-07-16 14:09:53 +0200 (Thu, 16 Jul 2015) | 5 lines
Special option for Dice 1D case with a new value of iflag_radia=2.
It allows to force LWUP computation via ground temperature
MPL
------------------------------------------------------------------------
r2327 | jyg | 2015-07-11 11:19:49 +0200 (Sat, 11 Jul 2015) | 5 lines
Changes in output levels of some variables and
Fixes of units for Ale and Alp variables.
Cleaning of variable declarations in cv3_yield.
------------------------------------------------------------------------
r2326 | emillour | 2015-07-10 14:24:29 +0200 (Fri, 10 Jul 2015) | 3 lines
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
------------------------------------------------------------------------
r2325 | oboucher | 2015-07-08 22:17:34 +0200 (Wed, 08 Jul 2015) | 3 lines
Enabling aerosols for the case RRTM + NSW=2
In this case we use the old aerosol properties
------------------------------------------------------------------------
r2324 | oboucher | 2015-07-08 17:20:22 +0200 (Wed, 08 Jul 2015) | 3 lines
Tidying up the old aerosol optical properties
especially use naero_tot dimension wherever needed
------------------------------------------------------------------------
r2323 | acozic | 2015-07-08 10:02:43 +0200 (Wed, 08 Jul 2015) | 5 lines
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.)
------------------------------------------------------------------------
r2322 | lguez | 2015-07-07 18:01:01 +0200 (Tue, 07 Jul 2015) | 7 lines
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.
------------------------------------------------------------------------
r2321 | emillour | 2015-07-06 18:22:14 +0200 (Mon, 06 Jul 2015) | 4 lines
Create the "obsolete" directory where old and unused stuff should go. And move some obsolete routine
s there.
Minor correction in "print_control_mod": use "opened" argument to inquire instead of "exist".
EM
------------------------------------------------------------------------
r2320 | emillour | 2015-07-01 15:57:32 +0200 (Wed, 01 Jul 2015) | 3 lines
Physics/dynamics separation: make an infotrac_phy module, which should be used from within the physi
cs, and is initialized from infotrac (dynamics) via iniphysiq.
EM
------------------------------------------------------------------------
r2319 | fairhead | 2015-06-30 13:49:22 +0200 (Tue, 30 Jun 2015) | 2 lines
Some initialisation missing
------------------------------------------------------------------------
r2318 | jyg | 2015-06-26 09:41:13 +0200 (Fri, 26 Jun 2015) | 4 lines
Bug fixing concerning the adiabatic adjustment of
the wake region (svn2309).
------------------------------------------------------------------------
r2317 | jyg | 2015-06-26 09:33:22 +0200 (Fri, 26 Jun 2015) | 5 lines
Bug fixing after physics/dynamic separation
enforcing in svn 2311 (corrections in
readaerosol.F90 and in dyn1d/lmdz1d.F90).
------------------------------------------------------------------------
r2315 | emillour | 2015-06-25 17:25:50 +0200 (Thu, 25 Jun 2015) | 3 lines
More on physics/dynamics separation: make a vertical_layers_mod module to contain information on the
vertical discretization. This module should be used from within the physics (instead of including c
omvert.h from dynamics).
EM
------------------------------------------------------------------------
r2313 | musat | 2015-06-25 14:59:12 +0200 (Thu, 25 Jun 2015) | 3 lines
Suite corrections missing_value pour compilation en mode debug
JG/IM
------------------------------------------------------------------------
r2312 | musat | 2015-06-25 14:56:33 +0200 (Thu, 25 Jun 2015) | 1 line
------------------------------------------------------------------------
r2311 | emillour | 2015-06-25 09:45:24 +0200 (Thu, 25 Jun 2015) | 8 lines
Further modifications to enforce physics/dynamics separation:
- moved iniprint.h and misc_mod back to dyn3d_common, as these should only be used by dynamics.
- created print_control_mod in the physics to store flags prt_level, lunout, debug to be local to ph
ysics (should be used rather than iniprint.h)
- created abort_physic.F90 , which does the same job as abort_gcm() did, but should be used instead
when in physics.
- reactivated inifis (turned it into a module, inifis_mod.F90) to initialize physical constants and
print_control_mod flags.
EM
------------------------------------------------------------------------
r2310 | fairhead | 2015-06-24 12:04:31 +0200 (Wed, 24 Jun 2015) | 3 lines
Modifications for 1d model
MPL
------------------------------------------------------------------------
r2309 | jyg | 2015-06-23 18:08:54 +0200 (Tue, 23 Jun 2015) | 7 lines
Modification of physiq.F90: when wakes are
present, an adiabatic adjustment is applied to the
wake region before calling the Emanuel convection
scheme. The switch 'ok_adjwk', read from .def
files, controls this adjustment.
------------------------------------------------------------------------
r2308 | jyg | 2015-06-23 13:23:49 +0200 (Tue, 23 Jun 2015) | 4 lines
Change to wake.F90 : when there are no wake,
the wake area and the wake height are set to zero
when leaving wake.F90.
------------------------------------------------------------------------
r2307 | fairhead | 2015-06-23 13:09:31 +0200 (Tue, 23 Jun 2015) | 3 lines
Modifications to run 1D Bomex case
MPL
------------------------------------------------------------------------
r2306 | jyg | 2015-06-23 11:33:09 +0200 (Tue, 23 Jun 2015) | 9 lines
Improved diagnostics: pmflxr and pmflxs are now
the true vertical profiles of liquid and solid
convective precipitation (previously they where
merely diagnosed through a test on temperature).
The convective scheme internal variable for ice
precipitation is Vprecipi.
------------------------------------------------------------------------
r2305 | jyg | 2015-06-22 23:25:33 +0200 (Mon, 22 Jun 2015) | 3 lines
bug fixing in conf_phys_m.F90 concerning the
initialization of iflag_albedo and ok_chlorophyll.
------------------------------------------------------------------------
r2302 | emillour | 2015-06-17 14:46:36 +0200 (Wed, 17 Jun 2015) | 4 lines
Move etat0phys_netcdf.F90 to "dynlonlat_phylonlat/phylmd" as it relies on "phylmd" routines.
Some cleanup to remove obsolete and unecessary CPP_EARTH preprocessing condition.
EM
------------------------------------------------------------------------
r2301 | emillour | 2015-06-17 09:33:27 +0200 (Wed, 17 Jun 2015) | 5 lines
Moved ioipsl_getin_p_mod.F90 to "phylmd" beacause files in "misc" should rely on neither physics nor
dynamics modules (or include) files.
Moved etat0dyn_netcdf.F90 to "dynlonlat_phylonlat/phylmd" as it is only used by ce0l and moreover re
quires "grid_atob_m" which is in "dynlonlat_phylonlat".
With these adjustements, compiling with '-nophys' works again.
EM
------------------------------------------------------------------------
r2300 | fhourdin | 2015-06-16 12:20:35 +0200 (Tue, 16 Jun 2015) | 3 lines
Deplacement de fichiers sources pour la compilation 1D.
------------------------------------------------------------------------
r2299 | dcugnet | 2015-06-15 20:48:31 +0200 (Mon, 15 Jun 2015) | 33 lines
In dyn3d/:
etat0dyn_netcdf.F90: "startget_dyn3d" syntax slightly simplified.
dynredem.F90: Shortcut routines (put_var*, cre_var,
dynredem_write_*, dynredem_read_u)
modified to match dyn3dmem version and put in
module dyredem_mod.F90.
dynetat0.F90 -> *.f90: Few simplifications (no usage of NC_DOUBLE
needed => no precompilation)
Add tracers initialization in the isotope case
suppressed by accident.
dynredem_mod.F90: Created to mimic dyn3dmem equivalent.
In dyn3dmem/:
dynetat0_loc.F -> *.f90: Converted into fortran 90 to match the dyn3d
version.
dynredem_loc.F -> *.F90: Converted into fortran 90.
dynredem_mod.F90: Add some shortcut routines to match the dyn3d
version.
In phylmd/:
phyredem.F90: Bug fix: nsw instead of nsoilmx was used as
Tsoil second maximum index.
Bug fix: fevap instead of snow was saved for
"SNOW".
etat0phys_netcdf.F90: "filtreg_mod" module usage suppressed.
Local variable rugo computation removed (not
used).
In dynlonlat_phylonlat/:
grid_atob_m.F90 -> *.f90 DOUBLE PRECISION variables usage removed.
Precompilation o longer needed => .F90 extension.
------------------------------------------------------------------------
r2297 | musat | 2015-06-12 23:38:57 +0200 (Fri, 12 Jun 2015) | 10 lines
RRTM :
Correction bug interface avec rrtm (radlwsw_m.F90)
Ajouter la possibilite d'utilsation d'anciennes proprites optique dans le cas RRTM 2bandes (readaero
sol_optic_rrtm.F90)
Correction rrtm (rrtm_rtrn1a_140gp.F90)
COSP :
Modifications pour traiter les valeurs indefines :
mise a 0 en cas d'utilisation de IOIPSL
mise a la valeur recuperee dans .xml en cas d'utilisation de XIOS
------------------------------------------------------------------------
r2295 | fhourdin | 2015-06-12 06:44:36 +0200 (Fri, 12 Jun 2015) | 4 lines
Correction du seuil sur la probabilite de non declenchement.
Changement de priorite de sortie sur deux variables.
------------------------------------------------------------------------
r2294 | fhourdin | 2015-06-11 20:06:49 +0200 (Thu, 11 Jun 2015) | 12 lines
Introduction d'un seuil sur la probabilité de non déclenchement
stochastique random_notrig_max
Si le tirage aleatoire uniforme entre 0 et 1 est > random_notrig_max
on ne déclenche pas.
random_notrig_max=1-epsilon avec epsilon petit.
Corrections dans grid_atob_m.F90 pour la compilation Linux.
A vérifier.
------------------------------------------------------------------------
r2293 | dcugnet | 2015-06-05 21:16:07 +0200 (Fri, 05 Jun 2015) | 35 lines
Initial states creation routines have been reorganized and simplified.
As far as possible, dynamics and physics related routines have been
separated.
Some routines have been converted to fortran 90 and repeated codes sections
have been "factorized".
Array/vector arguments have become implicit in some routines to avoid usage
of "dimensions.h" ; possible for routines with explicit interfaces and if
iim and jjm can be deduced from arguments sizes.
* dynlonlat_phylonlat/ce0l.F90 calls now phylmd/etat0phys_netcdf.F90 and
dyn3d/etat0dyn_netcdf.F90 that replace phylmd/etat0_netcdf.F90.
start.nc and startphy.nc creations are now independant.
* startvar.F90 has been suppressed ; corresponding operations have been
simplified and embedded in etat0*_netcdf.F90 routines as internal
procedures.
* Routines converted to fortran 90 and "factorized":
- dyn3d_common/conf_dat_m.F90 (replaces dyn3d_common/conf_dat2d.F
and dyn3d_common/conf_dat3d.F)
- dyn3d/dynredem.F90 (replaces dyn3d/dynredem.F)
- dyn3d/dynetat0.F90 (replaces dyn3d/dynetat0.F)
- phylmd/grid_noro_m.F90 (replaces dyn3d_common/grid_noro.F)
- dynlonlat_phylonlat/grid_atob_m.F90 (replaces dyn3d_common/grid_atob.F)
- dyn3d_common/caldyn0.F90 (replaces dyn3d_common/caldyn0.F)
- dyn3d_common/covcont.F90 (replaces dyn3d_common/covcont.F)
- dyn3d_common/pression.F90 (replaces dyn3d_common/pression.F)
* phylmd/phyredem.F90 and phylmd/limit_netcdf.F90 have been slightly
factorized.
TO DO:
* little fix needed in grid_noro_m.F90 ; untouched yet to ensure results
are exactly the same as before. Unsmoothed orography is used to compute
"zphi", but smoothed (should be unsmoothed) one is used at poles.
* add the dyn3dmem versions of dynredem.F90 and dynetat0.F90
(dynredem_loc.F90 and dynetat0_loc.F90, untested yet).
* test compilation in parallel mode for a single processor.
------------------------------------------------------------------------
r2291 | emillour | 2015-06-03 12:59:18 +0200 (Wed, 03 Jun 2015) | 3 lines