← Back to branch summary

~fluidity-core/fluidity/sediment_restructure

« back to all changes in this revision

Viewing changes to femtools/VTK_interfaces.F90

merge with branch including remap_surface_to_field but then removed this change and approached it a different way which seems to work

Show diffs side-by-side

added added

removed removed

Lines of Context:
femtools/VTK_interfaces.F90
443
443
        ewrite(-1,*) "Please switch to a non-periodic output_mesh."
444
444
        FLExit("Just remapped from an unperiodic to a periodic continuous field!")
445
445
      end if
 
446
    else if ((lstat/=0).and. &
 
447
             (lstat/=REMAP_ERR_BUBBLE_LAGRANGE).and. &
 
448
             (lstat/=REMAP_ERR_HIGHER_LOWER_CONTINUOUS)) then
 
449
      if(present(stat)) then
 
450
        stat = lstat
 
451
      else
 
452
        FLAbort("Unknown error when remapping coordinates while outputting to vtu.")
 
453
      end if
446
454
    end if
447
455
    ! we've just allowed remapping from a higher order to a lower order continuous field
448
456
 
483
491
              else
484
492
                FLAbort("Just remapped from an unperiodic to a periodic continuous field!")
485
493
              end if
 
494
            else if ((lstat/=0).and. &
 
495
                     (lstat/=REMAP_ERR_BUBBLE_LAGRANGE).and. &
 
496
                     (lstat/=REMAP_ERR_HIGHER_LOWER_CONTINUOUS)) then
 
497
              if(present(stat)) then
 
498
                stat = lstat
 
499
              else
 
500
                FLAbort("Unknown error when remapping field.")
 
501
              end if
486
502
            end if
487
503
            ! we've just allowed remapping from a higher order to a lower order continuous field
488
504
            
574
590
              else
575
591
                FLAbort("Just remapped from an unperiodic to a periodic continuous field!")
576
592
              end if
 
593
            else if ((lstat/=0).and. &
 
594
                     (lstat/=REMAP_ERR_BUBBLE_LAGRANGE).and. &
 
595
                     (lstat/=REMAP_ERR_HIGHER_LOWER_CONTINUOUS)) then
 
596
              if(present(stat)) then
 
597
                stat = lstat
 
598
              else
 
599
                FLAbort("Unknown error when remapping field.")
 
600
              end if
577
601
            end if
578
602
            ! we've just allowed remapping from a higher order to a lower order continuous field
579
603
          
655
679
              else
656
680
                FLAbort("Just remapped from an unperiodic to a periodic continuous field!")
657
681
              end if
 
682
            else if ((lstat/=0).and. &
 
683
                     (lstat/=REMAP_ERR_BUBBLE_LAGRANGE).and. &
 
684
                     (lstat/=REMAP_ERR_HIGHER_LOWER_CONTINUOUS)) then
 
685
              if(present(stat)) then
 
686
                stat = lstat
 
687
              else
 
688
                FLAbort("Unknown error when remapping field.")
 
689
              end if
658
690
            end if
659
691
            ! we've just allowed remapping from a higher order to a lower order continuous field
660
692