Guide/Read_reanalyse/read_reanalyse.f

module read_reanalyse_m

  IMPLICIT NONE

contains

  subroutine read_reanalyse(psi, u, v, t, q)

    ! From LMDZ4/libf/dyn3d/read_reanalyse.F, version 1.3, 2005/04/15 12:31:21

    USE conf_guide_m, ONLY: guide_q, guide_t, guide_u, guide_v
    USE dimens_m, ONLY: iim, jjm, llm
    use nat2gcm_m, only: nat2gcm
    USE netcdf, ONLY: nf90_nowrite
    USE netcdf95, ONLY: nf95_get_var, nf95_inq_dimid, nf95_inq_varid, &
         nf95_inquire_dimension, nf95_open, find_coord
    USE paramet_m, ONLY: iip1, jjp1
    use reanalyse2nat_m, only: reanalyse2nat

    real, intent(in):: psi(:, :) ! (iip1, jjp1)
    real, intent(out):: u(:, :, :) ! (iip1, jjp1, llm)
    real, intent(out):: v(:, :, :) ! (iip1, jjm, llm)
    real, intent(out):: t(:, :, :), q(:, :, :) ! (iip1, jjp1, llm)

    ! Local:
    integer nlevnc
    integer:: timestep = 0
    real pk(iip1, jjp1, llm)
    integer, save:: ncidu, varidu, ncidv, varidv, ncidt, varidt, ncidQ, varidQ
    integer ncid, varid, dimid
    real, allocatable, save:: unc(:, :, :) ! (iip1, jjp1, nlevnc)
    real, allocatable, save:: vnc(:, :, :) ! (iip1, jjm, nlevnc)
    real, allocatable, save:: tnc(:, :, :), Qnc(:, :, :) ! (iip1, jjp1, nlevnc)
    real, allocatable, save:: pl(:) ! (nlevnc)
    real latitude(jjm + 1)
    logical:: first = .true.
    logical, save:: invert_y

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

    ! Initialisation de la lecture des fichiers

    if (first) then
       print *, 'Intitialisation de read reanalyse'

       ! Vent zonal
       if (guide_u) then
          call nf95_open('u.nc', nf90_nowrite, ncidu)
          call nf95_inq_varid(ncidu, 'UWND', varidu)
       endif

       ! Vent meridien
       if (guide_v) then
          call nf95_open('v.nc', nf90_nowrite, ncidv)
          call nf95_inq_varid(ncidv, 'VWND', varidv)
       endif

       ! Temperature
       if (guide_T) then
          call nf95_open('T.nc', nf90_nowrite, ncidt)
          call nf95_inq_varid(ncidt, 'AIR', varidt)
       endif

       ! Humidite
       if (guide_Q) then
          call nf95_open('hur.nc', nf90_nowrite, ncidQ)
          call nf95_inq_varid(ncidQ, 'RH', varidQ)
       endif

       ! Coordonn\'ee verticale :

       if (guide_u) then
          ncid = ncidu
       else if (guide_v) then
          ncid = ncidv
       else if (guide_T) then
          ncid = ncidt
       else
          ncid = ncidq
       end if

       call find_coord(ncid, dimid = dimid, varid = varid, std_name = "plev")
       call nf95_inquire_dimension(ncid, dimid, nclen = nlevnc)
       PRINT *, 'nlevnc = ', nlevnc
       allocate(unc(iip1, jjp1, nlevnc), vnc(iip1, jjm, nlevnc))
       allocate(tnc(iip1, jjp1, nlevnc), Qnc(iip1, jjp1, nlevnc), pl(nlevnc))
       call NF95_GET_VAR(ncid, varid, pl)
       pl = 100. * pl ! passage en pascal

       ! Read latitude values just to know their order:
       call find_coord(ncid, varid = varid, std_name = "latitude")
       call nf95_get_var(ncid, varid, latitude)
       invert_y = latitude(1) < latitude(2)

       first = .false.
    endif

    ! lecture des champs u, v, T, q

    timestep = timestep + 1

    ! Vent zonal
    if (guide_u) then
       call NF95_GET_VAR(ncidu, varidu, unc, start = (/1, 1, 1, timestep/))
    else
       unc = 0.
    end if

    ! Temperature
    if (guide_T) then
       call NF95_GET_VAR(ncidt, varidt, tnc, start = (/1, 1, 1, timestep/))
    else
       tnc = 0.
    end if

    ! Humidite
    if (guide_Q) then
       call NF95_GET_VAR(ncidQ, varidQ, Qnc, start = (/1, 1, 1, timestep/))
    else
       Qnc = 0.
    end if

    ! Vent meridien
    if (guide_v) then
       call NF95_GET_VAR(ncidv, varidv, vnc, start = (/1, 1, 1, timestep/))
    else
       vnc = 0.
    end if

    call reanalyse2nat(invert_y, psi, unc, vnc, tnc, Qnc, pl, u, v, t, Q, pk)
    call nat2gcm(pk, u, v, t)

  end subroutine read_reanalyse

end module read_reanalyse_m
1	module read_reanalyse_m
2
3	IMPLICIT NONE
4
5	contains
6
7	subroutine read_reanalyse(psi, u, v, t, q)
8
9	! From LMDZ4/libf/dyn3d/read_reanalyse.F, version 1.3, 2005/04/15 12:31:21
10
11	USE conf_guide_m, ONLY: guide_q, guide_t, guide_u, guide_v
12	USE dimens_m, ONLY: iim, jjm, llm
13	use nat2gcm_m, only: nat2gcm
14	USE netcdf, ONLY: nf90_nowrite
15	USE netcdf95, ONLY: nf95_get_var, nf95_inq_dimid, nf95_inq_varid, &
16	nf95_inquire_dimension, nf95_open, find_coord
17	USE paramet_m, ONLY: iip1, jjp1
18	use reanalyse2nat_m, only: reanalyse2nat
19
20	real, intent(in):: psi(:, :) ! (iip1, jjp1)
21	real, intent(out):: u(:, :, :) ! (iip1, jjp1, llm)
22	real, intent(out):: v(:, :, :) ! (iip1, jjm, llm)
23	real, intent(out):: t(:, :, :), q(:, :, :) ! (iip1, jjp1, llm)
24
25	! Local:
26	integer nlevnc
27	integer:: timestep = 0
28	real pk(iip1, jjp1, llm)
29	integer, save:: ncidu, varidu, ncidv, varidv, ncidt, varidt, ncidQ, varidQ
30	integer ncid, varid, dimid
31	real, allocatable, save:: unc(:, :, :) ! (iip1, jjp1, nlevnc)
32	real, allocatable, save:: vnc(:, :, :) ! (iip1, jjm, nlevnc)
33	real, allocatable, save:: tnc(:, :, :), Qnc(:, :, :) ! (iip1, jjp1, nlevnc)
34	real, allocatable, save:: pl(:) ! (nlevnc)
35	real latitude(jjm + 1)
36	logical:: first = .true.
37	logical, save:: invert_y
38
39	! -----------------------------------------------------------------
40
41	! Initialisation de la lecture des fichiers
42
43	if (first) then
44	print *, 'Intitialisation de read reanalyse'
45
46	! Vent zonal
47	if (guide_u) then
48	call nf95_open('u.nc', nf90_nowrite, ncidu)
49	call nf95_inq_varid(ncidu, 'UWND', varidu)
50	endif
51
52	! Vent meridien
53	if (guide_v) then
54	call nf95_open('v.nc', nf90_nowrite, ncidv)
55	call nf95_inq_varid(ncidv, 'VWND', varidv)
56	endif
57
58	! Temperature
59	if (guide_T) then
60	call nf95_open('T.nc', nf90_nowrite, ncidt)
61	call nf95_inq_varid(ncidt, 'AIR', varidt)
62	endif
63
64	! Humidite
65	if (guide_Q) then
66	call nf95_open('hur.nc', nf90_nowrite, ncidQ)
67	call nf95_inq_varid(ncidQ, 'RH', varidQ)
68	endif
69
70	! Coordonn\'ee verticale :
71
72	if (guide_u) then
73	ncid = ncidu
74	else if (guide_v) then
75	ncid = ncidv
76	else if (guide_T) then
77	ncid = ncidt
78	else
79	ncid = ncidq
80	end if
81
82	call find_coord(ncid, dimid = dimid, varid = varid, std_name = "plev")
83	call nf95_inquire_dimension(ncid, dimid, nclen = nlevnc)
84	PRINT *, 'nlevnc = ', nlevnc
85	allocate(unc(iip1, jjp1, nlevnc), vnc(iip1, jjm, nlevnc))
86	allocate(tnc(iip1, jjp1, nlevnc), Qnc(iip1, jjp1, nlevnc), pl(nlevnc))
87	call NF95_GET_VAR(ncid, varid, pl)
88	pl = 100. * pl ! passage en pascal
89
90	! Read latitude values just to know their order:
91	call find_coord(ncid, varid = varid, std_name = "latitude")
92	call nf95_get_var(ncid, varid, latitude)
93	invert_y = latitude(1) < latitude(2)
94
95	first = .false.
96	endif
97
98	! lecture des champs u, v, T, q
99
100	timestep = timestep + 1
101
102	! Vent zonal
103	if (guide_u) then
104	call NF95_GET_VAR(ncidu, varidu, unc, start = (/1, 1, 1, timestep/))
105	else
106	unc = 0.
107	end if
108
109	! Temperature
110	if (guide_T) then
111	call NF95_GET_VAR(ncidt, varidt, tnc, start = (/1, 1, 1, timestep/))
112	else
113	tnc = 0.
114	end if
115
116	! Humidite
117	if (guide_Q) then
118	call NF95_GET_VAR(ncidQ, varidQ, Qnc, start = (/1, 1, 1, timestep/))
119	else
120	Qnc = 0.
121	end if
122
123	! Vent meridien
124	if (guide_v) then
125	call NF95_GET_VAR(ncidv, varidv, vnc, start = (/1, 1, 1, timestep/))
126	else
127	vnc = 0.
128	end if
129
130	call reanalyse2nat(invert_y, psi, unc, vnc, tnc, Qnc, pl, u, v, t, Q, pk)
131	call nat2gcm(pk, u, v, t)
132
133	end subroutine read_reanalyse
134
135	end module read_reanalyse_m