r1921-r1997

Différences entre les releases 1921 et 1997 svn log -r1997:1921

------------------------------------------------------------------------
r1997 | acaubel | 2014-03-18 13:16:18 +0100 (Tue, 18 Mar 2014) | 3 lines

SAVE attribute was missing to run on OpenMP mode.


------------------------------------------------------------------------
r1996 | emillour | 2014-03-13 09:53:16 +0100 (Thu, 13 Mar 2014) | 4 lines

Remove old unused file, which has nothing to do there.
EM


------------------------------------------------------------------------
r1995 | fairhead | 2014-03-10 17:00:43 +0100 (Mon, 10 Mar 2014) | 6 lines

Utile uniquement pour les tests

-------------------

Needed only for test purposes

------------------------------------------------------------------------
r1994 | emillour | 2014-03-05 16:42:42 +0100 (Wed, 05 Mar 2014) | 4 lines

Follow the trend; convert all fixed form sources in phydev to free form.
EM


------------------------------------------------------------------------
r1993 | emillour | 2014-03-05 15:38:40 +0100 (Wed, 05 Mar 2014) | 3 lines

 OpenMP instruction syntax must be free form if source file is free form Fortran.
EM

------------------------------------------------------------------------
r1992 | lguez | 2014-03-05 14:19:12 +0100 (Wed, 05 Mar 2014) | 20 lines

Converted to free source form files in libf/phylmd which were still in
fixed source form. The conversion was done using the polish mode of
the NAG Fortran Compiler.

In addition to converting to free source form, the processing of the
files also:

-- indented the code (including comments); 

-- set Fortran keywords to uppercase, and set all other identifiers
to lower case;

-- added qualifiers to end statements (for example "end subroutine
conflx", instead of "end");

-- changed the terminating statements of all DO loops so that each
loop ends with an ENDDO statement (instead of a labeled continue).

-- replaced #include by include.

------------------------------------------------------------------------
r1991 | fhourdin | 2014-03-05 12:55:30 +0100 (Wed, 05 Mar 2014) | 4 lines

Correction cosmétique pour RRTM
RRTM cosmetics


------------------------------------------------------------------------
r1990 | fairhead | 2014-03-05 11:33:03 +0100 (Wed, 05 Mar 2014) | 11 lines

Corrections  à la version r1989 pour permettre la compilation avec RRTM
Inclusion de la licence CeCILL_V2 pour RRTM
--------

Changes to revision r1989 to enable RRTM code compilation
RRTM part put under CeCILL_V2 licence



------------------------------------------------------------------------
r1989 | fairhead | 2014-03-03 10:57:40 +0100 (Mon, 03 Mar 2014) | 9 lines

Inclusion du code RRTM

------------------------

Adding RRTM code

MPL


------------------------------------------------------------------------
r1988 | emillour | 2014-02-27 08:43:23 +0100 (Thu, 27 Feb 2014) | 3 lines

Bug fix: in dyn3dpar seting and updating of surface pressure hallo was missing in physics step (only matters for cases
 where physics modify surface pressure).
EM

------------------------------------------------------------------------
r1987 | emillour | 2014-02-24 16:05:47 +0100 (Mon, 24 Feb 2014) | 4 lines

Add updating pressure, mass and Exner function (ie: all variables which depend on surface pressure) after adding physi
cs tendencies (which include a surface pressure tendency).
Note that this change induces slight changes in GCM results with respect to previous svn version of the code, even if 
surface pressure tendency is zero (because of recomputation of polar values as an average over polar points on the dyn
amics grid).
EM

------------------------------------------------------------------------
r1986 | jghattas | 2014-02-19 15:51:00 +0100 (Wed, 19 Feb 2014) | 4 lines

cosp : corrections for compiling with XIOS. Runnig with cosp not tested.  
xml : corrected error in field ozone_daylight. Error seen when read_climoz=2.


------------------------------------------------------------------------
r1985 | fhourdin | 2014-02-19 06:53:24 +0100 (Wed, 19 Feb 2014) | 6 lines

Modification de du schéma implicite des thermique pour que la concentration
    soit inchangée en cas de flux de masse uniformément nul.
Implicit scheme for thermal plume modified to keep unchanged the concentration
    when the mass flux is zero.


------------------------------------------------------------------------
r1984 | emillour | 2014-02-18 10:59:29 +0100 (Tue, 18 Feb 2014) | 3 lines

Cleanup and update of "grilles_gcm_netcdf_sub": removed old version, removed flag "grilles_gcm_netcdf" used to trigger
 output of grilles_gcm.nc file (file will now always be created) and fixed writing vertical levels presnivs in the fil
e.
JBM+EM

------------------------------------------------------------------------
r1983 | fhourdin | 2014-02-17 19:30:42 +0100 (Mon, 17 Feb 2014) | 4 lines

Changement de place de pres2lev_mod.F90
Changing location for pres2lev_mod.F90

------------------------------------------------------------------------
r1982 | fhourdin | 2014-02-16 23:30:48 +0100 (Sun, 16 Feb 2014) | 3 lines

Suppression d'impressions


------------------------------------------------------------------------
r1981 | fhourdin | 2014-02-14 12:22:31 +0100 (Fri, 14 Feb 2014) | 3 lines

Correction
Bug fixing

------------------------------------------------------------------------
r1980 | fhourdin | 2014-02-14 12:03:45 +0100 (Fri, 14 Feb 2014) | 4 lines

Interface legere de sorties netcdf. Peut etre utilisee hors de LMDZ
Light interface for netcdf. Can be used outside LMDZ


------------------------------------------------------------------------
r1979 | fhourdin | 2014-02-14 11:58:58 +0100 (Fri, 14 Feb 2014) | 3 lines

Suprresion de inigrads de 1DUTILS.h


------------------------------------------------------------------------
r1978 | fhourdin | 2014-02-14 11:56:27 +0100 (Fri, 14 Feb 2014) | 3 lines

Suppression de l'appel à wrgradsfi
Removing calls to wrgradsfi

------------------------------------------------------------------------
r1977 | fhourdin | 2014-02-14 10:57:32 +0100 (Fri, 14 Feb 2014) | 3 lines

La modification de phys_output_mod était devenue obsolète
The phys_output_mod change was obsolete

------------------------------------------------------------------------
r1976 | fhourdin | 2014-02-14 10:44:20 +0100 (Fri, 14 Feb 2014) | 3 lines

Suppression du traitement de wrgradsfi
Removing the line concrening wrgradsfi

------------------------------------------------------------------------
r1975 | fhourdin | 2014-02-14 10:42:27 +0100 (Fri, 14 Feb 2014) | 4 lines

On enleve wrgradsfi de 1DUTILS.h. La routine est deja dans phylmd.
Removing wrgradsfi from 1DUTILS.h (already in phylmd)


------------------------------------------------------------------------
r1974 | fhourdin | 2014-02-14 10:37:34 +0100 (Fri, 14 Feb 2014) | 3 lines

Plus beson de modifier physiq.F90 lors de l'installation du 1D
No more need to modify physiq.F90 when isntalling the SCM

------------------------------------------------------------------------
r1973 | fhourdin | 2014-02-14 10:35:14 +0100 (Fri, 14 Feb 2014) | 5 lines

Le add_pbl_tend est passé de phy1d à phylmd et appelé dans phylmd/physiq.F90
si klon_glo=1
add_pbl_tend is moved from phy1d to phylmd and called in phylmd/physiq.F90
if klon_glo=1

------------------------------------------------------------------------
r1972 | fhourdin | 2014-02-14 00:10:10 +0100 (Fri, 14 Feb 2014) | 4 lines

Suppression de fichiers inutiles dans phy1d
Removing useless files from phy1d


------------------------------------------------------------------------
r1971 | idelkadi | 2014-02-13 15:45:29 +0100 (Thu, 13 Feb 2014) | 11 lines

Corrections pour les sorties XIOS :

1. Dans phys_output_mod.F90 :
   Calcul de l'heure de départ sorti en dehors de la boucle sur les fichiers
   Definition des axes verticaux pour XIOS mis sous la cles de controle des sorties 
   Definition de l'axe vertical pour les fichiers NMC

2. Corrections divers pour les sorties NMC (unite, description, ...)
   phys_output_ctrlout_mod.F90 et phys_output_write_mod.F90


------------------------------------------------------------------------
r1970 | fhourdin | 2014-02-12 09:17:24 +0100 (Wed, 12 Feb 2014) | 3 lines

Correction de la commission précédente
Correction of previous commit

------------------------------------------------------------------------
r1969 | fhourdin | 2014-02-11 19:11:45 +0100 (Tue, 11 Feb 2014) | 4 lines

Correction. Declarations manquantes.
Bug fixing (mixing declarations).


------------------------------------------------------------------------
r1968 | fhourdin | 2014-02-11 16:01:12 +0100 (Tue, 11 Feb 2014) | 4 lines

Nouvelle version de thermcell_plume, permettant d'avoir des stratocumulus
    avec 40 couches dans  le cas fire.
New version of thermcell_plume, giving stratocumulus in the fire case

------------------------------------------------------------------------
r1967 | jghattas | 2014-02-11 13:54:07 +0100 (Tue, 11 Feb 2014) | 3 lines

Added path to XIOS for curie arch.


------------------------------------------------------------------------
r1966 | fairhead | 2014-02-11 11:52:47 +0100 (Tue, 11 Feb 2014) | 10 lines

Correction de bug pour le cas Sandu.
Nudge au dessus de 700hpa à prendre en compte. 
MPL

 --This line, and those below, will be ignored--

M    phy1d/1D_decl_cases.h
A    phy1d/1D_nudge_sandu_astex.h
M    phy1d/lmdz1d.F

------------------------------------------------------------------------
r1965 | acaubel | 2014-02-11 11:42:52 +0100 (Tue, 11 Feb 2014) | 3 lines

AC : Modified to run with Oasis-MCT as coupler in IPSLCM6 coupled configuration


------------------------------------------------------------------------
r1964 | fhourdin | 2014-02-10 13:54:50 +0100 (Mon, 10 Feb 2014) | 4 lines

Correction d'une erreur dans makelmdz.
Bug fixing in makelmdz

------------------------------------------------------------------------
r1963 | fhourdin | 2014-02-08 01:40:31 +0100 (Sat, 08 Feb 2014) | 1 line

Supprimé le fichier phy1d/ocean_forced_mod.F90
------------------------------------------------------------------------
r1962 | fhourdin | 2014-02-08 01:38:55 +0100 (Sat, 08 Feb 2014) | 3 lines

Le create_make_gcm manquait beaucoup de dépendances.


------------------------------------------------------------------------
r1961 | fhourdin | 2014-02-08 01:37:48 +0100 (Sat, 08 Feb 2014) | 6 lines

ocean_force_mod est mis à jour pour inclure en option le cas 1D, et éviter
    de dupliquer la routine dans phy1d.
clesphys.h est modifié pour placer les réels avant les entiers et logiques
    pour éviter des warnings à la compilation.


------------------------------------------------------------------------
r1960 | fhourdin | 2014-02-07 17:20:25 +0100 (Fri, 07 Feb 2014) | 7 lines

Nettoyage du code 1D pour limiter les warning en mode debug.
Essentiellement des declarations inutiles.

Cleaning of the 1D code in order to limit warnings in debug mode.
Essentially supressing unused declared variables.


------------------------------------------------------------------------
r1959 | fhourdin | 2014-02-07 10:37:38 +0100 (Fri, 07 Feb 2014) | 26 lines

Nouvelles options dans disvert.F90 :
=------------------------------------
1. avec vert_sampling="sigma" -> niveaux sigma plutôt qu'hybride.
    P_i,j,k= Bhyb_k * Ps  , Ahyb_k=0.
2. disgmin, l'épaisseur approximative du premier niveau en coordonnée
    sigma peut être défini dans les fichiers .def.
    Cette option est valable pour vert_sampling=sigma, tropo ou strato.
    Valeur par défaut : 1, sauf pour   (llm==39 .and. vert_sampling=="strato")
               ce qui correspond aux configurations CMIP5
=======================================================================
New options in disvert.F90 :
=---------------------------
1. with vert_samplig="sigma" -> sigma levels instead of hybrid 
    P_i,j,k= Bhyb_k * Ps  , Ahyb_k=0.
2. disgmin, the approximate thickness of the first layer in sigma
    coordinate can be specified from the .def files.
    Active for vert_sampling=sigma, tropo or strato.
    Defaut value : 1, except for  (llm==39 .and. vert_sampling=="strato")
                   which corresponds to the CMIP5 configurations. 



--Cette ligne, et les suivantes ci-dessous, seront ignorées-

M    dyn3d_common/disvert.F90

------------------------------------------------------------------------
r1958 | fairhead | 2014-02-03 15:29:35 +0100 (Mon, 03 Feb 2014) | 8 lines


Empêche l'exécution en mode debug

----------------

Crashes in debugging mode
------------------------------------------------------------------------
r1957 | fairhead | 2014-02-03 15:27:04 +0100 (Mon, 03 Feb 2014) | 6 lines

L'édition des liens avec la librairie idrmemmpi de l'idris est supprimée par
défaut pour éviter que les jobs ne se bloquent sur Ada. En effet, dans le cas
d'une simulation couplée, il faut que tout les exécutables aient été liés avec
la librairie pour éviter un blocage et comme ce n'est pas le cas par défaut ...


------------------------------------------------------------------------
r1956 | lguez | 2014-01-31 13:42:35 +0100 (Fri, 31 Jan 2014) | 1 line

Correcting bad spacing left over from fixed source form.
------------------------------------------------------------------------
r1955 | fhourdin | 2014-01-31 09:17:49 +0100 (Fri, 31 Jan 2014) | 7 lines

Options for idealized ozone (read_climoz=0) :
Imposing ozone be constant in case were solarlong0> -999 
(defaut value is -999.999 when seasonal cycle is compiuted) and
to be symetric with respect to equator when read_climoz=-1
FH


------------------------------------------------------------------------
r1954 | lguez | 2014-01-29 16:36:27 +0100 (Wed, 29 Jan 2014) | 7 lines

Files in dyn3d_common use files in bibio and vice-versa. So, for
compilation with makelmdz, we add a reference to dyn3d_comm at link
time.

Removed dyn3dmem/nxgrarot_p.F which was identical to
/dyn3dpar/nxgrarot_p.F but not used in dyn3dmem configuration.

------------------------------------------------------------------------
r1953 | fhourdin | 2014-01-29 11:40:39 +0100 (Wed, 29 Jan 2014) | 2 lines

Correction in for readprofiles_fire

------------------------------------------------------------------------
r1952 | lguez | 2014-01-28 14:05:47 +0100 (Tue, 28 Jan 2014) | 1 line

Other files with unneeded differences in dyn3d, dyn3dpar and dyn3dmem.
------------------------------------------------------------------------
r1951 | fairhead | 2014-01-28 13:38:55 +0100 (Tue, 28 Jan 2014) | 3 lines

Pour eviter que ça ne plante pour rien


------------------------------------------------------------------------
r1950 | emillour | 2014-01-27 17:16:37 +0100 (Mon, 27 Jan 2014) | 3 lines

Correction from last commit for iflag_thermals=0
FH

------------------------------------------------------------------------
r1949 | emillour | 2014-01-27 16:28:09 +0100 (Mon, 27 Jan 2014) | 3 lines

Due to a "bug" with ifort compiler version 13 (compilation failure at linking stage when using OpenMP), revert to usin
g static version of HDF5.
EM 

------------------------------------------------------------------------
r1948 | idelkadi | 2014-01-27 11:29:29 +0100 (Mon, 27 Jan 2014) | 6 lines

- Option pour tourner sans thermiques ni ajustement sec si ifalg_thermals<0
  (phylmd/physiq.F90)

- Modifications permettant d'imposer des profils initiaux de traceurs
lus dans un tracer.inp.001 (phylmd/1DUTILS.h, 1D_read_forc_cases.h et lmdz1d.F)
------------------------------------------------------------------------
r1947 | idelkadi | 2014-01-27 11:25:51 +0100 (Mon, 27 Jan 2014) | 2 lines

Changement de liens dans phy1d suite au regroupement de sous programmes de la dynamique dans dyn3d_common 

------------------------------------------------------------------------
r1946 | lguez | 2014-01-24 18:21:46 +0100 (Fri, 24 Jan 2014) | 1 line

A procedure in dyn3d_common calls a procedure in phylmd so we need an additional -lphylmd, for compilation with makelm
dz.
------------------------------------------------------------------------
r1945 | lguez | 2014-01-23 15:11:43 +0100 (Thu, 23 Jan 2014) | 1 line

Duplicated files moved to dyn3d_common.
------------------------------------------------------------------------
r1944 | lguez | 2014-01-22 18:39:15 +0100 (Wed, 22 Jan 2014) | 5 lines

Created directory dyn3d_common. Began moving to dyn3d_common files
that were duplicated in dyn3d, dyn3dpar and dyn3dmem. Adapted
compilation files for makegcm, makelmdz and makelmdz_fcm. Adapted
install_1d_src.sh.

------------------------------------------------------------------------
r1943 | idelkadi | 2014-01-22 10:51:36 +0100 (Wed, 22 Jan 2014) | 51 lines

Changes concerinng the Hourdin et al. 2002 version of the thermals and 
Ayotte cases

calltherm.F90
-------------
call to thermcell_2002 modified :
1) iflag_thermals changed from 1 to  >= 1000 ; iflag_thermals-1000 
controls sub options.
2) thermals w and fraction in output files.

hbtm.F
------
Singularity in the 1/heat dependency of the Monin Obukov length L when 
heat=0. Since 1/L is used rather than L, it is preferable to compute 
directly L. There is a dependency in 1/u* then which is treated with a 
threshold value.
(+ some cleaning in the syntax).

lmdz1d.F
--------
If nday<0, -nday is used as the total number of time steps of the 
simulation.
The option with imposed wtsurf and wqsurf read in lmdz1d.def was not 
active anymore.
<       IF(.NOT.ok_flux_surf) THEN
changed to
>       IF(.NOT.ok_flux_surf.or.max(abs(wtsurf),abs(wqsurf))>0.) THEN
before
        fsens=-wtsurf*rcpd*rho(1)
        flat=-wqsurf*rlvtt*rho(1)

phys_output_write.F90 et phys_local_var_mod.F90
-----------------------------------------------
Removing the d_u_ajsb contribution to duthe (same for dv).
Those tendencies are not computed by the model ...
<           zx_tmp_fi3d(1:klon,1:klev)=d_u_ajs(1:klon,1:klev)/pdtphys - &
<                d_u_ajsb(1:klon,1:klev)/pdtphys
---
>           zx_tmp_fi3d(1:klon,1:klev)=d_u_ajs(1:klon,1:klev)/pdtphys
> !              d_u_ajsb(1:klon,1:klev)/pdtphys

thermcell_dq.F90, thermcell_main.F90
-----------------------------------
Some cleaning

thermcell_old.F
---------------
Control by iflag_thermals >= 1000
wa_moy and fraca in outputs
+ cleaning

------------------------------------------------------------------------
r1939 | lguez | 2014-01-21 15:23:17 +0100 (Tue, 21 Jan 2014) | 5 lines

Same as revision 1930: replaced abort by abort_gcm.

Also replaced real*8 by real(kind=8) (was done way back in revision 1220 for
dyn3d/fxhyp.F but not dyn3dpar/fxhyp.F).

------------------------------------------------------------------------
r1938 | lguez | 2014-01-20 17:20:05 +0100 (Mon, 20 Jan 2014) | 8 lines

New run-time parameters: ok_gwd_rando, ruwmax and sat. Logical
variable ok_gwd_rando activates the new stochastic parameterization of
the drag of non-orographic gravity waves, by F. Lott. Default value of
ok_gwd_rando is false. A run with ok_gwd_rando creates new variables
in restartphy.nc: du_gwd_rando and dv_gwd_rando. Possible output in
history files of new variables du_gwd_rando, dv_gwd_rando and
vstr_gwd_rando.

------------------------------------------------------------------------
r1937 | lguez | 2014-01-20 16:31:20 +0100 (Mon, 20 Jan 2014) | 1 line

Added test on ipos. Program could crash because of a bad value of phys_out_filetimesteps, with an unclear message.
------------------------------------------------------------------------
r1933 | lguez | 2014-01-17 19:06:50 +0100 (Fri, 17 Jan 2014) | 3 lines

real*8 is not in the Fortran standard.
(There was a compilation error with gfortran and debugging options.)

------------------------------------------------------------------------
r1932 | lguez | 2014-01-17 19:03:16 +0100 (Fri, 17 Jan 2014) | 8 lines

Declaration of variables appearing in array bounds must be specified
before the declaration of the array.

Variables used in FCTTRE.h are declared in YOMCST.h and YOETHF.h so
YOMCST.h and YOETHF.h must be included before FCTTRE.h.

(There were compilation errors with gfortran and debugging options.)

------------------------------------------------------------------------
r1931 | lguez | 2014-01-17 18:49:48 +0100 (Fri, 17 Jan 2014) | 2 lines

abort and exit are not in the Fortran standard. Replaced them by abort_gcm.

------------------------------------------------------------------------
r1930 | lguez | 2014-01-17 17:45:09 +0100 (Fri, 17 Jan 2014) | 4 lines

abort, dfloat and pause are not in the Fortran standard. Replaced
abort by abort_gcm and dfloat by dble. Note: I modified dyn3dpar files
that were identical to dyn3d modified files.

------------------------------------------------------------------------
r1929 | fairhead | 2014-01-17 17:34:47 +0100 (Fri, 17 Jan 2014) | 3 lines

Correction d'une correction bugguée


------------------------------------------------------------------------
r1927 | fairhead | 2014-01-10 12:49:47 +0100 (Fri, 10 Jan 2014) | 6 lines

Ménage pour éviter une sortie de script intempestive

----------------

Cleanup to prevent spurrious script leaving

------------------------------------------------------------------------
r1926 | idelkadi | 2014-01-08 22:47:56 +0100 (Wed, 08 Jan 2014) | 19 lines

- Re-ecriture des sorties du simulateur Cosp :
Les 3 fichiers mensuel, journalier et haute frequence ne sont plus geres de facon classique par IOIPSL.
Les fichiers "includes" ini_hist*COSP.h et write_hist*COSP.h sont supprimes et remplaces par 2 module :
1. cosp_output_mod.F90 : ou sont crees les fichiers et ou sont definis les dimensions et les differents axes.
2. cosp_output_write_mod.F90 : ou sont definis les variables diagnostiques a stocker dans ces fichiers et ou est ge
ree leur ecriture.

L'utilisation d'XIO est prevu (A tester)

- Rajouts d'autre variables diagnostiques :
sunlit (=0 si nuit et =1 si jour)
parasol_crefl (reflectance integree)

- Correction du bug dans l'interface avec la physique :
On ne veut pas que la distinction entre les nuages convectifs et stratiformes soient prise en compte dans les calcul
s Cosp.
Dans l'interface avec la physique, sont passees en entree pour Cosp, les quantites totales (stratiforme + convective
) du contenus en eau et d'autres variables.
La fraction nuageuse convective calculee dans la physique est passee en entree pour Cosp dans la version buggee. Ell
e est remise a 0 dans cette version corrigee.

- Mise a jour pour ISCCP : la fraction d'ensoleillement calculee par LMDZ est passee en argument pour Isccp

------------------------------------------------------------------------
r1925 | idelkadi | 2014-01-08 22:38:56 +0100 (Wed, 08 Jan 2014) | 20 lines

L'utilisation d'XIO est prevu (A tester)

- Rajouts d'autre variables diagnostiques :
sunlit (=0 si nuit et =1 si jour)
parasol_crefl (reflectance integree)

- Mise a jour pour ISCCP : la fraction d'ensoleillement calculee par LMDZ est passee en argument pour Isccp

------------------------------------------------------------------------
r1924 | idelkadi | 2014-01-08 15:41:30 +0100 (Wed, 08 Jan 2014) | 7 lines

- Rajout de nouveaux diagnostiques pour la comparaison entre les fractions nuageuses sorties derictement dans LMDZ et
celles sorties via le simulateur ISCCP.
Routines concernees : physiq.F90, phys_output_ctrlout_mod.F90, phys_local_var_mod.F90 et phys_output_write_mod.F90
- Rajout de la fraction d'ensoleillement par jour aux arguments de l'interface avec Cosp (entree pour le simulateur IS
CCP).
Routine concernée : physiq.F90

------------------------------------------------------------------------
r1923 | lguez | 2014-01-06 17:07:02 +0100 (Mon, 06 Jan 2014) | 1 line

Replace test on level index (which was probably tailored for a vertical grid with 19 levels) by test on sigma = pressu
re / surface_pressure.
------------------------------------------------------------------------
r1922 | emillour | 2014-01-06 10:00:25 +0100 (Mon, 06 Jan 2014) | 4 lines

Revert Curie arch file to previous version (i.e. remove the large memory options) to remain compatible with other IPSL
softwares.
Once the appropriate testing has been done, all should then simultaneously add the large memory (-mcmodel=large at com
pilation and -shared-intel when linking) to enable running at high resolution.
EM

------------------------------------------------------------------------