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, intent(in):: itap |
integer, intent(in):: itap |
122 |
REAL t(klon,klev), q(klon,klev) |
REAL t(klon,klev), q(klon,klev) |
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) |
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 |
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 |
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 |
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 |