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

Timestamp:

Aug 6, 2003, 4:50:49 PM (21 years ago)

Author:

lmdzadmin

Message:

Modifs sur les seuils (cdrag etc...), inclusion des diagnostics ISCCP par Ionela
LF

File:

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

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

@@ -7 +7 @@
      .                  jour, rmu0,
      .                  ok_veget, ocean, npas, nexca, ts,
-     .                  soil_model,ftsoil,qsol,
+     .                  soil_model,cdmmax, cdhmax, ftsoil,qsol,
      .                  paprs,pplay,radsol,snow,qsurf,evap,albe,alblw,
      .                  fluxlat,
@@ -140 +140 @@
 c$$$ PB ajout pour soil
       LOGICAL soil_model
+cIM ajout seuils cdrm, cdrh
+      REAL cdmmax, cdhmax
       REAL ftsoil(klon,nsoilmx,nbsrf)
       REAL ytsoil(klon,nsoilmx)
@@ -481 +483 @@
       ENDDO
 c
-c
+cIM cf JLD : on seuille ycoefm et ycoefh
+      if (nsrf.eq.is_oce) then
+         do j=1,knon
+c           ycoefm(j,1)=min(ycoefm(j,1),1.1E-3)
+            ycoefm(j,1)=min(ycoefm(j,1),cdmmax)
+c           ycoefh(j,1)=min(ycoefh(j,1),1.1E-3)
+            ycoefh(j,1)=min(ycoefh(j,1),cdhmax)
+         enddo
+      endif
+
 c calculer la diffusion des vitesses "u" et "v"
       CALL clvent(knon,dtime,yu1,yv1,ycoefm,yt,yu,ypaprs,ypplay,ydelp,
@@ -493 +504 @@
 
 c FH modif sur le cdrag temperature
-      do i=1,knon
-         ycoefh(i,1)=ycoefm(i,1)*0.8
-      enddo
+c$$$PB : déplace dans clcdrag
+c$$$      do i=1,knon
+c$$$         ycoefh(i,1)=ycoefm(i,1)*0.8
+c$$$      enddo
 
 c calculer la diffusion de "q" et de "h"
@@ -520 +532 @@
       IF (nsrf.EQ.is_oce) THEN
       DO j = 1, knon
-         yrugm(j) = 0.018*ycoefm(j,1) * (yu1(j)**2+yv1(j)**2)/RG
+         yrugm(j) = 0.018*ycoefm(j,1) * (yu1(j)**2+yv1(j)**2)/RG 
+     $      +  0.11*14e-6 / sqrt(ycoefm(j,1) * (yu1(j)**2+yv1(j)**2))
          yrugm(j) = MAX(1.5e-05,yrugm(j))
       ENDDO
@@ -1286 +1299 @@
       REAL t_coup
       PARAMETER (t_coup=273.15)
+cIM
+      LOGICAL check
+      PARAMETER (check=.false.)
 c
 c contre-gradient pour la chaleur sensible: Kelvin/metre
@@ -1388 +1404 @@
       ENDDO
 
-c$$$      PRINT*,' isommet=',isommet,' knon=',knon
+      IF (check) THEN
+      PRINT*,' isommet=',isommet,' knon=',knon
+      ENDIF
 
       DO k = 2, isommet