Updates Version 1.02

2002/10/06

by Jean-Daniel GRIL

CNRM/GMAP/COOPE

E-mail : Jean-Daniel.GRIL++at++Meteo.fr

Tel : +33 (0) 561 078 416



News, Changes


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.



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


EGG_TOOLS_MOD.F90 : modication for lat-lon domain in old eggx files


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





Home Home Top Top