trunk/dyn3d/startdyn.f

MODULE startdyn

  ! From startvar.F, version 1.4
  ! January 27th, 2006 15:14:22

  INTEGER iml_dyn, jml_dyn, llm_dyn

  REAL, allocatable:: lon_ini(:), lat_ini(:)
  ! longitude and latitude from the input file, converted to rad

  real, allocatable:: levdyn_ini(:)

CONTAINS

  SUBROUTINE start_init_dyn(tsol_2d, phis, ps)

    use comgeom, only: aire_2d, apoln, apols
    use conf_dat2d_m, only: conf_dat2d
    use dimens_m, only: iim, jjm
    use dynetat0_m, only: rlonu, rlatv
    use gr_int_dyn_m, only: gr_int_dyn
    use inter_barxy_m, only: inter_barxy
    use netcdf, only: nf90_nowrite
    use netcdf95, only: nf95_open, nf95_close, nf95_get_var, nf95_inq_varid, &
         nf95_gw_var, find_coord
    use nr_util, only: assert, pi

    REAL, intent(in):: tsol_2d(:, :) ! (iim + 1, jjm + 1)

    REAL, intent(in):: phis(:, :) ! (iim + 1, jjm + 1)
    ! surface geopotential, in m2 s-2

    REAL, intent(out):: ps(:, :) ! (iim + 1, jjm + 1) surface pressure, in Pa

    ! Local:

    INTEGER ncid, varid
    REAL, ALLOCATABLE:: lon_rad(:), lat_rad(:)
    REAL, ALLOCATABLE:: var_ana(:, :)
    real z(iim + 1, jjm + 1)
    real tmp_var(iim, jjm + 1)

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

    print *, "Call sequence information: start_init_dyn"
    call assert((/size(tsol_2d, 1), size(phis, 1), size(ps, 1)/) == iim + 1, &
         "start_init_dyn size 1")
    call assert((/size(tsol_2d, 2), size(phis, 2), size(ps, 2)/) == jjm + 1, &
         "start_init_dyn size 2")

    call nf95_open('ECDYN.nc', nf90_nowrite, ncid)

    call find_coord(ncid, varid=varid, std_name="longitude")
    call nf95_gw_var(ncid, varid, lon_ini)
    lon_ini = lon_ini * pi / 180. ! convert to rad
    iml_dyn = size(lon_ini)
    print *, "iml_dyn = ", iml_dyn

    call find_coord(ncid, varid=varid, std_name="latitude")
    call nf95_gw_var(ncid, varid, lat_ini)
    lat_ini = lat_ini * pi / 180. ! convert to rad
    jml_dyn = size(lat_ini)
    print *, "jml_dyn = ", jml_dyn

    call nf95_inq_varid(ncid, "level", varid)
    call nf95_gw_var(ncid, varid, levdyn_ini)
    llm_dyn = size(levdyn_ini)
    print *, "llm_dyn = ", llm_dyn

    ALLOCATE(var_ana(iml_dyn, jml_dyn), lon_rad(iml_dyn), lat_rad(jml_dyn))

    ! 'Z': Surface geopotential
    call nf95_inq_varid(ncid, 'Z', varid)
    call nf95_get_var(ncid, varid, var_ana)
    CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, var_ana)
    CALL inter_barxy(lon_rad, lat_rad(:jml_dyn -1), var_ana, rlonu(:iim), &
         rlatv, tmp_var)
    z = gr_int_dyn(tmp_var)

    ! 'SP': Surface pressure
    call nf95_inq_varid(ncid, 'SP', varid)
    call nf95_get_var(ncid, varid, var_ana)
    CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, var_ana)
    CALL inter_barxy(lon_rad, lat_rad(:jml_dyn -1), var_ana, rlonu(:iim), &
         rlatv, tmp_var) 
    ps = gr_int_dyn(tmp_var)

    call nf95_close(ncid)

    ! Adapt the surface pressure to "phis", with the hypsometric equation:
    ps(:iim, :) = ps(:iim, :) &
         * (1. + (z(:iim, :) - phis(:iim, :)) / 287. / tsol_2d(:iim, :))

    ps(iim + 1, :) = ps(1, :)
    ps(:, 1) = SUM(aire_2d(:iim, 1) * ps(:iim, 1)) / apoln
    ps(:, jjm + 1) = SUM(aire_2d(:iim, jjm + 1) * ps(:iim, jjm + 1)) / apols

  END SUBROUTINE start_init_dyn

END MODULE startdyn
1	guez	3	MODULE startdyn
2
3			! From startvar.F, version 1.4
4	guez	73	! January 27th, 2006 15:14:22
5	guez	3
6	guez	54	INTEGER iml_dyn, jml_dyn, llm_dyn
7	guez	3
8	guez	225	REAL, allocatable:: lon_ini(:), lat_ini(:)
9	guez	73	! longitude and latitude from the input file, converted to rad
10	guez	3
11	guez	225	real, allocatable:: levdyn_ini(:)
12	guez	3
13			CONTAINS
14
15	guez	78	SUBROUTINE start_init_dyn(tsol_2d, phis, ps)
16	guez	3
17	guez	139	use comgeom, only: aire_2d, apoln, apols
18	guez	3	use conf_dat2d_m, only: conf_dat2d
19			use dimens_m, only: iim, jjm
20	guez	139	use dynetat0_m, only: rlonu, rlatv
21	guez	3	use gr_int_dyn_m, only: gr_int_dyn
22	guez	73	use inter_barxy_m, only: inter_barxy
23			use netcdf, only: nf90_nowrite
24			use netcdf95, only: nf95_open, nf95_close, nf95_get_var, nf95_inq_varid, &
25			nf95_gw_var, find_coord
26			use nr_util, only: assert, pi
27	guez	3
28	guez	43	REAL, intent(in):: tsol_2d(:, :) ! (iim + 1, jjm + 1)
29	guez	3
30	guez	78	REAL, intent(in):: phis(:, :) ! (iim + 1, jjm + 1)
31			! surface geopotential, in m2 s-2
32
33			REAL, intent(out):: ps(:, :) ! (iim + 1, jjm + 1) surface pressure, in Pa
34
35	guez	3	! Local:
36
37	guez	73	INTEGER ncid, varid
38	guez	3	REAL, ALLOCATABLE:: lon_rad(:), lat_rad(:)
39	guez	36	REAL, ALLOCATABLE:: var_ana(:, :)
40			real z(iim + 1, jjm + 1)
41	guez	3	real tmp_var(iim, jjm + 1)
42
43			!--------------------------
44
45			print *, "Call sequence information: start_init_dyn"
46	guez	78	call assert((/size(tsol_2d, 1), size(phis, 1), size(ps, 1)/) == iim + 1, &
47	guez	42	"start_init_dyn size 1")
48	guez	78	call assert((/size(tsol_2d, 2), size(phis, 2), size(ps, 2)/) == jjm + 1, &
49	guez	42	"start_init_dyn size 2")
50	guez	54
51	guez	73	call nf95_open('ECDYN.nc', nf90_nowrite, ncid)
52	guez	3
53	guez	73	call find_coord(ncid, varid=varid, std_name="longitude")
54			call nf95_gw_var(ncid, varid, lon_ini)
55			lon_ini = lon_ini * pi / 180. ! convert to rad
56			iml_dyn = size(lon_ini)
57			print *, "iml_dyn = ", iml_dyn
58	guez	3
59	guez	73	call find_coord(ncid, varid=varid, std_name="latitude")
60			call nf95_gw_var(ncid, varid, lat_ini)
61			lat_ini = lat_ini * pi / 180. ! convert to rad
62			jml_dyn = size(lat_ini)
63			print *, "jml_dyn = ", jml_dyn
64	guez	3
65	guez	73	call nf95_inq_varid(ncid, "level", varid)
66			call nf95_gw_var(ncid, varid, levdyn_ini)
67			llm_dyn = size(levdyn_ini)
68			print *, "llm_dyn = ", llm_dyn
69	guez	3
70	guez	73	ALLOCATE(var_ana(iml_dyn, jml_dyn), lon_rad(iml_dyn), lat_rad(jml_dyn))
71	guez	3
72			! 'Z': Surface geopotential
73	guez	48	call nf95_inq_varid(ncid, 'Z', varid)
74			call nf95_get_var(ncid, varid, var_ana)
75	guez	3	CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, var_ana)
76			CALL inter_barxy(lon_rad, lat_rad(:jml_dyn -1), var_ana, rlonu(:iim), &
77			rlatv, tmp_var)
78	guez	36	z = gr_int_dyn(tmp_var)
79	guez	3
80			! 'SP': Surface pressure
81	guez	48	call nf95_inq_varid(ncid, 'SP', varid)
82			call nf95_get_var(ncid, varid, var_ana)
83	guez	3	CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, var_ana)
84			CALL inter_barxy(lon_rad, lat_rad(:jml_dyn -1), var_ana, rlonu(:iim), &
85			rlatv, tmp_var)
86	guez	78	ps = gr_int_dyn(tmp_var)
87	guez	3
88	guez	48	call nf95_close(ncid)
89
90	guez	78	! Adapt the surface pressure to "phis", with the hypsometric equation:
91			ps(:iim, :) = ps(:iim, :) &
92	guez	36	* (1. + (z(:iim, :) - phis(:iim, :)) / 287. / tsol_2d(:iim, :))
93	guez	3
94	guez	78	ps(iim + 1, :) = ps(1, :)
95			ps(:, 1) = SUM(aire_2d(:iim, 1) * ps(:iim, 1)) / apoln
96			ps(:, jjm + 1) = SUM(aire_2d(:iim, jjm + 1) * ps(:iim, jjm + 1)) / apols
97	guez	3
98			END SUBROUTINE start_init_dyn
99
100			END MODULE startdyn