trunk/phylmd/phyredem.f

module phyredem_m

  IMPLICIT NONE

contains

  SUBROUTINE phyredem(fichnom, rlat, rlon, pctsrf, tsol, tsoil, tslab, &
       seaice, qsurf, qsol, snow, albedo, 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)

    ! From phylmd/phyredem.F, v 1.3 2005/05/25 13:10:09
    ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
    ! Objet: Ecriture de l'etat de démarrage ou redémarrage pour la physique

    USE indicesol, ONLY : is_lic, is_oce, is_sic, is_ter, nbsrf
    USE dimphy, ONLY : klev, klon, zmasq
    USE dimsoil, ONLY : nsoilmx
    USE temps, ONLY : itau_phy
    USE netcdf95, ONLY : nf95_create, nf95_put_att
    USE netcdf, ONLY : nf90_clobber, nf90_global

    INCLUDE 'netcdf.inc'

    CHARACTER(len=*) fichnom
    REAL, INTENT (IN) :: rlat(klon), rlon(klon)
    REAL :: tsol(klon, nbsrf)
    REAL :: tsoil(klon, nsoilmx, nbsrf)

    REAL :: tslab(klon), seaice(klon) !IM "slab" ocean
    REAL :: qsurf(klon, nbsrf)
    REAL :: qsol(klon)
    REAL :: snow(klon, nbsrf)
    REAL :: albedo(klon, nbsrf)

    REAL :: alblw(klon, nbsrf)

    REAL :: evap(klon, nbsrf)
    REAL :: rain_fall(klon)
    REAL :: snow_fall(klon)
    REAL :: solsw(klon)
    REAL :: sollw(klon)
    REAL :: fder(klon)
    REAL :: radsol(klon)
    REAL :: frugs(klon, nbsrf)
    REAL :: agesno(klon, nbsrf)
    REAL :: zmea(klon)
    REAL, intent(in):: zstd(klon)
    REAL, intent(in):: zsig(klon)
    REAL :: zgam(klon)
    REAL :: zthe(klon)
    REAL :: zpic(klon)
    REAL :: zval(klon)
    REAL :: pctsrf(klon, nbsrf)
    REAL :: t_ancien(klon, klev), q_ancien(klon, klev)
    REAL :: clwcon(klon, klev), rnebcon(klon, klev), ratqs(klon, klev)
    REAL :: run_off_lic_0(klon)

    INTEGER :: nid, nvarid, idim2, idim3
    INTEGER :: ierr

    INTEGER :: isoil, nsrf
    CHARACTER (7) :: str7
    CHARACTER (2) :: str2

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

    PRINT *, 'Call sequence information: phyredem'
    CALL nf95_create(fichnom, nf90_clobber, nid)

    call nf95_put_att(nid, nf90_global, 'title', &
         'Fichier redémarrage physique')
    call nf95_put_att(nid, nf90_global, "itau_phy", itau_phy)

    ierr = nf_def_dim(nid, 'points_physiques', klon, idim2)
    ierr = nf_def_dim(nid, 'horizon_vertical', klon*klev, idim3)

    ierr = nf_def_var(nid, 'longitude', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 32, &
         'Longitudes de la grille physique')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, rlon)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'latitude', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 31, &
         'Latitudes de la grille physique')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, rlat)

    ! PB ajout du masque terre/mer

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'masque', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 16, 'masque terre mer')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, zmasq)
    ! BP ajout des fraction de chaque sous-surface

    ! 1. fraction de terre

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'FTER', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 21, 'fraction de continent')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, pctsrf(1:klon, is_ter))

    ! 2. Fraction de glace de terre

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'FLIC', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 24, 'fraction glace de terre')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, pctsrf(1:klon, is_lic))

    ! 3. fraction ocean

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'FOCE', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 14, 'fraction ocean')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, pctsrf(1:klon, is_oce))

    ! 4. Fraction glace de mer

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'FSIC', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 18, 'fraction glace mer')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, pctsrf(1:klon, is_sic))


    DO nsrf = 1, nbsrf
       IF (nsrf<=99) THEN
          WRITE (str2, '(i2.2)') nsrf
          ierr = nf_redef(nid)
          ierr = nf_def_var(nid, 'TS'//str2, nf_float, 1, idim2, nvarid)
          ierr = nf_put_att_text(nid, nvarid, 'title', 28, &
               'Temperature de surface No.'//str2)
          ierr = nf_enddef(nid)
       ELSE
          PRINT *, 'Trop de sous-mailles'
          STOP 1
       END IF
       ierr = nf_put_var_real(nid, nvarid, tsol(1, nsrf))
    END DO

    DO nsrf = 1, nbsrf
       DO isoil = 1, nsoilmx
          IF (isoil<=99 .AND. nsrf<=99) THEN
             WRITE (str7, '(i2.2, "srf", i2.2)') isoil, nsrf
             ierr = nf_redef(nid)
             ierr = nf_def_var(nid, 'Tsoil'//str7, nf_float, 1, idim2, nvarid)
             ierr = nf_put_att_text(nid, nvarid, 'title', 29, &
                  'Temperature du sol No.'//str7)
             ierr = nf_enddef(nid)
          ELSE
             PRINT *, 'Trop de couches'
             STOP 1
          END IF
          ierr = nf_put_var_real(nid, nvarid, tsoil(1, isoil, nsrf))
       END DO
    END DO

    !IM "slab" ocean
    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'TSLAB', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 33, &
         'Ecart de la SST (pour slab-ocean)')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, tslab)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'SEAICE', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 33, &
         'Glace de mer kg/m2 (pour slab-ocean)')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, seaice)

    DO nsrf = 1, nbsrf
       IF (nsrf<=99) THEN
          WRITE (str2, '(i2.2)') nsrf
          ierr = nf_redef(nid)
          ierr = nf_def_var(nid, 'QS'//str2, nf_float, 1, idim2, nvarid)
          ierr = nf_put_att_text(nid, nvarid, 'title', 25, &
               'Humidite de surface No.'//str2)
          ierr = nf_enddef(nid)
       ELSE
          PRINT *, 'Trop de sous-mailles'
          STOP 1
       END IF
       ierr = nf_put_var_real(nid, nvarid, qsurf(1, nsrf))
    END DO

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'QSOL', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 20, 'Eau dans le sol (mm)')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, qsol)

    DO nsrf = 1, nbsrf
       IF (nsrf<=99) THEN
          WRITE (str2, '(i2.2)') nsrf
          ierr = nf_redef(nid)
          ierr = nf_def_var(nid, 'ALBE'//str2, nf_float, 1, idim2, nvarid)
          ierr = nf_put_att_text(nid, nvarid, 'title', 23, &
               'albedo de surface No.'//str2)
          ierr = nf_enddef(nid)
       ELSE
          PRINT *, 'Trop de sous-mailles'
          STOP 1
       END IF
       ierr = nf_put_var_real(nid, nvarid, albedo(1, nsrf))
    END DO

    !IM BEG albedo LW
    DO nsrf = 1, nbsrf
       IF (nsrf<=99) THEN
          WRITE (str2, '(i2.2)') nsrf
          ierr = nf_redef(nid)
          ierr = nf_def_var(nid, 'ALBLW'//str2, nf_float, 1, idim2, nvarid)
          ierr = nf_put_att_text(nid, nvarid, 'title', 23, &
               'albedo LW de surface No.'//str2)
          ierr = nf_enddef(nid)
       ELSE
          PRINT *, 'Trop de sous-mailles'
          STOP 1
       END IF
       ierr = nf_put_var_real(nid, nvarid, alblw(1, nsrf))
    END DO
    !IM END albedo LW

    DO nsrf = 1, nbsrf
       IF (nsrf<=99) THEN
          WRITE (str2, '(i2.2)') nsrf
          ierr = nf_redef(nid)
          ierr = nf_def_var(nid, 'EVAP'//str2, nf_float, 1, idim2, nvarid)
          ierr = nf_put_att_text(nid, nvarid, 'title', 28, &
               'Evaporation de surface No.'//str2)
          ierr = nf_enddef(nid)
       ELSE
          PRINT *, 'Trop de sous-mailles'
          STOP 1
       END IF
       ierr = nf_put_var_real(nid, nvarid, evap(1, nsrf))
    END DO


    DO nsrf = 1, nbsrf
       IF (nsrf<=99) THEN
          WRITE (str2, '(i2.2)') nsrf
          ierr = nf_redef(nid)
          ierr = nf_def_var(nid, 'SNOW'//str2, nf_float, 1, idim2, nvarid)
          ierr = nf_put_att_text(nid, nvarid, 'title', 22, &
               'Neige de surface No.'//str2)
          ierr = nf_enddef(nid)
       ELSE
          PRINT *, 'Trop de sous-mailles'
          STOP 1
       END IF
       ierr = nf_put_var_real(nid, nvarid, snow(1, nsrf))
    END DO


    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'RADS', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 28, &
         'Rayonnement net a la surface')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, radsol)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'solsw', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 32, &
         'Rayonnement solaire a la surface')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, solsw)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'sollw', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 27, &
         'Rayonnement IF a la surface')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, sollw)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'fder', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 14, 'Derive de flux')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, fder)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'rain_f', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 21, 'precipitation liquide')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, rain_fall)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'snow_f', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 20, 'precipitation solide')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, snow_fall)

    DO nsrf = 1, nbsrf
       IF (nsrf<=99) THEN
          WRITE (str2, '(i2.2)') nsrf
          ierr = nf_redef(nid)
          ierr = nf_def_var(nid, 'RUG'//str2, nf_float, 1, idim2, nvarid)
          ierr = nf_put_att_text(nid, nvarid, 'title', 23, &
               'rugosite de surface No.'//str2)
          ierr = nf_enddef(nid)
       ELSE
          PRINT *, 'Trop de sous-mailles'
          STOP 1
       END IF
       ierr = nf_put_var_real(nid, nvarid, frugs(1, nsrf))
    END DO

    DO nsrf = 1, nbsrf
       IF (nsrf<=99) THEN
          WRITE (str2, '(i2.2)') nsrf
          ierr = nf_redef(nid)
          ierr = nf_def_var(nid, 'AGESNO'//str2, nf_float, 1, idim2, nvarid)
          ierr = nf_put_att_text(nid, nvarid, 'title', 15, &
               'Age de la neige surface No.'//str2)
          ierr = nf_enddef(nid)
       ELSE
          PRINT *, 'Trop de sous-mailles'
          STOP 1
       END IF
       ierr = nf_put_var_real(nid, nvarid, agesno(1, nsrf))
    END DO

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'ZMEA', nf_float, 1, idim2, nvarid)
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, zmea)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'ZSTD', nf_float, 1, idim2, nvarid)
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, zstd)
    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'ZSIG', nf_float, 1, idim2, nvarid)
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, zsig)
    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'ZGAM', nf_float, 1, idim2, nvarid)
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, zgam)
    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'ZTHE', nf_float, 1, idim2, nvarid)
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, zthe)
    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'ZPIC', nf_float, 1, idim2, nvarid)
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, zpic)
    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'ZVAL', nf_float, 1, idim2, nvarid)
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, zval)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'TANCIEN', nf_float, 1, idim3, nvarid)
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, t_ancien)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'QANCIEN', nf_float, 1, idim3, nvarid)
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, q_ancien)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'RUGMER', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 28, &
         'Longueur de rugosite sur mer')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, frugs(1, is_oce))

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'CLWCON', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 28, 'Eau liquide convective')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, clwcon)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'RNEBCON', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 28, 'Nebulosite convective')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, rnebcon)

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'RATQS', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 5, 'Ratqs')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, ratqs)

    ! run_off_lic_0

    ierr = nf_redef(nid)
    ierr = nf_def_var(nid, 'RUNOFFLIC0', nf_float, 1, idim2, nvarid)
    ierr = nf_put_att_text(nid, nvarid, 'title', 10, 'Runofflic0')
    ierr = nf_enddef(nid)
    ierr = nf_put_var_real(nid, nvarid, run_off_lic_0)


    ierr = nf_close(nid)

  END SUBROUTINE phyredem

end module phyredem_m
1	guez	15	module phyredem_m
2	guez	3
3	guez	15	IMPLICIT NONE
4	guez	3
5	guez	15	contains
6	guez	3
7	guez	15	SUBROUTINE phyredem(fichnom, rlat, rlon, pctsrf, tsol, tsoil, tslab, &
8			seaice, qsurf, qsol, snow, albedo, alblw, evap, rain_fall,&
9			snow_fall, solsw, sollw, fder, radsol, frugs, agesno, zmea,&
10			zstd, zsig, zgam, zthe, zpic, zval, t_ancien, q_ancien,&
11			rnebcon, ratqs, clwcon, run_off_lic_0)
12	guez	3
13	guez	15	! From phylmd/phyredem.F, v 1.3 2005/05/25 13:10:09
14			! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
15			! Objet: Ecriture de l'etat de démarrage ou redémarrage pour la physique
16	guez	12
17	guez	15	USE indicesol, ONLY : is_lic, is_oce, is_sic, is_ter, nbsrf
18			USE dimphy, ONLY : klev, klon, zmasq
19			USE dimsoil, ONLY : nsoilmx
20			USE temps, ONLY : itau_phy
21			USE netcdf95, ONLY : nf95_create, nf95_put_att
22			USE netcdf, ONLY : nf90_clobber, nf90_global
23	guez	12
24	guez	15	INCLUDE 'netcdf.inc'
25	guez	12
26	guez	15	CHARACTER(len=*) fichnom
27			REAL, INTENT (IN) :: rlat(klon), rlon(klon)
28			REAL :: tsol(klon, nbsrf)
29			REAL :: tsoil(klon, nsoilmx, nbsrf)
30	guez	12
31	guez	15	REAL :: tslab(klon), seaice(klon) !IM "slab" ocean
32			REAL :: qsurf(klon, nbsrf)
33			REAL :: qsol(klon)
34			REAL :: snow(klon, nbsrf)
35			REAL :: albedo(klon, nbsrf)
36	guez	12
37	guez	15	REAL :: alblw(klon, nbsrf)
38	guez	12
39	guez	15	REAL :: evap(klon, nbsrf)
40			REAL :: rain_fall(klon)
41			REAL :: snow_fall(klon)
42			REAL :: solsw(klon)
43			REAL :: sollw(klon)
44			REAL :: fder(klon)
45			REAL :: radsol(klon)
46			REAL :: frugs(klon, nbsrf)
47			REAL :: agesno(klon, nbsrf)
48			REAL :: zmea(klon)
49			REAL, intent(in):: zstd(klon)
50			REAL, intent(in):: zsig(klon)
51			REAL :: zgam(klon)
52			REAL :: zthe(klon)
53			REAL :: zpic(klon)
54			REAL :: zval(klon)
55			REAL :: pctsrf(klon, nbsrf)
56			REAL :: t_ancien(klon, klev), q_ancien(klon, klev)
57			REAL :: clwcon(klon, klev), rnebcon(klon, klev), ratqs(klon, klev)
58			REAL :: run_off_lic_0(klon)
59	guez	12
60	guez	15	INTEGER :: nid, nvarid, idim2, idim3
61			INTEGER :: ierr
62	guez	12
63	guez	15	INTEGER :: isoil, nsrf
64			CHARACTER (7) :: str7
65			CHARACTER (2) :: str2
66	guez	12
67	guez	15	!------------------------------------------------------------
68	guez	12
69	guez	15	PRINT *, 'Call sequence information: phyredem'
70			CALL nf95_create(fichnom, nf90_clobber, nid)
71	guez	12
72	guez	19	call nf95_put_att(nid, nf90_global, 'title', &
73			'Fichier redémarrage physique')
74	guez	15	call nf95_put_att(nid, nf90_global, "itau_phy", itau_phy)
75	guez	12
76	guez	15	ierr = nf_def_dim(nid, 'points_physiques', klon, idim2)
77			ierr = nf_def_dim(nid, 'horizon_vertical', klon*klev, idim3)
78	guez	12
79	guez	15	ierr = nf_def_var(nid, 'longitude', nf_float, 1, idim2, nvarid)
80			ierr = nf_put_att_text(nid, nvarid, 'title', 32, &
81			'Longitudes de la grille physique')
82			ierr = nf_enddef(nid)
83			ierr = nf_put_var_real(nid, nvarid, rlon)
84	guez	12
85	guez	15	ierr = nf_redef(nid)
86			ierr = nf_def_var(nid, 'latitude', nf_float, 1, idim2, nvarid)
87			ierr = nf_put_att_text(nid, nvarid, 'title', 31, &
88			'Latitudes de la grille physique')
89			ierr = nf_enddef(nid)
90			ierr = nf_put_var_real(nid, nvarid, rlat)
91	guez	12
92	guez	15	! PB ajout du masque terre/mer
93	guez	12
94	guez	15	ierr = nf_redef(nid)
95			ierr = nf_def_var(nid, 'masque', nf_float, 1, idim2, nvarid)
96			ierr = nf_put_att_text(nid, nvarid, 'title', 16, 'masque terre mer')
97			ierr = nf_enddef(nid)
98			ierr = nf_put_var_real(nid, nvarid, zmasq)
99			! BP ajout des fraction de chaque sous-surface
100	guez	12
101	guez	15	! 1. fraction de terre
102	guez	12
103	guez	15	ierr = nf_redef(nid)
104			ierr = nf_def_var(nid, 'FTER', nf_float, 1, idim2, nvarid)
105			ierr = nf_put_att_text(nid, nvarid, 'title', 21, 'fraction de continent')
106			ierr = nf_enddef(nid)
107			ierr = nf_put_var_real(nid, nvarid, pctsrf(1:klon, is_ter))
108	guez	12
109	guez	15	! 2. Fraction de glace de terre
110	guez	12
111	guez	15	ierr = nf_redef(nid)
112			ierr = nf_def_var(nid, 'FLIC', nf_float, 1, idim2, nvarid)
113			ierr = nf_put_att_text(nid, nvarid, 'title', 24, 'fraction glace de terre')
114			ierr = nf_enddef(nid)
115			ierr = nf_put_var_real(nid, nvarid, pctsrf(1:klon, is_lic))
116	guez	12
117	guez	15	! 3. fraction ocean
118	guez	12
119	guez	15	ierr = nf_redef(nid)
120			ierr = nf_def_var(nid, 'FOCE', nf_float, 1, idim2, nvarid)
121			ierr = nf_put_att_text(nid, nvarid, 'title', 14, 'fraction ocean')
122			ierr = nf_enddef(nid)
123			ierr = nf_put_var_real(nid, nvarid, pctsrf(1:klon, is_oce))
124	guez	12
125	guez	15	! 4. Fraction glace de mer
126	guez	3
127	guez	15	ierr = nf_redef(nid)
128			ierr = nf_def_var(nid, 'FSIC', nf_float, 1, idim2, nvarid)
129			ierr = nf_put_att_text(nid, nvarid, 'title', 18, 'fraction glace mer')
130			ierr = nf_enddef(nid)
131			ierr = nf_put_var_real(nid, nvarid, pctsrf(1:klon, is_sic))
132	guez	12
133
134
135	guez	15	DO nsrf = 1, nbsrf
136			IF (nsrf<=99) THEN
137			WRITE (str2, '(i2.2)') nsrf
138			ierr = nf_redef(nid)
139			ierr = nf_def_var(nid, 'TS'//str2, nf_float, 1, idim2, nvarid)
140			ierr = nf_put_att_text(nid, nvarid, 'title', 28, &
141			'Temperature de surface No.'//str2)
142			ierr = nf_enddef(nid)
143			ELSE
144			PRINT *, 'Trop de sous-mailles'
145			STOP 1
146			END IF
147			ierr = nf_put_var_real(nid, nvarid, tsol(1, nsrf))
148			END DO
149	guez	12
150	guez	15	DO nsrf = 1, nbsrf
151			DO isoil = 1, nsoilmx
152			IF (isoil<=99 .AND. nsrf<=99) THEN
153			WRITE (str7, '(i2.2, "srf", i2.2)') isoil, nsrf
154			ierr = nf_redef(nid)
155			ierr = nf_def_var(nid, 'Tsoil'//str7, nf_float, 1, idim2, nvarid)
156			ierr = nf_put_att_text(nid, nvarid, 'title', 29, &
157			'Temperature du sol No.'//str7)
158			ierr = nf_enddef(nid)
159			ELSE
160			PRINT *, 'Trop de couches'
161			STOP 1
162			END IF
163			ierr = nf_put_var_real(nid, nvarid, tsoil(1, isoil, nsrf))
164			END DO
165			END DO
166	guez	12
167	guez	15	!IM "slab" ocean
168			ierr = nf_redef(nid)
169			ierr = nf_def_var(nid, 'TSLAB', nf_float, 1, idim2, nvarid)
170			ierr = nf_put_att_text(nid, nvarid, 'title', 33, &
171			'Ecart de la SST (pour slab-ocean)')
172			ierr = nf_enddef(nid)
173			ierr = nf_put_var_real(nid, nvarid, tslab)
174	guez	12
175	guez	15	ierr = nf_redef(nid)
176			ierr = nf_def_var(nid, 'SEAICE', nf_float, 1, idim2, nvarid)
177			ierr = nf_put_att_text(nid, nvarid, 'title', 33, &
178			'Glace de mer kg/m2 (pour slab-ocean)')
179			ierr = nf_enddef(nid)
180			ierr = nf_put_var_real(nid, nvarid, seaice)
181	guez	12
182	guez	15	DO nsrf = 1, nbsrf
183			IF (nsrf<=99) THEN
184			WRITE (str2, '(i2.2)') nsrf
185			ierr = nf_redef(nid)
186			ierr = nf_def_var(nid, 'QS'//str2, nf_float, 1, idim2, nvarid)
187			ierr = nf_put_att_text(nid, nvarid, 'title', 25, &
188			'Humidite de surface No.'//str2)
189			ierr = nf_enddef(nid)
190			ELSE
191			PRINT *, 'Trop de sous-mailles'
192			STOP 1
193			END IF
194			ierr = nf_put_var_real(nid, nvarid, qsurf(1, nsrf))
195			END DO
196	guez	12
197	guez	15	ierr = nf_redef(nid)
198			ierr = nf_def_var(nid, 'QSOL', nf_float, 1, idim2, nvarid)
199			ierr = nf_put_att_text(nid, nvarid, 'title', 20, 'Eau dans le sol (mm)')
200			ierr = nf_enddef(nid)
201			ierr = nf_put_var_real(nid, nvarid, qsol)
202	guez	12
203	guez	15	DO nsrf = 1, nbsrf
204			IF (nsrf<=99) THEN
205			WRITE (str2, '(i2.2)') nsrf
206			ierr = nf_redef(nid)
207			ierr = nf_def_var(nid, 'ALBE'//str2, nf_float, 1, idim2, nvarid)
208			ierr = nf_put_att_text(nid, nvarid, 'title', 23, &
209			'albedo de surface No.'//str2)
210			ierr = nf_enddef(nid)
211			ELSE
212			PRINT *, 'Trop de sous-mailles'
213			STOP 1
214			END IF
215			ierr = nf_put_var_real(nid, nvarid, albedo(1, nsrf))
216			END DO
217	guez	12
218	guez	15	!IM BEG albedo LW
219			DO nsrf = 1, nbsrf
220			IF (nsrf<=99) THEN
221			WRITE (str2, '(i2.2)') nsrf
222			ierr = nf_redef(nid)
223			ierr = nf_def_var(nid, 'ALBLW'//str2, nf_float, 1, idim2, nvarid)
224			ierr = nf_put_att_text(nid, nvarid, 'title', 23, &
225			'albedo LW de surface No.'//str2)
226			ierr = nf_enddef(nid)
227			ELSE
228			PRINT *, 'Trop de sous-mailles'
229			STOP 1
230			END IF
231			ierr = nf_put_var_real(nid, nvarid, alblw(1, nsrf))
232			END DO
233			!IM END albedo LW
234	guez	12
235	guez	15	DO nsrf = 1, nbsrf
236			IF (nsrf<=99) THEN
237			WRITE (str2, '(i2.2)') nsrf
238			ierr = nf_redef(nid)
239			ierr = nf_def_var(nid, 'EVAP'//str2, nf_float, 1, idim2, nvarid)
240			ierr = nf_put_att_text(nid, nvarid, 'title', 28, &
241			'Evaporation de surface No.'//str2)
242			ierr = nf_enddef(nid)
243			ELSE
244			PRINT *, 'Trop de sous-mailles'
245			STOP 1
246			END IF
247			ierr = nf_put_var_real(nid, nvarid, evap(1, nsrf))
248			END DO
249	guez	12
250
251	guez	15	DO nsrf = 1, nbsrf
252			IF (nsrf<=99) THEN
253			WRITE (str2, '(i2.2)') nsrf
254			ierr = nf_redef(nid)
255			ierr = nf_def_var(nid, 'SNOW'//str2, nf_float, 1, idim2, nvarid)
256			ierr = nf_put_att_text(nid, nvarid, 'title', 22, &
257			'Neige de surface No.'//str2)
258			ierr = nf_enddef(nid)
259			ELSE
260			PRINT *, 'Trop de sous-mailles'
261			STOP 1
262			END IF
263			ierr = nf_put_var_real(nid, nvarid, snow(1, nsrf))
264			END DO
265	guez	12
266
267	guez	15	ierr = nf_redef(nid)
268			ierr = nf_def_var(nid, 'RADS', nf_float, 1, idim2, nvarid)
269			ierr = nf_put_att_text(nid, nvarid, 'title', 28, &
270			'Rayonnement net a la surface')
271			ierr = nf_enddef(nid)
272			ierr = nf_put_var_real(nid, nvarid, radsol)
273	guez	12
274	guez	15	ierr = nf_redef(nid)
275			ierr = nf_def_var(nid, 'solsw', nf_float, 1, idim2, nvarid)
276			ierr = nf_put_att_text(nid, nvarid, 'title', 32, &
277			'Rayonnement solaire a la surface')
278			ierr = nf_enddef(nid)
279			ierr = nf_put_var_real(nid, nvarid, solsw)
280	guez	12
281	guez	15	ierr = nf_redef(nid)
282			ierr = nf_def_var(nid, 'sollw', nf_float, 1, idim2, nvarid)
283			ierr = nf_put_att_text(nid, nvarid, 'title', 27, &
284			'Rayonnement IF a la surface')
285			ierr = nf_enddef(nid)
286			ierr = nf_put_var_real(nid, nvarid, sollw)
287	guez	12
288	guez	15	ierr = nf_redef(nid)
289			ierr = nf_def_var(nid, 'fder', nf_float, 1, idim2, nvarid)
290			ierr = nf_put_att_text(nid, nvarid, 'title', 14, 'Derive de flux')
291			ierr = nf_enddef(nid)
292			ierr = nf_put_var_real(nid, nvarid, fder)
293	guez	12
294	guez	15	ierr = nf_redef(nid)
295			ierr = nf_def_var(nid, 'rain_f', nf_float, 1, idim2, nvarid)
296			ierr = nf_put_att_text(nid, nvarid, 'title', 21, 'precipitation liquide')
297			ierr = nf_enddef(nid)
298			ierr = nf_put_var_real(nid, nvarid, rain_fall)
299	guez	12
300	guez	15	ierr = nf_redef(nid)
301			ierr = nf_def_var(nid, 'snow_f', nf_float, 1, idim2, nvarid)
302			ierr = nf_put_att_text(nid, nvarid, 'title', 20, 'precipitation solide')
303			ierr = nf_enddef(nid)
304			ierr = nf_put_var_real(nid, nvarid, snow_fall)
305	guez	12
306	guez	15	DO nsrf = 1, nbsrf
307			IF (nsrf<=99) THEN
308			WRITE (str2, '(i2.2)') nsrf
309			ierr = nf_redef(nid)
310			ierr = nf_def_var(nid, 'RUG'//str2, nf_float, 1, idim2, nvarid)
311			ierr = nf_put_att_text(nid, nvarid, 'title', 23, &
312			'rugosite de surface No.'//str2)
313			ierr = nf_enddef(nid)
314			ELSE
315			PRINT *, 'Trop de sous-mailles'
316			STOP 1
317			END IF
318			ierr = nf_put_var_real(nid, nvarid, frugs(1, nsrf))
319			END DO
320	guez	12
321	guez	15	DO nsrf = 1, nbsrf
322			IF (nsrf<=99) THEN
323			WRITE (str2, '(i2.2)') nsrf
324			ierr = nf_redef(nid)
325			ierr = nf_def_var(nid, 'AGESNO'//str2, nf_float, 1, idim2, nvarid)
326			ierr = nf_put_att_text(nid, nvarid, 'title', 15, &
327			'Age de la neige surface No.'//str2)
328			ierr = nf_enddef(nid)
329			ELSE
330			PRINT *, 'Trop de sous-mailles'
331			STOP 1
332			END IF
333			ierr = nf_put_var_real(nid, nvarid, agesno(1, nsrf))
334			END DO
335	guez	12
336	guez	15	ierr = nf_redef(nid)
337			ierr = nf_def_var(nid, 'ZMEA', nf_float, 1, idim2, nvarid)
338			ierr = nf_enddef(nid)
339			ierr = nf_put_var_real(nid, nvarid, zmea)
340	guez	12
341	guez	15	ierr = nf_redef(nid)
342			ierr = nf_def_var(nid, 'ZSTD', nf_float, 1, idim2, nvarid)
343			ierr = nf_enddef(nid)
344			ierr = nf_put_var_real(nid, nvarid, zstd)
345			ierr = nf_redef(nid)
346			ierr = nf_def_var(nid, 'ZSIG', nf_float, 1, idim2, nvarid)
347			ierr = nf_enddef(nid)
348			ierr = nf_put_var_real(nid, nvarid, zsig)
349			ierr = nf_redef(nid)
350			ierr = nf_def_var(nid, 'ZGAM', nf_float, 1, idim2, nvarid)
351			ierr = nf_enddef(nid)
352			ierr = nf_put_var_real(nid, nvarid, zgam)
353			ierr = nf_redef(nid)
354			ierr = nf_def_var(nid, 'ZTHE', nf_float, 1, idim2, nvarid)
355			ierr = nf_enddef(nid)
356			ierr = nf_put_var_real(nid, nvarid, zthe)
357			ierr = nf_redef(nid)
358			ierr = nf_def_var(nid, 'ZPIC', nf_float, 1, idim2, nvarid)
359			ierr = nf_enddef(nid)
360			ierr = nf_put_var_real(nid, nvarid, zpic)
361			ierr = nf_redef(nid)
362			ierr = nf_def_var(nid, 'ZVAL', nf_float, 1, idim2, nvarid)
363			ierr = nf_enddef(nid)
364			ierr = nf_put_var_real(nid, nvarid, zval)
365	guez	12
366	guez	15	ierr = nf_redef(nid)
367			ierr = nf_def_var(nid, 'TANCIEN', nf_float, 1, idim3, nvarid)
368			ierr = nf_enddef(nid)
369			ierr = nf_put_var_real(nid, nvarid, t_ancien)
370	guez	12
371	guez	15	ierr = nf_redef(nid)
372			ierr = nf_def_var(nid, 'QANCIEN', nf_float, 1, idim3, nvarid)
373			ierr = nf_enddef(nid)
374			ierr = nf_put_var_real(nid, nvarid, q_ancien)
375	guez	12
376	guez	15	ierr = nf_redef(nid)
377			ierr = nf_def_var(nid, 'RUGMER', nf_float, 1, idim2, nvarid)
378			ierr = nf_put_att_text(nid, nvarid, 'title', 28, &
379			'Longueur de rugosite sur mer')
380			ierr = nf_enddef(nid)
381			ierr = nf_put_var_real(nid, nvarid, frugs(1, is_oce))
382	guez	12
383	guez	15	ierr = nf_redef(nid)
384			ierr = nf_def_var(nid, 'CLWCON', nf_float, 1, idim2, nvarid)
385			ierr = nf_put_att_text(nid, nvarid, 'title', 28, 'Eau liquide convective')
386			ierr = nf_enddef(nid)
387			ierr = nf_put_var_real(nid, nvarid, clwcon)
388	guez	12
389	guez	15	ierr = nf_redef(nid)
390			ierr = nf_def_var(nid, 'RNEBCON', nf_float, 1, idim2, nvarid)
391			ierr = nf_put_att_text(nid, nvarid, 'title', 28, 'Nebulosite convective')
392			ierr = nf_enddef(nid)
393			ierr = nf_put_var_real(nid, nvarid, rnebcon)
394	guez	12
395	guez	15	ierr = nf_redef(nid)
396			ierr = nf_def_var(nid, 'RATQS', nf_float, 1, idim2, nvarid)
397			ierr = nf_put_att_text(nid, nvarid, 'title', 5, 'Ratqs')
398			ierr = nf_enddef(nid)
399			ierr = nf_put_var_real(nid, nvarid, ratqs)
400
401			! run_off_lic_0
402
403			ierr = nf_redef(nid)
404			ierr = nf_def_var(nid, 'RUNOFFLIC0', nf_float, 1, idim2, nvarid)
405			ierr = nf_put_att_text(nid, nvarid, 'title', 10, 'Runofflic0')
406			ierr = nf_enddef(nid)
407			ierr = nf_put_var_real(nid, nvarid, run_off_lic_0)
408
409
410			ierr = nf_close(nid)
411
412			END SUBROUTINE phyredem
413
414			end module phyredem_m