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

Diff of /trunk/phylmd/clmain.f

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

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 12 by guez, Mon Jul 21 16:05:07 2008 UTC
# Line 116  c pblT------- T au nveau HCL Line 116  c pblT------- T au nveau HCL
116  c======================================================================  c======================================================================
117  c$$$ PB ajout pour soil  c$$$ PB ajout pour soil
118  c  c
119        REAL dtime        REAL, intent(in):: dtime
120        real date0        real date0
121        integer itap        integer, intent(in):: itap
122        REAL t(klon,klev), q(klon,klev)        REAL t(klon,klev), q(klon,klev)
123        REAL u(klon,klev), v(klon,klev)        REAL u(klon,klev), v(klon,klev)
124  cIM 230604 BAD  REAL radsol(klon) ???  cIM 230604 BAD  REAL radsol(klon) ???
125        REAL, intent(in):: paprs(klon,klev+1)        REAL, intent(in):: paprs(klon,klev+1)
126        real pplay(klon,klev)        real, intent(in):: pplay(klon,klev)
127        REAL, intent(in):: rlon(klon), rlat(klon)        REAL, intent(in):: rlon(klon), rlat(klon)
128        real cufi(klon), cvfi(klon)        real cufi(klon), cvfi(klon)
129        REAL d_t(klon, klev), d_q(klon, klev)        REAL d_t(klon, klev), d_q(klon, klev)
# Line 149  cIM cf JLD Line 149  cIM cf JLD
149        LOGICAL, intent(in):: debut        LOGICAL, intent(in):: debut
150        logical, intent(in):: lafin        logical, intent(in):: lafin
151        logical ok_veget        logical ok_veget
152        character*6 ocean        character(len=*), intent(IN):: ocean
153        integer npas, nexca        integer npas, nexca
154  c  c
155        REAL pctsrf(klon,nbsrf)        REAL pctsrf(klon,nbsrf)
# Line 176  cAA Line 176  cAA
176        REAL zv1(klon)        REAL zv1(klon)
177  cAA  cAA
178  c$$$ PB ajout pour soil  c$$$ PB ajout pour soil
179        LOGICAL soil_model        LOGICAL, intent(in):: soil_model
180  cIM ajout seuils cdrm, cdrh  cIM ajout seuils cdrm, cdrh
181        REAL cdmmax, cdhmax        REAL cdmmax, cdhmax
182  cIM: 261103  cIM: 261103
# Line 356  CC        call flush(6) Line 356  CC        call flush(6)
356        ENDDO        ENDDO
357        ENDDO        ENDDO
358        DO i = 1, klon  ! vent de la premiere couche        DO i = 1, klon  ! vent de la premiere couche
 ccc         zx_alf1 = (paprs(i,1)-pplay(i,2))/(pplay(i,1)-pplay(i,2))  
359           zx_alf1 = 1.0           zx_alf1 = 1.0
360           zx_alf2 = 1.0 - zx_alf1           zx_alf2 = 1.0 - zx_alf1
361           u1lay(i) = u(i,1)*zx_alf1 + u(i,2)*zx_alf2           u1lay(i) = u(i,1)*zx_alf1 + u(i,2)*zx_alf2
# Line 658  c  du papier sur les thermiques. Line 657  c  du papier sur les thermiques.
657           call ustarhb(knon,yu,yv,y_cd_m, yustar)           call ustarhb(knon,yu,yv,y_cd_m, yustar)
658    
659          if (prt_level > 9) THEN          if (prt_level > 9) THEN
660            WRITE(lunout,*)'USTAR = ',yustar            print *,'USTAR = ',yustar
661          ENDIF          ENDIF
662    
663  c   iflag_pbl peut etre utilise comme longuer de melange  c   iflag_pbl peut etre utilise comme longuer de melange
# Line 912  c Line 911  c
911           enddo           enddo
912        enddo        enddo
913  cIM "slab" ocean  cIM "slab" ocean
       IF(OCEAN.EQ.'slab  '.OR.OCEAN.EQ.'force ') THEN  
914         IF (nsrf.EQ.is_oce) THEN         IF (nsrf.EQ.is_oce) THEN
915          DO j = 1, knon          DO j = 1, knon
916  c on projette sur la grille globale  c on projette sur la grille globale
# Line 936  cIM 230604 on pondere lorsque l'on fait Line 934  cIM 230604 on pondere lorsque l'on fait
934           ENDIF           ENDIF
935          ENDDO          ENDDO
936         ENDIF !nsrf.EQ.is_sic         ENDIF !nsrf.EQ.is_sic
       ENDIF !OCEAN  
937  c  c
938        IF(OCEAN.EQ.'slab  ') THEN        IF(OCEAN.EQ.'slab  ') THEN
939         IF(nsrf.EQ.is_oce) then         IF(nsrf.EQ.is_oce) then
Legend:
Removed from v.3  
changed lines
  Added in v.12
  ViewVC Help
Powered by ViewVC 1.1.21