Sources/phylmd/phyetat0.f

module phyetat0_m

  use dimphy, only: klon

  IMPLICIT none

  REAL, save:: rlat(klon), rlon(klon)
  ! latitude and longitude of a point of the scalar grid identified
  ! by a simple index, in degrees

  private klon

contains

  SUBROUTINE phyetat0(pctsrf, tsol, tsoil, tslab, seaice, qsurf, qsol, &
       snow, albe, evap, rain_fall, snow_fall, solsw, sollw, fder, &
       radsol, frugs, agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, &
       t_ancien, q_ancien, ancien_ok, rnebcon, ratqs, clwcon, run_off_lic_0, &
       sig1, w01, ncid_startphy)

    ! From phylmd/phyetat0.F, version 1.4 2005/06/03 10:03:07
    ! Author: Z.X. Li (LMD/CNRS)
    ! Date: 1993/08/18
    ! Objet : lecture de l'état initial pour la physique

    use dimphy, only: zmasq, klev
    USE dimsoil, ONLY : nsoilmx
    USE indicesol, ONLY : epsfra, is_lic, is_oce, is_sic, is_ter, nbsrf
    use netcdf, only: nf90_global, nf90_inq_varid, NF90_NOERR, &
         NF90_NOWRITE
    use netcdf95, only: nf95_close, nf95_get_att, nf95_get_var, &
         nf95_inq_varid, nf95_inquire_variable, NF95_OPEN
    USE temps, ONLY : itau_phy

    REAL pctsrf(klon, nbsrf)
    REAL tsol(klon, nbsrf)
    REAL tsoil(klon, nsoilmx, nbsrf)
    REAL tslab(klon), seaice(klon)
    REAL qsurf(klon, nbsrf)
    REAL, intent(out):: qsol(:) ! (klon)
    REAL snow(klon, nbsrf)
    REAL albe(klon, nbsrf)
    REAL evap(klon, nbsrf)
    REAL, intent(out):: rain_fall(klon)
    REAL snow_fall(klon)
    real solsw(klon)
    REAL, intent(out):: sollw(klon)
    real fder(klon)
    REAL radsol(klon)
    REAL frugs(klon, nbsrf)
    REAL agesno(klon, nbsrf)
    REAL, intent(out):: zmea(klon)
    REAL, intent(out):: zstd(klon)
    REAL, intent(out):: zsig(klon)
    REAL zgam(klon)
    REAL zthe(klon)
    REAL zpic(klon)
    REAL zval(klon)
    REAL t_ancien(klon, klev), q_ancien(klon, klev)
    LOGICAL, intent(out):: ancien_ok
    real rnebcon(klon, klev), ratqs(klon, klev), clwcon(klon, klev)
    REAL run_off_lic_0(klon)
    real, intent(out):: sig1(klon, klev) ! section adiabatic updraft

    real, intent(out):: w01(klon, klev) 
    ! vertical velocity within adiabatic updraft

    integer, intent(out):: ncid_startphy

    ! Local:
    REAL fractint(klon)
    INTEGER varid, ndims
    INTEGER ierr, i

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

    print *, "Call sequence information: phyetat0"

    ! Fichier contenant l'état initial :
    call NF95_OPEN("startphy.nc", NF90_NOWRITE, ncid_startphy)

    call nf95_get_att(ncid_startphy, nf90_global, "itau_phy", itau_phy)

    ! Lecture des latitudes (coordonnees):

    call NF95_INQ_VARID(ncid_startphy, "latitude", varid)
    call NF95_GET_VAR(ncid_startphy, varid, rlat)

    ! Lecture des longitudes (coordonnees):

    call NF95_INQ_VARID(ncid_startphy, "longitude", varid)
    call NF95_GET_VAR(ncid_startphy, varid, rlon)

    ! Lecture du masque terre mer

    call NF95_INQ_VARID(ncid_startphy, "masque", varid)
    call nf95_get_var(ncid_startphy, varid, zmasq)

    ! Lecture des fractions pour chaque sous-surface

    ! initialisation des sous-surfaces

    pctsrf = 0.

    ! fraction de terre

    ierr = NF90_INQ_VARID(ncid_startphy, "FTER", varid)
    IF (ierr == NF90_NOERR) THEN
       call nf95_get_var(ncid_startphy, varid, pctsrf(:, is_ter))
    else
       PRINT *, 'phyetat0: Le champ <FTER> est absent'
    ENDIF

    ! fraction de glace de terre

    ierr = NF90_INQ_VARID(ncid_startphy, "FLIC", varid)
    IF (ierr == NF90_NOERR) THEN
       call nf95_get_var(ncid_startphy, varid, pctsrf(:, is_lic))
    else
       PRINT *, 'phyetat0: Le champ <FLIC> est absent'
    ENDIF

    ! fraction d'ocean

    ierr = NF90_INQ_VARID(ncid_startphy, "FOCE", varid)
    IF (ierr == NF90_NOERR) THEN
       call nf95_get_var(ncid_startphy, varid, pctsrf(:, is_oce))
    else
       PRINT *, 'phyetat0: Le champ <FOCE> est absent'
    ENDIF

    ! fraction glace de mer

    ierr = NF90_INQ_VARID(ncid_startphy, "FSIC", varid)
    IF (ierr == NF90_NOERR) THEN
       call nf95_get_var(ncid_startphy, varid, pctsrf(:, is_sic))
    else
       PRINT *, 'phyetat0: Le champ <FSIC> est absent'
    ENDIF

    ! Verification de l'adequation entre le masque et les sous-surfaces

    fractint = pctsrf(:, is_ter) + pctsrf(:, is_lic)
    DO i = 1 , klon
       IF ( abs(fractint(i) - zmasq(i) ) > EPSFRA ) THEN
          WRITE(*, *) 'phyetat0: attention fraction terre pas ', &
               'coherente ', i, zmasq(i), pctsrf(i, is_ter) &
               , pctsrf(i, is_lic)
       ENDIF
    END DO
    fractint = pctsrf(:, is_oce) + pctsrf(:, is_sic)
    DO i = 1 , klon
       IF ( abs( fractint(i) - (1. - zmasq(i))) > EPSFRA ) THEN
          WRITE(*, *) 'phyetat0 attention fraction ocean pas ', &
               'coherente ', i, zmasq(i) , pctsrf(i, is_oce) &
               , pctsrf(i, is_sic)
       ENDIF
    END DO

    ! Lecture des temperatures du sol:
    call NF95_INQ_VARID(ncid_startphy, "TS", varid)
    call nf95_inquire_variable(ncid_startphy, varid, ndims = ndims)
    if (ndims == 2) then
       call NF95_GET_VAR(ncid_startphy, varid, tsol)
    else
       print *, "Found only one surface type for soil temperature."
       call nf95_get_var(ncid_startphy, varid, tsol(:, 1))
       tsol(:, 2:nbsrf) = spread(tsol(:, 1), dim = 2, ncopies = nbsrf - 1)
    end if

    ! Lecture des temperatures du sol profond:

    call NF95_INQ_VARID(ncid_startphy, 'Tsoil', varid)
    call NF95_GET_VAR(ncid_startphy, varid, tsoil)

    !IM "slab" ocean 
    ! Lecture de tslab (pour slab ocean seulement): 
    tslab = 0.
    seaice = 0.

    ! Lecture de l'humidite de l'air juste au dessus du sol:

    call NF95_INQ_VARID(ncid_startphy, "QS", varid)
    call nf95_get_var(ncid_startphy, varid, qsurf)

    ! Eau dans le sol (pour le modele de sol "bucket")

    ierr = NF90_INQ_VARID(ncid_startphy, "QSOL", varid)
    IF (ierr == NF90_NOERR) THEN
       call nf95_get_var(ncid_startphy, varid, qsol)
    else
       PRINT *, 'phyetat0: Le champ <QSOL> est absent'
       PRINT *, ' Valeur par defaut nulle'
       qsol = 0.
    ENDIF

    ! Lecture de neige au sol:

    call NF95_INQ_VARID(ncid_startphy, "SNOW", varid)
    call nf95_get_var(ncid_startphy, varid, snow)

    ! Lecture de albedo au sol:

    call NF95_INQ_VARID(ncid_startphy, "ALBE", varid)
    call nf95_get_var(ncid_startphy, varid, albe)

    ! Lecture de evaporation: 

    call NF95_INQ_VARID(ncid_startphy, "EVAP", varid)
    call nf95_get_var(ncid_startphy, varid, evap)

    ! Lecture precipitation liquide:

    call NF95_INQ_VARID(ncid_startphy, "rain_f", varid)
    call NF95_GET_VAR(ncid_startphy, varid, rain_fall)

    ! Lecture precipitation solide:

    call NF95_INQ_VARID(ncid_startphy, "snow_f", varid)
    call NF95_GET_VAR(ncid_startphy, varid, snow_fall)

    ! Lecture rayonnement solaire au sol:

    ierr = NF90_INQ_VARID(ncid_startphy, "solsw", varid)
    IF (ierr /= NF90_NOERR) THEN
       PRINT *, 'phyetat0: Le champ <solsw> est absent'
       PRINT *, 'mis a zero'
       solsw = 0.
    ELSE
       call nf95_get_var(ncid_startphy, varid, solsw)
    ENDIF

    ! Lecture rayonnement IF au sol:

    ierr = NF90_INQ_VARID(ncid_startphy, "sollw", varid)
    IF (ierr /= NF90_NOERR) THEN
       PRINT *, 'phyetat0: Le champ <sollw> est absent'
       PRINT *, 'mis a zero'
       sollw = 0.
    ELSE
       call nf95_get_var(ncid_startphy, varid, sollw)
    ENDIF

    ! Lecture derive des flux:

    ierr = NF90_INQ_VARID(ncid_startphy, "fder", varid)
    IF (ierr /= NF90_NOERR) THEN
       PRINT *, 'phyetat0: Le champ <fder> est absent'
       PRINT *, 'mis a zero'
       fder = 0.
    ELSE
       call nf95_get_var(ncid_startphy, varid, fder)
    ENDIF

    ! Lecture du rayonnement net au sol:

    call NF95_INQ_VARID(ncid_startphy, "RADS", varid)
    call NF95_GET_VAR(ncid_startphy, varid, radsol)

    ! Lecture de la longueur de rugosite 

    call NF95_INQ_VARID(ncid_startphy, "RUG", varid)
    call nf95_get_var(ncid_startphy, varid, frugs)

    ! Lecture de l'age de la neige:

    call NF95_INQ_VARID(ncid_startphy, "AGESNO", varid)
    call nf95_get_var(ncid_startphy, varid, agesno)

    call NF95_INQ_VARID(ncid_startphy, "ZMEA", varid)
    call NF95_GET_VAR(ncid_startphy, varid, zmea)

    call NF95_INQ_VARID(ncid_startphy, "ZSTD", varid)
    call NF95_GET_VAR(ncid_startphy, varid, zstd)

    call NF95_INQ_VARID(ncid_startphy, "ZSIG", varid)
    call NF95_GET_VAR(ncid_startphy, varid, zsig)

    call NF95_INQ_VARID(ncid_startphy, "ZGAM", varid)
    call NF95_GET_VAR(ncid_startphy, varid, zgam)

    call NF95_INQ_VARID(ncid_startphy, "ZTHE", varid)
    call NF95_GET_VAR(ncid_startphy, varid, zthe)

    call NF95_INQ_VARID(ncid_startphy, "ZPIC", varid)
    call NF95_GET_VAR(ncid_startphy, varid, zpic)

    call NF95_INQ_VARID(ncid_startphy, "ZVAL", varid)
    call NF95_GET_VAR(ncid_startphy, varid, zval)

    ancien_ok = .TRUE.

    ierr = NF90_INQ_VARID(ncid_startphy, "TANCIEN", varid)
    IF (ierr /= NF90_NOERR) THEN
       PRINT *, "phyetat0: Le champ <TANCIEN> est absent"
       PRINT *, "Depart legerement fausse. Mais je continue"
       ancien_ok = .FALSE.
    ELSE
       call nf95_get_var(ncid_startphy, varid, t_ancien)
    ENDIF

    ierr = NF90_INQ_VARID(ncid_startphy, "QANCIEN", varid)
    IF (ierr /= NF90_NOERR) THEN
       PRINT *, "phyetat0: Le champ <QANCIEN> est absent"
       PRINT *, "Depart legerement fausse. Mais je continue"
       ancien_ok = .FALSE.
    ELSE
       call nf95_get_var(ncid_startphy, varid, q_ancien)
    ENDIF

    ierr = NF90_INQ_VARID(ncid_startphy, "CLWCON", varid)
    IF (ierr /= NF90_NOERR) THEN
       PRINT *, "phyetat0: Le champ CLWCON est absent"
       PRINT *, "Depart legerement fausse. Mais je continue"
       clwcon = 0.
    ELSE
       call nf95_get_var(ncid_startphy, varid, clwcon(:, 1))
       clwcon(:, 2:) = 0.
    ENDIF

    ierr = NF90_INQ_VARID(ncid_startphy, "RNEBCON", varid)
    IF (ierr /= NF90_NOERR) THEN
       PRINT *, "phyetat0: Le champ RNEBCON est absent"
       PRINT *, "Depart legerement fausse. Mais je continue"
       rnebcon = 0.
    ELSE
       call nf95_get_var(ncid_startphy, varid, rnebcon(:, 1))
       rnebcon(:, 2:) = 0.
    ENDIF

    ! Lecture ratqs

    ierr = NF90_INQ_VARID(ncid_startphy, "RATQS", varid)
    IF (ierr /= NF90_NOERR) THEN
       PRINT *, "phyetat0: Le champ <RATQS> est absent"
       PRINT *, "Depart legerement fausse. Mais je continue"
       ratqs = 0.
    ELSE
       call nf95_get_var(ncid_startphy, varid, ratqs(:, 1))
       ratqs(:, 2:) = 0.
    ENDIF

    ! Lecture run_off_lic_0

    ierr = NF90_INQ_VARID(ncid_startphy, "RUNOFFLIC0", varid)
    IF (ierr /= NF90_NOERR) THEN
       PRINT *, "phyetat0: Le champ <RUNOFFLIC0> est absent"
       PRINT *, "Depart legerement fausse. Mais je continue"
       run_off_lic_0 = 0.
    ELSE
       call nf95_get_var(ncid_startphy, varid, run_off_lic_0)
    ENDIF

    call nf95_inq_varid(ncid_startphy, "sig1", varid)
    call nf95_get_var(ncid_startphy, varid, sig1)

    call nf95_inq_varid(ncid_startphy, "w01", varid)
    call nf95_get_var(ncid_startphy, varid, w01)

  END SUBROUTINE phyetat0

end module phyetat0_m
1	guez	3	module phyetat0_m
2
3	guez	15	use dimphy, only: klon
4	guez	12
5	guez	3	IMPLICIT none
6
7	guez	138	REAL, save:: rlat(klon), rlon(klon)
8			! latitude and longitude of a point of the scalar grid identified
9			! by a simple index, in degrees
10	guez	3
11	guez	15	private klon
12	guez	3
13			contains
14
15	guez	99	SUBROUTINE phyetat0(pctsrf, tsol, tsoil, tslab, seaice, qsurf, qsol, &
16	guez	155	snow, albe, evap, rain_fall, snow_fall, solsw, sollw, fder, &
17	guez	99	radsol, frugs, agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, &
18			t_ancien, q_ancien, ancien_ok, rnebcon, ratqs, clwcon, run_off_lic_0, &
19	guez	157	sig1, w01, ncid_startphy)
20	guez	3
21	guez	49	! From phylmd/phyetat0.F, version 1.4 2005/06/03 10:03:07
22			! Author: Z.X. Li (LMD/CNRS)
23	guez	50	! Date: 1993/08/18
24	guez	69	! Objet : lecture de l'état initial pour la physique
25	guez	3
26	guez	69	use dimphy, only: zmasq, klev
27			USE dimsoil, ONLY : nsoilmx
28	guez	12	USE indicesol, ONLY : epsfra, is_lic, is_oce, is_sic, is_ter, nbsrf
29	guez	101	use netcdf, only: nf90_global, nf90_inq_varid, NF90_NOERR, &
30	guez	50	NF90_NOWRITE
31	guez	101	use netcdf95, only: nf95_close, nf95_get_att, nf95_get_var, &
32			nf95_inq_varid, nf95_inquire_variable, NF95_OPEN
33	guez	69	USE temps, ONLY : itau_phy
34	guez	12
35	guez	72	REAL pctsrf(klon, nbsrf)
36	guez	49	REAL tsol(klon, nbsrf)
37			REAL tsoil(klon, nsoilmx, nbsrf)
38	guez	3	REAL tslab(klon), seaice(klon)
39	guez	49	REAL qsurf(klon, nbsrf)
40	guez	99	REAL, intent(out):: qsol(:) ! (klon)
41	guez	49	REAL snow(klon, nbsrf)
42			REAL albe(klon, nbsrf)
43			REAL evap(klon, nbsrf)
44	guez	62	REAL, intent(out):: rain_fall(klon)
45	guez	3	REAL snow_fall(klon)
46			real solsw(klon)
47	guez	72	REAL, intent(out):: sollw(klon)
48	guez	3	real fder(klon)
49	guez	72	REAL radsol(klon)
50	guez	49	REAL frugs(klon, nbsrf)
51			REAL agesno(klon, nbsrf)
52	guez	174	REAL, intent(out):: zmea(klon)
53	guez	13	REAL, intent(out):: zstd(klon)
54			REAL, intent(out):: zsig(klon)
55	guez	3	REAL zgam(klon)
56			REAL zthe(klon)
57			REAL zpic(klon)
58			REAL zval(klon)
59	guez	49	REAL t_ancien(klon, klev), q_ancien(klon, klev)
60			LOGICAL, intent(out):: ancien_ok
61	guez	72	real rnebcon(klon, klev), ratqs(klon, klev), clwcon(klon, klev)
62			REAL run_off_lic_0(klon)
63			real, intent(out):: sig1(klon, klev) ! section adiabatic updraft
64	guez	3
65	guez	72	real, intent(out):: w01(klon, klev)
66			! vertical velocity within adiabatic updraft
67	guez	3
68	guez	157	integer, intent(out):: ncid_startphy
69
70	guez	72	! Local:
71			REAL fractint(klon)
72	guez	157	INTEGER varid, ndims
73	guez	156	INTEGER ierr, i
74	guez	3
75			!---------------------------------------------------------------
76
77			print *, "Call sequence information: phyetat0"
78
79	guez	72	! Fichier contenant l'état initial :
80	guez	157	call NF95_OPEN("startphy.nc", NF90_NOWRITE, ncid_startphy)
81	guez	3
82	guez	157	call nf95_get_att(ncid_startphy, nf90_global, "itau_phy", itau_phy)
83	guez	3
84			! Lecture des latitudes (coordonnees):
85
86	guez	157	call NF95_INQ_VARID(ncid_startphy, "latitude", varid)
87			call NF95_GET_VAR(ncid_startphy, varid, rlat)
88	guez	3
89			! Lecture des longitudes (coordonnees):
90
91	guez	157	call NF95_INQ_VARID(ncid_startphy, "longitude", varid)
92			call NF95_GET_VAR(ncid_startphy, varid, rlon)
93	guez	3
94			! Lecture du masque terre mer
95
96	guez	157	call NF95_INQ_VARID(ncid_startphy, "masque", varid)
97			call nf95_get_var(ncid_startphy, varid, zmasq)
98	guez	101
99	guez	3	! Lecture des fractions pour chaque sous-surface
100
101			! initialisation des sous-surfaces
102
103			pctsrf = 0.
104
105			! fraction de terre
106
107	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "FTER", varid)
108	guez	50	IF (ierr == NF90_NOERR) THEN
109	guez	157	call nf95_get_var(ncid_startphy, varid, pctsrf(:, is_ter))
110	guez	3	else
111	guez	43	PRINT *, 'phyetat0: Le champ <FTER> est absent'
112	guez	3	ENDIF
113
114			! fraction de glace de terre
115
116	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "FLIC", varid)
117	guez	50	IF (ierr == NF90_NOERR) THEN
118	guez	157	call nf95_get_var(ncid_startphy, varid, pctsrf(:, is_lic))
119	guez	3	else
120	guez	43	PRINT *, 'phyetat0: Le champ <FLIC> est absent'
121	guez	3	ENDIF
122
123			! fraction d'ocean
124
125	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "FOCE", varid)
126	guez	50	IF (ierr == NF90_NOERR) THEN
127	guez	157	call nf95_get_var(ncid_startphy, varid, pctsrf(:, is_oce))
128	guez	3	else
129	guez	43	PRINT *, 'phyetat0: Le champ <FOCE> est absent'
130	guez	3	ENDIF
131
132			! fraction glace de mer
133
134	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "FSIC", varid)
135	guez	50	IF (ierr == NF90_NOERR) THEN
136	guez	157	call nf95_get_var(ncid_startphy, varid, pctsrf(:, is_sic))
137	guez	3	else
138	guez	43	PRINT *, 'phyetat0: Le champ <FSIC> est absent'
139	guez	3	ENDIF
140
141	guez	50	! Verification de l'adequation entre le masque et les sous-surfaces
142	guez	3
143	guez	50	fractint = pctsrf(:, is_ter) + pctsrf(:, is_lic)
144	guez	3	DO i = 1 , klon
145	guez	50	IF ( abs(fractint(i) - zmasq(i) ) > EPSFRA ) THEN
146			WRITE(, ) 'phyetat0: attention fraction terre pas ', &
147	guez	3	'coherente ', i, zmasq(i), pctsrf(i, is_ter) &
148	guez	49	, pctsrf(i, is_lic)
149	guez	3	ENDIF
150			END DO
151	guez	50	fractint = pctsrf(:, is_oce) + pctsrf(:, is_sic)
152	guez	3	DO i = 1 , klon
153	guez	50	IF ( abs( fractint(i) - (1. - zmasq(i))) > EPSFRA ) THEN
154			WRITE(, ) 'phyetat0 attention fraction ocean pas ', &
155	guez	3	'coherente ', i, zmasq(i) , pctsrf(i, is_oce) &
156	guez	49	, pctsrf(i, is_sic)
157	guez	3	ENDIF
158			END DO
159
160			! Lecture des temperatures du sol:
161	guez	157	call NF95_INQ_VARID(ncid_startphy, "TS", varid)
162			call nf95_inquire_variable(ncid_startphy, varid, ndims = ndims)
163	guez	101	if (ndims == 2) then
164	guez	157	call NF95_GET_VAR(ncid_startphy, varid, tsol)
165	guez	101	else
166			print *, "Found only one surface type for soil temperature."
167	guez	157	call nf95_get_var(ncid_startphy, varid, tsol(:, 1))
168	guez	101	tsol(:, 2:nbsrf) = spread(tsol(:, 1), dim = 2, ncopies = nbsrf - 1)
169	guez	156	end if
170	guez	3
171	guez	156	! Lecture des temperatures du sol profond:
172	guez	3
173	guez	157	call NF95_INQ_VARID(ncid_startphy, 'Tsoil', varid)
174			call NF95_GET_VAR(ncid_startphy, varid, tsoil)
175	guez	3
176			!IM "slab" ocean
177	guez	50	! Lecture de tslab (pour slab ocean seulement):
178	guez	99	tslab = 0.
179			seaice = 0.
180	guez	3
181			! Lecture de l'humidite de l'air juste au dessus du sol:
182
183	guez	157	call NF95_INQ_VARID(ncid_startphy, "QS", varid)
184			call nf95_get_var(ncid_startphy, varid, qsurf)
185	guez	3
186			! Eau dans le sol (pour le modele de sol "bucket")
187
188	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "QSOL", varid)
189	guez	50	IF (ierr == NF90_NOERR) THEN
190	guez	157	call nf95_get_var(ncid_startphy, varid, qsol)
191	guez	3	else
192	guez	43	PRINT *, 'phyetat0: Le champ <QSOL> est absent'
193	guez	50	PRINT *, ' Valeur par defaut nulle'
194	guez	43	qsol = 0.
195	guez	3	ENDIF
196
197			! Lecture de neige au sol:
198
199	guez	157	call NF95_INQ_VARID(ncid_startphy, "SNOW", varid)
200			call nf95_get_var(ncid_startphy, varid, snow)
201	guez	3
202			! Lecture de albedo au sol:
203
204	guez	157	call NF95_INQ_VARID(ncid_startphy, "ALBE", varid)
205			call nf95_get_var(ncid_startphy, varid, albe)
206	guez	3
207	guez	50	! Lecture de evaporation:
208	guez	3
209	guez	157	call NF95_INQ_VARID(ncid_startphy, "EVAP", varid)
210			call nf95_get_var(ncid_startphy, varid, evap)
211	guez	3
212			! Lecture precipitation liquide:
213
214	guez	157	call NF95_INQ_VARID(ncid_startphy, "rain_f", varid)
215			call NF95_GET_VAR(ncid_startphy, varid, rain_fall)
216	guez	3
217			! Lecture precipitation solide:
218
219	guez	157	call NF95_INQ_VARID(ncid_startphy, "snow_f", varid)
220			call NF95_GET_VAR(ncid_startphy, varid, snow_fall)
221	guez	3
222			! Lecture rayonnement solaire au sol:
223
224	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "solsw", varid)
225	guez	49	IF (ierr /= NF90_NOERR) THEN
226	guez	43	PRINT *, 'phyetat0: Le champ <solsw> est absent'
227			PRINT *, 'mis a zero'
228	guez	3	solsw = 0.
229			ELSE
230	guez	157	call nf95_get_var(ncid_startphy, varid, solsw)
231	guez	3	ENDIF
232
233			! Lecture rayonnement IF au sol:
234
235	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "sollw", varid)
236	guez	49	IF (ierr /= NF90_NOERR) THEN
237	guez	43	PRINT *, 'phyetat0: Le champ <sollw> est absent'
238			PRINT *, 'mis a zero'
239	guez	3	sollw = 0.
240			ELSE
241	guez	157	call nf95_get_var(ncid_startphy, varid, sollw)
242	guez	3	ENDIF
243
244			! Lecture derive des flux:
245
246	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "fder", varid)
247	guez	49	IF (ierr /= NF90_NOERR) THEN
248	guez	43	PRINT *, 'phyetat0: Le champ <fder> est absent'
249			PRINT *, 'mis a zero'
250	guez	3	fder = 0.
251			ELSE
252	guez	157	call nf95_get_var(ncid_startphy, varid, fder)
253	guez	3	ENDIF
254
255			! Lecture du rayonnement net au sol:
256
257	guez	157	call NF95_INQ_VARID(ncid_startphy, "RADS", varid)
258			call NF95_GET_VAR(ncid_startphy, varid, radsol)
259	guez	3
260			! Lecture de la longueur de rugosite
261
262	guez	157	call NF95_INQ_VARID(ncid_startphy, "RUG", varid)
263			call nf95_get_var(ncid_startphy, varid, frugs)
264	guez	3
265			! Lecture de l'age de la neige:
266
267	guez	157	call NF95_INQ_VARID(ncid_startphy, "AGESNO", varid)
268			call nf95_get_var(ncid_startphy, varid, agesno)
269	guez	3
270	guez	157	call NF95_INQ_VARID(ncid_startphy, "ZMEA", varid)
271			call NF95_GET_VAR(ncid_startphy, varid, zmea)
272	guez	3
273	guez	157	call NF95_INQ_VARID(ncid_startphy, "ZSTD", varid)
274			call NF95_GET_VAR(ncid_startphy, varid, zstd)
275	guez	3
276	guez	157	call NF95_INQ_VARID(ncid_startphy, "ZSIG", varid)
277			call NF95_GET_VAR(ncid_startphy, varid, zsig)
278	guez	3
279	guez	157	call NF95_INQ_VARID(ncid_startphy, "ZGAM", varid)
280			call NF95_GET_VAR(ncid_startphy, varid, zgam)
281	guez	3
282	guez	157	call NF95_INQ_VARID(ncid_startphy, "ZTHE", varid)
283			call NF95_GET_VAR(ncid_startphy, varid, zthe)
284	guez	3
285	guez	157	call NF95_INQ_VARID(ncid_startphy, "ZPIC", varid)
286			call NF95_GET_VAR(ncid_startphy, varid, zpic)
287	guez	3
288	guez	157	call NF95_INQ_VARID(ncid_startphy, "ZVAL", varid)
289			call NF95_GET_VAR(ncid_startphy, varid, zval)
290	guez	3
291			ancien_ok = .TRUE.
292
293	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "TANCIEN", varid)
294	guez	49	IF (ierr /= NF90_NOERR) THEN
295	guez	43	PRINT *, "phyetat0: Le champ <TANCIEN> est absent"
296			PRINT *, "Depart legerement fausse. Mais je continue"
297	guez	3	ancien_ok = .FALSE.
298			ELSE
299	guez	157	call nf95_get_var(ncid_startphy, varid, t_ancien)
300	guez	3	ENDIF
301
302	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "QANCIEN", varid)
303	guez	49	IF (ierr /= NF90_NOERR) THEN
304	guez	43	PRINT *, "phyetat0: Le champ <QANCIEN> est absent"
305			PRINT *, "Depart legerement fausse. Mais je continue"
306	guez	3	ancien_ok = .FALSE.
307			ELSE
308	guez	157	call nf95_get_var(ncid_startphy, varid, q_ancien)
309	guez	3	ENDIF
310
311	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "CLWCON", varid)
312	guez	49	IF (ierr /= NF90_NOERR) THEN
313	guez	43	PRINT *, "phyetat0: Le champ CLWCON est absent"
314			PRINT *, "Depart legerement fausse. Mais je continue"
315	guez	3	clwcon = 0.
316			ELSE
317	guez	157	call nf95_get_var(ncid_startphy, varid, clwcon(:, 1))
318	guez	72	clwcon(:, 2:) = 0.
319	guez	3	ENDIF
320
321	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "RNEBCON", varid)
322	guez	49	IF (ierr /= NF90_NOERR) THEN
323	guez	43	PRINT *, "phyetat0: Le champ RNEBCON est absent"
324			PRINT *, "Depart legerement fausse. Mais je continue"
325	guez	3	rnebcon = 0.
326			ELSE
327	guez	157	call nf95_get_var(ncid_startphy, varid, rnebcon(:, 1))
328	guez	72	rnebcon(:, 2:) = 0.
329	guez	3	ENDIF
330
331			! Lecture ratqs
332
333	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "RATQS", varid)
334	guez	49	IF (ierr /= NF90_NOERR) THEN
335	guez	43	PRINT *, "phyetat0: Le champ <RATQS> est absent"
336			PRINT *, "Depart legerement fausse. Mais je continue"
337	guez	3	ratqs = 0.
338			ELSE
339	guez	157	call nf95_get_var(ncid_startphy, varid, ratqs(:, 1))
340	guez	72	ratqs(:, 2:) = 0.
341	guez	3	ENDIF
342
343			! Lecture run_off_lic_0
344
345	guez	157	ierr = NF90_INQ_VARID(ncid_startphy, "RUNOFFLIC0", varid)
346	guez	49	IF (ierr /= NF90_NOERR) THEN
347	guez	43	PRINT *, "phyetat0: Le champ <RUNOFFLIC0> est absent"
348			PRINT *, "Depart legerement fausse. Mais je continue"
349	guez	3	run_off_lic_0 = 0.
350			ELSE
351	guez	157	call nf95_get_var(ncid_startphy, varid, run_off_lic_0)
352	guez	3	ENDIF
353
354	guez	157	call nf95_inq_varid(ncid_startphy, "sig1", varid)
355			call nf95_get_var(ncid_startphy, varid, sig1)
356	guez	72
357	guez	157	call nf95_inq_varid(ncid_startphy, "w01", varid)
358			call nf95_get_var(ncid_startphy, varid, w01)
359	guez	72
360	guez	3	END SUBROUTINE phyetat0
361
362			end module phyetat0_m