--- trunk/dyn3d/etat0.f 2014/04/25 14:58:31 97 +++ trunk/Sources/dyn3d/etat0.f 2015/07/08 17:03:45 155 @@ -13,26 +13,27 @@ contains - SUBROUTINE etat0 + SUBROUTINE etat0(phis) ! From "etat0_netcdf.F", version 1.3, 2005/05/25 13:10:09 use caldyn0_m, only: caldyn0 use comconst, only: cpp, kappa, iniconst - use comgeom, only: rlatu, rlonv, rlonu, rlatv, aire_2d, apoln, apols, & - cu_2d, cv_2d, inigeom - use conf_gcm_m, only: dayref, anneeref + use comgeom, only: aire_2d, apoln, apols, cu_2d, cv_2d, inigeom use dimens_m, only: iim, jjm, llm, nqmx use dimphy, only: zmasq use dimsoil, only: nsoilmx use disvert_m, only: ap, bp, preff, pa, disvert + use dynetat0_m, only: day_ref, annee_ref, xprimp025, xprimm025, rlatu1, & + rlatu2, rlatu, rlatv, yprimu1, yprimu2, rlonu, rlonv, xprimu, xprimv use dynredem0_m, only: dynredem0 use dynredem1_m, only: dynredem1 use exner_hyb_m, only: exner_hyb + use fxhyp_m, only: fxhyp + use fyhyp_m, only: fyhyp use geopot_m, only: geopot use grid_atob, only: grille_m use grid_change, only: init_dyn_phy, dyn_phy - use histclo_m, only: histclo use indicesol, only: is_oce, is_sic, is_ter, is_lic, epsfra use iniadvtrac_m, only: iniadvtrac use inifilr_m, only: inifilr @@ -42,22 +43,23 @@ nf95_inq_varid, nf95_open use nr_util, only: pi, assert use paramet_m, only: ip1jm, ip1jmp1 + use phyetat0_m, only: rlat, rlon use phyredem_m, only: phyredem use q_sat_m, only: q_sat use regr_lat_time_coefoz_m, only: regr_lat_time_coefoz use regr_pr_o3_m, only: regr_pr_o3 - use serre, only: alphax use startdyn, only: start_init_dyn USE start_init_orog_m, only: start_init_orog, mask use start_init_phys_m, only: start_init_phys use start_inter_3d_m, only: start_inter_3d - use temps, only: itau_phy, annee_ref, day_ref + use temps, only: itau_phy + use test_disvert_m, only: test_disvert + use unit_nml_m, only: unit_nml - ! Variables local to the procedure: + REAL, intent(out):: phis(:, :) ! (iim + 1, jjm + 1) + ! surface geopotential, in m2 s-2 - REAL latfi(klon), lonfi(klon) - ! (latitude and longitude of a point of the scalar grid identified - ! by a simple index, in degrees) + ! Local: REAL, dimension(iim + 1, jjm + 1, llm):: ucov, t3d, teta REAL vcov(iim + 1, jjm, llm) @@ -68,19 +70,17 @@ ! and pressure level "pls(i, j, l)".) real qsat(iim + 1, jjm + 1, llm) ! mass fraction of saturating water vapor - REAL tsol(klon), qsol(klon), sn(klon) - REAL tsolsrf(klon, nbsrf), qsolsrf(klon, nbsrf), snsrf(klon, nbsrf) + REAL sn(klon) + REAL qsolsrf(klon, nbsrf), snsrf(klon, nbsrf) REAL albe(klon, nbsrf), evap(klon, nbsrf) REAL alblw(klon, nbsrf) REAL tsoil(klon, nsoilmx, nbsrf) REAL radsol(klon), rain_fall(klon), snow_fall(klon) REAL solsw(klon), sollw(klon), fder(klon) !IM "slab" ocean - REAL tslab(klon) real seaice(klon) ! kg m-2 REAL frugs(klon, nbsrf), agesno(klon, nbsrf) REAL rugmer(klon) - REAL phis(iim + 1, jjm + 1) ! surface geopotential, in m2 s-2 real, dimension(iim + 1, jjm + 1):: zmea_2d, zstd_2d, zsig_2d, zgam_2d real, dimension(iim + 1, jjm + 1):: zthe_2d, zpic_2d, zval_2d real, dimension(iim + 1, jjm + 1):: tsol_2d, qsol_2d, ps @@ -121,28 +121,42 @@ ! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)", ! for interface "l") + namelist /etat0_nml/ day_ref, annee_ref + !--------------------------------- print *, "Call sequence information: etat0" + print *, "Enter namelist 'etat0_nml'." + read(unit=*, nml=etat0_nml) + write(unit_nml, nml=etat0_nml) + CALL iniconst ! Construct a grid: pa = 5e4 CALL disvert + call test_disvert + + CALL fyhyp(rlatu, rlatv, rlatu2, yprimu2, rlatu1, yprimu1) + CALL fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025) + + rlatu(1) = pi / 2. + rlatu(jjm + 1) = -rlatu(1) + CALL inigeom CALL inifilr - latfi(1) = 90. - latfi(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.) * 180. / pi + rlat(1) = 90. + rlat(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.) * 180. / pi ! (with conversion to degrees) - latfi(klon) = - 90. + rlat(klon) = - 90. - lonfi(1) = 0. - lonfi(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.) * 180. / pi + rlon(1) = 0. + rlon(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.) * 180. / pi ! (with conversion to degrees) - lonfi(klon) = 0. + rlon(klon) = 0. call start_init_orog(phis, zmea_2d, zstd_2d, zsig_2d, zgam_2d, zthe_2d, & zpic_2d, zval_2d) ! also compute "mask" @@ -208,11 +222,8 @@ q(:, :, :, 5) = q(:, :, :, 5) * 48. / 29. end if - tsol = pack(tsol_2d, dyn_phy) - qsol = pack(qsol_2d, dyn_phy) sn = 0. ! snow radsol = 0. - tslab = 0. ! IM "slab" ocean seaice = 0. rugmer = 0.001 zmea = pack(zmea_2d, dyn_phy) @@ -296,8 +307,6 @@ ! Calcul interm\'ediaire : CALL massdair(p3d, masse) - print *, 'ALPHAX = ', alphax - forall (l = 1:llm) masse(:, 1, l) = SUM(aire_2d(:iim, 1) * masse(:iim, 1, l)) / apoln masse(:, jjm + 1, l) = & @@ -307,20 +316,14 @@ ! Initialisation pour traceurs: call iniadvtrac itau_phy = 0 - day_ref = dayref - annee_ref = anneeref CALL geopot(teta, pk , pks, phis, phi) CALL caldyn0(ucov, vcov, teta, ps, masse, pk, phis, phi, w, pbaru, & pbarv) - CALL dynredem0("start.nc", dayref, phis) + CALL dynredem0("start.nc", day_ref, phis) CALL dynredem1("start.nc", vcov, ucov, teta, q, masse, ps, itau=0) ! Initialisations : - tsolsrf(:, is_ter) = tsol - tsolsrf(:, is_lic) = tsol - tsolsrf(:, is_oce) = tsol - tsolsrf(:, is_sic) = tsol snsrf(:, is_ter) = sn snsrf(:, is_lic) = sn snsrf(:, is_oce) = sn @@ -331,11 +334,8 @@ albe(:, is_sic) = 0.6 alblw = albe evap = 0. - qsolsrf(:, is_ter) = 150. - qsolsrf(:, is_lic) = 150. - qsolsrf(:, is_oce) = 150. - qsolsrf(:, is_sic) = 150. - tsoil = spread(spread(tsol, 2, nsoilmx), 3, nbsrf) + qsolsrf = 150. + tsoil = spread(spread(pack(tsol_2d, dyn_phy), 2, nsoilmx), 3, nbsrf) rain_fall = 0. snow_fall = 0. solsw = 165. @@ -343,8 +343,6 @@ t_ancien = 273.15 q_ancien = 0. agesno = 0. - !IM "slab" ocean - tslab = tsolsrf(:, is_oce) seaice = 0. frugs(:, is_oce) = rugmer @@ -359,12 +357,11 @@ sig1 = 0. w01 = 0. - call phyredem("startphy.nc", latfi, lonfi, pctsrf, & - tsolsrf, tsoil, tslab, seaice, qsolsrf, qsol, snsrf, albe, alblw, & - evap, rain_fall, snow_fall, solsw, sollw, fder, radsol, frugs, & - agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, & - t_ancien, q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01) - CALL histclo + call phyredem("startphy.nc", pctsrf, tsoil(:, 1, :), tsoil, & + tsoil(:, 1, is_oce), seaice, qsolsrf, pack(qsol_2d, dyn_phy), snsrf, & + albe, evap, rain_fall, snow_fall, solsw, sollw, fder, radsol, frugs, & + agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, t_ancien, & + q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01) END SUBROUTINE etat0