Advertisement
Guest User

Untitled

a guest
Sep 20th, 2017
59
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.45 KB | None | 0 0
  1. MODULE m
  2. IMPLICIT NONE
  3.  
  4. REAL, ALLOCATABLE :: &
  5. pa2c (:)
  6. ! REAL, ALLOCATABLE :: &
  7. ! g_arr(:), &
  8. ! g_arrinp1(:, :), &
  9. ! g_arrinp2(:, :)
  10. ! g_arr_const_5_5(5, 5)
  11.  
  12. CONTAINS
  13. SUBROUTINE inv_so(ki1sd, ki1ed, &
  14. ki3sd, ki3ed, &
  15. ki1sc, ki1ec, &
  16. ki3sc, ki3ec, &
  17. pflfu)
  18. INTEGER, INTENT(IN) :: &
  19. ki1sd, ki1ed, &
  20. ki3sd, ki3ed, &
  21. ki1sc, ki1ec, &
  22. ki3sc, ki3ec
  23. REAL, INTENT(OUT) :: &
  24. ! pflfu (ki1sd:ki1ed,ki3sd:ki3ed+1) ! flux up cloud-free
  25. pflfu (ki1sd:ki1ed,ki3sd:ki3ed) ! flux up cloud-free
  26. INTEGER :: &
  27. j1, j3
  28.  
  29. DO j1 = ki1sc, ki1ec
  30. pa2c(j1) = 10
  31. END DO
  32.  
  33. do j3 = ki3sc+1, ki3ec
  34. DO j1 = ki1sc, ki1ec
  35. pa2c(j1) = 10
  36. END DO
  37. END DO
  38.  
  39. DO j1 = ki1sc, ki3ec
  40. pflfu(3, ki3ec+1) = j1
  41. END DO
  42.  
  43. END SUBROUTINE inv_so
  44. END MODULE m
  45.  
  46. PROGRAM p
  47. USE m
  48. INTEGER :: si, ei
  49. REAL, ALLOCATABLE :: &
  50. pflfu(:,:), pflcu(:,:)
  51.  
  52. si=1
  53. ei=10
  54.  
  55. ALLOCATE(pa2c (si-4:ei+1))
  56. ALLOCATE(pflfu (si-4:ei+1,si-4:ei+1))
  57. ALLOCATE(pflcu (si-4:ei+1, si-4:ei+1))
  58.  
  59. CALL inv_so (si, ei, &
  60. si, ei, &
  61. si, ei, &
  62. si, ei, &
  63. pflfu)
  64. END PROGRAM p
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement