[MITgcm-support] output files from the FLT package to netcdf
Camille Mazoyer
mazoyer at univ-tln.fr
Tue May 2 05:59:10 EDT 2017
Dear all,
I would like to compare read drifters trajectories (surface buoys) to my
simulation. For that, I see that I can use flt package. So I start with
the flt example and try to understand how it works.
I would like to transform outputs float_trajectories.* from the example,
to netcdf files. I see that they are two ways to do that:
- cvfloat.F90 from the extra directory of flt example
- a python toolbox named floater
I'm trying both manners, and I encounter problems with both ... For
floater toolbox, I'm asking for help in they github url. For
cvfloat.f90, I would like to know if the script should work with the
example without doing nothing? Is it the last version? flt output are in
float32 (according to the meta files), so I had to change some double
precision into single precision in cvfloat.F90, and I change a flag for
the grid:
For the flt example, theses lines:
logical :: usingSphericalPolarGrid=.true.
real*8, dimension(:), allocatable :: tmp
ilen=IMAX*8
become :
logical :: usingSphericalPolarGrid=.false.
real*4, dimension(:), allocatable :: tmp
ilen=IMAX*4
Now, the error for me is when I try to read the first record tmp. If I
write tmp at line
tmp: 7.1838967E-41 3.4893733E-41 1.7646863E-38 3.4893733E-41
2.3048557E-41
4.3986759E-42 2.8789677E-41 0.0000000E+00 0.0000000E+00 0.0000000E+00
0.0000000E+00 0.0000000E+00 0.0000000E+00
Where, according to the matlab script read_flt_traj.m, I should read
100, 3600, 10800, 3600, 3, 140, 13, 0, 0, 0, 0, 0, 0
I am stuck here, so any advices is welcome,
Thank you very much for your help,
Camille
--
------------------------------------------
Camille Mazoyer
Phd Student
Mediterranean Institute of Oceanography (MIO)
Institut de Mathématiques de Toulon (IMATH)
Université de TOULON
Bat X - CS 60584
83041 TOULON cedex 9
France
tel: +33.4.94.14.24.50
Email: mazoyer at univ-tln.fr
http://mio.pytheas.univ-amu.fr/
http://imath.fr/
More information about the MITgcm-support
mailing list