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, pointer:: lon_ini(:), lat_ini(:)
  ! longitude and latitude from the input file, converted to rad

  real, pointer:: levdyn_ini(:)

CONTAINS

  SUBROUTINE start_init_dyn(tsol_2d, psol)

    use comgeom, only: rlonu, rlatv, aire_2d, apoln, apols
    use conf_dat2d_m, only: conf_dat2d
    use dimens_m, only: iim, jjm
    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
    use start_init_orog_m, only: phis

    REAL, intent(in):: tsol_2d(:, :) ! (iim + 1, jjm + 1)
    REAL, intent(out):: psol(:, :) ! (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(psol, 1)/) == iim + 1, &
         "start_init_dyn size 1")
    call assert((/size(tsol_2d, 2), size(psol, 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) 
    psol = gr_int_dyn(tmp_var)

    call nf95_close(ncid)

    psol(:iim, :) = psol(:iim, :) &
         * (1. + (z(:iim, :) - phis(:iim, :)) / 287. / tsol_2d(:iim, :))
    psol(iim + 1, :) = psol(1, :)

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