search

IMAU-paleo / UFEMISM2.0

Version 2.0 of the Utrecht FinitE voluMe Ice-Sheet Model
6 stars 0 forks source link

Bug restart thermodynamics #22

Open erwinlambert opened 1 month ago

erwinlambert commented 1 month ago

I run into an error when trying to initialize a new run with temperatures read_from_file from a restart_thermodynamics_xxxxx.nc file:

Error termination. Backtrace:
At line 756 of file mesh/mesh_remapping.f90
Fortran runtime error: Array bound mismatch for dimension 1 of array 'd_dst_partial'
erwinlambert commented 1 month ago

In addition, also tested restarting from files either main_output_ANT_grid.nc with error message:

 ERROR: data fields are the wrong size! in UFEMISM_program/initialise_model_region/initialise_thermodynamics_model/initialise_ice_temperature/initialise_ice_temperature_from_file/read_field_from_file_3D/map_from_xy_grid_to_mesh_3D/apply_map_xy_grid_to_mesh_3D on process 43/64 (0 = master)
ERROR STOP

and main_output_ANT_xxxxx.nc with the same error message as for restart_thermodynamics_xxxxx.nc

TijnBerends commented 1 month ago

It might be that I set up the thermodynamics restart to only work when restarting with the same mesh (i.e. choice_initial_mesh = 'read_from_file'), you could take a look at that.

erwinlambert commented 1 month ago

Nope, that's what I did. In fact, it crashes after recognising the meshes are identical:

    ! If the two meshes are identical, the remapping operation is trivial
    CALL check_if_meshes_are_identical( mesh_src, mesh_dst, are_identical)
    IF (are_identical) THEN
      d_dst_partial = d_src_partial
      CALL finalise_routine( routine_name)
      RETURN
    END IF

in map_from_mesh_to_mesh_3D called by read_field_from_file_3D