From yunheng at ou.edu  Tue Apr  1 16:42:17 2008
From: yunheng at ou.edu (Yunheng Wang)
Date: Tue, 01 Apr 2008 16:42:17 -0500
Subject: [ARPSSUPPORT] test2
Message-ID: <47F2AC39.2010601@ou.edu>

This is to test the link from arpssupport at ou.edu to arpssupport at caps.ou.edu.

Yunheng.

From yunheng at ou.edu  Thu Apr  3 09:19:43 2008
From: yunheng at ou.edu (Yunheng Wang)
Date: Thu, 03 Apr 2008 09:19:43 -0500
Subject: [ARPSSUPPORT] arps5.2.8
In-Reply-To: <20080403020010.GA13324@lr6.msi.umn.edu>
References: <20080318202809.GE5662@lr6.msi.umn.edu> <47E031F2.2070707@ou.edu>
	<20080318223323.GG5662@lr6.msi.umn.edu> <47E046F0.9050402@ou.edu>
	<20080403020010.GA13324@lr6.msi.umn.edu>
Message-ID: <47F4E77F.3050109@ou.edu>




----- Original Message -----
*From:* "Robert W. Numrich" <rwn at msi.umn.edu>
*To:* "Wang, Yunheng" <yunheng at ou.edu>
*Sent:* 04/02/2008 9:00:11 PM -0500
*Subject:* arps5.2.8



> Can you please help me with these two problems?
>
> 1) First Problem
>
> make[1]: Entering directory `/home/xe1/numrich/arps/arps5.2.8/src/arpsens'
> ifort    -O3  -o /home/xe1/numrich/arps/arps5.2.8/bin/bin2gem bin2gem.o arpsenslib.o postcore.o arpspostlib.o gemio.o          \
>                                /home/xe1/numrich/arps/arps5.2.8/lib/libadas.a     \
>                                /home/xe1/numrich/arps/arps5.2.8/lib/libarps.a     \
>
> bin2gem.o: In function `MAIN__':
> bin2gem.f90:(.text+0x662f): undefined reference to `in_bdta_'
> bin2gem.f90:(.text+0x6a86): undefined reference to `gr_mnav_'
> gemio.o: In function `enswrtgem_':
> gemio.f90:(.text+0x8fa): undefined reference to `gd_cref_'
> gemio.f90:(.text+0x963): undefined reference to `gd_swrt_'
> gemio.f90:(.text+0x2b10): undefined reference to `gd_clos_'
> gemio.o: In function `initgemio_':
> gemio.f90:(.text+0x2d32): undefined reference to `in_bdta_'
> gemio.f90:(.text+0x3011): undefined reference to `gr_mnav_'
> gemio.o: In function `gemdump2d_':
> gemio.f90:(.text+0x35b1): undefined reference to `gd_wpgd_'
> gemio.f90:(.text+0x36ca): undefined reference to `gd_wpgd_'
> make[1]: *** [/home/xe1/numrich/arps/arps5.2.8/bin/bin2gem] Error 1
> make[1]: Leaving directory `/home/xe1/numrich/arps/arps5.2.8/src/arpsens'
> make: *** [bin2gem] Error 2
>
> ------------------------------------------------------------------------
>
>   
GEMPAK is not linked. Since your post is not complete, I do not know 
which library is missed. Please do the following

1. install GEMPAK
2. source $GEMPAK_DIR/Gemenviron
3. setenv NAWIPS $GEMPAK_DIR
4 compile the program again.


> 2)  Second Problem
>
> ARPS Makefile configuration utility
> Updated: 03/14/2007 (ARPS version 5.2.8)
>
> Making all for linux
>
> Compilation options: 3
> Support for HDF I/O.
>
> WARNING: HDF library is required for 'all' and it is not found in the default
>          location '/soft/hdf5/hdf5-1.8.0/HDF5'. Please install HDF library, define environment
>          variable 'HDFPATH' to be the root of the installed HDF package
>          and try again.
>
>          If you do not want to use HDF data format, please try option '-io nohdf'
>          For example:
>                        ./makearps -io nohdf all
>
> The libraries have been installed as follows:
>
> calhoun> echo $HDFPATH
> /soft/hdf5/hdf5-1.8.0/HDF5
>
> Directory: /soft/hdf5/hdf5-1.8.0/HDF5/lib
> calhoun> ls -lrt
> total 12912
> -rwxr-xr-x 1 swinst swinst  554508 2008-03-21 18:01 libhdf5_cpp.so.0.0.0*
> lrwxrwxrwx 1 swinst swinst      20 2008-03-21 18:01 libhdf5_cpp.so.0 -> libhdf5_cpp.so.0.0.0*
> lrwxrwxrwx 1 swinst swinst      20 2008-03-21 18:01 libhdf5_cpp.so -> libhdf5_cpp.so.0.0.0*
> -rwxr-xr-x 1 swinst swinst     905 2008-03-21 18:01 libhdf5_cpp.la*
> -rw-r--r-- 1 swinst swinst 1141110 2008-03-21 18:01 libhdf5_cpp.a
> -rwxr-xr-x 1 swinst swinst   16630 2008-03-21 18:01 libhdf5_hl_cpp.so.0.0.0*
> lrwxrwxrwx 1 swinst swinst      23 2008-03-21 18:01 libhdf5_hl_cpp.so.0 -> libhdf5_hl_cpp.so.0.0.0*
> lrwxrwxrwx 1 swinst swinst      23 2008-03-21 18:01 libhdf5_hl_cpp.so -> libhdf5_hl_cpp.so.0.0.0*
> -rwxr-xr-x 1 swinst swinst     926 2008-03-21 18:01 libhdf5_hl_cpp.la*
> -rw-r--r-- 1 swinst swinst   11098 2008-03-21 18:01 libhdf5_hl_cpp.a
> -rwxr-xr-x 1 swinst swinst 3299584 2008-03-21 21:24 libhdf5.so.5.0.0*
> lrwxrwxrwx 1 swinst swinst      16 2008-03-21 21:24 libhdf5.so.5 -> libhdf5.so.5.0.0*
> lrwxrwxrwx 1 swinst swinst      16 2008-03-21 21:24 libhdf5.so -> libhdf5.so.5.0.0*
> -rwxr-xr-x 1 swinst swinst     829 2008-03-21 21:24 libhdf5.la*
> -rw-r--r-- 1 swinst swinst 7266194 2008-03-21 21:24 libhdf5.a
> -rw-r--r-- 1 swinst swinst    2405 2008-03-21 21:24 libhdf5.settings
> -rwxr-xr-x 1 swinst swinst     801 2008-03-21 21:24 libhdf5_fortran.la*
> -rw-r--r-- 1 swinst swinst  314508 2008-03-21 21:24 libhdf5_fortran.a
> -rw-r--r-- 1 swinst swinst     670 2008-03-21 21:24 libhdf5_fortran.settings
> -rwxr-xr-x 1 swinst swinst  147232 2008-03-21 21:24 libhdf5_hl.so.0.0.0*
> lrwxrwxrwx 1 swinst swinst      19 2008-03-21 21:24 libhdf5_hl.so.0 -> libhdf5_hl.so.0.0.0*
> lrwxrwxrwx 1 swinst swinst      19 2008-03-21 21:24 libhdf5_hl.so -> libhdf5_hl.so.0.0.0*
> -rwxr-xr-x 1 swinst swinst     850 2008-03-21 21:24 libhdf5_hl.la*
> -rw-r--r-- 1 swinst swinst  275964 2008-03-21 21:24 libhdf5_hl.a
> -rwxr-xr-x 1 swinst swinst     807 2008-03-21 21:24 libhdf5hl_fortran.la*
> -rw-r--r-- 1 swinst swinst   96038 2008-03-21 21:24 libhdf5hl_fortran.a
>
>
>
> --
>
>   
The ARPS package uses HDF library, not HDF5 library. Please visit NCSA 
web about HDF for details. HDF and HDF5 are different.

Yunheng.

>                         -bob numrich
>
> --
> Robert W. Numrich, PhD
> Senior Research Associate
> Minnesota Supercomputing Institute
> 465 Walter Library
> 117 Pleasant Street SE
> Minneapolis, MN 55455
>
> Tel:  612-624-4341
>   

From mdaniels at ce.berkeley.edu  Thu Apr  3 18:13:48 2008
From: mdaniels at ce.berkeley.edu (Megan Daniels)
Date: Thu, 03 Apr 2008 16:13:48 -0700
Subject: [ARPSSUPPORT] plot clouds
Message-ID: <9c44ef6a0804031613i5d7e3b40jbef67cdd9a29b671@mail.gmail.com>

Dear ARPS Support,

What is the simplest way to plot clouds if I am interested in whether or not
a cloud is affecting the incoming solar radiation at the surface?  For
example, it would be desirable to plot the presence (or absence) of clouds
in the horizontal (x and y).  Which arps variable(s) do you recommend
plotting?

Thank you,
Megan

-- 
Megan Daniels
Ph.D. Candidate
Department of Civil and Environmental Engineering
760 Davis Hall
University of California, Berkeley
94720-1710

Tel. (510) 643-0018
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.caps.ou.edu/pipermail/arpssupport/attachments/20080403/0204a3ad/attachment.html 

From mka94 at student.canterbury.ac.nz  Wed Apr  9 17:01:31 2008
From: mka94 at student.canterbury.ac.nz (Marwan Katurji)
Date: Thu, 10 Apr 2008 10:01:31 +1200
Subject: [ARPSSUPPORT] ARPS Help:trouble shooting
Message-ID: <47FD3CBB.6080500@student.canterbury.ac.nz>

Hi there,

I am a registered user for arps and have downloaded version 5.2.8 and am 
using it successfully on my PC under cygwin "but not in parallel" for 
boundary layer modeling and LES. Recently arps was installed and 
compiled on a high performance computer, IBM P5-575 containing16 IBM 
Power 5, 1.9 Giga Hertz (GHz) processors and 32 GigaBytes (GB) of 
memory. The purpose is to run it in parallel mode under one node and 
using 8 processors. The arps_mpi and all other utilities I need compiled 
fine with no errors an I am using the serial version of arps on one 
processor successfully. For testing the arps_mpi I tried to run the demo 
may20 simulation and it appeared to work fine too with 2 and 4 
processors using ascii dump files. The parallel job is submitted to the 
hpc via a load leveler script.

I have some issues that I need your help with, can't really figure out 
if they are code related or platform:

1) The may20 simulation was not successful using 8 processors (nxproc=2, 
nyproc=4 and vice versa) using ascii dump, and the only way the 
simulation ran was using the combination (nxproc=1 and nyproc=8, the 
opposite did not work). Unfortunately I don't get a comprehensive error 
message but the program stops directly after when it starts. Is this 
something related to the platform or the code?

2) Using netcdf history dump format the combination of (nxproc=2, 
nyproc=4 and vice versa) worked fine but only when I dumped the output 
from each processor individually with out joining. The join option is 
not working with netcdf and I get a segmentation fault in the load 
leveler error output and an error in DTADUMP in arps output.

3) When using a different simulation that involves terrain file, surface 
physics and other parameters that you can see in the attached log file, 
the simulation progressed normally using 4 processors (nxproc=2, 
nyproc=2) but suddenly at time =7200s when it was dumping the history 
joined file an error occurred stating "error occurred when opening file 
.//xxxx.asc007200 using fortran unit 39" and when looking at some 
formatted dumped data in the arps.output file they appeared to be 
incomplete, that is some missing lines in the gridded output matrix and 
some redundant ones. (I suspect that the dumping from processors is not 
working correctly).

When decreasing the processors to 2 (nxproc=1,nxproc=2) the same happend 
but occured at a later time, 46800s. You can see a sample of the output 
file attached.

Thanks for your time and would appreciate your help. This work is part 
of my PhD study at University of Canterbury, New Zealand.

-- 
Marwan Katurji

PhD Candidate
Center for Atmospheric Research
Department of Geography, 2nd floor
University of Canterbury
Private Bag 4800
Christchurch 8020, New Zealand
-----------------
Work:  +64 3 364 2987  ext: 7925
Fax:   +64 3 364 2907

Email: marwan.houda at gmail.com
       mka94 at student.canterbury.ac.nz

Home Address:
1/92 Office Road
Merivale, Christchurch 8014
New Zealand
Tel: +64 3 3555760   (land line)
     +64 3 212049534 (mobile)


From mka94 at student.canterbury.ac.nz  Wed Apr  9 18:37:32 2008
From: mka94 at student.canterbury.ac.nz (Marwan Katurji)
Date: Thu, 10 Apr 2008 11:37:32 +1200
Subject: [ARPSSUPPORT] ARPS Help:trouble shooting
In-Reply-To: <47FD4901.7080100@caps.ou.edu>
References: <47FD3CBB.6080500@student.canterbury.ac.nz>
	<47FD4901.7080100@caps.ou.edu>
Message-ID: <47FD533C.5010809@student.canterbury.ac.nz>

Thanks Yunheng for your prompt reply,

You mean that ascii dump is not supported in the mpi mode? because the 
serial arps works fine with ascii dump and it is listed in the name list 
as one of the dump options.

Is there a way I can join splitted netcdf files? The joinfiles utility I 
have works on hdf4 and binary output, or do you recommend in this case 
dumping binary in joined mode thus using less memory and consequently I 
can use arpspltpost for plotting binary history files?

Many Thanks

Marwan Katurji

PhD Candidate
Center for Atmospheric Research
Department of Geography, 2nd floor
University of Canterbury
Private Bag 4800
Christchurch 8020, New Zealand
-----------------
Work:  +64 3 364 2987  ext: 7925
Fax:   +64 3 364 2907

Email: marwan.houda at gmail.com
       mka94 at student.canterbury.ac.nz

Home Address:
1/92 Office Road
Merivale, Christchurch 8014
New Zealand
Tel: +64 3 3555760   (land line)
     +64 3 212049534 (mobile)



ARPSSUPPORT wrote:
> No surprise. ASCII format is NOT supported any more. For the netCDF 
> format, it may be that the platform does not have enough memory to 
> hold joined arrays.
>
> The supported format are
>
> !          = 0, no data dump is produced;
> !          = 1, unformatted binary data dump;
> !          = 3, NCSA HDF4 format data dump;
> !          = 7, NetCDF format;
> !          = 8, NetCDF format with data of all time levels in one file;
> !               It is suitable for small domain and is provided for LEAD.
> !          = 9, GrADS data dump;
> !          = 11, Vis5D data dump.
>
> Thanks,
> Yunheng Wang.
>
>
> ----- Original Message -----
> *From:* Marwan Katurji <mka94 at student.canterbury.ac.nz>
> *To:* arpssupport at ou.edu
> *Sent:* 04/09/2008 5:01:31 PM -0500
> *Subject:* [ARPSSUPPORT] ARPS Help:trouble shooting
>
>
>
>> Hi there,
>>
>> I am a registered user for arps and have downloaded version 5.2.8 and 
>> am using it successfully on my PC under cygwin "but not in parallel" 
>> for boundary layer modeling and LES. Recently arps was installed and 
>> compiled on a high performance computer, IBM P5-575 containing16 IBM 
>> Power 5, 1.9 Giga Hertz (GHz) processors and 32 GigaBytes (GB) of 
>> memory. The purpose is to run it in parallel mode under one node and 
>> using 8 processors. The arps_mpi and all other utilities I need 
>> compiled fine with no errors an I am using the serial version of arps 
>> on one processor successfully. For testing the arps_mpi I tried to 
>> run the demo may20 simulation and it appeared to work fine too with 2 
>> and 4 processors using ascii dump files. The parallel job is 
>> submitted to the hpc via a load leveler script.
>>
>> I have some issues that I need your help with, can't really figure 
>> out if they are code related or platform:
>>
>> 1) The may20 simulation was not successful using 8 processors 
>> (nxproc=2, nyproc=4 and vice versa) using ascii dump, and the only 
>> way the simulation ran was using the combination (nxproc=1 and 
>> nyproc=8, the opposite did not work). Unfortunately I don't get a 
>> comprehensive error message but the program stops directly after when 
>> it starts. Is this something related to the platform or the code?
>>
>> 2) Using netcdf history dump format the combination of (nxproc=2, 
>> nyproc=4 and vice versa) worked fine but only when I dumped the 
>> output from each processor individually with out joining. The join 
>> option is not working with netcdf and I get a segmentation fault in 
>> the load leveler error output and an error in DTADUMP in arps output.
>>
>> 3) When using a different simulation that involves terrain file, 
>> surface physics and other parameters that you can see in the attached 
>> log file, the simulation progressed normally using 4 processors 
>> (nxproc=2, nyproc=2) but suddenly at time =7200s when it was dumping 
>> the history joined file an error occurred stating "error occurred 
>> when opening file .//xxxx.asc007200 using fortran unit 39" and when 
>> looking at some formatted dumped data in the arps.output file they 
>> appeared to be incomplete, that is some missing lines in the gridded 
>> output matrix and some redundant ones. (I suspect that the dumping 
>> from processors is not working correctly).
>>
>> When decreasing the processors to 2 (nxproc=1,nxproc=2) the same 
>> happend but occured at a later time, 46800s. You can see a sample of 
>> the output file attached.
>>
>> Thanks for your time and would appreciate your help. This work is 
>> part of my PhD study at University of Canterbury, New Zealand.
>>
>>   
>

From arpssupport at ou.edu  Thu Apr 10 07:47:42 2008
From: arpssupport at ou.edu (arpssupport at ou.edu)
Date: Thu, 10 Apr 2008 07:47:42 -0500 (CDT)
Subject: [ARPSSUPPORT] April 77% OFF
Message-ID: <20080410111427.6707.qmail@ally12>


From leanderson at cefet-rj.br  Thu Apr 10 09:19:56 2008
From: leanderson at cefet-rj.br (Leanderson Marcos)
Date: Thu, 10 Apr 2008 11:19:56 -0300
Subject: [ARPSSUPPORT] terrain 3sec and GFS 0.5 degree
Message-ID: <HHEIJOJABICFEFAGHOHFCEFOCKAA.leanderson@cefet-rj.br>

Dear friend,
1)	I have 3 sec terrain data file that was processed by arpstrn with these
grid parameters:
nx = 4501
ny = 3501
dx = 90.
dy = 90.

2)	In another way, I have 0.5 degree NCEP/GFS data file in GRIB2, that I
should processes by ext2arps program to create boundary conditions for run
arps program. On the arps.input file processed to ext2arps I have:
nx= 7
ny = 5
dx = 50000.
dy = 50000.

If I try processes ext2arps reading the terrain file created by arpstrn I
generate a ERROR due the not matching grids. So, what can I use high
resolution terrain on 0.5 degree GFS processed by ext2arps? Is it correct to
set dx = dy = 50000.00 on arpstrn.input to process high resolution terrain
that merges with ext2arps grid?
The mergetrn program has anything about this question? Should I use mergetrn
program in this case? In future, I intend to create nested grid compared
with the resolution grid of the terrain. Can you help me?

Thanks a lot.
Leanderson.


----------------------------------------------------
Prof. Leanderson Marcos da Silva Paiva
Centro Federal de Educa??o Tecnol?gica
Celso Suckow da Fonseca - CEFET/RJ
Coordenadoria de Meteorologia - COMET
Av. Maracan?, n? 229, bloco A, Torre, Maracan?
CEP: 20271-110, Rio de Janeiro, RJ, Brasil
tel.:(55)(0xx21)2566-3163
telefax:(55)(0xx21)2567-0844


From leanderson at cefet-rj.br  Thu Apr 10 09:19:56 2008
From: leanderson at cefet-rj.br (Leanderson Marcos)
Date: Thu, 10 Apr 2008 11:19:56 -0300
Subject: [ARPSSUPPORT] terrain 3sec and GFS 0.5 degree
Message-ID: <HHEIJOJABICFEFAGHOHFCEFOCKAA.leanderson@cefet-rj.br>

Dear friend,
1)	I have 3 sec terrain data file that was processed by arpstrn with these
grid parameters:
nx = 4501
ny = 3501
dx = 90.
dy = 90.

2)	In another way, I have 0.5 degree NCEP/GFS data file in GRIB2, that I
should processes by ext2arps program to create boundary conditions for run
arps program. On the arps.input file processed to ext2arps I have:
nx= 7
ny = 5
dx = 50000.
dy = 50000.

If I try processes ext2arps reading the terrain file created by arpstrn I
generate a ERROR due the not matching grids. So, what can I use high
resolution terrain on 0.5 degree GFS processed by ext2arps? Is it correct to
set dx = dy = 50000.00 on arpstrn.input to process high resolution terrain
that merges with ext2arps grid?
The mergetrn program has anything about this question? Should I use mergetrn
program in this case? In future, I intend to create nested grid compared
with the resolution grid of the terrain. Can you help me?

Thanks a lot.
Leanderson.


----------------------------------------------------
Prof. Leanderson Marcos da Silva Paiva
Centro Federal de Educa??o Tecnol?gica
Celso Suckow da Fonseca - CEFET/RJ
Coordenadoria de Meteorologia - COMET
Av. Maracan?, n? 229, bloco A, Torre, Maracan?
CEP: 20271-110, Rio de Janeiro, RJ, Brasil
tel.:(55)(0xx21)2566-3163
telefax:(55)(0xx21)2567-0844


From mxue at ou.edu  Thu Apr 10 14:15:42 2008
From: mxue at ou.edu (Ming Xue)
Date: Thu, 10 Apr 2008 14:15:42 -0500
Subject: [ARPSSUPPORT] ARPS Help:trouble shooting
In-Reply-To: <47FD3CBB.6080500@student.canterbury.ac.nz>
References: <47FD3CBB.6080500@student.canterbury.ac.nz>
Message-ID: <11c5901c89b3f$45732c80$d0598580$@edu>



Thanks for reporting these problems (seem to be all related to ASCII and
NetCDF I/O
with MPI jobs). We hardly use ASCII format at all because there is no reason
to. 
NetCDF format support is relatively new.  We are in the middle of realtime
experiments so 
we may not be able to address these issues immediately.

We use HDF format almost exclusively and you are strongly encouraged to do
so - the chance to encountering problems is much smaller. Binary is another
option - we enforce bigendian binary across all platforms
so that ARPS binary data are portable.


Ming Xue

> -----Original Message-----
> From: arpssupport-bounces at caps.ou.edu [mailto:arpssupport-
> bounces at caps.ou.edu] On Behalf Of Marwan Katurji
> Sent: Wednesday, April 09, 2008 5:02 PM
> To: arpssupport at ou.edu
> Subject: [ARPSSUPPORT] ARPS Help:trouble shooting
> 
> Hi there,
> 
> I am a registered user for arps and have downloaded version 5.2.8 and
> am
> using it successfully on my PC under cygwin "but not in parallel" for
> boundary layer modeling and LES. Recently arps was installed and
> compiled on a high performance computer, IBM P5-575 containing16 IBM
> Power 5, 1.9 Giga Hertz (GHz) processors and 32 GigaBytes (GB) of
> memory. The purpose is to run it in parallel mode under one node and
> using 8 processors. The arps_mpi and all other utilities I need
> compiled
> fine with no errors an I am using the serial version of arps on one
> processor successfully. For testing the arps_mpi I tried to run the
> demo
> may20 simulation and it appeared to work fine too with 2 and 4
> processors using ascii dump files. The parallel job is submitted to the
> hpc via a load leveler script.
> 
> I have some issues that I need your help with, can't really figure out
> if they are code related or platform:
> 
> 1) The may20 simulation was not successful using 8 processors (nxproc=2,
> nyproc=4 and vice versa) using ascii dump, and the only way the
> simulation ran was using the combination (nxproc=1 and nyproc=8, the
> opposite did not work). Unfortunately I don't get a comprehensive error
> message but the program stops directly after when it starts. Is this
> something related to the platform or the code?
> 
> 2) Using netcdf history dump format the combination of (nxproc=2,
> nyproc=4 and vice versa) worked fine but only when I dumped the output
> from each processor individually with out joining. The join option is
> not working with netcdf and I get a segmentation fault in the load
> leveler error output and an error in DTADUMP in arps output.
> 
> 3) When using a different simulation that involves terrain file,
> surface
> physics and other parameters that you can see in the attached log file,
> the simulation progressed normally using 4 processors (nxproc=2,
> nyproc=2) but suddenly at time =7200s when it was dumping the history
> joined file an error occurred stating "error occurred when opening file
> .//xxxx.asc007200 using fortran unit 39" and when looking at some
> formatted dumped data in the arps.output file they appeared to be
> incomplete, that is some missing lines in the gridded output matrix and
> some redundant ones. (I suspect that the dumping from processors is not
> working correctly).
> 
> When decreasing the processors to 2 (nxproc=1,nxproc=2) the same
> happend
> but occured at a later time, 46800s. You can see a sample of the output
> file attached.
> 
> Thanks for your time and would appreciate your help. This work is part
> of my PhD study at University of Canterbury, New Zealand.
> 
> --
> Marwan Katurji
> 
> PhD Candidate
> Center for Atmospheric Research
> Department of Geography, 2nd floor
> University of Canterbury
> Private Bag 4800
> Christchurch 8020, New Zealand
> -----------------
> Work:  +64 3 364 2987  ext: 7925
> Fax:   +64 3 364 2907
> 
> Email: marwan.houda at gmail.com
>        mka94 at student.canterbury.ac.nz
> 
> Home Address:
> 1/92 Office Road
> Merivale, Christchurch 8014
> New Zealand
> Tel: +64 3 3555760   (land line)
>      +64 3 212049534 (mobile)
> 
> _______________________________________________
> ARPSSUPPORT mailing list
> ARPSSUPPORT at caps.ou.edu
> http://www.caps.ou.edu/mailman/listinfo/arpssupport


From leanderson at cefet-rj.br  Thu Apr 10 14:21:43 2008
From: leanderson at cefet-rj.br (Leanderson Marcos)
Date: Thu, 10 Apr 2008 16:21:43 -0300
Subject: [ARPSSUPPORT] RES:  terrain 3sec and GFS 0.5 degree
In-Reply-To: <47FE2892.8040902@ou.edu>
Message-ID: <HHEIJOJABICFEFAGHOHFOEFPCKAA.leanderson@cefet-rj.br>

That?s ok ... but, How can I nest grids on the run, to reach high resolution
(100 m approximately), using 3 sec terrain data, leaving from 0.5 degree
GFS? How is the numerical fraction beteween the grids that it is recomended?
Have I make runs with differents terrain database for diferents resolutions
grids or Can I use only the 3 sec terrain data for any numerical grid?

Sorry, but I am trying to make this things alone and I have many doubts.

Thanks in advance ...
Leanderson.

-----Mensagem original-----
De: ARPSSUPPORT [mailto:arpssupport at ou.edu]
Enviada em: quinta-feira, 10 de abril de 2008 11:48
Para: leanderson at cefet-rj.br
Assunto: Re: [ARPSSUPPORT] terrain 3sec and GFS 0.5 degree


It is obvious that the grid is not matched. Those 4 parameters should be
always the same values for both programs.

Yunheng.



----- Original Message -----
*From:* Leanderson Marcos <leanderson at cefet-rj.br>
*To:* arpssupport at ou.edu, ARPSSUPPORT at lists.ou.edu, ywang at gcn.ou.edu
*Sent:* 04/10/2008 9:19:56 AM -0500
*Subject:* [ARPSSUPPORT] terrain 3sec and GFS 0.5 degree



> Dear friend,
> 1)	I have 3 sec terrain data file that was processed by arpstrn with these
> grid parameters:
> nx = 4501
> ny = 3501
> dx = 90.
> dy = 90.
>
> 2)	In another way, I have 0.5 degree NCEP/GFS data file in GRIB2, that I
> should processes by ext2arps program to create boundary conditions for run
> arps program. On the arps.input file processed to ext2arps I have:
> nx= 7
> ny = 5
> dx = 50000.
> dy = 50000.
>
> If I try processes ext2arps reading the terrain file created by arpstrn I
> generate a ERROR due the not matching grids. So, what can I use high
> resolution terrain on 0.5 degree GFS processed by ext2arps? Is it correct
to
> set dx = dy = 50000.00 on arpstrn.input to process high resolution terrain
> that merges with ext2arps grid?
> The mergetrn program has anything about this question? Should I use
mergetrn
> program in this case? In future, I intend to create nested grid compared
> with the resolution grid of the terrain. Can you help me?
>
> Thanks a lot.
> Leanderson.
>
>
> ----------------------------------------------------
> Prof. Leanderson Marcos da Silva Paiva
> Centro Federal de Educa??o Tecnol?gica
> Celso Suckow da Fonseca - CEFET/RJ
> Coordenadoria de Meteorologia - COMET
> Av. Maracan?, n? 229, bloco A, Torre, Maracan?
> CEP: 20271-110, Rio de Janeiro, RJ, Brasil
> tel.:(55)(0xx21)2566-3163
> telefax:(55)(0xx21)2567-0844
>
> _______________________________________________
> ARPSSUPPORT mailing list
> ARPSSUPPORT at caps.ou.edu
> http://www.caps.ou.edu/mailman/listinfo/arpssupport
>
>


From mka94 at student.canterbury.ac.nz  Thu Apr 10 15:41:23 2008
From: mka94 at student.canterbury.ac.nz (Marwan Katurji)
Date: Fri, 11 Apr 2008 08:41:23 +1200
Subject: [ARPSSUPPORT] ARPS Help:trouble shooting
In-Reply-To: <11c5901c89b3f$45732c80$d0598580$@edu>
References: <47FD3CBB.6080500@student.canterbury.ac.nz>
	<11c5901c89b3f$45732c80$d0598580$@edu>
Message-ID: <47FE7B73.5060008@student.canterbury.ac.nz>

Hi Ming,

I see the issue now, below is my email to Yunheng about my trouble 
shooting with netcdf output dump maybe it helps in your experiments or 
maybe you can give me a thought on it.
****************
Hi Yunheng,

Until now all ascii dumps under cygwin in serial arps are working 
perfectly fine, good physical results, except that they are large files. 
I am surprised why you would suspect the results. Actually I can't use 
except ascii dump under cygwin in the meant time. And I am still trouble 
shooting on our hpc for the other dump options. The netcdf seems to work 
fine on two processors except that when I go to 4 or 8 I get this 
segmentation error that I told you about previously, I need to check 
your idea if it is a memory problem. With two processors the formatted 
output in the arps.output file are complete and continuous but other 
than that it is discontinuous and incomplete.

***************

I will compile with HDF-4 now and use it.

Many Thanks
 

Marwan Katurji

PhD Candidate
Center for Atmospheric Research
Department of Geography, 2nd floor
University of Canterbury
Private Bag 4800
Christchurch 8020, New Zealand
-----------------
Work:  +64 3 364 2987  ext: 7925
Fax:   +64 3 364 2907

Email: marwan.houda at gmail.com
       mka94 at student.canterbury.ac.nz

Home Address:
1/92 Office Road
Merivale, Christchurch 8014
New Zealand
Tel: +64 3 3555760   (land line)
     +64 3 212049534 (mobile)



Ming Xue wrote:
> Thanks for reporting these problems (seem to be all related to ASCII and
> NetCDF I/O
> with MPI jobs). We hardly use ASCII format at all because there is no reason
> to. 
> NetCDF format support is relatively new.  We are in the middle of realtime
> experiments so 
> we may not be able to address these issues immediately.
>
> We use HDF format almost exclusively and you are strongly encouraged to do
> so - the chance to encountering problems is much smaller. Binary is another
> option - we enforce bigendian binary across all platforms
> so that ARPS binary data are portable.
>
>
> Ming Xue
>
>   
>> -----Original Message-----
>> From: arpssupport-bounces at caps.ou.edu [mailto:arpssupport-
>> bounces at caps.ou.edu] On Behalf Of Marwan Katurji
>> Sent: Wednesday, April 09, 2008 5:02 PM
>> To: arpssupport at ou.edu
>> Subject: [ARPSSUPPORT] ARPS Help:trouble shooting
>>
>> Hi there,
>>
>> I am a registered user for arps and have downloaded version 5.2.8 and
>> am
>> using it successfully on my PC under cygwin "but not in parallel" for
>> boundary layer modeling and LES. Recently arps was installed and
>> compiled on a high performance computer, IBM P5-575 containing16 IBM
>> Power 5, 1.9 Giga Hertz (GHz) processors and 32 GigaBytes (GB) of
>> memory. The purpose is to run it in parallel mode under one node and
>> using 8 processors. The arps_mpi and all other utilities I need
>> compiled
>> fine with no errors an I am using the serial version of arps on one
>> processor successfully. For testing the arps_mpi I tried to run the
>> demo
>> may20 simulation and it appeared to work fine too with 2 and 4
>> processors using ascii dump files. The parallel job is submitted to the
>> hpc via a load leveler script.
>>
>> I have some issues that I need your help with, can't really figure out
>> if they are code related or platform:
>>
>> 1) The may20 simulation was not successful using 8 processors (nxproc=2,
>> nyproc=4 and vice versa) using ascii dump, and the only way the
>> simulation ran was using the combination (nxproc=1 and nyproc=8, the
>> opposite did not work). Unfortunately I don't get a comprehensive error
>> message but the program stops directly after when it starts. Is this
>> something related to the platform or the code?
>>
>> 2) Using netcdf history dump format the combination of (nxproc=2,
>> nyproc=4 and vice versa) worked fine but only when I dumped the output
>> from each processor individually with out joining. The join option is
>> not working with netcdf and I get a segmentation fault in the load
>> leveler error output and an error in DTADUMP in arps output.
>>
>> 3) When using a different simulation that involves terrain file,
>> surface
>> physics and other parameters that you can see in the attached log file,
>> the simulation progressed normally using 4 processors (nxproc=2,
>> nyproc=2) but suddenly at time =7200s when it was dumping the history
>> joined file an error occurred stating "error occurred when opening file
>> .//xxxx.asc007200 using fortran unit 39" and when looking at some
>> formatted dumped data in the arps.output file they appeared to be
>> incomplete, that is some missing lines in the gridded output matrix and
>> some redundant ones. (I suspect that the dumping from processors is not
>> working correctly).
>>
>> When decreasing the processors to 2 (nxproc=1,nxproc=2) the same
>> happend
>> but occured at a later time, 46800s. You can see a sample of the output
>> file attached.
>>
>> Thanks for your time and would appreciate your help. This work is part
>> of my PhD study at University of Canterbury, New Zealand.
>>
>> --
>> Marwan Katurji
>>
>> PhD Candidate
>> Center for Atmospheric Research
>> Department of Geography, 2nd floor
>> University of Canterbury
>> Private Bag 4800
>> Christchurch 8020, New Zealand
>> -----------------
>> Work:  +64 3 364 2987  ext: 7925
>> Fax:   +64 3 364 2907
>>
>> Email: marwan.houda at gmail.com
>>        mka94 at student.canterbury.ac.nz
>>
>> Home Address:
>> 1/92 Office Road
>> Merivale, Christchurch 8014
>> New Zealand
>> Tel: +64 3 3555760   (land line)
>>      +64 3 212049534 (mobile)
>>
>> _______________________________________________
>> ARPSSUPPORT mailing list
>> ARPSSUPPORT at caps.ou.edu
>> http://www.caps.ou.edu/mailman/listinfo/arpssupport
>>     
>
>
>   

From ywang at gcn.ou.edu  Sun Apr 13 11:24:18 2008
From: ywang at gcn.ou.edu (Yunheng Wang)
Date: Sun, 13 Apr 2008 11:24:18 -0500
Subject: [ARPSSUPPORT] ARPS Simulation
In-Reply-To: <000c01c89d20$afba8510$9c87c3de@liziliangouc>
References: <000c01c89d20$afba8510$9c87c3de@liziliangouc>
Message-ID: <480233B2.2070000@gcn.ou.edu>

Sorry, I have no idea about what /fnl_080208_12_00 /is/.

If it is an external forecast data set or analysis data, please check
whether it is NAM or GFS data and what is the NCEP grid number it is
defined on. Then you can check the instruction for variable (extdopt) in
file input/arps.input.

If you did not find an instruction for this data set, then it is not
support by the ARPS system. You can, if you want to, read the ARPS
User's Guide (Chapter 8) about how to write a decoder for a specific
data set yourself.

Yunheng.
/
-------- Original Message --------
Subject: ARPS Simulation
From: liziliang <liziliang at ouc.edu.cn>
To: Wang, Yunheng <yunheng at ou.edu>
Date: 4/12/2008 11:41 PM
> Dear Pro. /wang,/
> //
> / I am a student of Ocean University of Chnia, my name is Zi-Liang Li,
> Would you like to tell me how to convert the fnl_080208_12_00 (from
> the ftpprd.ncep.noaa.gov) into fnl.2008020812f00 (EXT2ARPS format)? /
> / Thanks you for your considering and your time./
> //
> /Sincerely yours,/
> /Zi-Liang, Li /
> /2008.4.12/


From arpssupport at ou.edu  Sun Apr 13 12:01:08 2008
From: arpssupport at ou.edu (arpssupport at ou.edu)
Date: Sun, 13 Apr 2008 12:01:08 -0500 (CDT)
Subject: [ARPSSUPPORT] MensHelp 07220
Message-ID: <20080406071215.4888.qmail@ador47.neoplus.adsl.tpnet.pl>


From mxue at ou.edu  Sun Apr 13 13:14:11 2008
From: mxue at ou.edu (Ming Xue)
Date: Sun, 13 Apr 2008 13:14:11 -0500
Subject: [ARPSSUPPORT] ARPS Simulation
In-Reply-To: <480233B2.2070000@gcn.ou.edu>
References: <000c01c89d20$afba8510$9c87c3de@liziliangouc>
	<480233B2.2070000@gcn.ou.edu>
Message-ID: <BAY115-W499830DC8EC79C80E7E5E7FBE90@phx.gbl>

FNL should be the final analysis of NCEP GFS model, the analysis done after all observations have come in. 
 
They should be 1 x 1 degree data and in the same format as the 1 x 1 degree AVN data set.
 
So the files should be renamed to something like avn_grb.2008041500f03 (not sure of the naming convention).
 
Ming Xue



> Date: Sun, 13 Apr 2008 11:24:18 -0500> From: ywang at gcn.ou.edu> Subject: Re: [ARPSSUPPORT] ARPS Simulation> To: liziliang at ouc.edu.cn; arpssupport at ou.edu> > Sorry, I have no idea about what /fnl_080208_12_00 /is/.> > If it is an external forecast data set or analysis data, please check> whether it is NAM or GFS data and what is the NCEP grid number it is> defined on. Then you can check the instruction for variable (extdopt) in> file input/arps.input.> > If you did not find an instruction for this data set, then it is not> support by the ARPS system. You can, if you want to, read the ARPS> User's Guide (Chapter 8) about how to write a decoder for a specific> data set yourself.> > Yunheng.> /> -------- Original Message --------> Subject: ARPS Simulation> From: liziliang <liziliang at ouc.edu.cn>> To: Wang, Yunheng <yunheng at ou.edu>> Date: 4/12/2008 11:41 PM> > Dear Pro. /wang,/> > //> > / I am a student of Ocean University of Chnia, my name is Zi-Liang Li,> > Would you like to tell me how to convert the fnl_080208_12_00 (from> > the ftpprd.ncep.noaa.gov) into fnl.2008020812f00 (EXT2ARPS format)? /> > / Thanks you for your considering and your time./> > //> > /Sincerely yours,/> > /Zi-Liang, Li /> > /2008.4.12/> 
_________________________________________________________________
Use video conversation to talk face-to-face with Windows Live Messenger.
http://www.windowslive.com/messenger/connect_your_way.html?ocid=TXT_TAGLM_WL_Refresh_messenger_video_042008

From mka94 at student.canterbury.ac.nz  Mon Apr 14 02:10:37 2008
From: mka94 at student.canterbury.ac.nz (Marwan Katurji)
Date: Mon, 14 Apr 2008 19:10:37 +1200
Subject: [ARPSSUPPORT] ARPS Help:trouble shooting
In-Reply-To: <11c5901c89b3f$45732c80$d0598580$@edu>
References: <47FD3CBB.6080500@student.canterbury.ac.nz>
	<11c5901c89b3f$45732c80$d0598580$@edu>
Message-ID: <4803036D.7030404@student.canterbury.ac.nz>

Hi,

Thanks for your help, I am dumping split binary files from 8 processors 
and joining them successfully with joinfiles. But I still can't 
understand how the ascii dump is working fine on a single processor with 
meaningful results whilst you mentioned it's not implemented in version 
5.2.8 anymore!!

Many thanks.

Marwan Katurji

PhD Candidate
Center for Atmospheric Research
Department of Geography, 2nd floor
University of Canterbury
Private Bag 4800
Christchurch 8020, New Zealand
-----------------
Work:  +64 3 364 2987  ext: 7925
Fax:   +64 3 364 2907

Email: marwan.houda at gmail.com
       mka94 at student.canterbury.ac.nz

Home Address:
1/92 Office Road
Merivale, Christchurch 8014
New Zealand
Tel: +64 3 3555760   (land line)
     +64 3 212049534 (mobile)



Ming Xue wrote:
> Thanks for reporting these problems (seem to be all related to ASCII and
> NetCDF I/O
> with MPI jobs). We hardly use ASCII format at all because there is no reason
> to. 
> NetCDF format support is relatively new.  We are in the middle of realtime
> experiments so 
> we may not be able to address these issues immediately.
>
> We use HDF format almost exclusively and you are strongly encouraged to do
> so - the chance to encountering problems is much smaller. Binary is another
> option - we enforce bigendian binary across all platforms
> so that ARPS binary data are portable.
>
>
> Ming Xue
>
>   
>> -----Original Message-----
>> From: arpssupport-bounces at caps.ou.edu [mailto:arpssupport-
>> bounces at caps.ou.edu] On Behalf Of Marwan Katurji
>> Sent: Wednesday, April 09, 2008 5:02 PM
>> To: arpssupport at ou.edu
>> Subject: [ARPSSUPPORT] ARPS Help:trouble shooting
>>
>> Hi there,
>>
>> I am a registered user for arps and have downloaded version 5.2.8 and
>> am
>> using it successfully on my PC under cygwin "but not in parallel" for
>> boundary layer modeling and LES. Recently arps was installed and
>> compiled on a high performance computer, IBM P5-575 containing16 IBM
>> Power 5, 1.9 Giga Hertz (GHz) processors and 32 GigaBytes (GB) of
>> memory. The purpose is to run it in parallel mode under one node and
>> using 8 processors. The arps_mpi and all other utilities I need
>> compiled
>> fine with no errors an I am using the serial version of arps on one
>> processor successfully. For testing the arps_mpi I tried to run the
>> demo
>> may20 simulation and it appeared to work fine too with 2 and 4
>> processors using ascii dump files. The parallel job is submitted to the
>> hpc via a load leveler script.
>>
>> I have some issues that I need your help with, can't really figure out
>> if they are code related or platform:
>>
>> 1) The may20 simulation was not successful using 8 processors (nxproc=2,
>> nyproc=4 and vice versa) using ascii dump, and the only way the
>> simulation ran was using the combination (nxproc=1 and nyproc=8, the
>> opposite did not work). Unfortunately I don't get a comprehensive error
>> message but the program stops directly after when it starts. Is this
>> something related to the platform or the code?
>>
>> 2) Using netcdf history dump format the combination of (nxproc=2,
>> nyproc=4 and vice versa) worked fine but only when I dumped the output
>> from each processor individually with out joining. The join option is
>> not working with netcdf and I get a segmentation fault in the load
>> leveler error output and an error in DTADUMP in arps output.
>>
>> 3) When using a different simulation that involves terrain file,
>> surface
>> physics and other parameters that you can see in the attached log file,
>> the simulation progressed normally using 4 processors (nxproc=2,
>> nyproc=2) but suddenly at time =7200s when it was dumping the history
>> joined file an error occurred stating "error occurred when opening file
>> .//xxxx.asc007200 using fortran unit 39" and when looking at some
>> formatted dumped data in the arps.output file they appeared to be
>> incomplete, that is some missing lines in the gridded output matrix and
>> some redundant ones. (I suspect that the dumping from processors is not
>> working correctly).
>>
>> When decreasing the processors to 2 (nxproc=1,nxproc=2) the same
>> happend
>> but occured at a later time, 46800s. You can see a sample of the output
>> file attached.
>>
>> Thanks for your time and would appreciate your help. This work is part
>> of my PhD study at University of Canterbury, New Zealand.
>>
>> --
>> Marwan Katurji
>>
>> PhD Candidate
>> Center for Atmospheric Research
>> Department of Geography, 2nd floor
>> University of Canterbury
>> Private Bag 4800
>> Christchurch 8020, New Zealand
>> -----------------
>> Work:  +64 3 364 2987  ext: 7925
>> Fax:   +64 3 364 2907
>>
>> Email: marwan.houda at gmail.com
>>        mka94 at student.canterbury.ac.nz
>>
>> Home Address:
>> 1/92 Office Road
>> Merivale, Christchurch 8014
>> New Zealand
>> Tel: +64 3 3555760   (land line)
>>      +64 3 212049534 (mobile)
>>
>> _______________________________________________
>> ARPSSUPPORT mailing list
>> ARPSSUPPORT at caps.ou.edu
>> http://www.caps.ou.edu/mailman/listinfo/arpssupport
>>     
>
>
>   

From ywang at gcn.ou.edu  Mon Apr 14 09:03:57 2008
From: ywang at gcn.ou.edu (Yunheng Wang)
Date: Mon, 14 Apr 2008 09:03:57 -0500
Subject: [ARPSSUPPORT] ARPS Simulation
In-Reply-To: <BAY115-W499830DC8EC79C80E7E5E7FBE90@phx.gbl>
References: <000c01c89d20$afba8510$9c87c3de@liziliangouc>
	<480233B2.2070000@gcn.ou.edu>
	<BAY115-W499830DC8EC79C80E7E5E7FBE90@phx.gbl>
Message-ID: <4803644D.2070200@gcn.ou.edu>

Hi Ziliang,

If your data is what Dr Xue described below, please do the following to 
ingest the data into ext2arps.

1. Download the ARPS package and compile ext2arps;
2. Rename your file to something like,  gnl.YYYYMMDDHHfhh. For example, 
fnl.2008041500f03.
3. set the following parameters in the namelist file
    extdname = 'fnl',
    nextdfil = 1,
    extdtime(1) = '2008-04-15.00:00:00+003:00:00',
    ....
    extdopt = 13,

    and other parameters.
4. Run ext2arps.

Yunheng.



----- Original Message -----
*From:* Ming Xue <mxue at ou.edu>
*To:* Yunheng Wang <ywang at gcn.ou.edu>, liziliang <liziliang at ouc.edu.cn>, 
arpssupport at ou.edu
*Sent:* 04/13/2008 1:14:11 PM -0500
*Subject:* [ARPSSUPPORT] ARPS Simulation



>
> FNL should be the final analysis of NCEP GFS model, the analysis done 
> after all observations have come in.
>  
> They should be 1 x 1 degree data and in the same format as the 1 x 1 
> degree AVN data set.
>  
> So the files should be renamed to something like avn_grb.2008041500f03 
> (not sure of the naming convention).
>  
> Ming Xue
>
>
>
> ------------------------------------------------------------------------
>
> > Date: Sun, 13 Apr 2008 11:24:18 -0500
> > From: ywang at gcn.ou.edu
> > Subject: Re: [ARPSSUPPORT] ARPS Simulation
> > To: liziliang at ouc.edu.cn; arpssupport at ou.edu
> >
> > Sorry, I have no idea about what /fnl_080208_12_00 /is/.
> >
> > If it is an external forecast data set or analysis data, please check
> > whether it is NAM or GFS data and what is the NCEP grid number it is
> > defined on. Then you can check the instruction for variable (extdopt) in
> > file input/arps.input.
> >
> > If you did not find an instruction for this data set, then it is not
> > support by the ARPS system. You can, if you want to, read the ARPS
> > User's Guide (Chapter 8) about how to write a decoder for a specific
> > data set yourself.
> >
> > Yunheng.
> > /
> > -------- Original Message --------
> > Subject: ARPS Simulation
> > From: liziliang <liziliang at ouc.edu.cn>
> > To: Wang, Yunheng <yunheng at ou.edu>
> > Date: 4/12/2008 11:41 PM
> > > Dear Pro. /wang,/
> > > //
> > > / I am a student of Ocean University of Chnia, my name is Zi-Liang Li,
> > > Would you like to tell me how to convert the fnl_080208_12_00 (from
> > > the ftpprd.ncep.noaa.gov) into fnl.2008020812f00 (EXT2ARPS format)? /
> > > / Thanks you for your considering and your time./
> > > //
> > > /Sincerely yours,/
> > > /Zi-Liang, Li /
> > > /2008.4.12/
> >
>
>
> ------------------------------------------------------------------------
> Use video conversation to talk face-to-face with Windows Live 
> Messenger. Get started! 
> <http://www.windowslive.com/messenger/connect_your_way.html?ocid=TXT_TAGLM_WL_Refresh_messenger_video_042008>

From ywang at gcn.ou.edu  Mon Apr 14 08:43:59 2008
From: ywang at gcn.ou.edu (Yunheng Wang)
Date: Mon, 14 Apr 2008 08:43:59 -0500
Subject: [ARPSSUPPORT] ARPS Help:trouble shooting
In-Reply-To: <4803036D.7030404@student.canterbury.ac.nz>
References: <47FD3CBB.6080500@student.canterbury.ac.nz>
	<11c5901c89b3f$45732c80$d0598580$@edu>
	<4803036D.7030404@student.canterbury.ac.nz>
Message-ID: <48035F9F.7070101@gcn.ou.edu>

We have changed the ARPS history format since version 5.0.0IHOP_3. 
However, the ASCII format was not updated since lack of usage. 
Furthermore, the ASCII format is suitable for small grid only, which is 
usually a serial job.

So even you got results with serial runs, file contents are not 
guaranteed compatible with the current ARPS system. It can be used for 
diagnostic purpose, but cannot be ingested to the ARPS system anymore. 
ACII format is never tested in parallel mode and it is also not intended 
to be used in parallel mode.

Hope this information clarify things.

Yunheng.


----- Original Message -----
*From:* Marwan Katurji <mka94 at student.canterbury.ac.nz>
*To:* "Xue, Ming" <mxue at ou.edu>, "arpssupport at ou.edu" 
<arpssupport at ou.edu>, "Wang, Yunheng" <yunheng at ou.edu>
*CC:* Peyman Zawar-Reza <peyman.zawar-reza at canterbury.ac.nz>
*Sent:* 04/14/2008 2:10:37 AM -0500
*Subject:* [ARPSSUPPORT] ARPS Help:trouble shooting



> Hi,
>
> Thanks for your help, I am dumping split binary files from 8 processors
> and joining them successfully with joinfiles. But I still can't
> understand how the ascii dump is working fine on a single processor with
> meaningful results whilst you mentioned it's not implemented in version
> 5.2.8 anymore!!
>
> Many thanks.
>
> Marwan Katurji
>
> PhD Candidate
> Center for Atmospheric Research
> Department of Geography, 2nd floor
> University of Canterbury
> Private Bag 4800
> Christchurch 8020, New Zealand
> -----------------
> Work:  +64 3 364 2987  ext: 7925
> Fax:   +64 3 364 2907
>
> Email: marwan.houda at gmail.com
>        mka94 at student.canterbury.ac.nz
>
> Home Address:
> 1/92 Office Road
> Merivale, Christchurch 8014
> New Zealand
> Tel: +64 3 3555760   (land line)
>      +64 3 212049534 (mobile)
>
>
>
> Ming Xue wrote:
>   
>> Thanks for reporting these problems (seem to be all related to ASCII and
>> NetCDF I/O
>> with MPI jobs). We hardly use ASCII format at all because there is no reason
>> to.
>> NetCDF format support is relatively new.  We are in the middle of realtime
>> experiments so
>> we may not be able to address these issues immediately.
>>
>> We use HDF format almost exclusively and you are strongly encouraged to do
>> so - the chance to encountering problems is much smaller. Binary is another
>> option - we enforce bigendian binary across all platforms
>> so that ARPS binary data are portable.
>>
>>
>> Ming Xue
>>
>>
>>     
>>> -----Original Message-----
>>> From: arpssupport-bounces at caps.ou.edu [mailto:arpssupport-
>>> bounces at caps.ou.edu] On Behalf Of Marwan Katurji
>>> Sent: Wednesday, April 09, 2008 5:02 PM
>>> To: arpssupport at ou.edu
>>> Subject: [ARPSSUPPORT] ARPS Help:trouble shooting
>>>
>>> Hi there,
>>>
>>> I am a registered user for arps and have downloaded version 5.2.8 and
>>> am
>>> using it successfully on my PC under cygwin "but not in parallel" for
>>> boundary layer modeling and LES. Recently arps was installed and
>>> compiled on a high performance computer, IBM P5-575 containing16 IBM
>>> Power 5, 1.9 Giga Hertz (GHz) processors and 32 GigaBytes (GB) of
>>> memory. The purpose is to run it in parallel mode under one node and
>>> using 8 processors. The arps_mpi and all other utilities I need
>>> compiled
>>> fine with no errors an I am using the serial version of arps on one
>>> processor successfully. For testing the arps_mpi I tried to run the
>>> demo
>>> may20 simulation and it appeared to work fine too with 2 and 4
>>> processors using ascii dump files. The parallel job is submitted to the
>>> hpc via a load leveler script.
>>>
>>> I have some issues that I need your help with, can't really figure out
>>> if they are code related or platform:
>>>
>>> 1) The may20 simulation was not successful using 8 processors (nxproc=2,
>>> nyproc=4 and vice versa) using ascii dump, and the only way the
>>> simulation ran was using the combination (nxproc=1 and nyproc=8, the
>>> opposite did not work). Unfortunately I don't get a comprehensive error
>>> message but the program stops directly after when it starts. Is this
>>> something related to the platform or the code?
>>>
>>> 2) Using netcdf history dump format the combination of (nxproc=2,
>>> nyproc=4 and vice versa) worked fine but only when I dumped the output
>>> from each processor individually with out joining. The join option is
>>> not working with netcdf and I get a segmentation fault in the load
>>> leveler error output and an error in DTADUMP in arps output.
>>>
>>> 3) When using a different simulation that involves terrain file,
>>> surface
>>> physics and other parameters that you can see in the attached log file,
>>> the simulation progressed normally using 4 processors (nxproc=2,
>>> nyproc=2) but suddenly at time =7200s when it was dumping the history
>>> joined file an error occurred stating "error occurred when opening file
>>> .//xxxx.asc007200 using fortran unit 39" and when looking at some
>>> formatted dumped data in the arps.output file they appeared to be
>>> incomplete, that is some missing lines in the gridded output matrix and
>>> some redundant ones. (I suspect that the dumping from processors is not
>>> working correctly).
>>>
>>> When decreasing the processors to 2 (nxproc=1,nxproc=2) the same
>>> happend
>>> but occured at a later time, 46800s. You can see a sample of the output
>>> file attached.
>>>
>>> Thanks for your time and would appreciate your help. This work is part
>>> of my PhD study at University of Canterbury, New Zealand.
>>>
>>> --
>>> Marwan Katurji
>>>
>>> PhD Candidate
>>> Center for Atmospheric Research
>>> Department of Geography, 2nd floor
>>> University of Canterbury
>>> Private Bag 4800
>>> Christchurch 8020, New Zealand
>>> -----------------
>>> Work:  +64 3 364 2987  ext: 7925
>>> Fax:   +64 3 364 2907
>>>
>>> Email: marwan.houda at gmail.com
>>>        mka94 at student.canterbury.ac.nz
>>>
>>> Home Address:
>>> 1/92 Office Road
>>> Merivale, Christchurch 8014
>>> New Zealand
>>> Tel: +64 3 3555760   (land line)
>>>      +64 3 212049534 (mobile)
>>>
>>> _______________________________________________
>>> ARPSSUPPORT mailing list
>>> ARPSSUPPORT at caps.ou.edu
>>> http://www.caps.ou.edu/mailman/listinfo/arpssupport
>>>
>>>       
>>
>>     

From mka94 at student.canterbury.ac.nz  Tue Apr 15 02:27:31 2008
From: mka94 at student.canterbury.ac.nz (Marwan Katurji)
Date: Tue, 15 Apr 2008 19:27:31 +1200
Subject: [ARPSSUPPORT] simulation question
Message-ID: <480458E3.5060609@student.canterbury.ac.nz>

Hi there,

I was trying a simulation in arps over a flat terrain with sounding 
initialization. The sounding file is a stable atmospheric profile and 
the atmosphere is initially at rest. I started from the may20 simulation 
and removed the potential temperature perturbation option and added user 
input soil model initializations and soil data, the model simulation 
finished ok. The next step was introducing a real surface file through 
arpssfc and that was successful with the creation of the may20.sfcdata, 
but when running the model with this surface file I got the following 
message and it kept repeating on all grid points:

 Warning: The PBL is found to extend the entire model domain at i,j=  29  66
 The atmosphere is either neutral or unstable throughout the domain depth.
 In this case, PBL parameterization, if used, will operate on the entire 
domain depth.

 Warning: The PBL is found to extend the entire model domain at i,j=  30  66
 The atmosphere is either neutral or unstable throughout the domain depth.
 In this case, PBL parameterization, if used, will operate on the entire 
domain depth.

 Warning: The PBL is found to extend the entire model domain at i,j=  31  66
 The atmosphere is either neutral or unstable throughout the domain depth.
 In this case, PBL parameterization, if used, will operate on the entire 
domain depth.

After that I changed the option for PBL depth determination from 
pbldopt=2 to 1 and then the simulation started and went on fine. Do you 
have an idea on why this is happening? I want my planetary boundary 
layer depth diagnosed and not set initially.

Many thanks for your help, again!

-- 
Marwan Katurji

PhD Candidate
Center for Atmospheric Research
Department of Geography, 2nd floor
University of Canterbury
Private Bag 4800
Christchurch 8020, New Zealand
-----------------
Work:  +64 3 364 2987  ext: 7925
Fax:   +64 3 364 2907

Email: marwan.houda at gmail.com
       mka94 at student.canterbury.ac.nz

Home Address:
1/92 Office Road
Merivale, Christchurch 8014
New Zealand
Tel: +64 3 3555760   (land line)
     +64 3 212049534 (mobile)


From mxue at ou.edu  Tue Apr 15 09:18:48 2008
From: mxue at ou.edu (Ming Xue)
Date: Tue, 15 Apr 2008 09:18:48 -0500
Subject: [ARPSSUPPORT] simulation question
In-Reply-To: <480458E3.5060609@student.canterbury.ac.nz>
References: <480458E3.5060609@student.canterbury.ac.nz>
Message-ID: <027a01c89f03$9f522c20$ddf68460$@edu>


That usually means that the soil model skin temperature went bad.
You should carefully check the soil/veg type data, as well as their 
the soil model initial condition. 


Ming Xue

> -----Original Message-----
> From: arpssupport-bounces at caps.ou.edu [mailto:arpssupport-
> bounces at caps.ou.edu] On Behalf Of Marwan Katurji
> Sent: Tuesday, April 15, 2008 2:28 AM
> To: arpssupport at ou.edu
> Cc: Peyman Zawar-Reza
> Subject: [ARPSSUPPORT] simulation question
> 
> Hi there,
> 
> I was trying a simulation in arps over a flat terrain with sounding
> initialization. The sounding file is a stable atmospheric profile and
> the atmosphere is initially at rest. I started from the may20
> simulation
> and removed the potential temperature perturbation option and added
> user
> input soil model initializations and soil data, the model simulation
> finished ok. The next step was introducing a real surface file through
> arpssfc and that was successful with the creation of the may20.sfcdata,
> but when running the model with this surface file I got the following
> message and it kept repeating on all grid points:
> 
>  Warning: The PBL is found to extend the entire model domain at i,j=
> 29  66
>  The atmosphere is either neutral or unstable throughout the domain
> depth.
>  In this case, PBL parameterization, if used, will operate on the
> entire
> domain depth.
> 
>  Warning: The PBL is found to extend the entire model domain at i,j=
> 30  66
>  The atmosphere is either neutral or unstable throughout the domain
> depth.
>  In this case, PBL parameterization, if used, will operate on the
> entire
> domain depth.
> 
>  Warning: The PBL is found to extend the entire model domain at i,j=
> 31  66
>  The atmosphere is either neutral or unstable throughout the domain
> depth.
>  In this case, PBL parameterization, if used, will operate on the
> entire
> domain depth.
> 
> After that I changed the option for PBL depth determination from
> pbldopt=2 to 1 and then the simulation started and went on fine. Do you
> have an idea on why this is happening? I want my planetary boundary
> layer depth diagnosed and not set initially.
> 
> Many thanks for your help, again!
> 
> --
> Marwan Katurji
> 
> PhD Candidate
> Center for Atmospheric Research
> Department of Geography, 2nd floor
> University of Canterbury
> Private Bag 4800
> Christchurch 8020, New Zealand
> -----------------
> Work:  +64 3 364 2987  ext: 7925
> Fax:   +64 3 364 2907
> 
> Email: marwan.houda at gmail.com
>        mka94 at student.canterbury.ac.nz
> 
> Home Address:
> 1/92 Office Road
> Merivale, Christchurch 8014
> New Zealand
> Tel: +64 3 3555760   (land line)
>      +64 3 212049534 (mobile)
> 
> _______________________________________________
> ARPSSUPPORT mailing list
> ARPSSUPPORT at caps.ou.edu
> http://www.caps.ou.edu/mailman/listinfo/arpssupport


From bernard.dudon at univ-ag.fr  Tue Apr 15 11:41:11 2008
From: bernard.dudon at univ-ag.fr (Bernard Dudon)
Date: Tue, 15 Apr 2008 12:41:11 -0400
Subject: [ARPSSUPPORT] using ext2arps with a subregion of NCEP GFS data
Message-ID: <8D0FFF65-3CCD-4806-B8ED-02CECA13FFF6@univ-ag.fr>

Dear sir

I'd like to initialize arps with a subregion of NCEP GFS forecast (1  
deg) to reduce the amount of data to download. (available at http://nomad5.ncep.noaa.gov/cgi-bin/ftp2u_gfs_dir.sh)

ext2arps doesn't work with these data , I think because of the  
dimension of external data grid nx_ext and ny_ext  in extract_avn.f90  
which values are 360 and 180.

Is there an issue to use a subregion of GFS forecast ?

I have another question : we'd like to print the track of a hurricane  
(one of the goals of our laboratory in Guadeloupe French West Indies).  
Have you developped a program to do that ?

The third things is a suggestion : I have modified arpstrn.f90 and  
arpstrn.input to use SRTM (Shuttle Radar Topography Mission) 3" data  
available for all the world from 60N to 56S .

Are you interested to introduce in your next release of arps and do  
you want I send you a copy to test ?

Have a look at our site (http://calamar.univ-ag.fr/uag/physique/lpat/html/index.php 
  : "Simulation ARPS" ) we have a four times daily forecast with arps  
in atlantic area and our region. For the moment only initialized with  
the analyse data of GFS and not with also the forecast files. It's the  
reason why we'd like to use subregion in place of global file.

Thank you for your help.

Best regards

Bernard


Bernard Dudon
Universite Antilles Guyane
UFR Sciences Exactes et Naturelles
Laboratoire de Physique de l'Atmosphere Tropicale
BP 592
97157 Pointe a Pitre Cedex
Guadeloupe (FWI)
Fax 05 90  48 31 05
Phone 05 90 48 31 07
International :
Fax +590 590 48 31 05
Phone +590 590 48 31 07
http://calamar.univ-ag.fr/uag/physique/lpat/html/index.php



From ywang at gcn.ou.edu  Tue Apr 15 11:56:27 2008
From: ywang at gcn.ou.edu (Yunheng Wang)
Date: Tue, 15 Apr 2008 11:56:27 -0500
Subject: [ARPSSUPPORT] test from GCN account
Message-ID: <4804DE3B.5060304@gcn.ou.edu>

test

From ywang at gcn.ou.edu  Tue Apr 15 22:14:57 2008
From: ywang at gcn.ou.edu (Yunheng Wang)
Date: Tue, 15 Apr 2008 22:14:57 -0500
Subject: [ARPSSUPPORT] [Fwd: Re: [Fwd: Questions about zxplot in ARPS]]
Message-ID: <48056F31.9080305@gcn.ou.edu>

Gordon,

I am forwarding your message to ARPSSUPPORT at caps.ou.edu to seek for an 
answer.

Yunheng.

-------- Original Message --------
Subject: 	Re: [Fwd: Questions about zxplot in ARPS]
Date: 	Wed, 16 Apr 2008 08:51:00 +0800
From: 	Gordon <gordonpann at gmail.com>
To: 	Yunheng Wang <ywang at gcn.ou.edu>
References: 	<45B3A10D.3070400 at gcn.ou.edu> <479623FE.2000205 at gcn.ou.edu>



 
Hi, Dr. Wang,
 
I used ARPS5.2.6 to simulate the wind fields for South China and the tke 
output is very important in my research, but the results seem a little 
weird. The attachment is a graph for the tke profile on 2004.10.12. The 
profile has two maximums and a minimum unexpectedly appears at 200m. 
This problem came out in each day result I simulated.  I don't know why, 
could you give me some suggestions? Is there anything wrong in my input 
file?
 
Thank you for your help.
 
Regards,
 
Gordon
 


 
2008/1/23, Yunheng Wang <ywang at gcn.ou.edu <mailto:ywang at gcn.ou.edu>>:

    Hi Gordon,

    There are several problems with your map data.

    First, zxplot expects data in format 'F8.3', but some of your data
    is in 'F9.3' format. That is why it messed up with zxplot.

    Second, there are many strange lines in your file, for example,

         0       2   0.0008000.000   0.0008000.000

    They should be removed or reformatted.

    Finally, some headers for each block has format 'F7.3'. It also
    confuses zxplot and they should all be 'F8.3'.

    I have a corrected copy attached for your convenience. BTW, I will
    include this data file in the official release later if you do not
    object.

    Yunheng.




    ----- Original Message -----
    *From:* Yunheng Wang <ywang at gcn.ou.edu <mailto:ywang at gcn.ou.edu>>
    *To:* ARPSSUPPORT at lists.ou.edu <mailto:ARPSSUPPORT at lists.ou.edu>
    *Sent:* 01/21/2007 11:21:17 AM -0600
    *Subject:* [Fwd: Questions about zxplot in ARPS]






        -------- Original Message --------
        Subject: Questions about zxplot in ARPS
        Date: Fri, 18 Jan 2008 14:48:52 +0800
        From: Gordon <gordonpann at gmail.com <mailto:gordonpann at gmail.com>>
        To: Yunheng Wang <ywang at gcn.ou.edu <mailto:ywang at gcn.ou.edu>>

        Dear Dr. Wang,

        The version of ARPS in my system is 5.2.6. <http://5.2.6./>
        Because the world mapdata
        in ARPS is not good enough to show Chinese political boundary, so I
        make a china.mapdata, which is changed from the .bln data I used in
        Surfer with the same format as your mapdata. But when I used
        arpspltpost to plot the map, I get this error "Dead loop encountered
        in XTSTWD, job stopped". Is there any additional format in your
        mapdata or do you have any ideas to solve this problem? Thank
        you for
        your help.

        The accessory is my china.mapdata, which includes the provincial
        boundary. Is there anything wrong in it?

        Thank you for your help.

        Gordon


    CHINA POLITICAL BOUNDARIES.
    Header: Number_of_lat_lon_data_values, group_id, max, min lat, max,
    min lon of this segment
        200       2  34.646  34.150 117.072 116.369
     34.644 116.369  34.643 116.383  34.643 116.395  34.646 116.406
     34.646 116.418
     34.640 116.428  34.632 116.432  34.623 116.436  34.616 116.444
     34.613 116.455
     34.610 116.466  34.603 116.473  34.593 116.475  34.585 116.480
     34.576 116.483
     34.567 116.488  34.559 116.495  34.553 116.504  34.550 116.515
     34.547 116.526
     34.544 116.537  34.539 116.546  34.531 116.554  34.522 116.558
     34.517 116.568
     34.512 116.578  34.504 116.582  34.500 116.572  34.496 116.561
     34.489 116.554
     34.483 116.561  34.480 116.572  34.478 116.583  34.479 116.595
     34.480 116.607
     34.479 116.619  34.480 116.630  34.479 116.642  34.478 116.654
     34.478 116.666
     34.475 116.677  34.473 116.688  34.471 116.700  34.472 116.712
     34.472 116.724
     34.470 116.735  34.466 116.746  34.459 116.754  34.451 116.760
     34.443 116.766
     34.436 116.774  34.428 116.781  34.420 116.790  34.414 116.799
     34.408 116.810
     34.400 116.816  34.395 116.826  34.393 116.837  34.392 116.849
     34.394 116.860
     34.398 116.871  34.401 116.882  34.403 116.893  34.405 116.905
     34.404 116.917
     34.403 116.928  34.401 116.940  34.394 116.949  34.387 116.956
     34.377 116.958
     34.368 116.960  34.358 116.962  34.349 116.966  34.340 116.963
     34.331 116.958
     34.321 116.959  34.312 116.962  34.305 116.968  34.295 116.969
     34.286 116.969
     34.278 116.975  34.271 116.984  34.265 116.993  34.260 117.003
     34.255 117.013
     34.251 117.023  34.245 117.032  34.238 117.041  34.230 117.041
     34.221 117.035
     34.213 117.029  34.205 117.024  34.196 117.020  34.186 117.020
     34.178 117.024
     34.171 117.033  34.165 117.043  34.161 117.053  34.156 117.062
     34.150 117.072
        200       2  34.143  33.769 117.749 117.080
     34.143 117.080  34.137 117.090  34.131 117.099  34.127 117.109
     34.120 117.118
     34.110 117.123  34.105 117.132  34.097 117.138  34.088 117.143
     34.082 117.152
     34.078 117.162  34.075 117.173  34.074 117.185  34.072 117.196
     34.071 117.208
     34.070 117.220  34.069 117.232  34.072 117.242  34.077 117.252
     34.078 117.264
     34.074 117.274  34.070 117.285  34.068 117.295  34.074 117.305
     34.082 117.311
     34.089 117.319  34.092 117.330  34.090 117.342  34.085 117.352
     34.081 117.362
     34.074 117.370  34.066 117.378  34.058 117.382  34.049 117.386
     34.040 117.390
     34.030 117.395  34.027 117.406  34.029 117.417  34.031 117.428
     34.035 117.438
     34.038 117.449  34.041 117.460  34.044 117.471  34.047 117.482
     34.052 117.491
     34.057 117.501  34.061 117.511  34.053 117.519  34.046 117.528
     34.038 117.534
     34.028 117.537  34.019 117.534  34.010 117.537  34.002 117.544
     33.999 117.555
     33.990 117.557  33.981 117.562  33.983 117.572  33.987 117.583
     33.989 117.594
     33.997 117.599  34.004 117.603  34.013 117.602  34.024 117.602
     34.031 117.608
     34.029 117.619  34.025 117.630  34.019 117.639  34.011 117.646
     34.003 117.652
     34.000 117.659  33.991 117.662  33.981 117.662  33.972 117.659
     33.962 117.657
     33.953 117.660  33.944 117.663  33.935 117.662  33.927 117.668
     33.921 117.678
     33.915 117.686  33.906 117.690  33.897 117.693  33.889 117.701
     33.884 117.710
     33.884 117.721  33.887 117.732  33.886 117.744  33.877 117.748
     33.867 117.749
     33.857 117.747  33.847 117.745  33.837 117.744  33.828 117.743
     33.817 117.743
     33.807 117.743  33.798 117.745  33.788 117.744  33.778 117.743
     33.769 117.739
        200       2  33.762  33.349 118.164 117.710
     33.760 117.735  33.752 117.729  33.746 117.720  33.739 117.710
     33.730 117.713
     33.723 117.722  33.718 117.731  33.717 117.743  33.719 117.754
     33.722 117.765
     33.726 117.775  33.730 117.785  33.734 117.796  33.735 117.807
     33.735 117.819
     33.734 117.831  33.732 117.842  33.729 117.853  33.726 117.863
     33.723 117.874
     33.719 117.885  33.726 117.892  33.732 117.902  33.734 117.913
     33.732 117.924
     33.730 117.936  33.736 117.946  33.745 117.948  33.755 117.950
     33.760 117.960
     33.759 117.972  33.756 117.982  33.754 117.993  33.748 118.000
     33.746 118.010
     33.748 118.021  33.749 118.033  33.753 118.043  33.756 118.054
     33.760 118.064
     33.762 118.075  33.762 118.087  33.760 118.098  33.755 118.108
     33.748 118.116
     33.743 118.126  33.739 118.136  33.739 118.148  33.736 118.155
     33.729 118.164
     33.723 118.155  33.714 118.152  33.704 118.152  33.695 118.157
     33.686 118.159
     33.676 118.160  33.667 118.157  33.658 118.153  33.649 118.147
     33.643 118.139
     33.636 118.130  33.629 118.122  33.620 118.117  33.612 118.113
     33.602 118.110
     33.593 118.109  33.583 118.107  33.573 118.106  33.563 118.106
     33.553 118.105
     33.543 118.103  33.534 118.100  33.525 118.097  33.516 118.094
     33.506 118.093
     33.496 118.094  33.486 118.092  33.479 118.085  33.475 118.074
     33.479 118.064
     33.484 118.054  33.482 118.043  33.473 118.039  33.464 118.036
     33.455 118.033
     33.446 118.029  33.437 118.026  33.427 118.024  33.417 118.023
     33.407 118.024
     33.397 118.024  33.387 118.025  33.377 118.023  33.369 118.019
     33.360 118.015
     33.351 118.012  33.349 118.001  33.351 117.989  33.352 117.977
     33.349 117.966
        200       2  33.341  32.723 118.395 117.929
     33.341 117.962  33.333 117.965  33.326 117.972  33.315 117.972
     33.306 117.969
     33.298 117.963  33.289 117.958  33.281 117.955  33.274 117.946
     33.268 117.938
     33.258 117.937  33.248 117.935  33.240 117.929  33.231 117.930
     33.226 117.940
     33.223 117.950  33.221 117.962  33.213 117.969  33.203 117.968
     33.192 117.968
     33.183 117.970  33.173 117.972  33.165 117.977  33.157 117.984
     33.152 117.994
     33.146 118.003  33.140 118.012  33.135 118.022  33.134 118.034
     33.135 118.045
     33.138 118.056  33.141 118.069  33.150 118.080  33.152 118.091
     33.152 118.104
     33.151 118.116  33.153 118.127  33.157 118.137  33.161 118.147
     33.168 118.155
     33.177 118.158  33.187 118.160  33.195 118.165  33.202 118.173
     33.203 118.185
     33.200 118.196  33.197 118.206  33.194 118.217  33.191 118.228
     33.188 118.238
     33.158 118.227  33.121 118.228  33.069 118.233  33.031 118.239
     32.995 118.270
     32.975 118.290  32.949 118.278  32.946 118.268  32.940 118.258
     32.936 118.248
     32.930 118.239  32.922 118.235  32.912 118.233  32.902 118.233
     32.892 118.235
     32.883 118.240  32.874 118.241  32.864 118.239  32.854 118.239
     32.849 118.249
     32.852 118.259  32.853 118.270  32.846 118.279  32.839 118.286
     32.836 118.297
     32.835 118.309  32.831 118.319  32.825 118.328  32.819 118.337
     32.816 118.348
     32.816 118.359  32.813 118.370  32.804 118.366  32.799 118.356
     32.795 118.346
     32.792 118.335  32.792 118.324  32.787 118.314  32.778 118.313
     32.773 118.322
     32.771 118.333  32.769 118.344  32.765 118.354  32.760 118.364
     32.752 118.370
     32.742 118.371  32.732 118.371  32.723 118.374  32.726 118.384
     32.728 118.395
        200       2  32.959  32.708 119.112 118.396
     32.722 118.400  32.713 118.396  32.708 118.403  32.709 118.414
     32.714 118.418
     32.721 118.425  32.726 118.435  32.725 118.447  32.722 118.458
     32.722 118.469
     32.720 118.481  32.719 118.492  32.720 118.503  32.721 118.515
     32.721 118.527
     32.719 118.539  32.719 118.551  32.718 118.563  32.719 118.574
     32.720 118.586
     32.725 118.595  32.731 118.603  32.735 118.614  32.737 118.625
     32.740 118.635
     32.736 118.646  32.731 118.656  32.733 118.667  32.731 118.678
     32.727 118.688
     32.719 118.695  32.711 118.703  32.719 118.716  32.726 118.724
     32.731 118.734
     32.737 118.742  32.746 118.746  32.756 118.743  32.764 118.737
     32.772 118.733
     32.782 118.730  32.792 118.730  32.801 118.734  32.809 118.739
     32.819 118.737
     32.828 118.734  32.838 118.733  32.845 118.740  32.850 118.749
     32.852 118.761
     32.849 118.772  32.850 118.783  32.853 118.794  32.857 118.804
     32.866 118.803
     32.875 118.802  32.884 118.807  32.892 118.809  32.902 118.811
     32.911 118.813
     32.916 118.822  32.916 118.834  32.923 118.840  32.933 118.842
     32.943 118.844
     32.952 118.846  32.958 118.853  32.958 118.864  32.958 118.875
     32.959 118.887
     32.954 118.896  32.945 118.899  32.941 118.906  32.939 118.919
     32.940 118.931
     32.941 118.942  32.942 118.953  32.944 118.965  32.947 118.975
     32.950 118.986
     32.952 118.997  32.953 119.009  32.948 119.018  32.939 119.022
     32.929 119.023
     32.920 119.020  32.910 119.018  32.907 119.029  32.903 119.040
     32.894 119.044
     32.886 119.050  32.879 119.057  32.871 119.064  32.862 119.068
     32.854 119.075
     32.846 119.083  32.840 119.092  32.832 119.099  32.823 119.103
     32.817 119.112
        200       2  32.827  32.473 119.222 118.697
     32.819 119.122  32.822 119.133  32.823 119.144  32.825 119.156
     32.827 119.167
     32.825 119.178  32.821 119.188  32.811 119.186  32.801 119.186
     32.792 119.188
     32.782 119.191  32.773 119.193  32.764 119.197  32.755 119.199
     32.745 119.201
     32.735 119.202  32.726 119.203  32.716 119.205  32.706 119.205
     32.697 119.202
     32.688 119.197  32.680 119.194  32.669 119.193  32.660 119.195
     32.652 119.201
     32.644 119.207  32.635 119.212  32.626 119.216  32.617 119.218
     32.607 119.221
     32.598 119.222  32.589 119.221  32.582 119.213  32.580 119.202
     32.581 119.190
     32.581 119.178  32.576 119.169  32.569 119.162  32.560 119.157
     32.551 119.157
     32.543 119.164  32.534 119.162  32.526 119.157  32.517 119.151
     32.508 119.149
     32.500 119.143  32.493 119.134  32.487 119.125  32.481 119.116
     32.475 119.107
     32.473 119.096  32.475 119.085  32.481 119.075  32.483 119.063
     32.489 119.054
     32.497 119.049  32.506 119.045  32.514 119.039  32.517 119.029
     32.516 119.017
     32.511 119.008  32.505 118.998  32.506 118.987  32.508 118.976
     32.512 118.966
     32.520 118.960  32.528 118.954  32.536 118.946  32.541 118.937
     32.545 118.926
     32.546 118.915  32.550 118.905  32.555 118.896  32.564 118.896
     32.574 118.898
     32.584 118.897  32.588 118.887  32.583 118.878  32.574 118.874
     32.568 118.865
     32.569 118.854  32.571 118.843  32.572 118.831  32.573 118.820
     32.583 118.818
     32.593 118.818  32.595 118.807  32.589 118.799  32.583 118.791
     32.580 118.780
     32.588 118.773  32.597 118.768  32.596 118.757  32.591 118.748
     32.587 118.737
     32.591 118.727  32.598 118.718  32.603 118.709  32.607 118.699
     32.597 118.697
        200       2  32.593  32.118 118.709 118.377
     32.589 118.692  32.583 118.683  32.582 118.671  32.583 118.660
     32.583 118.648
     32.584 118.636  32.586 118.625  32.589 118.614  32.592 118.603
     32.593 118.592
     32.593 118.581  32.587 118.571  32.580 118.564  32.571 118.559
     32.561 118.560
     32.554 118.568  32.548 118.577  32.541 118.585  32.533 118.592
     32.526 118.600
     32.518 118.607  32.509 118.605  32.503 118.595  32.495 118.588
     32.487 118.585
     32.474 118.591  32.472 118.601  32.469 118.612  32.469 118.624
     32.469 118.636
     32.470 118.647  32.469 118.659  32.467 118.670  32.461 118.678
     32.452 118.682
     32.442 118.682  32.431 118.681  32.422 118.683  32.412 118.685
     32.402 118.685
     32.392 118.683  32.383 118.684  32.381 118.695  32.373 118.701
     32.363 118.698
     32.358 118.688  32.348 118.689  32.340 118.694  32.331 118.698
     32.323 118.704
     32.314 118.709  32.305 118.709  32.297 118.702  32.292 118.692
     32.288 118.682
     32.286 118.671  32.279 118.663  32.269 118.664  32.260 118.665
     32.249 118.674
     32.241 118.668  32.237 118.658  32.232 118.649  32.226 118.640
     32.220 118.630
     32.214 118.625  32.211 118.615  32.207 118.605  32.202 118.595
     32.197 118.585
     32.193 118.576  32.186 118.567  32.179 118.559  32.175 118.549
     32.176 118.537
     32.177 118.526  32.181 118.516  32.187 118.507  32.191 118.497
     32.192 118.485
     32.191 118.474  32.193 118.463  32.195 118.452  32.198 118.441
     32.202 118.431
     32.207 118.422  32.212 118.412  32.214 118.404  32.207 118.396
     32.200 118.388
     32.192 118.382  32.186 118.377  32.179 118.385  32.170 118.388
     32.162 118.382
     32.153 118.379  32.143 118.379  32.134 118.383  32.126 118.389
     32.118 118.396
        200       2  32.110  31.639 118.778 118.345
     32.110 118.403  32.102 118.408  32.093 118.406  32.086 118.399
     32.080 118.390
     32.074 118.382  32.065 118.378  32.056 118.380  32.049 118.380
     32.041 118.386
     32.032 118.390  32.022 118.391  32.012 118.390  32.003 118.385
     31.996 118.377
     31.990 118.369  31.983 118.361  31.975 118.355  31.966 118.350
     31.958 118.345
     31.948 118.345  31.939 118.348  31.930 118.353  31.923 118.361
     31.918 118.370
     31.912 118.380  31.904 118.385  31.895 118.388  31.888 118.397
     31.888 118.408
     31.891 118.418  31.891 118.429  31.887 118.439  31.879 118.447
     31.871 118.454
     31.864 118.461  31.856 118.467  31.849 118.476  31.845 118.487
     31.836 118.491
     31.827 118.493  31.818 118.489  31.810 118.482  31.804 118.473
     31.795 118.468
     31.782 118.469  31.779 118.480  31.779 118.492  31.778 118.504
     31.776 118.515
     31.767 118.514  31.759 118.509  31.749 118.506  31.741 118.510
     31.737 118.520
     31.737 118.532  31.739 118.543  31.744 118.553  31.749 118.562
     31.755 118.570
     31.762 118.579  31.766 118.589  31.766 118.601  31.762 118.611
     31.754 118.617
     31.750 118.627  31.748 118.638  31.746 118.649  31.746 118.661
     31.738 118.666
     31.728 118.663  31.719 118.661  31.711 118.668  31.703 118.661
     31.699 118.651
     31.697 118.640  31.691 118.632  31.682 118.630  31.672 118.629
     31.662 118.632
     31.653 118.635  31.647 118.636  31.644 118.650  31.644 118.662
     31.645 118.674
     31.644 118.685  31.641 118.695  31.639 118.707  31.643 118.717
     31.650 118.724
     31.659 118.727  31.669 118.729  31.677 118.734  31.683 118.742
     31.686 118.753
     31.687 118.764  31.681 118.773  31.673 118.778  31.664 118.773
     31.657 118.766
        200       2  31.648  31.234 119.145 118.692
     31.648 118.768  31.639 118.771  31.630 118.776  31.624 118.785
     31.624 118.796
     31.629 118.805  31.634 118.815  31.627 118.823  31.625 118.834
     31.624 118.846
     31.618 118.854  31.609 118.857  31.599 118.855  31.590 118.853
     31.580 118.854
     31.571 118.858  31.561 118.859  31.552 118.856  31.542 118.855
     31.533 118.859
     31.525 118.867  31.517 118.872  31.507 118.873  31.497 118.872
     31.487 118.871
     31.478 118.869  31.469 118.869  31.460 118.866  31.450 118.865
     31.441 118.863
     31.431 118.861  31.423 118.856  31.415 118.850  31.407 118.844
     31.399 118.837
     31.393 118.829  31.387 118.819  31.381 118.810  31.375 118.801
     31.370 118.792
     31.366 118.782  31.363 118.772  31.360 118.761  31.359 118.750
     31.356 118.739
     31.351 118.730  31.342 118.725  31.334 118.719  31.325 118.715
     31.318 118.708
     31.310 118.701  31.300 118.692  31.294 118.699  31.290 118.709
     31.285 118.719
     31.282 118.729  31.280 118.740  31.281 118.752  31.280 118.763
     31.271 118.767
     31.261 118.768  31.252 118.770  31.243 118.773  31.235 118.780
     31.235 118.791
     31.237 118.802  31.239 118.813  31.241 118.824  31.242 118.835
     31.242 118.847
     31.241 118.858  31.239 118.869  31.237 118.880  31.236 118.891
     31.240 118.900
     31.243 118.911  31.246 118.921  31.245 118.933  31.242 118.944
     31.240 118.955
     31.238 118.966  31.237 118.977  31.240 118.988  31.245 118.997
     31.246 119.009
     31.241 119.019  31.234 119.026  31.237 119.037  31.238 119.048
     31.235 119.059
     31.235 119.070  31.240 119.079  31.245 119.088  31.250 119.098
     31.256 119.107
     31.263 119.115  31.270 119.124  31.277 119.131  31.285 119.137
     31.292 119.145
        200       2  31.301  30.820 119.627 119.156
     31.296 119.156  31.298 119.166  31.300 119.177  31.292 119.183
     31.283 119.187
     31.273 119.187  31.269 119.196  31.268 119.207  31.263 119.217
     31.258 119.226
     31.256 119.236  31.263 119.243  31.270 119.249  31.267 119.260
     31.258 119.263
     31.259 119.273  31.262 119.283  31.260 119.294  31.259 119.305
     31.263 119.315
     31.268 119.324  31.277 119.329  31.286 119.331  31.295 119.336
     31.301 119.344
     31.295 119.351  31.285 119.352  31.277 119.359  31.271 119.368
     31.261 119.370
     31.253 119.366  31.244 119.361  31.236 119.356  31.227 119.351
     31.217 119.350
     31.208 119.349  31.200 119.356  31.202 119.373  31.205 119.384
     31.198 119.388
     31.190 119.381  31.183 119.385  31.182 119.397  31.183 119.408
     31.183 119.420
     31.181 119.431  31.175 119.440  31.170 119.449  31.167 119.460
     31.163 119.470
     31.160 119.480  31.159 119.491  31.160 119.503  31.160 119.515
     31.161 119.526
     31.155 119.537  31.150 119.546  31.144 119.555  31.139 119.564
     31.135 119.575
     31.135 119.586  31.136 119.597  31.139 119.607  31.143 119.617
     31.135 119.622
     31.125 119.624  31.116 119.625  31.106 119.627  31.098 119.622
     31.088 119.622
     31.077 119.622  31.068 119.620  31.058 119.619  31.048 119.619
     31.038 119.620
     31.028 119.623  31.019 119.621  31.010 119.617  31.003 119.609
     30.995 119.601
     30.989 119.592  30.982 119.584  30.973 119.580  30.964 119.577
     30.954 119.575
     30.945 119.572  30.937 119.567  30.929 119.561  30.921 119.555
     30.912 119.551
     30.902 119.551  30.892 119.552  30.882 119.553  30.873 119.557
     30.863 119.565
     30.853 119.564  30.844 119.561  30.836 119.555  30.828 119.548
     30.820 119.542
        200       2  30.813  30.342 119.535 119.233
     30.813 119.535  30.805 119.530  30.794 119.530  30.785 119.528
     30.781 119.518
     30.785 119.509  30.789 119.499  30.793 119.489  30.790 119.478
     30.781 119.474
     30.772 119.472  30.762 119.472  30.751 119.472  30.742 119.474
     30.732 119.475
     30.722 119.474  30.713 119.472  30.704 119.466  30.697 119.458
     30.690 119.450
     30.684 119.442  30.678 119.433  30.669 119.427  30.661 119.423
     30.654 119.415
     30.650 119.406  30.658 119.401  30.667 119.397  30.676 119.395
     30.685 119.391
     30.693 119.386  30.694 119.375  30.693 119.364  30.688 119.355
     30.682 119.346
     30.676 119.338  30.669 119.329  30.661 119.322  30.653 119.318
     30.645 119.313
     30.637 119.305  30.632 119.296  30.628 119.286  30.626 119.275
     30.623 119.265
     30.622 119.253  30.620 119.243  30.608 119.239  30.600 119.244
     30.591 119.249
     30.582 119.250  30.573 119.246  30.566 119.240  30.558 119.233
     30.548 119.233
     30.541 119.241  30.539 119.251  30.536 119.262  30.529 119.269
     30.524 119.278
     30.525 119.289  30.530 119.299  30.535 119.308  30.533 119.319
     30.524 119.324
     30.515 119.327  30.505 119.327  30.495 119.326  30.486 119.328
     30.479 119.336
     30.470 119.340  30.461 119.338  30.453 119.343  30.447 119.352
     30.438 119.355
     30.429 119.349  30.420 119.350  30.414 119.359  30.409 119.369
     30.403 119.377
     30.396 119.385  30.393 119.395  30.388 119.398  30.380 119.391
     30.375 119.382
     30.371 119.372  30.367 119.362  30.363 119.352  30.361 119.342
     30.367 119.333
     30.375 119.325  30.376 119.314  30.375 119.303  30.372 119.293
     30.366 119.284
     30.361 119.274  30.357 119.265  30.354 119.254  30.350 119.244
     30.342 119.238
        200       2  30.359  29.896 119.233 118.821
     30.334 119.233  30.325 119.229  30.315 119.228  30.306 119.225
     30.300 119.217
     30.299 119.205  30.301 119.194  30.305 119.185  30.309 119.174
     30.312 119.164
     30.314 119.154  30.316 119.143  30.317 119.132  30.318 119.122
     30.322 119.112
     30.325 119.102  30.329 119.091  30.331 119.081  30.332 119.070
     30.326 119.062
     30.318 119.056  30.311 119.048  30.315 119.038  30.321 119.029
     30.328 119.022
     30.335 119.014  30.340 119.005  30.345 118.995  30.350 118.986
     30.356 118.977
     30.359 118.967  30.359 118.956  30.359 118.944  30.357 118.934
     30.353 118.924
     30.347 118.915  30.341 118.906  30.334 118.897  30.329 118.888
     30.321 118.877
     30.311 118.876  30.301 118.875  30.291 118.875  30.282 118.878
     30.274 118.884
     30.265 118.888  30.255 118.889  30.245 118.889  30.237 118.894
     30.232 118.903
     30.229 118.914  30.223 118.921  30.215 118.916  30.210 118.907
     30.204 118.898
     30.197 118.890  30.191 118.881  30.187 118.871  30.183 118.861
     30.177 118.852
     30.170 118.845  30.161 118.846  30.158 118.857  30.160 118.868
     30.163 118.879
     30.163 118.890  30.160 118.901  30.151 118.904  30.143 118.899
     30.137 118.890
     30.133 118.880  30.128 118.871  30.118 118.861  30.109 118.865
     30.100 118.867
     30.091 118.871  30.083 118.877  30.075 118.883  30.066 118.887
     30.057 118.889
     30.047 118.891  30.037 118.892  30.025 118.897  30.016 118.892
     30.006 118.886
     29.997 118.882  29.988 118.884  29.978 118.883  29.969 118.882
     29.958 118.882
     29.950 118.877  29.943 118.870  29.941 118.859  29.938 118.849
     29.939 118.838
     29.935 118.828  29.926 118.825  29.916 118.825  29.906 118.824
     29.896 118.821
        200       2  29.888  29.393 118.816 118.136
     29.888 118.816  29.881 118.809  29.873 118.802  29.866 118.794
     29.858 118.787
     29.851 118.779  29.845 118.771  29.852 118.762  29.853 118.752
     29.845 118.747
     29.836 118.743  29.826 118.745  29.819 118.738  29.815 118.728
     29.806 118.724
     29.796 118.725  29.787 118.728  29.778 118.731  29.768 118.731
     29.758 118.729
     29.749 118.726  29.740 118.723  29.731 118.719  29.723 118.713
     29.716 118.706
     29.710 118.698  29.703 118.689  29.696 118.682  29.688 118.675
     29.680 118.669
     29.673 118.661  29.666 118.655  29.658 118.647  29.651 118.640
     29.645 118.631
     29.645 118.622  29.649 118.612  29.650 118.601  29.648 118.590
     29.641 118.582
     29.636 118.572  29.628 118.568  29.619 118.565  29.614 118.556
     29.615 118.545
     29.608 118.540  29.599 118.536  29.592 118.530  29.586 118.521
     29.581 118.512
     29.573 118.503  29.566 118.497  29.557 118.491  29.548 118.490
     29.538 118.490
     29.528 118.489  29.518 118.487  29.511 118.482  29.507 118.472
     29.505 118.461
     29.504 118.451  29.508 118.441  29.511 118.436  29.504 118.427
     29.502 118.417
     29.500 118.406  29.492 118.401  29.482 118.402  29.475 118.395
     29.475 118.384
     29.480 118.375  29.481 118.363  29.478 118.353  29.474 118.339
     29.480 118.330
     29.487 118.322  29.488 118.311  29.482 118.305  29.473 118.303
     29.462 118.303
     29.454 118.307  29.445 118.310  29.435 118.308  29.428 118.302
     29.427 118.292
     29.427 118.280  29.427 118.269  29.427 118.258  29.426 118.246
     29.425 118.235
     29.421 118.226  29.416 118.216  29.408 118.209  29.401 118.201
     29.393 118.187
     29.401 118.176  29.405 118.166  29.409 118.156  29.413 118.146
     29.416 118.136
        200       2  29.786  29.421 118.133 117.417
     29.421 118.127  29.427 118.118  29.435 118.122  29.442 118.128
     29.452 118.126
     29.460 118.122  29.470 118.119  29.477 118.124  29.484 118.133
     29.493 118.133
     29.502 118.129  29.510 118.124  29.518 118.117  29.521 118.107
     29.525 118.097
     29.528 118.086  29.530 118.076  29.535 118.067  29.541 118.058
     29.548 118.049
     29.555 118.041  29.562 118.033  29.568 118.024  29.570 118.013
     29.569 118.002
     29.564 117.994  29.557 117.987  29.560 117.977  29.562 117.966
     29.562 117.955
     29.559 117.945  29.554 117.936  29.548 117.927  29.545 117.917
     29.545 117.905
     29.546 117.894  29.545 117.883  29.544 117.872  29.549 117.863
     29.557 117.854
     29.560 117.844  29.561 117.833  29.564 117.823  29.571 117.814
     29.573 117.804
     29.572 117.793  29.568 117.783  29.563 117.774  29.558 117.765
     29.555 117.754
     29.555 117.743  29.553 117.732  29.552 117.721  29.551 117.710
     29.549 117.699
     29.560 117.690  29.570 117.688  29.579 117.685  29.587 117.679
     29.594 117.672
     29.601 117.663  29.607 117.658  29.615 117.652  29.609 117.643
     29.604 117.634
     29.598 117.626  29.593 117.616  29.588 117.607  29.584 117.597
     29.583 117.586
     29.586 117.577  29.589 117.566  29.593 117.556  29.596 117.546
     29.600 117.537
     29.608 117.531  29.617 117.527  29.624 117.519  29.634 117.521
     29.642 117.526
     29.651 117.530  29.660 117.525  29.663 117.514  29.664 117.503
     29.664 117.493
     29.670 117.484  29.678 117.477  29.685 117.469  29.691 117.461
     29.698 117.453
     29.708 117.453  29.718 117.453  29.728 117.454  29.738 117.452
     29.747 117.448
     29.755 117.443  29.764 117.440  29.772 117.433  29.779 117.426
     29.786 117.417
        200       2  29.934  29.605 117.426 116.873
     29.794 117.410  29.808 117.415  29.817 117.416  29.828 117.416
     29.837 117.420
     29.845 117.426  29.852 117.421  29.851 117.409  29.851 117.398
     29.848 117.387
     29.841 117.381  29.831 117.380  29.822 117.378  29.817 117.368
     29.818 117.358
     29.825 117.350  29.834 117.346  29.842 117.342  29.849 117.334
     29.852 117.324
     29.848 117.314  29.839 117.310  29.836 117.299  29.832 117.291
     29.824 117.286
     29.826 117.275  29.827 117.264  29.832 117.255  29.842 117.253
     29.851 117.257
     29.860 117.260  29.869 117.260  29.878 117.256  29.886 117.252
     29.896 117.249
     29.904 117.245  29.913 117.240  29.919 117.232  29.923 117.222
     29.930 117.214
     29.934 117.207  29.927 117.198  29.922 117.189  29.919 117.179
     29.916 117.169
     29.912 117.159  29.910 117.148  29.906 117.138  29.899 117.132
     29.890 117.128
     29.880 117.128  29.870 117.127  29.861 117.124  29.856 117.114
     29.851 117.105
     29.845 117.097  29.838 117.088  29.831 117.080  29.823 117.075
     29.815 117.080
     29.809 117.088  29.805 117.098  29.801 117.108  29.798 117.118
     29.793 117.128
     29.784 117.131  29.774 117.130  29.766 117.123  29.760 117.114
     29.751 117.110
     29.742 117.109  29.732 117.112  29.724 117.116  29.715 117.113
     29.708 117.105
     29.700 117.098  29.694 117.089  29.694 117.078  29.692 117.067
     29.688 117.057
     29.684 117.047  29.681 117.038  29.678 117.027  29.682 117.018
     29.683 117.007
     29.680 116.997  29.675 116.988  29.667 116.982  29.659 116.976
     29.653 116.967
     29.649 116.957  29.646 116.947  29.644 116.936  29.641 116.926
     29.636 116.917
     29.631 116.907  29.625 116.899  29.619 116.890  29.612 116.881
     29.605 116.873
        200       2  30.057  29.559 116.907 116.658
     29.599 116.864  29.592 116.856  29.584 116.848  29.576 116.844
     29.568 116.838
     29.563 116.829  29.559 116.818  29.561 116.807  29.565 116.798
     29.570 116.788
     29.578 116.783  29.587 116.780  29.593 116.772  29.594 116.760
     29.590 116.751
     29.582 116.746  29.574 116.740  29.566 116.734  29.568 116.724
     29.573 116.715
     29.582 116.710  29.591 116.708  29.598 116.700  29.602 116.694
     29.609 116.694
     29.618 116.691  29.626 116.685  29.632 116.676  29.636 116.667
     29.644 116.662
     29.654 116.662  29.662 116.666  29.669 116.675  29.675 116.683
     29.682 116.690
     29.689 116.687  29.688 116.675  29.689 116.665  29.696 116.658
     29.705 116.662
     29.712 116.669  29.719 116.678  29.726 116.687  29.733 116.694
     29.740 116.702
     29.747 116.709  29.754 116.718  29.761 116.725  29.769 116.732
     29.777 116.737
     29.786 116.742  29.794 116.749  29.800 116.757  29.803 116.768
     29.802 116.778
     29.800 116.789  29.802 116.800  29.807 116.810  29.814 116.817
     29.823 116.822
     29.831 116.827  29.839 116.833  29.847 116.840  29.853 116.849
     29.857 116.858
     29.864 116.867  29.871 116.874  29.879 116.880  29.888 116.885
     29.897 116.887
     29.906 116.891  29.915 116.894  29.923 116.899  29.932 116.904
     29.941 116.907
     29.951 116.906  29.960 116.904  29.968 116.898  29.976 116.891
     29.981 116.882
     29.981 116.870  29.978 116.860  29.974 116.850  29.973 116.840
     29.982 116.836
     29.994 116.829  30.001 116.822  30.011 116.821  30.015 116.812
     30.015 116.801
     30.016 116.790  30.023 116.781  30.031 116.776  30.040 116.774
     30.046 116.765
     30.051 116.755  30.057 116.747  30.053 116.737  30.048 116.727
     30.049 116.717
        200       2  30.073  29.778 116.706 116.079
     30.052 116.706  30.055 116.696  30.059 116.687  30.064 116.677
     30.068 116.667
     30.071 116.657  30.073 116.646  30.072 116.635  30.069 116.625
     30.066 116.614
     30.059 116.606  30.051 116.599  30.047 116.593  30.039 116.588
     30.030 116.583
     30.022 116.578  30.013 116.574  30.004 116.571  29.995 116.572
     29.985 116.573
     29.974 116.574  29.964 116.575  29.955 116.573  29.946 116.568
     29.938 116.564
     29.929 116.561  29.920 116.558  29.911 116.553  29.905 116.545
     29.902 116.534
     29.902 116.523  29.904 116.512  29.905 116.501  29.904 116.490
     29.902 116.480
     29.898 116.469  29.895 116.459  29.891 116.450  29.885 116.441
     29.879 116.432
     29.873 116.423  29.867 116.414  29.861 116.406  29.854 116.398
     29.849 116.394
     29.850 116.384  29.848 116.373  29.844 116.363  29.839 116.354
     29.832 116.346
     29.827 116.337  29.823 116.327  29.818 116.318  29.813 116.309
     29.808 116.300
     29.802 116.291  29.795 116.283  29.789 116.274  29.782 116.266
     29.778 116.256
     29.781 116.247  29.789 116.242  29.798 116.239  29.806 116.232
     29.814 116.226
     29.821 116.219  29.825 116.209  29.829 116.199  29.828 116.188
     29.827 116.177
     29.825 116.166  29.824 116.155  29.822 116.145  29.825 116.131
     29.838 116.131
     29.848 116.131  29.857 116.129  29.865 116.129  29.875 116.130
     29.885 116.129
     29.894 116.125  29.902 116.121  29.910 116.113  29.918 116.107
     29.926 116.100
     29.934 116.095  29.942 116.091  29.951 116.087  29.960 116.082
     29.969 116.079
     29.979 116.079  29.989 116.081  29.999 116.083  30.009 116.083
     30.017 116.087
     30.026 116.091  30.036 116.093  30.044 116.087  30.052 116.081
     30.062 116.080
        200       2  30.762  30.072 116.084 115.758
     30.072 116.078  30.081 116.079  30.091 116.080  30.101 116.082
     30.110 116.084
     30.120 116.084  30.127 116.076  30.133 116.068  30.141 116.061
     30.150 116.056
     30.159 116.055  30.169 116.056  30.179 116.058  30.188 116.060
     30.198 116.061
     30.208 116.061  30.215 116.054  30.223 116.046  30.228 116.037
     30.233 116.027
     30.240 116.019  30.246 116.010  30.252 116.001  30.258 115.992
     30.266 115.987
     30.276 115.986  30.285 115.983  30.293 115.978  30.301 115.971
     30.307 115.961
     30.308 115.950  30.308 115.939  30.308 115.927  30.306 115.916
     30.304 115.906
     30.314 115.902  30.324 115.900  30.333 115.902  30.343 115.902
     30.351 115.897
     30.359 115.891  30.368 115.887  30.378 115.886  30.387 115.888
     30.396 115.891
     30.405 115.896  30.411 115.905  30.415 115.915  30.420 115.924
     30.425 115.933
     30.429 115.928  30.435 115.919  30.442 115.911  30.448 115.902
     30.455 115.894
     30.464 115.893  30.473 115.896  30.482 115.901  30.491 115.904
     30.500 115.906
     30.509 115.909  30.519 115.908  30.527 115.902  30.534 115.894
     30.542 115.889
     30.550 115.882  30.558 115.877  30.568 115.877  30.577 115.873
     30.586 115.869
     30.593 115.862  30.599 115.852  30.601 115.842  30.601 115.830
     30.598 115.820
     30.597 115.810  30.606 115.804  30.615 115.802  30.625 115.801
     30.634 115.799
     30.643 115.795  30.651 115.788  30.658 115.780  30.664 115.771
     30.668 115.761
     30.674 115.760  30.683 115.758  30.693 115.759  30.703 115.759
     30.713 115.760
     30.723 115.762  30.732 115.764  30.741 115.767  30.750 115.771
     30.758 115.777
     30.762 115.787  30.760 115.799  30.757 115.809  30.759 115.820
     30.760 115.831
        200       2  31.213  30.766 116.063 115.589
     30.766 115.837  30.773 115.844  30.778 115.854  30.783 115.863
     30.792 115.861
     30.800 115.855  30.809 115.850  30.818 115.846  30.826 115.843
     30.835 115.839
     30.845 115.840  30.854 115.845  30.862 115.852  30.868 115.860
     30.874 115.870
     30.879 115.879  30.883 115.889  30.887 115.899  30.890 115.909
     30.892 115.920
     30.896 115.930  30.904 115.938  30.911 115.945  30.919 115.951
     30.926 115.959
     30.931 115.969  30.936 115.978  30.940 115.988  30.942 115.999
     30.943 116.010
     30.944 116.021  30.944 116.033  30.947 116.043  30.954 116.051
     30.961 116.058
     30.970 116.063  30.977 116.055  30.986 116.051  30.996 116.051
     31.006 116.049
     31.009 116.038  31.012 116.028  31.010 116.017  31.009 116.005
     31.013 115.995
     31.019 115.987  31.029 115.985  31.039 115.984  31.044 115.974
     31.043 115.963
     31.041 115.952  31.038 115.942  31.042 115.932  31.051 115.929
     31.062 115.929
     31.071 115.931  31.080 115.926  31.086 115.917  31.090 115.907
     31.093 115.897
     31.099 115.888  31.109 115.886  31.119 115.887  31.131 115.880
     31.139 115.874
     31.144 115.865  31.147 115.855  31.146 115.843  31.140 115.834
     31.134 115.826
     31.129 115.816  31.127 115.805  31.127 115.793  31.125 115.782
     31.121 115.772
     31.121 115.762  31.129 115.757  31.137 115.750  31.143 115.742
     31.151 115.737
     31.159 115.731  31.167 115.725  31.175 115.718  31.183 115.712
     31.191 115.706
     31.197 115.697  31.201 115.687  31.204 115.677  31.208 115.667
     31.212 115.657
     31.213 115.645  31.208 115.637  31.200 115.631  31.191 115.627
     31.183 115.622
     31.175 115.615  31.169 115.607  31.162 115.598  31.154 115.594
     31.145 115.589
        200       2  31.811  31.145 115.620 115.365
     31.145 115.578  31.151 115.569  31.156 115.560  31.161 115.550
     31.169 115.543
     31.177 115.538  31.185 115.533  31.194 115.529  31.204 115.527
     31.214 115.527
     31.224 115.528  31.233 115.528  31.243 115.526  31.251 115.520
     31.257 115.511
     31.263 115.502  31.266 115.492  31.268 115.481  31.267 115.469
     31.278 115.463
     31.287 115.459  31.296 115.457  31.305 115.452  31.312 115.446
     31.321 115.443
     31.330 115.446  31.340 115.447  31.349 115.442  31.354 115.433
     31.353 115.422
     31.348 115.413  31.345 115.402  31.343 115.391  31.342 115.381
     31.348 115.372
     31.356 115.367  31.366 115.365  31.376 115.365  31.385 115.367
     31.395 115.369
     31.405 115.368  31.416 115.369  31.425 115.372  31.434 115.374
     31.444 115.375
     31.454 115.374  31.463 115.370  31.472 115.366  31.482 115.366
     31.491 115.369
     31.499 115.375  31.507 115.383  31.514 115.391  31.520 115.400
     31.528 115.408
     31.535 115.415  31.544 115.418  31.554 115.420  31.563 115.424
     31.571 115.428
     31.580 115.432  31.588 115.439  31.593 115.448  31.599 115.457
     31.605 115.466
     31.612 115.474  31.621 115.478  31.630 115.475  31.640 115.475
     31.649 115.477
     31.659 115.479  31.667 115.485  31.674 115.492  31.680 115.501
     31.686 115.510
     31.694 115.518  31.701 115.527  31.708 115.534  31.715 115.543
     31.717 115.553
     31.717 115.565  31.721 115.576  31.724 115.586  31.729 115.596
     31.734 115.605
     31.740 115.614  31.748 115.620  31.757 115.617  31.765 115.609
     31.772 115.601
     31.777 115.592  31.784 115.583  31.790 115.574  31.797 115.566
     31.811 115.560
     31.808 115.569  31.809 115.581  31.808 115.592  31.807 115.604
     31.805 115.615
        200       2  32.443  31.772 115.928 115.625
     31.802 115.625  31.800 115.636  31.796 115.646  31.792 115.656
     31.786 115.665
     31.782 115.675  31.777 115.685  31.774 115.696  31.773 115.707
     31.774 115.719
     31.780 115.728  31.788 115.734  31.793 115.743  31.792 115.754
     31.789 115.765
     31.785 115.775  31.782 115.785  31.779 115.796  31.778 115.807
     31.776 115.818
     31.773 115.829  31.772 115.840  31.773 115.851  31.774 115.863
     31.777 115.873
     31.781 115.883  31.787 115.893  31.795 115.898  31.802 115.906
     31.807 115.915
     31.818 115.922  31.827 115.917  31.836 115.914  31.846 115.913
     31.854 115.906
     31.859 115.896  31.867 115.891  31.877 115.892  31.885 115.897
     31.893 115.904
     31.901 115.910  31.909 115.914  31.919 115.915  31.929 115.914
     31.938 115.910
     31.947 115.908  31.957 115.905  31.966 115.908  31.974 115.914
     31.982 115.919
     31.991 115.923  31.999 115.925  32.009 115.921  32.019 115.919
     32.029 115.919
     32.039 115.920  32.048 115.922  32.057 115.926  32.066 115.927
     32.075 115.924
     32.084 115.918  32.093 115.916  32.103 115.917  32.113 115.919
     32.122 115.921
     32.131 115.923  32.141 115.925  32.150 115.927  32.160 115.928
     32.170 115.926
     32.179 115.921  32.187 115.916  32.195 115.910  32.204 115.906
     32.213 115.903
     32.223 115.903  32.233 115.904  32.243 115.905  32.253 115.905
     32.264 115.905
     32.273 115.904  32.282 115.899  32.287 115.889  32.295 115.885
     32.305 115.888
     32.312 115.893  32.321 115.898  32.331 115.899  32.341 115.897
     32.351 115.896
     32.360 115.894  32.370 115.893  32.380 115.891  32.389 115.889
     32.399 115.887
     32.408 115.884  32.417 115.880  32.425 115.876  32.434 115.873
     32.443 115.869
        200       2  32.725  32.411 115.869 115.169
     32.451 115.863  32.461 115.860  32.470 115.863  32.478 115.869
     32.487 115.865
     32.492 115.855  32.497 115.845  32.499 115.835  32.506 115.826
     32.503 115.815
     32.496 115.806  32.489 115.798  32.483 115.790  32.479 115.779
     32.486 115.772
     32.494 115.767  32.495 115.757  32.489 115.747  32.483 115.740
     32.478 115.730
     32.476 115.719  32.477 115.708  32.483 115.699  32.487 115.689
     32.483 115.679
     32.473 115.680  32.464 115.680  32.459 115.670  32.452 115.662
     32.443 115.660
     32.434 115.656  32.426 115.650  32.418 115.642  32.414 115.632
     32.416 115.621
     32.417 115.610  32.422 115.600  32.422 115.592  32.422 115.582
     32.416 115.573
     32.412 115.562  32.411 115.551  32.419 115.546  32.427 115.539
     32.435 115.532
     32.440 115.522  32.444 115.513  32.452 115.505  32.461 115.501
     32.470 115.499
     32.480 115.499  32.490 115.497  32.499 115.494  32.509 115.492
     32.517 115.486
     32.525 115.479  32.532 115.472  32.538 115.462  32.545 115.453
     32.550 115.444
     32.556 115.435  32.564 115.429  32.572 115.423  32.577 115.413
     32.582 115.403
     32.584 115.392  32.584 115.380  32.582 115.369  32.578 115.359
     32.574 115.349
     32.572 115.338  32.569 115.327  32.568 115.315  32.570 115.304
     32.575 115.294
     32.582 115.286  32.590 115.271  32.587 115.260  32.585 115.249
     32.586 115.238
     32.593 115.230  32.596 115.220  32.594 115.209  32.593 115.197
     32.595 115.186
     32.599 115.175  32.603 115.169  32.610 115.170  32.619 115.174
     32.628 115.176
     32.638 115.178  32.648 115.179  32.657 115.178  32.666 115.173
     32.676 115.173
     32.686 115.173  32.695 115.174  32.705 115.175  32.716 115.175
     32.725 115.177
        200       2  33.158  32.734 115.214 114.876
     32.734 115.180  32.744 115.181  32.753 115.180  32.762 115.177
     32.772 115.175
     32.782 115.176  32.786 115.184  32.787 115.195  32.796 115.196
     32.803 115.187
     32.812 115.183  32.822 115.185  32.830 115.189  32.840 115.190
     32.849 115.184
     32.856 115.176  32.857 115.165  32.859 115.153  32.864 115.144
     32.874 115.142
     32.882 115.135  32.889 115.126  32.893 115.116  32.892 115.105
     32.890 115.094
     32.889 115.083  32.895 115.074  32.902 115.065  32.906 115.055
     32.907 115.043
     32.907 115.031  32.909 115.020  32.919 115.020  32.928 115.017
     32.934 115.008
     32.937 114.997  32.939 114.986  32.941 114.975  32.941 114.963
     32.940 114.951
     32.939 114.939  32.939 114.928  32.942 114.917  32.950 114.911
     32.959 114.906
     32.967 114.901  32.972 114.891  32.975 114.881  32.983 114.876
     32.993 114.879
     33.001 114.885  33.007 114.894  33.014 114.903  33.022 114.910
     33.030 114.915
     33.040 114.912  33.049 114.910  33.058 114.905  33.066 114.900
     33.075 114.895
     33.083 114.891  33.095 114.893  33.103 114.896  33.112 114.901
     33.121 114.905
     33.129 114.910  33.138 114.914  33.146 114.919  33.155 114.926
     33.158 114.937
     33.155 114.948  33.151 114.958  33.148 114.969  33.134 114.968
     33.125 114.966
     33.116 114.970  33.109 114.978  33.103 114.988  33.098 114.997
     33.091 115.006
     33.087 115.016  33.087 115.029  33.090 115.040  33.086 115.050
     33.082 115.060
     33.082 115.072  33.082 115.084  33.082 115.095  33.081 115.107
     33.080 115.119
     33.082 115.129  33.086 115.139  33.092 115.148  33.100 115.155
     33.108 115.162
     33.114 115.172  33.117 115.182  33.119 115.193  33.123 115.204
     33.127 115.214
        200       2  33.777  33.128 115.625 115.224
     33.131 115.224  33.131 115.236  33.132 115.250  33.129 115.260
     33.128 115.272
     33.136 115.278  33.144 115.283  33.154 115.284  33.164 115.285
     33.174 115.285
     33.184 115.287  33.193 115.290  33.201 115.295  33.209 115.303
     33.214 115.312
     33.222 115.320  33.230 115.324  33.240 115.327  33.249 115.329
     33.259 115.329
     33.269 115.330  33.279 115.331  33.288 115.334  33.297 115.337
     33.306 115.341
     33.314 115.347  33.323 115.352  33.332 115.351  33.341 115.346
     33.349 115.341
     33.359 115.338  33.367 115.334  33.376 115.330  33.385 115.327
     33.394 115.323
     33.403 115.318  33.411 115.313  33.419 115.307  33.428 115.303
     33.438 115.307
     33.446 115.314  33.451 115.324  33.454 115.334  33.463 115.341
     33.472 115.340
     33.482 115.339  33.491 115.336  33.501 115.339  33.509 115.345
     33.516 115.352
     33.515 115.363  33.510 115.373  33.503 115.381  33.502 115.392
     33.512 115.396
     33.521 115.399  33.530 115.403  33.539 115.408  33.547 115.414
     33.553 115.423
     33.556 115.434  33.557 115.446  33.557 115.458  33.555 115.469
     33.553 115.481
     33.553 115.493  33.555 115.504  33.560 115.514  33.565 115.524
     33.569 115.534
     33.572 115.545  33.575 115.555  33.577 115.567  33.577 115.579
     33.578 115.591
     33.577 115.603  33.577 115.614  33.580 115.625  33.594 115.621
     33.603 115.618
     33.612 115.615  33.622 115.614  33.631 115.610  33.639 115.605
     33.647 115.599
     33.656 115.595  33.666 115.597  33.675 115.599  33.685 115.601
     33.695 115.603
     33.704 115.602  33.713 115.596  33.720 115.588  33.725 115.588
     33.734 115.583
     33.741 115.576  33.749 115.570  33.758 115.565  33.767 115.562
     33.777 115.561
        200       2  34.078  33.771 116.028 115.548
     33.781 115.571  33.781 115.583  33.775 115.593  33.771 115.603
     33.778 115.610
     33.787 115.616  33.795 115.620  33.805 115.623  33.814 115.625
     33.823 115.628
     33.833 115.628  33.843 115.630  33.853 115.630  33.863 115.629
     33.872 115.626
     33.874 115.615  33.873 115.603  33.872 115.591  33.872 115.579
     33.872 115.567
     33.873 115.555  33.880 115.548  33.890 115.550  33.899 115.553
     33.908 115.555
     33.917 115.559  33.926 115.563  33.935 115.567  33.944 115.571
     33.953 115.575
     33.961 115.579  33.971 115.582  33.980 115.584  33.989 115.586
     33.999 115.588
     34.001 115.599  34.008 115.594  34.019 115.594  34.029 115.595
     34.036 115.602
     34.040 115.612  34.040 115.624  34.045 115.634  34.053 115.640
     34.062 115.644
     34.062 115.656  34.066 115.667  34.071 115.677  34.074 115.687
     34.076 115.699
     34.078 115.710  34.076 115.722  34.073 115.732  34.069 115.743
     34.064 115.752
     34.065 115.764  34.073 115.769  34.072 115.785  34.074 115.797
     34.070 115.806
     34.060 115.808  34.051 115.810  34.046 115.819  34.041 115.829
     34.034 115.837
     34.026 115.842  34.016 115.844  34.008 115.849  34.014 115.859
     34.021 115.866
     34.030 115.872  34.032 115.882  34.024 115.889  34.017 115.896
     34.017 115.907
     34.018 115.919  34.013 115.929  34.010 115.940  34.005 115.949
     34.004 115.961
     34.001 115.972  33.999 115.983  33.997 115.995  33.996 116.004
     33.985 116.002
     33.975 116.001  33.965 116.000  33.955 115.999  33.946 115.994
     33.937 115.990
     33.929 115.984  33.921 115.978  33.911 115.978  33.902 115.983
     33.897 115.993
     33.890 115.998  33.881 116.001  33.876 116.010  33.868 116.018
     33.862 116.028
        200       2  34.103  33.703 116.642 116.038
     33.858 116.038  33.850 116.044  33.840 116.046  33.830 116.048
     33.821 116.051
     33.812 116.055  33.803 116.058  33.794 116.060  33.784 116.062
     33.780 116.071
     33.774 116.080  33.766 116.087  33.757 116.092  33.750 116.100
     33.746 116.110
     33.741 116.120  33.735 116.129  33.728 116.137  33.719 116.141
     33.710 116.145
     33.703 116.154  33.707 116.164  33.711 116.174  33.715 116.185
     33.720 116.194
     33.726 116.204  33.729 116.215  33.729 116.226  33.729 116.238
     33.731 116.250
     33.734 116.260  33.741 116.268  33.746 116.278  33.749 116.288
     33.752 116.299
     33.760 116.305  33.767 116.311  33.772 116.321  33.776 116.331
     33.780 116.341
     33.783 116.352  33.785 116.363  33.788 116.374  33.791 116.385
     33.797 116.394
     33.795 116.405  33.799 116.413  33.803 116.423  33.809 116.430
     33.818 116.426
     33.828 116.426  33.837 116.430  33.846 116.435  33.852 116.444
     33.855 116.455
     33.858 116.466  33.861 116.476  33.866 116.486  33.870 116.496
     33.871 116.508
     33.878 116.516  33.883 116.526  33.884 116.538  33.886 116.549
     33.895 116.553
     33.903 116.559  33.908 116.569  33.909 116.580  33.906 116.591
     33.898 116.599
     33.891 116.607  33.887 116.618  33.888 116.629  33.893 116.638
     33.902 116.641
     33.912 116.640  33.922 116.639  33.931 116.636  33.941 116.634
     33.950 116.639
     33.961 116.642  33.971 116.641  33.979 116.634  33.984 116.625
     33.991 116.616
     34.001 116.605  34.010 116.601  34.018 116.597  34.027 116.593
     34.035 116.589
     34.044 116.584  34.050 116.584  34.058 116.581  34.066 116.574
     34.071 116.564
     34.080 116.559  34.090 116.558  34.098 116.553  34.101 116.542
     34.103 116.530
        200       2  34.622  34.110 116.576 116.156
     34.110 116.522  34.120 116.522  34.129 116.525  34.138 116.529
     34.147 116.534
     34.154 116.542  34.161 116.551  34.169 116.559  34.177 116.565
     34.187 116.567
     34.196 116.564  34.203 116.556  34.208 116.546  34.215 116.537
     34.224 116.534
     34.233 116.536  34.241 116.543  34.249 116.551  34.254 116.561
     34.260 116.570
     34.273 116.576  34.281 116.568  34.286 116.559  34.289 116.548
     34.292 116.537
     34.294 116.526  34.295 116.514  34.296 116.502  34.296 116.489
     34.293 116.478
     34.288 116.469  34.281 116.460  34.270 116.451  34.267 116.440
     34.267 116.428
     34.268 116.417  34.269 116.405  34.269 116.393  34.268 116.381
     34.266 116.370
     34.268 116.358  34.274 116.350  34.283 116.356  34.288 116.365
     34.297 116.369
     34.310 116.363  34.315 116.353  34.321 116.344  34.324 116.333
     34.328 116.322
     34.331 116.312  34.335 116.302  34.340 116.291  34.346 116.282
     34.354 116.276
     34.363 116.272  34.371 116.266  34.375 116.255  34.375 116.243
     34.375 116.231
     34.375 116.220  34.378 116.209  34.387 116.204  34.396 116.200
     34.405 116.197
     34.413 116.190  34.420 116.182  34.422 116.172  34.433 116.166
     34.441 116.161
     34.451 116.158  34.460 116.156  34.470 116.157  34.479 116.160
     34.488 116.165
     34.496 116.171  34.503 116.180  34.508 116.189  34.518 116.188
     34.528 116.186
     34.537 116.187  34.548 116.187  34.557 116.189  34.567 116.191
     34.569 116.198
     34.562 116.207  34.555 116.216  34.552 116.227  34.558 116.236
     34.566 116.243
     34.572 116.252  34.579 116.260  34.588 116.265  34.596 116.272
     34.603 116.280
     34.610 116.289  34.615 116.299  34.617 116.310  34.619 116.322
     34.622 116.333
          8       2  34.644  34.625 116.369 116.343
     34.625 116.343  34.629 116.354  34.636 116.363  34.644 116.369
        200       2  40.231  39.762 117.393 116.750
     40.231 117.379  40.220 117.379  40.211 117.383  40.204 117.393
     40.195 117.393
     40.191 117.381  40.188 117.370  40.181 117.362  40.172 117.357
     40.164 117.350
     40.159 117.339  40.160 117.327  40.163 117.315  40.166 117.304
     40.166 117.290
     40.162 117.279  40.154 117.271  40.144 117.271  40.135 117.277
     40.127 117.282
     40.117 117.279  40.114 117.267  40.117 117.255  40.116 117.242
     40.113 117.231
     40.106 117.221  40.099 117.212  40.087 117.202  40.086 117.189
     40.082 117.179
     40.076 117.168  40.072 117.157  40.068 117.146  40.066 117.134
     40.066 117.121
     40.069 117.109  40.072 117.098  40.073 117.085  40.069 117.074
     40.063 117.064
     40.057 117.054  40.052 117.043  40.046 117.033  40.042 117.022
     40.039 117.010
     40.038 116.998  40.040 116.985  40.044 116.975  40.047 116.963
     40.053 116.952
     40.050 116.940  40.049 116.928  40.048 116.915  40.047 116.903
     40.045 116.890
     40.043 116.878  40.042 116.866  40.043 116.854  40.047 116.843
     40.044 116.830
     40.041 116.819  40.037 116.808  40.034 116.796  40.032 116.784
     40.028 116.778
     40.019 116.776  40.010 116.772  40.000 116.767  39.992 116.760
     39.983 116.756
     39.974 116.751  39.964 116.750  39.956 116.757  39.950 116.767
     39.937 116.768
     39.924 116.769  39.913 116.773  39.905 116.777  39.900 116.782
     39.892 116.790
     39.889 116.798  39.887 116.810  39.882 116.819  39.877 116.830
     39.868 116.841
     39.860 116.849  39.854 116.866  39.845 116.879  39.837 116.899
     39.834 116.912
     39.832 116.919  39.825 116.919  39.816 116.923  39.807 116.929
     39.798 116.934
     39.789 116.938  39.780 116.934  39.772 116.928  39.766 116.918
     39.762 116.907
        200       2  39.756  39.447 116.901 116.180
     39.756 116.898  39.747 116.899  39.737 116.901  39.729 116.895
     39.721 116.887
     39.712 116.884  39.704 116.890  39.697 116.900  39.687 116.900
     39.678 116.896
     39.673 116.886  39.668 116.875  39.666 116.863  39.662 116.852
     39.654 116.845
     39.645 116.840  39.637 116.834  39.627 116.831  39.619 116.825
     39.613 116.816
     39.607 116.807  39.607 116.796  39.600 116.789  39.594 116.780
     39.595 116.770
     39.602 116.760  39.608 116.751  39.617 116.736  39.614 116.724
     39.608 116.715
     39.599 116.718  39.590 116.717  39.585 116.707  39.585 116.694
     39.590 116.684
     39.594 116.673  39.596 116.661  39.597 116.648  39.599 116.636
     39.604 116.625
     39.612 116.617  39.620 116.609  39.622 116.597  39.622 116.584
     39.621 116.572
     39.617 116.560  39.613 116.550  39.607 116.539  39.602 116.529
     39.595 116.519
     39.586 116.517  39.576 116.515  39.568 116.510  39.560 116.503
     39.552 116.494
     39.546 116.484  39.539 116.475  39.532 116.465  39.527 116.455
     39.524 116.443
     39.522 116.432  39.519 116.420  39.513 116.410  39.505 116.404
     39.495 116.404
     39.488 116.413  39.484 116.424  39.480 116.435  39.471 116.440
     39.462 116.443
     39.449 116.443  39.451 116.431  39.454 116.419  39.457 116.408
     39.458 116.395
     39.457 116.382  39.453 116.371  39.451 116.359  39.447 116.348
     39.450 116.336
     39.457 116.326  39.463 116.317  39.469 116.307  39.473 116.296
     39.477 116.285
     39.481 116.274  39.485 116.263  39.492 116.254  39.499 116.246
     39.511 116.240
     39.520 116.245  39.530 116.247  39.540 116.246  39.549 116.242
     39.557 116.236
     39.565 116.228  39.573 116.223  39.576 116.205  39.577 116.193
     39.576 116.180
        200       2  39.866  39.496 116.167 115.424
     39.576 116.167  39.574 116.154  39.568 116.144  39.563 116.134
     39.563 116.122
     39.563 116.109  39.565 116.097  39.566 116.084  39.567 116.072
     39.568 116.059
     39.570 116.047  39.572 116.035  39.577 116.024  39.584 116.015
     39.586 116.003
     39.591 115.994  39.590 115.982  39.583 115.977  39.574 115.971
     39.571 115.960
     39.573 115.948  39.580 115.939  39.588 115.932  39.593 115.921
     39.590 115.910
     39.581 115.906  39.572 115.902  39.564 115.896  39.557 115.887
     39.553 115.876
     39.550 115.865  39.547 115.853  39.544 115.842  39.539 115.831
     39.534 115.821
     39.525 115.814  39.516 115.811  39.508 115.804  39.503 115.794
     39.498 115.783
     39.496 115.772  39.500 115.760  39.504 115.753  39.511 115.744
     39.520 115.740
     39.529 115.736  39.537 115.730  39.544 115.721  39.549 115.710
     39.553 115.699
     39.561 115.691  39.570 115.687  39.578 115.681  39.586 115.673
     39.593 115.665
     39.595 115.653  39.594 115.640  39.593 115.627  39.592 115.615
     39.591 115.602
     39.592 115.590  39.595 115.578  39.597 115.566  39.601 115.555
     39.607 115.545
     39.611 115.534  39.615 115.523  39.623 115.516  39.631 115.508
     39.637 115.498
     39.643 115.488  39.650 115.479  39.660 115.477  39.670 115.477
     39.680 115.476
     39.690 115.476  39.700 115.475  39.709 115.473  39.719 115.470
     39.728 115.466
     39.736 115.461  39.743 115.452  39.747 115.441  39.754 115.431
     39.762 115.424
     39.771 115.427  39.777 115.437  39.781 115.448  39.786 115.458
     39.791 115.468
     39.795 115.479  39.798 115.491  39.806 115.497  39.815 115.494
     39.825 115.492
     39.835 115.494  39.843 115.501  39.849 115.510  39.858 115.515
     39.866 115.521
        200       2  40.416  39.875 115.951 115.418
     39.875 115.518  39.884 115.513  39.894 115.513  39.904 115.513
     39.912 115.507
     39.919 115.498  39.924 115.487  39.928 115.476  39.933 115.465
     39.938 115.455
     39.943 115.444  39.948 115.433  39.950 115.422  39.962 115.418
     39.972 115.420
     39.981 115.424  39.989 115.428  39.998 115.432  40.008 115.434
     40.017 115.439
     40.025 115.446  40.033 115.454  40.040 115.462  40.048 115.471
     40.056 115.478
     40.062 115.488  40.067 115.498  40.074 115.508  40.080 115.516
     40.084 115.529
     40.089 115.540  40.094 115.551  40.097 115.562  40.097 115.575
     40.099 115.588
     40.102 115.599  40.113 115.599  40.118 115.606  40.120 115.619
     40.122 115.631
     40.125 115.642  40.127 115.655  40.131 115.665  40.137 115.675
     40.135 115.687
     40.133 115.699  40.132 115.712  40.134 115.725  40.137 115.736
     40.143 115.746
     40.151 115.753  40.161 115.754  40.171 115.755  40.179 115.762
     40.179 115.774
     40.174 115.785  40.164 115.788  40.155 115.791  40.150 115.801
     40.148 115.814
     40.141 115.823  40.144 115.835  40.151 115.845  40.160 115.839
     40.170 115.838
     40.180 115.840  40.188 115.846  40.194 115.856  40.201 115.865
     40.210 115.868
     40.219 115.872  40.226 115.881  40.234 115.888  40.242 115.895
     40.251 115.900
     40.254 115.911  40.255 115.924  40.260 115.935  40.266 115.944
     40.275 115.951
     40.281 115.942  40.290 115.937  40.299 115.932  40.308 115.928
     40.316 115.923
     40.325 115.919  40.334 115.915  40.344 115.913  40.352 115.908
     40.358 115.899
     40.356 115.887  40.357 115.874  40.361 115.863  40.369 115.854
     40.376 115.845
     40.383 115.836  40.392 115.829  40.400 115.823  40.408 115.817
     40.416 115.808
        200       2  40.844  40.423 116.454 115.728
     40.423 115.799  40.432 115.793  40.439 115.785  40.447 115.776
     40.453 115.766
     40.462 115.761  40.471 115.765  40.481 115.767  40.491 115.767
     40.500 115.763
     40.503 115.751  40.504 115.738  40.511 115.728  40.519 115.730
     40.527 115.735
     40.535 115.745  40.542 115.752  40.551 115.757  40.559 115.764
     40.563 115.775
     40.564 115.788  40.567 115.800  40.572 115.810  40.580 115.818
     40.589 115.821
     40.595 115.830  40.595 115.843  40.598 115.855  40.602 115.866
     40.607 115.876
     40.613 115.887  40.617 115.898  40.619 115.910  40.618 115.922
     40.613 115.933
     40.608 115.943  40.601 115.952  40.593 115.960  40.587 115.970
     40.587 115.982
     40.591 115.993  40.598 116.000  40.601 116.011  40.605 116.022
     40.609 116.033
     40.613 116.045  40.617 116.056  40.621 116.067  40.624 116.078
     40.627 116.090
     40.632 116.101  40.642 116.104  40.651 116.100  40.659 116.106
     40.665 116.116
     40.668 116.127  40.672 116.138  40.678 116.149  40.686 116.154
     40.696 116.153
     40.706 116.154  40.715 116.156  40.721 116.166  40.728 116.176
     40.734 116.186
     40.739 116.197  40.747 116.205  40.755 116.212  40.764 116.216
     40.773 116.218
     40.782 116.222  40.790 116.229  40.790 116.241  40.785 116.252
     40.780 116.262
     40.773 116.272  40.768 116.282  40.773 116.294  40.778 116.304
     40.780 116.316
     40.780 116.329  40.778 116.342  40.777 116.354  40.780 116.365
     40.785 116.376
     40.788 116.388  40.790 116.401  40.783 116.408  40.774 116.412
     40.772 116.423
     40.776 116.434  40.782 116.445  40.788 116.454  40.801 116.451
     40.808 116.443
     40.816 116.434  40.824 116.427  40.831 116.417  40.838 116.408
     40.844 116.398
        200       2  41.066  40.725 116.956 116.329
     40.850 116.388  40.857 116.379  40.867 116.377  40.875 116.371
     40.883 116.362
     40.889 116.352  40.896 116.342  40.903 116.334  40.912 116.329
     40.922 116.330
     40.931 116.334  40.939 116.341  40.947 116.348  40.944 116.355
     40.939 116.366
     40.932 116.374  40.923 116.380  40.916 116.388  40.912 116.400
     40.913 116.412
     40.913 116.426  40.912 116.438  40.912 116.451  40.915 116.463
     40.922 116.467
     40.931 116.464  40.939 116.458  40.948 116.453  40.957 116.450
     40.967 116.448
     40.977 116.449  40.981 116.459  40.984 116.471  40.986 116.483
     40.987 116.496
     40.990 116.508  40.991 116.521  40.992 116.534  40.991 116.547
     40.988 116.558
     40.986 116.571  40.988 116.583  40.992 116.595  40.998 116.605
     41.007 116.607
     41.017 116.605  41.027 116.604  41.036 116.601  41.046 116.600
     41.056 116.602
     41.062 116.613  41.066 116.624  41.064 116.637  41.062 116.648
     41.058 116.660
     41.054 116.671  41.047 116.680  41.038 116.684  41.028 116.684
     41.018 116.686
     41.010 116.680  41.000 116.677  40.991 116.675  40.981 116.679
     40.973 116.684
     40.964 116.687  40.955 116.691  40.947 116.698  40.939 116.706
     40.931 116.713
     40.921 116.714  40.912 116.717  40.905 116.727  40.901 116.738
     40.897 116.749
     40.890 116.759  40.883 116.767  40.875 116.775  40.868 116.784
     40.859 116.790
     40.852 116.798  40.848 116.809  40.846 116.822  40.844 116.834
     40.839 116.845
     40.831 116.853  40.824 116.860  40.816 116.868  40.810 116.871
     40.805 116.886
     40.797 116.891  40.788 116.898  40.781 116.906  40.773 116.915
     40.765 116.920
     40.756 116.925  40.748 116.931  40.740 116.940  40.732 116.947
     40.725 116.956
        200       2  40.717  40.359 117.494 116.964
     40.717 116.964  40.711 116.974  40.709 116.987  40.709 117.000
     40.708 117.013
     40.708 117.026  40.710 117.039  40.711 117.051  40.714 117.063
     40.715 117.076
     40.716 117.088  40.716 117.102  40.717 117.115  40.715 117.127
     40.712 117.139
     40.709 117.151  40.705 117.162  40.704 117.175  40.702 117.187
     40.700 117.199
     40.696 117.211  40.693 117.222  40.691 117.234  40.689 117.247
     40.685 117.258
     40.681 117.269  40.677 117.281  40.671 117.291  40.667 117.302
     40.667 117.314
     40.671 117.325  40.678 117.331  40.686 117.338  40.690 117.349
     40.694 117.361
     40.697 117.372  40.699 117.385  40.698 117.398  40.697 117.411
     40.696 117.424
     40.693 117.435  40.689 117.447  40.688 117.460  40.685 117.472
     40.683 117.484
     40.678 117.494  40.646 117.489  40.649 117.477  40.652 117.466
     40.654 117.453
     40.653 117.440  40.648 117.430  40.638 117.428  40.629 117.428
     40.620 117.430
     40.611 117.434  40.601 117.436  40.592 117.438  40.583 117.432
     40.581 117.420
     40.581 117.408  40.580 117.395  40.576 117.384  40.572 117.372
     40.571 117.359
     40.578 117.351  40.583 117.341  40.586 117.329  40.585 117.316
     40.581 117.304
     40.579 117.292  40.576 117.280  40.572 117.269  40.568 117.258
     40.563 117.247
     40.556 117.239  40.547 117.236  40.539 117.242  40.530 117.245
     40.523 117.237
     40.520 117.226  40.520 117.213  40.512 117.206  40.502 117.206
     40.493 117.210
     40.484 117.214  40.475 117.218  40.466 117.223  40.457 117.227
     40.448 117.229
     40.439 117.232  40.428 117.232  40.419 117.229  40.410 117.224
     40.401 117.220
     40.390 117.213  40.383 117.223  40.375 117.231  40.367 117.239
     40.359 117.245
         36       2  40.351  40.231 117.379 117.252
     40.351 117.252  40.343 117.257  40.335 117.264  40.326 117.268
     40.317 117.274
     40.309 117.280  40.300 117.285  40.294 117.295  40.288 117.305
     40.282 117.315
     40.274 117.324  40.266 117.328  40.257 117.333  40.249 117.340
     40.242 117.348
     40.236 117.359  40.233 117.370  40.231 117.379
        200       2  27.251  27.083 120.426 120.213
     27.241 120.213  27.242 120.221  27.243 120.233  27.243 120.240
     27.248 120.244
     27.250 120.250  27.251 120.257  27.246 120.263  27.238 120.262
     27.232 120.272
     27.230 120.280  27.225 120.286  27.217 120.290  27.211 120.295
     27.209 120.304
     27.209 120.311  27.216 120.314  27.217 120.323  27.217 120.334
     27.219 120.346
     27.221 120.356  27.220 120.366  27.214 120.367  27.208 120.368
     27.201 120.371
     27.193 120.369  27.189 120.364  27.186 120.359  27.184 120.353
     27.182 120.347
     27.178 120.343  27.171 120.349  27.168 120.360  27.171 120.365
     27.176 120.368
     27.179 120.373  27.170 120.381  27.167 120.386  27.163 120.390
     27.158 120.393
     27.153 120.396  27.149 120.401  27.147 120.406  27.146 120.413
     27.145 120.419
     27.144 120.426  27.134 120.424  27.138 120.412  27.139 120.406
     27.137 120.400
     27.134 120.391  27.125 120.399  27.122 120.403  27.119 120.399
     27.120 120.392
     27.120 120.386  27.121 120.379  27.120 120.373  27.117 120.368
     27.111 120.374
     27.105 120.372  27.099 120.368  27.094 120.367  27.094 120.379
     27.089 120.386
     27.085 120.381  27.084 120.375  27.083 120.369  27.084 120.363
     27.086 120.356
     27.088 120.351  27.090 120.345  27.091 120.339  27.090 120.333
     27.090 120.326
     27.092 120.321  27.097 120.318  27.104 120.312  27.102 120.306
     27.097 120.301
     27.095 120.295  27.094 120.289  27.097 120.284  27.100 120.278
     27.102 120.272
     27.105 120.267  27.109 120.263  27.114 120.260  27.119 120.256
     27.122 120.252
     27.123 120.245  27.118 120.243  27.112 120.246  27.108 120.250
     27.103 120.249
     27.099 120.245  27.095 120.240  27.090 120.244  27.087 120.254
     27.089 120.260
        200       2  27.091  26.877 120.278 120.076
     27.091 120.266  27.085 120.271  27.079 120.277  27.072 120.277
     27.066 120.276
     27.061 120.278  27.056 120.275  27.045 120.278  27.046 120.270
     27.055 120.265
     27.058 120.254  27.058 120.248  27.058 120.241  27.052 120.241
     27.047 120.245
     27.045 120.252  27.039 120.251  27.040 120.245  27.037 120.239
     27.031 120.235
     27.029 120.241  27.029 120.248  27.028 120.255  27.025 120.260
     27.020 120.264
     27.015 120.262  27.010 120.258  27.005 120.257  27.000 120.259
     26.996 120.264
     26.991 120.263  26.991 120.254  26.985 120.245  26.981 120.240
     26.980 120.234
     26.979 120.227  26.979 120.221  26.979 120.215  26.980 120.209
     26.984 120.204
     26.984 120.197  26.981 120.192  26.983 120.186  26.984 120.180
     26.983 120.174
     26.980 120.169  26.974 120.168  26.972 120.180  26.969 120.185
     26.963 120.186
     26.962 120.198  26.964 120.203  26.967 120.209  26.970 120.214
     26.972 120.220
     26.971 120.226  26.968 120.232  26.957 120.232  26.951 120.230
     26.946 120.230
     26.941 120.234  26.937 120.238  26.932 120.235  26.927 120.232
     26.922 120.229
     26.917 120.225  26.915 120.220  26.920 120.217  26.922 120.208
     26.920 120.202
     26.918 120.196  26.917 120.189  26.918 120.183  26.921 120.178
     26.924 120.172
     26.926 120.166  26.928 120.161  26.929 120.154  26.927 120.149
     26.924 120.144
     26.926 120.137  26.922 120.132  26.924 120.123  26.926 120.118
     26.925 120.112
     26.920 120.108  26.915 120.110  26.911 120.115  26.906 120.119
     26.901 120.119
     26.897 120.115  26.893 120.109  26.890 120.104  26.885 120.101
     26.880 120.098
     26.877 120.093  26.877 120.087  26.880 120.081  26.885 120.078
     26.890 120.076
        200       2  26.897  26.620 120.151 120.018
     26.894 120.071  26.893 120.065  26.897 120.054  26.896 120.047
     26.890 120.046
     26.888 120.051  26.885 120.057  26.880 120.062  26.876 120.052
     26.877 120.045
     26.877 120.039  26.875 120.033  26.873 120.027  26.869 120.023
     26.864 120.020
     26.858 120.018  26.853 120.019  26.847 120.021  26.843 120.024
     26.838 120.028
     26.832 120.029  26.826 120.036  26.826 120.043  26.826 120.049
     26.825 120.056
     26.824 120.062  26.823 120.068  26.821 120.074  26.819 120.080
     26.815 120.085
     26.810 120.086  26.804 120.088  26.799 120.090  26.797 120.096
     26.800 120.101
     26.803 120.107  26.803 120.113  26.798 120.116  26.801 120.126
     26.797 120.129
     26.792 120.124  26.789 120.119  26.784 120.116  26.779 120.113
     26.774 120.110
     26.770 120.106  26.766 120.102  26.763 120.096  26.758 120.094
     26.753 120.096
     26.748 120.100  26.745 120.105  26.742 120.111  26.742 120.117
     26.745 120.123
     26.748 120.128  26.742 120.134  26.737 120.135  26.736 120.142
     26.737 120.148
     26.731 120.151  26.727 120.148  26.722 120.144  26.719 120.139
     26.717 120.133
     26.715 120.127  26.713 120.121  26.709 120.116  26.707 120.111
     26.704 120.105
     26.700 120.101  26.692 120.103  26.687 120.106  26.681 120.105
     26.675 120.104
     26.671 120.108  26.667 120.113  26.657 120.110  26.654 120.116
     26.655 120.123
     26.649 120.125  26.644 120.121  26.640 120.116  26.637 120.111
     26.634 120.106
     26.630 120.101  26.626 120.097  26.622 120.092  26.620 120.086
     26.621 120.080
     26.627 120.078  26.632 120.080  26.636 120.084  26.642 120.087
     26.647 120.085
     26.650 120.080  26.648 120.074  26.644 120.069  26.641 120.064
     26.637 120.059
        200       2  26.718  26.521 120.054 119.823
     26.634 120.054  26.631 120.049  26.628 120.043  26.626 120.038
     26.623 120.032
     26.620 120.027  26.617 120.022  26.613 120.018  26.610 120.012
     26.608 120.006
     26.607 120.000  26.605 119.995  26.603 119.989  26.604 119.982
     26.605 119.976
     26.607 119.970  26.608 119.964  26.609 119.958  26.610 119.952
     26.607 119.947
     26.605 119.942  26.603 119.935  26.593 119.932  26.586 119.923
     26.577 119.928
     26.568 119.935  26.567 119.929  26.564 119.924  26.558 119.922
     26.556 119.911
     26.554 119.905  26.550 119.901  26.546 119.897  26.542 119.892
     26.537 119.889
     26.532 119.888  26.527 119.885  26.523 119.880  26.522 119.874
     26.522 119.868
     26.521 119.857  26.521 119.851  26.524 119.845  26.526 119.839
     26.527