Changeset 150 for LMDZ.3.3

Timestamp:

Oct 30, 2000, 11:53:53 AM (24 years ago)

Author:

lmdzadmin

Message:

Probleme d'initialisation f77/f90
LF

File:

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

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

@@ -226 +226 @@
 cAA      ENDDO
 cAA      ENDIF
+
 c
 c Boucler sur toutes les sous-fractions du sol:
@@ -332 +333 @@
       ytaux = y_flux_u(:,1)
       ytauy = y_flux_v(:,1)
+
 c calculer la diffusion de "q" et de "h"
       CALL clqh(dtime, itap, jour, debut,lafin,
@@ -679 +681 @@
 C Appel a interfsurf (appel generique) routine d'interface avec la surface
 
-      petAcoef=zx_ch(:,1)
-      peqAcoef=zx_cq(:,1)
-      petBcoef=zx_dh(:,1)
-      peqBcoef=zx_dq(:,1)
-      tq_cdrag=coef(:,1)
-      temp_air=t(:,1)
-      spechum=q(:,1)
-      p1lay = pplay(:,1)
-      zlev1 = delp(:,1)
-      swnet = swdown * (1 - albedo)
+      do i = 1, knon
+        petAcoef=zx_ch(i,1)
+        peqAcoef=zx_cq(i,1)
+        petBcoef=zx_dh(i,1)
+        peqBcoef=zx_dq(i,1)
+        tq_cdrag=coef(i,1)
+        temp_air=t(i,1)
+        spechum=q(i,1)
+        p1lay = pplay(i,1)
+        zlev1 = delp(i,1)
+        swnet(i) = swdown(i) * (1. - albedo(i))
+      enddo
 c En attendant mieux
       hum_air = 0.
@@ -707 +711 @@
      s pctsrf_new, agesno)
 
-      flux_t(:,1) = fluxsens
-      flux_q(:,1) = - evap
-      d_ts = tsurf_new - ts
-      albedo = alb_new
+      do i = 1, knon
+      flux_t(i,1) = fluxsens(i)
+      flux_q(i,1) = - evap(i)
+      d_ts(i) = tsurf_new(i) - ts(i)
+      albedo(i) = alb_new(i)
+      enddo
 
 c==== une fois on a zx_h_ts, on peut faire l'iteration ========