/[lmdze]/trunk/dyn3d/etat0.f
ViewVC logotype

Diff of /trunk/dyn3d/etat0.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 97 by guez, Fri Apr 25 14:58:31 2014 UTC revision 99 by guez, Wed Jul 2 18:39:15 2014 UTC
# Line 52  contains Line 52  contains
52      use start_init_phys_m, only: start_init_phys      use start_init_phys_m, only: start_init_phys
53      use start_inter_3d_m, only: start_inter_3d      use start_inter_3d_m, only: start_inter_3d
54      use temps, only: itau_phy, annee_ref, day_ref      use temps, only: itau_phy, annee_ref, day_ref
55        use test_disvert_m, only: test_disvert
56    
57      ! Variables local to the procedure:      ! Variables local to the procedure:
58    
# Line 68  contains Line 69  contains
69      ! and pressure level "pls(i, j, l)".)      ! and pressure level "pls(i, j, l)".)
70    
71      real qsat(iim + 1, jjm + 1, llm) ! mass fraction of saturating water vapor      real qsat(iim + 1, jjm + 1, llm) ! mass fraction of saturating water vapor
72      REAL tsol(klon), qsol(klon), sn(klon)      REAL sn(klon)
73      REAL tsolsrf(klon, nbsrf), qsolsrf(klon, nbsrf), snsrf(klon, nbsrf)      REAL qsolsrf(klon, nbsrf), snsrf(klon, nbsrf)
74      REAL albe(klon, nbsrf), evap(klon, nbsrf)      REAL albe(klon, nbsrf), evap(klon, nbsrf)
75      REAL alblw(klon, nbsrf)      REAL alblw(klon, nbsrf)
76      REAL tsoil(klon, nsoilmx, nbsrf)      REAL tsoil(klon, nsoilmx, nbsrf)
77      REAL radsol(klon), rain_fall(klon), snow_fall(klon)      REAL radsol(klon), rain_fall(klon), snow_fall(klon)
78      REAL solsw(klon), sollw(klon), fder(klon)      REAL solsw(klon), sollw(klon), fder(klon)
79      !IM "slab" ocean      !IM "slab" ocean
     REAL tslab(klon)  
80      real seaice(klon) ! kg m-2      real seaice(klon) ! kg m-2
81      REAL frugs(klon, nbsrf), agesno(klon, nbsrf)      REAL frugs(klon, nbsrf), agesno(klon, nbsrf)
82      REAL rugmer(klon)      REAL rugmer(klon)
# Line 131  contains Line 131  contains
131    
132      pa = 5e4      pa = 5e4
133      CALL disvert      CALL disvert
134        call test_disvert
135      CALL inigeom      CALL inigeom
136      CALL inifilr      CALL inifilr
137    
# Line 208  contains Line 209  contains
209         q(:, :, :, 5) = q(:, :, :, 5) * 48. / 29.         q(:, :, :, 5) = q(:, :, :, 5) * 48. / 29.
210      end if      end if
211    
     tsol = pack(tsol_2d, dyn_phy)  
     qsol = pack(qsol_2d, dyn_phy)  
212      sn = 0. ! snow      sn = 0. ! snow
213      radsol = 0.      radsol = 0.
     tslab = 0. ! IM "slab" ocean  
214      seaice = 0.      seaice = 0.
215      rugmer = 0.001      rugmer = 0.001
216      zmea = pack(zmea_2d, dyn_phy)      zmea = pack(zmea_2d, dyn_phy)
# Line 317  contains Line 315  contains
315      CALL dynredem1("start.nc", vcov, ucov, teta, q, masse, ps, itau=0)      CALL dynredem1("start.nc", vcov, ucov, teta, q, masse, ps, itau=0)
316    
317      ! Initialisations :      ! Initialisations :
     tsolsrf(:, is_ter) = tsol  
     tsolsrf(:, is_lic) = tsol  
     tsolsrf(:, is_oce) = tsol  
     tsolsrf(:, is_sic) = tsol  
318      snsrf(:, is_ter) = sn      snsrf(:, is_ter) = sn
319      snsrf(:, is_lic) = sn      snsrf(:, is_lic) = sn
320      snsrf(:, is_oce) = sn      snsrf(:, is_oce) = sn
# Line 331  contains Line 325  contains
325      albe(:, is_sic) = 0.6      albe(:, is_sic) = 0.6
326      alblw = albe      alblw = albe
327      evap = 0.      evap = 0.
328      qsolsrf(:, is_ter) = 150.      qsolsrf = 150.
329      qsolsrf(:, is_lic) = 150.      tsoil = spread(spread(pack(tsol_2d, dyn_phy), 2, nsoilmx), 3, nbsrf)
     qsolsrf(:, is_oce) = 150.  
     qsolsrf(:, is_sic) = 150.  
     tsoil = spread(spread(tsol, 2, nsoilmx), 3, nbsrf)  
330      rain_fall = 0.      rain_fall = 0.
331      snow_fall = 0.      snow_fall = 0.
332      solsw = 165.      solsw = 165.
# Line 343  contains Line 334  contains
334      t_ancien = 273.15      t_ancien = 273.15
335      q_ancien = 0.      q_ancien = 0.
336      agesno = 0.      agesno = 0.
     !IM "slab" ocean  
     tslab = tsolsrf(:, is_oce)  
337      seaice = 0.      seaice = 0.
338    
339      frugs(:, is_oce) = rugmer      frugs(:, is_oce) = rugmer
# Line 359  contains Line 348  contains
348      sig1 = 0.      sig1 = 0.
349      w01 = 0.      w01 = 0.
350    
351      call phyredem("startphy.nc", latfi, lonfi, pctsrf, &      call phyredem("startphy.nc", latfi, lonfi, pctsrf, tsoil(:, 1, :), tsoil, &
352           tsolsrf, tsoil, tslab, seaice, qsolsrf, qsol, snsrf, albe, alblw, &           tsoil(:, 1, is_oce), seaice, qsolsrf, pack(qsol_2d, dyn_phy), snsrf, &
353           evap, rain_fall, snow_fall, solsw, sollw, fder, radsol, frugs, &           albe, alblw, evap, rain_fall, snow_fall, solsw, sollw, fder, radsol, &
354           agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, &           frugs, agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, t_ancien, &
355           t_ancien, q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01)           q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01)
356      CALL histclo      CALL histclo
357    
358    END SUBROUTINE etat0    END SUBROUTINE etat0
Legend:
Removed from v.97  
changed lines
  Added in v.99
  ViewVC Help
Powered by ViewVC 1.1.21