Changeset 109 for LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

Timestamp:

Jul 26, 2000, 2:58:36 PM (24 years ago)

Author:

lmdzadmin

Message:

Deplacement des dernieres lignes concernant le sol hors de la physique
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F (modified) (4 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

@@ -1509 +1509 @@
          idayvrai = NINT(xjour)
          PRINT *,' PHYS cond  julien ',julien,idayvrai
-         CALL condsurf(julien,idayvrai, pctsrf ,
-     .                  lmt_sst,lmt_alb,lmt_rug,lmt_bils  )
+c         CALL condsurf(julien,idayvrai, pctsrf ,
+c     .                  lmt_sst,lmt_alb,lmt_rug,lmt_bils  )
          CALL ozonecm( FLOAT(julien), rlat, paprs, wo)
       ENDIF
@@ -1535 +1535 @@
 c
       DO i = 1, klon
-         frugs(i,is_ter) = SQRT(lmt_rug(i)**2+rugoro(i)**2)
+        if (.not. ok_veget) then
+         frugs(i,is_ter) = SQRT(frugs(i,is_ter)**2+rugoro(i)**2)
          frugs(i,is_lic) = rugoro(i)
          frugs(i,is_oce) = rugmer(i)
@@ -1552 +1553 @@
       ENDDO
 c
+C calculs necessaires au calcul de l'albedo dans l'interface
+c
+      CALL orbite(FLOAT(julien),zlongi,dist)
+      IF (cycle_diurne) THEN
+        zdtime=dtime*FLOAT(radpas) ! pas de temps du rayonnement (s)
+        CALL zenang(zlongi,gmtime,zdtime,rlat,rlon,rmu0,fract)
+      ELSE
+        rmu0 = -999.999
+      ENDIF
+
       CALL clmain(dtime,pctsrf,
-     e            t_seri,q_seri,u_seri,v_seri,ok_veget,
-     e            ftsol,paprs,pplay,radsol,
-     e            fsnow,fqsol,fevap,falbe,
+     e            t_seri,q_seri,u_seri,v_seri,
+     e            julien, rmu0, 
+     e            ok_veget, ftsol,
+     e            paprs,pplay,radsol, fsnow,fqsol,fevap,falbe,
      e            rain_fall, snow_fall, solsw, sollw, fder,
      e            rlon, rlat, frugs,
@@ -1858 +1870 @@
 c
       IF (MOD(itaprad,radpas).EQ.0) THEN
-      CALL orbite(FLOAT(julien),zlongi,dist)
-      IF (cycle_diurne) THEN
-        zdtime=dtime*FLOAT(radpas) ! pas de temps du rayonnement (s)
-        CALL zenang(zlongi,gmtime,zdtime,rlat,rlon,rmu0,fract)
+c      CALL orbite(FLOAT(julien),zlongi,dist)
+c      IF (cycle_diurne) THEN
+c        zdtime=dtime*FLOAT(radpas) ! pas de temps du rayonnement (s)
+c        CALL zenang(zlongi,gmtime,zdtime,rlat,rlon,rmu0,fract)
 c        CALL zenith(zlongi,gmtime,rlat,rlon,rmu0,fract) !va disparaitre
-        CALL alboc_cd(rmu0,alb_eau)
-      ELSE
-        CALL angle(zlongi,rlat,fract,rmu0)
-        CALL alboc(FLOAT(julien),rlat,alb_eau)
-      ENDIF
-      CALL albsno(veget,agesno,alb_neig)
-      DO i = 1, klon
-         falbe(i,is_oce) = alb_eau(i)
-         IF (pctsrf(i,is_oce).GT.epsfra .AND. ftsol(i,is_oce).LT.271.35)
-     .   falbe(i,is_oce) = 0.6 ! pour slab_ocean
-         zfra = MAX(0.0,MIN(1.0,fsnow(i,is_lic)/(fsnow(i,is_lic)+10.0)))
-         falbe(i,is_lic) = alb_neig(i)*zfra + 0.6*(1.0-zfra)
-         zfra = MAX(0.0,MIN(1.0,fsnow(i,is_ter)/(fsnow(i,is_ter)+10.0)))
-         falbe(i,is_ter) = alb_neig(i)*zfra + lmt_alb(i)*(1.0-zfra)
-         zfra = MAX(0.0,MIN(1.0,fsnow(i,is_sic)/(fsnow(i,is_sic)+10.0)))
+c        CALL alboc_cd(rmu0,alb_eau)
+c      ELSE
+c        CALL angle(zlongi,rlat,fract,rmu0)
+c        CALL alboc(FLOAT(julien),rlat,alb_eau)
+c      ENDIF
+c      CALL albsno(veget,agesno,alb_neig)
+      DO i = 1, klon
+c         falbe(i,is_oce) = alb_eau(i)
+c         IF (pctsrf(i,is_oce).GT.epsfra .AND. ftsol(i,is_oce).LT.271.35)
+c     .   falbe(i,is_oce) = 0.6 ! pour slab_ocean
+c         zfra = MAX(0.0,MIN(1.0,fsnow(i,is_lic)/(fsnow(i,is_lic)+10.0)))
+c         falbe(i,is_lic) = alb_neig(i)*zfra + 0.6*(1.0-zfra)
+c         zfra = MAX(0.0,MIN(1.0,fsnow(i,is_ter)/(fsnow(i,is_ter)+10.0)))
+c         falbe(i,is_ter) = alb_neig(i)*zfra + lmt_alb(i)*(1.0-zfra)
+c         zfra = MAX(0.0,MIN(1.0,fsnow(i,is_sic)/(fsnow(i,is_sic)+10.0)))
          falbe(i,is_sic) = alb_neig(i)*zfra + 0.6*(1.0-zfra)
          albsol(i) = falbe(i,is_oce) * pctsrf(i,is_oce)