/[lmdze]/trunk/Sources/phylmd/phyetat0.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/phyetat0.f

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

revision 69 by guez, Mon Feb 18 16:33:12 2013 UTC revision 72 by guez, Tue Jul 23 13:00:07 2013 UTC
# Line 14  contains Line 14  contains
14         qsurf, qsol, snow, albe, alblw, evap, rain_fall, snow_fall, solsw, &         qsurf, qsol, snow, albe, alblw, evap, rain_fall, snow_fall, solsw, &
15         sollw, fder, radsol, frugs, agesno, zmea, zstd, zsig, zgam, zthe, &         sollw, fder, radsol, frugs, agesno, zmea, zstd, zsig, zgam, zthe, &
16         zpic, zval, t_ancien, q_ancien, ancien_ok, rnebcon, ratqs, clwcon, &         zpic, zval, t_ancien, q_ancien, ancien_ok, rnebcon, ratqs, clwcon, &
17         run_off_lic_0)         run_off_lic_0, sig1, w01)
18    
19      ! From phylmd/phyetat0.F, version 1.4 2005/06/03 10:03:07      ! From phylmd/phyetat0.F, version 1.4 2005/06/03 10:03:07
20      ! Author: Z.X. Li (LMD/CNRS)      ! Author: Z.X. Li (LMD/CNRS)
# Line 31  contains Line 31  contains
31      USE temps, ONLY : itau_phy      USE temps, ONLY : itau_phy
32    
33      CHARACTER(len=*), intent(in):: fichnom      CHARACTER(len=*), intent(in):: fichnom
34        REAL pctsrf(klon, nbsrf)
35      REAL tsol(klon, nbsrf)      REAL tsol(klon, nbsrf)
36      REAL tsoil(klon, nsoilmx, nbsrf)      REAL tsoil(klon, nsoilmx, nbsrf)
37        CHARACTER(len=*), intent(in):: ocean
38      REAL tslab(klon), seaice(klon)      REAL tslab(klon), seaice(klon)
39      REAL qsurf(klon, nbsrf)      REAL qsurf(klon, nbsrf)
40      REAL qsol(klon)      REAL qsol(klon)
# Line 40  contains Line 42  contains
42      REAL albe(klon, nbsrf)      REAL albe(klon, nbsrf)
43      REAL alblw(klon, nbsrf)      REAL alblw(klon, nbsrf)
44      REAL evap(klon, nbsrf)      REAL evap(klon, nbsrf)
     REAL radsol(klon)  
45      REAL, intent(out):: rain_fall(klon)      REAL, intent(out):: rain_fall(klon)
46      REAL snow_fall(klon)      REAL snow_fall(klon)
     REAL sollw(klon)  
47      real solsw(klon)      real solsw(klon)
48        REAL, intent(out):: sollw(klon)
49      real fder(klon)      real fder(klon)
50        REAL radsol(klon)
51      REAL frugs(klon, nbsrf)      REAL frugs(klon, nbsrf)
52      REAL agesno(klon, nbsrf)      REAL agesno(klon, nbsrf)
53      REAL zmea(klon)      REAL zmea(klon)
# Line 55  contains Line 57  contains
57      REAL zthe(klon)      REAL zthe(klon)
58      REAL zpic(klon)      REAL zpic(klon)
59      REAL zval(klon)      REAL zval(klon)
     REAL pctsrf(klon, nbsrf)  
     REAL fractint(klon)  
     REAL run_off_lic_0(klon)  
   
60      REAL t_ancien(klon, klev), q_ancien(klon, klev)      REAL t_ancien(klon, klev), q_ancien(klon, klev)
     real rnebcon(klon, klev), clwcon(klon, klev), ratqs(klon, klev)  
61      LOGICAL, intent(out):: ancien_ok      LOGICAL, intent(out):: ancien_ok
62        real rnebcon(klon, klev), ratqs(klon, klev), clwcon(klon, klev)
63        REAL run_off_lic_0(klon)
64        real, intent(out):: sig1(klon, klev) ! section adiabatic updraft
65    
66      CHARACTER(len=*), intent(in):: ocean      real, intent(out):: w01(klon, klev)
67        ! vertical velocity within adiabatic updraft
68    
69        ! Local:
70        REAL fractint(klon)
71      REAL xmin, xmax      REAL xmin, xmax
   
72      INTEGER ncid, varid      INTEGER ncid, varid
73      INTEGER ierr, i, nsrf, isoil      INTEGER ierr, i, nsrf, isoil
74      CHARACTER*7 str7      CHARACTER(len=7) str7
75      CHARACTER*2 str2      CHARACTER(len=2) str2
76    
77      !---------------------------------------------------------------      !---------------------------------------------------------------
78    
79      print *, "Call sequence information: phyetat0"      print *, "Call sequence information: phyetat0"
80    
81      ! Ouvrir le fichier contenant l'etat initial:      ! Fichier contenant l'état initial :
82      print *, 'fichnom = ', fichnom      print *, 'fichnom = ', fichnom
83      call NF95_OPEN(fichnom, NF90_NOWRITE, ncid)      call NF95_OPEN(fichnom, NF90_NOWRITE, ncid)
84    
# Line 516  contains Line 518  contains
518      ELSE      ELSE
519         call nf95_get_var(ncid, varid, sollw)         call nf95_get_var(ncid, varid, sollw)
520      ENDIF      ENDIF
521      xmin = 1.0E+20      PRINT *, 'Rayonnement IF au sol sollw:', minval(sollw), maxval(sollw)
     xmax = -1.0E+20  
     DO i = 1, klon  
        xmin = MIN(sollw(i), xmin)  
        xmax = MAX(sollw(i), xmax)  
     ENDDO  
     PRINT *, 'Rayonnement IF au sol sollw:', xmin, xmax  
522    
523      ! Lecture derive des flux:      ! Lecture derive des flux:
524    
# Line 734  contains Line 730  contains
730         PRINT *, "Depart legerement fausse. Mais je continue"         PRINT *, "Depart legerement fausse. Mais je continue"
731         clwcon = 0.         clwcon = 0.
732      ELSE      ELSE
733         call nf95_get_var(ncid, varid, clwcon)         call nf95_get_var(ncid, varid, clwcon(:, 1))
734           clwcon(:, 2:) = 0.
735      ENDIF      ENDIF
736      xmin = 1.0E+20      xmin = 1.0E+20
737      xmax = -1.0E+20      xmax = -1.0E+20
# Line 748  contains Line 745  contains
745         PRINT *, "Depart legerement fausse. Mais je continue"         PRINT *, "Depart legerement fausse. Mais je continue"
746         rnebcon = 0.         rnebcon = 0.
747      ELSE      ELSE
748         call nf95_get_var(ncid, varid, rnebcon)         call nf95_get_var(ncid, varid, rnebcon(:, 1))
749           rnebcon(:, 2:) = 0.
750      ENDIF      ENDIF
751      xmin = 1.0E+20      xmin = 1.0E+20
752      xmax = -1.0E+20      xmax = -1.0E+20
# Line 764  contains Line 762  contains
762         PRINT *, "Depart legerement fausse. Mais je continue"         PRINT *, "Depart legerement fausse. Mais je continue"
763         ratqs = 0.         ratqs = 0.
764      ELSE      ELSE
765         call nf95_get_var(ncid, varid, ratqs)         call nf95_get_var(ncid, varid, ratqs(:, 1))
766           ratqs(:, 2:) = 0.
767      ENDIF      ENDIF
768      xmin = 1.0E+20      xmin = 1.0E+20
769      xmax = -1.0E+20      xmax = -1.0E+20
# Line 788  contains Line 787  contains
787      xmax = MAXval(run_off_lic_0)      xmax = MAXval(run_off_lic_0)
788      PRINT *, '(ecart-type) run_off_lic_0:', xmin, xmax      PRINT *, '(ecart-type) run_off_lic_0:', xmin, xmax
789    
790        call nf95_inq_varid(ncid, "sig1", varid)
791        call nf95_get_var(ncid, varid, sig1)
792    
793        call nf95_inq_varid(ncid, "w01", varid)
794        call nf95_get_var(ncid, varid, w01)
795    
796      call NF95_CLOSE(ncid)      call NF95_CLOSE(ncid)
797    
798    END SUBROUTINE phyetat0    END SUBROUTINE phyetat0
Legend:
Removed from v.69  
changed lines
  Added in v.72
  ViewVC Help
Powered by ViewVC 1.1.21