by Jean-Daniel GRIL
CNRM/GMAP/COOPE
E-mail : Jean-Daniel.GRIL++at++Meteo.fr
Tel : +33 (0) 561 078 416
For the size of reals and integers :
for files in *.F90 the size is selected by JPIM and JPRB (in tmsbkind.h and PARKIND1.F90) for all variables.
It's no necessary to add something like -r8 -i4 (the classic values of JPRB JPIM) in F90FLAGS (Makefile.in.OSMFC).
for files in *.F you must add this compiler flags in F77FLAGS.
From the V1.00beta version :
2002.07.08 : V1.00
In ./src/xrd/grib_mf/Makefile.loc with the IS_LIT_E key for Little Endian
In ./src/jdglib/module/EGG_TOOLS_MOD.F90, I add lat-lon case (ERPK < 0.0) in SUBROUTINE OUTEGGX
2002.09.06 : V1.01
EGG_TOOLS_MOD.F90 : modication for lat-lon domain in old eggx files
2002.10.06 : V1.02
eggpack.F90 : correction of bug in xy_to_rtheta : if y=0 then theta is not given by "atan(x/y)".
EGG_TOOLS_MOD.F90 : lines 78 & 79 too long for some compilers (breaking lines)
mykind.h : add a DEBUG define key to make outputs or not (see domolalo.F90 in src/jdgprg/domolalo)
domolalo.F90 : a new program. Define an Aladin domain (plus eight points on each sides containing a lat-lon domain).
Inputs are : SOUTH NORTH WEST EAST RESOLUTION (the same in X and Y).
the Reference point latitude is (NORTH+SOUTH)/2 | excepted in some cases
the Reference point longitude is (WEST+EAST)/2 | excepted in some cases
West and East may have the same value : in this case the domain is from West to East.
For MERCATOR you need to use SOUTH=-NORTH (if NORTH > 70° then a question asks to cut or stop)
For STEREOPOLAR : If NORTH = |SOUTH| = 90° then the domain is centered on pole ( Center domain
latitude = 90° )
On outputs the SW,NE,SE,NW have the same latitude
This case don't work with old wersion of EGGX
For LAMBERT : If |(NORTH+SOUTH)/2| > 70° then Reference point Latitude=90° and program make a STEREOPOLAR
If |(NORTH+SOUTH)/2| < 20° then program make a MERCATOR (except if NORTH or |SOUTH| > 70°)
If pole will appear in lambert output domain then program make a STEREOPOLAR
If WEST = EAST then LAMBERT is changed in STEREOPOLAR or MERCATOR and the added points are supress on X.
You can give input datas : from command line (five value : SOUTH,NORTH,WEST,EAST,RESOL)
from namelist (one cmd line value : name of namelist filethat contains
NAMELIST /NAML/SOUTH,NORTH,WEST,EAST,RESOL)
from interactive mode (no parameters in cmd line)
Outputs are in a namelist file named : namelist_domaine