Changeset 1055 for trunk/LMDZ.TITAN/libf

Timestamp:

Oct 7, 2013, 2:40:02 PM (11 years ago)

Author:

slebonnois

Message:

SL: small details...

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• interface_surf.F90 (modified) (4 diffs)
• readstartphy.F (modified) (4 diffs)

trunk/LMDZ.TITAN/libf/phytitan/interface_surf.F90

@@ -45 +45 @@
       & tsol_rad, tsurf_new, alb_new)
 
-      use write_field_phy
       use cpdet_mod, only: cpdet
 
@@ -160 +159 @@
            CALL soil(dtime, knon, tsurf, tsoil,soilcap, soilflux)
            cal(1:knon) = zcp(1:knon) / soilcap(1:knon)
-! for tests:
-!  call writefield_phy('interfsurf_hq_zcp',zcp,1)
-!  call writefield_phy('interfsurf_hq_cal',cal,1)
-!  call writefield_phy('interfsurf_hq_soilcap',soilcap,1)
 !       print*,"DIAGNOSTIC SOIL"
 !       print*,"soilcap=",soilcap
@@ -196 +191 @@
      & tsurf_new, fluxsens, dflux_s)
 
-  use write_field_phy
   use cpdet_mod, only: t2tpot, tpot2t
 
@@ -332 +326 @@
   ENDDO
 
-! for tests: write output fields...
-!  call writefield_phy('calcul_fluxs_d_ts',d_ts,1)
-!  call writefield_phy('calcul_fluxs_fluxsens',fluxsens,1)
-!  call writefield_phy('calcul_fluxs_dflux_s',dflux_s,1)
-
   END SUBROUTINE calcul_fluxs
 !

trunk/LMDZ.TITAN/libf/phytitan/readstartphy.F

@@ -8 +8 @@
      .           albe, solsw, sollw,
      .           fder,radsol,
+     .    zmea, zstd, zsig, zgam, zthe, zpic, zval,
      .           tabcntr0)
 c======================================================================
@@ -35 +36 @@
       real solsw(ngridmx)
       real fder(ngridmx)
+      REAL zmea(ngridmx), zstd(ngridmx)
+      REAL zsig(ngridmx), zgam(ngridmx), zthe(ngridmx)
+      REAL zpic(ngridmx), zval(ngridmx)
       INTEGER length
       PARAMETER (length=100)
@@ -294 +298 @@
       ENDDO
       PRINT*,'Rayonnement net au sol radsol:', xmin, xmax
+
+c
+c Lecture des parametres orographie sous-maille:
+c
+      ierr = NF_INQ_VARID (nid, "ZMEA", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Le champ <ZMEA> est absent'
+         PRINT*, 'mis a zero'
+         zmea = 0.
+      ELSE
+#ifdef NC_DOUBLE
+       ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmea)
+#else
+       ierr = NF_GET_VAR_REAL(nid, nvarid, zmea)
+#endif
+       IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Lecture echouee pour <ZMEA>'
+         CALL abort
+       ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, ngridmx
+         xmin = MIN(zmea(i),xmin)
+         xmax = MAX(zmea(i),xmax)
+      ENDDO
+      PRINT*,'zmea:', xmin, xmax
+c
+      ierr = NF_INQ_VARID (nid, "ZSTD", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Le champ <ZSTD> est absent'
+         PRINT*, 'mis a zero'
+         zstd = 0.
+      ELSE
+#ifdef NC_DOUBLE
+       ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zstd)
+#else
+       ierr = NF_GET_VAR_REAL(nid, nvarid, zstd)
+#endif
+       IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Lecture echouee pour <ZSTD>'
+         CALL abort
+       ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, ngridmx
+         xmin = MIN(zstd(i),xmin)
+         xmax = MAX(zstd(i),xmax)
+      ENDDO
+      PRINT*,'zstd:', xmin, xmax
+c
+      ierr = NF_INQ_VARID (nid, "ZSIG", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Le champ <ZSIG> est absent'
+         PRINT*, 'mis a zero'
+         zsig = 0.
+      ELSE
+#ifdef NC_DOUBLE
+       ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zsig)
+#else
+       ierr = NF_GET_VAR_REAL(nid, nvarid, zsig)
+#endif
+       IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Lecture echouee pour <ZSIG>'
+         CALL abort
+       ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, ngridmx
+         xmin = MIN(zsig(i),xmin)
+         xmax = MAX(zsig(i),xmax)
+      ENDDO
+      PRINT*,'zsig:', xmin, xmax
+c
+      ierr = NF_INQ_VARID (nid, "ZGAM", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Le champ <ZGAM> est absent'
+         PRINT*, 'mis a zero'
+         zgam = 0.
+      ELSE
+#ifdef NC_DOUBLE
+       ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zgam)
+#else
+       ierr = NF_GET_VAR_REAL(nid, nvarid, zgam)
+#endif
+       IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Lecture echouee pour <ZGAM>'
+         CALL abort
+       ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, ngridmx
+         xmin = MIN(zgam(i),xmin)
+         xmax = MAX(zgam(i),xmax)
+      ENDDO
+      PRINT*,'zgam:', xmin, xmax
+c
+      ierr = NF_INQ_VARID (nid, "ZTHE", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Le champ <ZTHE> est absent'
+         PRINT*, 'mis a zero'
+         zthe = 0.
+      ELSE
+#ifdef NC_DOUBLE
+       ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zthe)
+#else
+       ierr = NF_GET_VAR_REAL(nid, nvarid, zthe)
+#endif
+       IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Lecture echouee pour <ZTHE>'
+         CALL abort
+       ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, ngridmx
+         xmin = MIN(zthe(i),xmin)
+         xmax = MAX(zthe(i),xmax)
+      ENDDO
+      PRINT*,'zthe:', xmin, xmax
+c
+      ierr = NF_INQ_VARID (nid, "ZPIC", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Le champ <ZPIC> est absent'
+         PRINT*, 'mis a zero'
+         zpic = 0.
+      ELSE
+#ifdef NC_DOUBLE
+       ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zpic)
+#else
+       ierr = NF_GET_VAR_REAL(nid, nvarid, zpic)
+#endif
+       IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Lecture echouee pour <ZPIC>'
+         CALL abort
+       ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, ngridmx
+         xmin = MIN(zpic(i),xmin)
+         xmax = MAX(zpic(i),xmax)
+      ENDDO
+      PRINT*,'zpic:', xmin, xmax
+c
+      ierr = NF_INQ_VARID (nid, "ZVAL", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Le champ <ZVAL> est absent'
+         PRINT*, 'mis a zero'
+         zval = 0.
+      ELSE
+#ifdef NC_DOUBLE
+       ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zval)
+#else
+       ierr = NF_GET_VAR_REAL(nid, nvarid, zval)
+#endif
+       IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Lecture echouee pour <ZVAL>'
+         CALL abort
+       ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, ngridmx
+         xmin = MIN(zval(i),xmin)
+         xmax = MAX(zval(i),xmax)
+      ENDDO
+      PRINT*,'zval:', xmin, xmax
 c
 c Fermer le fichier:
@@ -301 +476 @@
       RETURN
       END
+