Changeset 852 for LMDZ4/trunk/libf/phytherm/phyetat0.F

Timestamp:

Oct 11, 2007, 3:43:42 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Mise a jour de la physique avec thermiques avec la version de FH d'aout 2007
LF

File:

• LMDZ4/trunk/libf/phytherm/phyetat0.F (modified) (7 diffs)

LMDZ4/trunk/libf/phytherm/phyetat0.F

@@ -13 +13 @@
      .           rugsrel_p,tabcntr0,
      .           t_ancien_p,q_ancien_p,ancien_ok_p, rnebcon_p, ratqs_p,
-     .           clwcon_p)
+     .           clwcon_p,pbl_tke_p)
 
       USE dimphy
@@ -37 +37 @@
 #include "clesphys.h"
 #include "temps.h"
+#include "thermcell.h"
+#include "compbl.h"
 c======================================================================
       CHARACTER*(*) fichnom
@@ -45 +47 @@
       REAL solaire_etat0
       REAL tsol_p(klon,nbsrf)
+      REAL pbl_tke_p(klon,klev,nbsrf)
       REAL tsoil_p(klon,nsoilmx,nbsrf)
       REAL tslab_p(klon), seaice_p(klon)
@@ -83 +86 @@
       REAL rlat(klon_glo), rlon(klon_glo)
       REAL tsol(klon_glo,nbsrf)
+      REAL pbl_tke(klon_glo,klev,nbsrf)
       REAL tsoil(klon_glo,nsoilmx,nbsrf)
 cIM "slab" ocean
@@ -127 +131 @@
 c
       INTEGER nid, nvarid
-      INTEGER ierr, i, nsrf, isoil 
+      INTEGER ierr, i, nsrf, isoil ,k
       INTEGER length
       PARAMETER (length=100)
@@ -1492 +1496 @@
       xmax = MAXval(run_off_lic_0)
       PRINT*,'(ecart-type) run_off_lic_0:', xmin, xmax
+
+
+c Lecture de l'energie cinetique turbulente
+c
+
+      IF (iflag_pbl>1) then
+         PRINT*, 'phyetat0: Le champ <TKE> est absent'
+         PRINT*, '          Mais je vais essayer de lire TKE**'
+         DO nsrf = 1, nbsrf
+           IF (nsrf.GT.99) THEN
+             PRINT*, "Trop de sous-mailles"
+             CALL abort
+           ENDIF
+           WRITE(str2,'(i2.2)') nsrf
+           ierr = NF_INQ_VARID (nid, "TKE"//str2, nvarid)
+           IF (ierr.NE.NF_NOERR) THEN
+              PRINT*, "WARNING phyetat0: <TKE"//str2//"> est absent"
+              pbl_tke(:,:,nsrf)=1.e-8
+           ELSE
+#ifdef NC_DOUBLE
+              ierr = NF_GET_VAR_DOUBLE(nid, nvarid, pbl_tke(1,1,nsrf))
+#else
+              ierr = NF_GET_VAR_REAL(nid, nvarid, pbl_tke(1,1,nsrf))
+#endif
+              IF (ierr.NE.NF_NOERR) THEN
+                PRINT*, "WARNING phyetat0: echec lect <TKE"//str2//">"
+                CALL abort
+              ENDIF
+           ENDIF
+
+           xmin = 1.0E+20
+           xmax = -1.0E+20
+           DO k = 1, klev
+           DO i = 1, klon_glo
+              xmin = MIN(pbl_tke(i,k,nsrf),xmin)
+              xmax = MAX(pbl_tke(i,k,nsrf),xmax)
+           ENDDO
+           ENDDO
+           PRINT*,'Temperature du sol TKE**:', nsrf, xmin, xmax
+         ENDDO
+      ENDIF
+
 c
 c Fermer le fichier:
@@ -1516 +1562 @@
       call Scatter( rlon,rlon_p)
       call Scatter( tsol,tsol_p)
+      IF (iflag_pbl>1) then
+         call Scatter( pbl_tke,pbl_tke_p)
+      endif
       call Scatter( tsoil,tsoil_p)
       call Scatter( tslab,tslab_p)