/[lmdze]/trunk/phylmd/phyredem.f
ViewVC logotype

Diff of /trunk/phylmd/phyredem.f

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

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