[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