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

Timestamp:

Dec 19, 2002, 5:46:39 PM (22 years ago)

Author:

lmdzadmin

Message:

Convergence avec la version de Ionela dec 2002

YOMCST.? : suppression RI0 (IM)
albedo.F : facteur 1.2 sur le nouveau calcul (IM)
clesphys.h : rajout de différentes ctes (concentration des gaz) (IM)
clmain.F : separation des flux LW, SW (JLD)

remplace qsurf par yqsol (IM)

conf_phys.F90 : rajout de différentes ctes (gaz + orbite) (IM)
convect3.F : DPINV+SIGD*0.5*(EVAP(1)+EVAP(2)) (SBL)
cv3_routines.F:
cvparam3.h : compatibilite avec conema3 TEMPORAIRE (FH)
phyetat0.F : lecture de co2_ppm et solaire pour tests de coherence
phyredem.F : co2_ppm et solaire passé en common
physiq.F : separation flux LW, SW

rajout diagnostiques (slp, w500)
suppression iflag_con = 4
clwcon0=qcondc (FH)
position dU "ENDIF ! ok_cvl"

radlwsw.F : passage des concentrations gaz dans un common (IM)

PEMIS(i) = 1.0 (JLD pour cohérence ORCHIDEE)

stdlevvar.F90 :
suphec.F : suppression init. des ctes orbitales (IM)

nouvelles E/S (ini_hist..., write_hist...)

File:

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

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

@@ -119 +119 @@
       real rain_f(klon), snow_f(klon)
       REAL fder(klon)
-      REAL sollw(klon), solsw(klon), sollwdown(klon)
+cIM cf. JLD   REAL sollw(klon), solsw(klon), sollwdown(klon)
+      REAL sollw(klon,nbsrf), solsw(klon,nbsrf), sollwdown(klon)
       REAL rugos(klon,nbsrf)
 C la nouvelle repartition des surfaces sortie de l'interface
@@ -407 +408 @@
         ytaux(j) = flux_u(i,1,nsrf)
         ytauy(j) = flux_v(i,1,nsrf)
-c$$$        ysolsw(j) = solsw(i)
-       ysolsw(j) = (1 - albe(i,nsrf)) 
-     $      /(1 - pctsrf(i,is_ter) * albe(i,is_ter) 
-     $      - pctsrf(i, is_lic) *albe(i,is_lic) 
-     $      - pctsrf(i, is_oce) *albe(i,is_oce)
-     $      - pctsrf(i, is_sic) *albe(i,is_sic) 
-     $      ) * solsw(i)
-        ysollw(j) = sollw(i)
+c$$$    ysolsw(j) = solsw(i)
+cIM cf. JLD
+        ysolsw(j) = solsw(i,nsrf)
+c       ysolsw(j) = (1 - albe(i,nsrf)) 
+c    $      /(1 - pctsrf(i,is_ter) * albe(i,is_ter) 
+c    $      - pctsrf(i, is_lic) *albe(i,is_lic) 
+c    $      - pctsrf(i, is_oce) *albe(i,is_oce)
+c    $      - pctsrf(i, is_sic) *albe(i,is_sic) 
+c    $      ) * solsw(i)
+cIM cf. JLD     ysollw(j) = sollw(i)
+        ysollw(j) = sollw(i,nsrf)
         ysollwdown(j) = sollwdown(i)
         yrugos(j) = rugos(i,nsrf)
@@ -421 +425 @@
         yv1(j) = v1lay(i)
 c$$$        yrads(j) = totalflu(i)
-        yrads(j) = (1 - albe(i,nsrf)) 
-     $      /(1 - pctsrf(i,is_ter) * albe(i,is_ter) 
-     $      - pctsrf(i, is_lic) *albe(i,is_lic) 
-     $      - pctsrf(i, is_oce) *albe(i,is_oce)
-     $      - pctsrf(i, is_sic) *albe(i,is_sic) 
-     $      ) * solsw(i) + sollw(i)
+cIM cf. JLD
+        yrads(j) =  ysolsw(j)+ ysollw(j)
+c       yrads(j) = (1 - albe(i,nsrf)) 
+c    $      /(1 - pctsrf(i,is_ter) * albe(i,is_ter) 
+c    $      - pctsrf(i, is_lic) *albe(i,is_lic) 
+c    $      - pctsrf(i, is_oce) *albe(i,is_oce)
+c    $      - pctsrf(i, is_sic) *albe(i,is_sic) 
+cIM cf. JLD  $      ) * solsw(i) + sollw(i)
+c    $      ) * solsw(i) + ysollw(j)
         ypaprs(j,klev+1) = paprs(i,klev+1)
       END DO
@@ -453 +460 @@
       CALL coefkz(nsrf, knon, ypaprs, ypplay,
      .            yts, yrugos, yu, yv, yt, yq,
+cIM remplace qsurf par yqsol
+     .            yqsol, 
      .            ycoefm, ycoefh)
       CALL coefkz2(nsrf, knon, ypaprs, ypplay,yt,
@@ -595 +604 @@
 c
 #undef T2m     
+#define T2m     
 #ifdef T2m
 ccc diagnostic t,q a 2m et u, v a 10m
@@ -1136 +1146 @@
      .                  ts, rugos,
      .                  u,v,t,q,
+cIM remplace qsurf par yqsol
+     .                  qsol, 
      .                  pcfm, pcfh)
       IMPLICIT none
@@ -1226 +1238 @@
       REAL gamt(2:klev)
 c essai qsurf
-      real qsurf(klon)
+cIM   real qsurf(klon)
+      real qsol(klon) 
 c
       LOGICAL appel1er
@@ -1264 +1277 @@
       ENDDO
 
+cIM remplace qsurf par qsol
+      IF(nsrf.NE.1) THEN
       do i = 1, knon
-        qsurf(i) = qsatl(ts(i))/paprs(i,1)
+cIM     qsurf(i) = qsatl(ts(i))/paprs(i,1)
+        qsol(i) = qsatl(ts(i))/paprs(i,1)
       enddo
+      ENDIF
 
 c
@@ -1298 +1315 @@
 c Calculer le frottement au sol (Cdrag)
 c
-c     DO i = 1, knon
-      DO i = 1, klon
+      DO i = 1, knon
        u1(i) = u(i,1)
        v1(i) = v(i,1)
@@ -1309 +1325 @@
       CALL clcdrag(klon, knon, nsrf, zxli, 
      $             u1, v1, t1, q1, z1,
-     $             ts, qsurf, rugos,
+     $             ts, qsol, rugos,
      $             pcfm1, pcfh1) 
+cIM  $             ts, qsurf, rugos,
 C
       DO i = 1, knon