Guide/Read_reanalyse/nat2gcm.f

module nat2gcm_m

  implicit none

contains

  subroutine nat2gcm(pk, u, v, t)

    ! Passage aux variables du mod\`ele (vents covariants,
    ! temp\'erature potentielle et humidit\'e sp\'ecifique).

    use comconst, only: cpp
    use comgeom, only: cu_2d, cv_2d
    use dimens_m, only: iim, jjm, llm
    use paramet_m, only: iip1, jjp1

    real, intent(in):: pk(iip1, jjp1, llm)
    real, intent(inout):: u(iip1, jjp1, llm), v(iip1, jjm, llm)
    real, intent(inout):: t(iip1, jjp1, llm)

    ! Local:
    integer i, j, l

    !----------------------------------------------------------------------

    print *, "Call sequence information: nat2gcm"

    ! calcul de ucov et de la temperature potentielle
    do l = 1, llm
       do j = 1, jjp1
          do i = 1, iim
             u(i, j, l) = u(i, j, l) * cu_2d(i, j)
             t(i, j, l) = t(i, j, l) * cpp / pk(i, j, l)
          enddo
          u(iip1, j, l) = u(1, j, l)
          t(iip1, j, l) = t(1, j, l)
       enddo
       do i = 1, iip1
          u(i, 1, l) = 0.
          u(i, jjp1, l) = 0.
          t(i, 1, l) = t(1, 1, l)
          t(i, jjp1, l) = t(1, jjp1, l)
       enddo
    enddo

    do l = 1, llm
       do j = 1, jjm
          do i = 1, iim
             v(i, j, l) = v(i, j, l) * cv_2d(i, j)
          enddo
          v(iip1, j, l) = v(1, j, l)
       enddo
    enddo

  end subroutine nat2gcm

end module nat2gcm_m
1	module nat2gcm_m
2
3	implicit none
4
5	contains
6
7	subroutine nat2gcm(pk, u, v, t)
8
9	! Passage aux variables du mod\`ele (vents covariants,
10	! temp\'erature potentielle et humidit\'e sp\'ecifique).
11
12	use comconst, only: cpp
13	use comgeom, only: cu_2d, cv_2d
14	use dimens_m, only: iim, jjm, llm
15	use paramet_m, only: iip1, jjp1
16
17	real, intent(in):: pk(iip1, jjp1, llm)
18	real, intent(inout):: u(iip1, jjp1, llm), v(iip1, jjm, llm)
19	real, intent(inout):: t(iip1, jjp1, llm)
20
21	! Local:
22	integer i, j, l
23
24	!----------------------------------------------------------------------
25
26	print *, "Call sequence information: nat2gcm"
27
28	! calcul de ucov et de la temperature potentielle
29	do l = 1, llm
30	do j = 1, jjp1
31	do i = 1, iim
32	u(i, j, l) = u(i, j, l) * cu_2d(i, j)
33	t(i, j, l) = t(i, j, l) * cpp / pk(i, j, l)
34	enddo
35	u(iip1, j, l) = u(1, j, l)
36	t(iip1, j, l) = t(1, j, l)
37	enddo
38	do i = 1, iip1
39	u(i, 1, l) = 0.
40	u(i, jjp1, l) = 0.
41	t(i, 1, l) = t(1, 1, l)
42	t(i, jjp1, l) = t(1, jjp1, l)
43	enddo
44	enddo
45
46	do l = 1, llm
47	do j = 1, jjm
48	do i = 1, iim
49	v(i, j, l) = v(i, j, l) * cv_2d(i, j)
50	enddo
51	v(iip1, j, l) = v(1, j, l)
52	enddo
53	enddo
54
55	end subroutine nat2gcm
56
57	end module nat2gcm_m