用Fortran编写以下程序，且不用trim函数：找出2~1000之间的全部同构数。所谓同构数是指这样一个数，它出现在它的平方数的右端，例如5的平方是25，且5出现在25的右端，则5是同构数:又如25 的平方是625，且25出现在625的右端，则25也是一个同构数。

在Fortran中编写一个程序来查找2到1000之间的所有同构数，我们可以采用循环遍历的方式，对每个数字进行操作，判断其平方后的字符串是否在原字符串的右侧。由于Fortran的标准库可能没有内置的`strfind`等类似功能来直接找到子串的位置，我们需要手动实现这个逻辑。这里是一个简单的示例：

program HomomorphicNumbers
    implicit none
    integer :: num, i, j, square, len_num
    character(len=4) :: num_str, square_str

    ! 遍历从2到1000
    do num = 2, 1000
        len_num = len_trim(num_ch_to_string(num)) ! 计算数字的长度
        square = num**2
        square_str = num_ch_to_string(square)

        ! 检查num在square_str的最后
        if (len_num == len(square_str) .or. &
            (len(square_str) > len_num & &
             square_str(len_num+1:) == num_ch_to_string(num))) then
            write(*,*) num, " is a homomorphic number."
        end if
    end do

contains
    function num_ch_to_string(n)
        integer, intent(in) :: n
        character(len=10) :: num_ch_to_string

        ! 将整数转换成字符数组
        num_ch_to_string = trim(int_to_char(n))
    end function num_ch_to_string
end program HomomorphicNumbers

! 注意：此代码片段假设Fortran编译器支持将整数转换为字符的功能（int_to_char），如果实际环境中没有这样的函数，你需要自行实现这个转换。

在这个程序中，我们首先定义了一个`num_ch_to_string`函数用于将整数转换为字符数组，然后在主循环中计算每个数及其平方的字符串，并检查它们是否满足同构数的条件。