hmc-dev/debug_2dVar.f90 at main · c-hydro/hmc-dev · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
!-----------------------------------------------------------------------------------------------------------------------------------
! Subroutine to analyze and plot 2d variable(s) used in Hydrological Model Continuum
! Version 2.0.0
! Date 2023-06-05
!
! PYTHON SIDE:
! USE THE "HMC_DEBUG_LAUNCHER.SH" TO INSTALL THE DEPS AND RUN THE DEBUGGER
!
! MATLAB SIDE:
! COMMAND LINE (SAME FOLDER):   matlab -nosplash -nodesktop -r "debug_2dVar()"
! COMMAND LINE (OTHER FOLDER):  matlab -nosplash -nodesktop -r "debug_2dVar('/home/fabio/NetBeansProjects/Fortran_Continuum_RegMarche/')"
!
! FORTRAN SIDE:
! COMMAND LINE (CODE SIDE GENERIC): call debug_2dVar(dble(var_2d_data), rows, cols, figure, interpreter)
! figure=[1 ... ]; interpreter=[1-matlab, 2-python]
!
! EXAMPLE(S):
! COMMAND LINE (CODE SIDE FOR MATLAB): call debug_2dVar(dble(var_2d_data), rows, cols, 1, 1)
! COMMAND LINE (CODE SIDE FOR PYTHON): call debug_2dVar(dble(var_2d_data), rows, cols, 1, 2)
!-----------------------------------------------------------------------------------------------------------------------------------

!-----------------------------------------------------------------------------------------------------------------------------------
! subroutine to write var 2d in format columns
subroutine debug_2dVar(z, irows, icols, input, interpreter)

    !implicit none
    real(kind=8), intent(in)            :: z(irows, icols)
    real(kind=8)			:: xyz(3,size(z(:,1)),size(z(1,:)))
    integer				:: nx, ny, nrow
    integer                             :: i, j, file_unit, icols, irows
    integer                             :: ios, ierror
    integer ,optional                   :: input, interpreter

    character(len=100)                  :: data_file_name, string_input

    if (interpreter == 1) then
        print *,'debug_2dVar - write data using matlab interpreter.'
    else if (interpreter == 2) then
        print *,'debug_2dVar - write data using python interpreter.'
    else
        print *,'debug_2dVar - fatal error! Interpreter code is not allowed.'
        stop
    end if

    nx=size(z(:,1))
    ny=size(z(1,:))
    nrow=nx*ny

    do i=1,nx
        do j=1,ny
            xyz(1,i,j)=dble(i)
            xyz(2,i,j)=dble(j)
            xyz(3,i,j)=z(i,j)
        end do
    end do

    write(string_input,'(I2.2)') input

    if (present(input)) then
        data_file_name='data_file_'//trim(string_input)//'.txt'
    else
        data_file_name='data_file.txt'
    end if

    ierror=0
    call get_unit(file_unit)
    if (file_unit==0) then
        ierror=1
        print *,'debug_2dVar - fatal error! Could not get a free FORTRAN unit.'
        stop
    end if
    open (unit=file_unit, file=data_file_name, status='replace', iostat=ios)
    if (ios/=0) then
        ierror=2
        print *,'debug_2dVar - fatal error! Could not open the terminal data file.'
        stop
    end if

    if (interpreter == 1) then

        do j=1,ny
            do i=1,nx
                write (file_unit,'(3E15.7)') xyz(1:3,i,j)
            end do
        end do

    else if (interpreter == 2) then

        do j=1,ny
            do i=1,nx
                write (file_unit, "(*(G0,:,','))") xyz(1:3,i,j)
            end do
        end do

    end if

    close (unit=file_unit)

end subroutine debug_2dVar
!-----------------------------------------------------------------------------------------------------------------------------------


!-----------------------------------------------------------------------------------------------------------------------------------
! subroutine to get unit number automatically
subroutine get_unit(iunit)

    !implicit none
    integer i
    integer ios
    integer iunit
    logical lopen

    iunit=0
    do i=1,99
            if (i/= 5 .and. i/=6) then
                    inquire (unit=i, opened=lopen, iostat=ios)
                    if (ios==0) then
                            if (.not.lopen) then
                                    iunit=i
                                    return
                            end if
                    end if

            end if
    end do
    return

end subroutine get_unit
!-----------------------------------------------------------------------------------------------------------------------------------