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

Timestamp:

Jul 18, 2001, 1:28:31 PM (23 years ago)

Author:

lmdzadmin

Message:

Phasage avec la version de PB pour le sol, dlw (juillet 2001)
LF

File:

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

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

@@ -278 +278 @@
       INTEGER igwd,idx(klon),itest(klon)
 c
-      REAL agesno(klon)
+      REAL agesno(klon,nbsrf)
       SAVE agesno                 ! age de la neige
 c
@@ -340 +340 @@
       REAL evap(klon), devap(klon) ! evaporation et sa derivee
       REAL sens(klon), dsens(klon) ! chaleur sensible et sa derivee
+      REAL dlw(klon)    ! derivee infra rouge
       REAL bils(klon) ! bilan de chaleur au sol
       REAL fder(klon) ! Derive de flux (sensible et latente) 
@@ -627 +628 @@
      .       rlat,rlon,pctsrf, ftsol,ftsoil,deltat,fqsol,fsnow,
      .       falbe, fevap, rain_fall,snow_fall,solsw, sollwdown,
-     .       fder,radsol,frugs,agesno,clesphy0,
+     .       dlw,radsol,frugs,agesno,clesphy0,
      .       zmea,zstd,zsig,zgam,zthe,zpic,zval,rugoro,tabcntr0,
      .       t_ancien, q_ancien, ancien_ok )
@@ -1022 +1023 @@
      .                "ave(X)", zsto,zout)
 c
-         CALL histdef(nid_mth, "ages", "Snow age", "day",
-     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
-     .                "ave(X)", zsto,zout)
-c
          CALL histdef(nid_mth, "evap", "Evaporation", "mm/day",
      .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
@@ -1127 +1124 @@
      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
      $         "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
+     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
+     .                "ave(X)", zsto,zout)
+
          END DO
 C
@@ -1440 +1442 @@
      .                "inst(X)", zsto,zout)
 c
+        CALL histdef(nid_ins, "qsol", "Surface humidity", "mm",
+     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
+     .                "inst(X)", zsto,zout)
+c
+         CALL histdef(nid_ins, "rain", "Precipitation", "mm/day",
+     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
+     .                "inst(X)", zsto,zout)
+c
+         CALL histdef(nid_ins, "snow", "Snow fall", "mm/day",
+     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
+     .                "inst(X)", zsto,zout)
+c
+         CALL histdef(nid_ins, "snow_cov", "Snow cover", "mm",
+     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
+     .                "inst(X)", zsto,zout)
+c
          CALL histdef(nid_ins, "topl", "OLR", "W/m2",
      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
@@ -1458 +1476 @@
          CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface", 
      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
-     .                "ave(X)", zsto,zout)
+     .                "inst(X)", zsto,zout)
 c
          CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2",
@@ -1521 +1539 @@
 c
            call histdef(nid_ins, "albe_"//clnsurf(nsrf), 
-     $         "Albedo surf. "//clnsurf(nsrf), "-",  
+     $         "Albedo "//clnsurf(nsrf), "-",  
      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
      $         "inst(X)", zsto,zout)
@@ -1540 +1558 @@
      .                "inst(X)", zsto,zout)
 c
-         CALL histdef(nid_ins, "snow_cov", "Snow cover", "mm",
-     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
-     .                "inst(X)", zsto,zout)
 c
 c Champs 3D:
@@ -1741 +1756 @@
       ENDIF
 
-      fder = 0.
-      date0 = day_ini
+      fder = dlw
 
       CALL clmain(dtime,itap,date0,pctsrf,
@@ -1784 +1798 @@
 c         evap(i) = - fluxq(i,1) ! flux d'evaporation au sol
          evap(i) = - zxfluxq(i,1) ! flux d'evaporation au sol
-C LF test signe flux
-         sens(i) = zxfluxt(i,1)
-         evap(i) = zxfluxq(i,1)
-         fder(i) = dsens(i) + devap(i)
+         fder(i) = dlw(i) + dsens(i) + devap(i)
       ENDDO
 
@@ -1831 +1842 @@
 c Calculer la derive du flux infrarouge
 c
-      DO nsrf = 1, nbsrf
-      DO i = 1, klon
-         fder(i) = fder(i) - 4.0*RSIGMA*zxtsol(i)**3 *
-     .                       (ftsol(i,nsrf)-zxtsol(i))
-     .                      *pctsrf(i,nsrf)
-      ENDDO
+c$$$      DO nsrf = 1, nbsrf
+      DO i = 1, klon
+c$$$        IF (pctsrf(i,nsrf) .GE. EPSFRA) THEN 
+            dlw(i) = - 4.0*RSIGMA*zxtsol(i)**3 
+c$$$     .          *(ftsol(i,nsrf)-zxtsol(i))
+c$$$     .          *pctsrf(i,nsrf)
+c$$$        ENDIF 
+c$$$      ENDDO
       ENDDO
 c
@@ -2536 +2549 @@
       CALL histwrite(nid_mth,"snow_cov",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
 c
-      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
-      CALL histwrite(nid_mth,"ages",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
-c
       CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
       CALL histwrite(nid_mth,"evap",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
@@ -2632 +2642 @@
         CALL histwrite(nid_mth,"rugs_"//clnsurf(nsrf),itap,
      $      zx_tmp_2d,iim*jjmp1,ndex2d)
+c
+      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
+      CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),itap
+     $    ,zx_tmp_2d,iim*jjmp1,ndex2d)
 
       END DO  
@@ -2908 +2923 @@
       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
       CALL histwrite(nid_ins,"tsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
+c
+      DO i = 1, klon
+         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
+      ENDDO
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
+      CALL histwrite(nid_ins,"rain",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
+c
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
+      CALL histwrite(nid_ins,"snow",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
+
 c
       CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
@@ -3116 +3141 @@
      .      rlat, rlon, pctsrf, ftsol, ftsoil, deltat, fqsol, fsnow,
      .      falbe, fevap, rain_fall, snow_fall,
-     .      solsw, sollwdown,fder,
+     .      solsw, sollwdown,dlw,
      .      radsol,frugs,agesno,
      .      zmea,zstd,zsig,zgam,zthe,zpic,zval,rugoro,