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	guez	173	module nat2gcm_m
2	guez	37
3	guez	172	implicit none
4	guez	37
5	guez	173	contains
6	guez	37
7	guez	173	subroutine nat2gcm(pk, u, v, t)
8	guez	37
9	guez	173	! Passage aux variables du mod\`ele (vents covariants,
10			! temp\'erature potentielle et humidit\'e sp\'ecifique).
11	guez	37
12	guez	178	use comconst, only: cpp
13	guez	173	use comgeom, only: cu_2d, cv_2d
14			use dimens_m, only: iim, jjm, llm
15			use paramet_m, only: iip1, jjp1
16	guez	37
17	guez	173	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	guez	37
21	guez	173	! Local:
22			integer i, j, l
23	guez	37
24	guez	173	!----------------------------------------------------------------------
25	guez	37
26	guez	173	print *, "Call sequence information: nat2gcm"
27	guez	37
28	guez	173	! 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	guez	37
46	guez	173	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	guez	37
55	guez	173	end subroutine nat2gcm
56	guez	37
57	guez	173	end module nat2gcm_m