bugfix in out_movie_file.f90

proper values for n_r_mov_tot

Author: tgastine

src/movie.f90

@@ -27,7 +27,7 @@ module movie_data
    character(len=80), public :: movie(n_movies_max)  ! Only for input
    character(len=72), public :: movie_file(n_movies_max)
 
-   logical, public :: lGeosField(n_movies_max)
+   logical, public :: lICField(n_movies_max), lGeosField(n_movies_max)
    logical :: lPhaseField(n_movies_max)
    integer, public :: n_movies
    integer, public :: n_movie_surface(n_movies_max)
@@ -1119,6 +1119,7 @@ subroutine get_movie_type
          !--- Now store the necessary information:
          !------ Increase number of movies:
          n_movies=n_movies+1
+         lICField(n_movies)=lIC
          lGeosField(n_movies)=lGeos
          lPhaseField(n_movies)=lPhase
 

src/out_movie_file.f90

@@ -13,7 +13,7 @@ module out_movie
    use movie_data, only: frames, n_movie_fields, n_movies, n_movie_surface, &
        &                 n_movie_const, n_movie_field_type, lGeosField,     &
        &                 n_movie_field_start,n_movie_field_stop,            &
-       &                 movie_const, movie_file, n_movie_file,             &
+       &                 movie_const, movie_file, n_movie_file, lICField,   &
        &                 n_movie_fields_ic
    use radial_data, only: n_r_icb, n_r_cmb
    use radial_functions, only: orho1, orho2, or1, or2, or3, or4, beta,  &
@@ -214,7 +214,7 @@ subroutine write_movie_frame(n_frame, time)
             n_r_mov_tot=n_r_max
             r_mov_tot(:n_r_max)=r(:)
             if ( n_r_ic_max > 0 ) then
-               n_r_mov_tot=n_r_mov_tot+n_r_ic_max
+               if ( lICField(n_movie) ) n_r_mov_tot=n_r_mov_tot+n_r_ic_max
                r_mov_tot(n_r_max+1:)=r_ic(:)
             end if