Changeset 98 for LMDZ.3.3/branches/rel-LF/libf/phylmd/condsurf.F

Timestamp:

Jul 5, 2000, 4:58:04 PM (24 years ago)

Author:

lmdzadmin

Message:

Interface avec les differentes surface, version de travail.LF

File:

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

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

@@ -31 +31 @@
 
       LOGICAL newlmt
-      PARAMETER (newlmt=.FALSE.)
+      PARAMETER (newlmt=.TRUE.)
 
       INTEGER     nannemax
@@ -114 +114 @@
       IF (newlmt) THEN
 c
-c Fraction "ocean":
-      ierr = NF_INQ_VARID (nid, "FOCE", nvarid)
-      IF (ierr .NE. NF_NOERR) THEN
-         PRINT*, "condsurf: Le champ <FOCE> est absent"
-         CALL abort
-      ENDIF
-#ifdef NC_DOUBLE
-      ierr = NF_GET_VARA_DOUBLE(nid,nvarid,debut,epais,pctsrf(1,is_oce))
-#else
-      ierr = NF_GET_VARA_REAL(nid,nvarid,debut,epais,pctsrf(1,is_oce))
-#endif
-      IF (ierr .NE. NF_NOERR) THEN
-         PRINT*, "condsurf: Lecture echouee pour <FOCE>"
-         CALL abort
-      ENDIF
+c$$$c Fraction "ocean":
+c$$$      ierr = NF_INQ_VARID (nid, "FOCE", nvarid)
+c$$$      IF (ierr .NE. NF_NOERR) THEN
+c$$$         PRINT*, "condsurf: Le champ <FOCE> est absent"
+c$$$         CALL abort
+c$$$      ENDIF
+c$$$#ifdef NC_DOUBLE
+c$$$      ierr = NF_GET_VARA_DOUBLE(nid,nvarid,debut,epais,pctsrf(1,is_oce))
+c$$$#else
+c$$$      ierr = NF_GET_VARA_REAL(nid,nvarid,debut,epais,pctsrf(1,is_oce))
+c$$$#endif
+c$$$      IF (ierr .NE. NF_NOERR) THEN
+c$$$         PRINT*, "condsurf: Lecture echouee pour <FOCE>"
+c$$$         CALL abort
+c$$$      ENDIF
 c
 c Fraction "glace de mer":
+c
+c
       ierr = NF_INQ_VARID (nid, "FSIC", nvarid)
       IF (ierr .NE. NF_NOERR) THEN
@@ -145 +147 @@
          CALL abort
       ENDIF
-c
-c Fraction "terre":
-      ierr = NF_INQ_VARID (nid, "FTER", nvarid)
-      IF (ierr .NE. NF_NOERR) THEN
-         PRINT*, "condsurf: Le champ <FTER> est absent"
-         CALL abort
-      ENDIF
-#ifdef NC_DOUBLE
-      ierr = NF_GET_VARA_DOUBLE(nid,nvarid,debut,epais,pctsrf(1,is_ter))
-#else
-      ierr = NF_GET_VARA_REAL(nid,nvarid,debut,epais,pctsrf(1,is_ter))
-#endif
-      IF (ierr .NE. NF_NOERR) THEN
-         PRINT*, "condsurf: Lecture echouee pour <FTER>"
-         CALL abort
-      ENDIF
-c
-c Fraction "glacier terre":
-      ierr = NF_INQ_VARID (nid, "FLIC", nvarid)
-      IF (ierr .NE. NF_NOERR) THEN
-         PRINT*, "condsurf: Le champ <FLIC> est absent"
-         CALL abort
-      ENDIF
-#ifdef NC_DOUBLE
-      ierr = NF_GET_VARA_DOUBLE(nid,nvarid,debut,epais,pctsrf(1,is_lic))
-#else
-      ierr = NF_GET_VARA_REAL(nid,nvarid,debut,epais,pctsrf(1,is_lic))
-#endif
-      IF (ierr .NE. 0) THEN
-         PRINT*, "condsurf: Lecture echouee pour <FLIC>"
-         CALL abort
-      ENDIF
+C 
+C positionnement % ocean libre et verification qu'il y a compatibilite des soussurfaces
+      pctsrf(1 : klon, is_oce) = (1. - zmasq(1 : klon)) 
+     $    - pctsrf(1 : klon, is_sic)
+      DO i = 1, klon
+        IF ( pctsrf(i, is_sic) .GT. (1. - zmasq(i)) ) THEN
+            WRITE(*,*) 'condsurf : sea-ice et masque pb en ', i,
+     $     pctsrf(i, is_sic), (1. - zmasq(i))     
+            pctsrf(i, is_sic) = (1. - zmasq(i))
+            pctsrf(i, is_oce) = 0.
+        ENDIF          
+        IF ( abs( pctsrf(i, is_ter) + pctsrf(i, is_lic) + 
+     $       pctsrf(i, is_oce) + pctsrf(i, is_sic)  - 1.) .GT. EPSFRA) 
+     $       THEN 
+             WRITE(*,*) 'physiq : pb sous surface au point ', i, 
+     $           pctsrf(i, 1 : nbsrf)
+         ENDIF 
+      END DO 
+c
+c$$$c Fraction "terre":
+c$$$      ierr = NF_INQ_VARID (nid, "FTER", nvarid)
+c$$$      IF (ierr .NE. NF_NOERR) THEN
+c$$$         PRINT*, "condsurf: Le champ <FTER> est absent"
+c$$$         CALL abort
+c$$$      ENDIF
+c$$$#ifdef NC_DOUBLE
+c$$$      ierr = NF_GET_VARA_DOUBLE(nid,nvarid,debut,epais,pctsrf(1,is_ter))
+c$$$#else
+c$$$      ierr = NF_GET_VARA_REAL(nid,nvarid,debut,epais,pctsrf(1,is_ter))
+c$$$#endif
+c$$$      IF (ierr .NE. NF_NOERR) THEN
+c$$$         PRINT*, "condsurf: Lecture echouee pour <FTER>"
+c$$$         CALL abort
+c$$$      ENDIF
+c$$$c
+c$$$c Fraction "glacier terre":
+c$$$      ierr = NF_INQ_VARID (nid, "FLIC", nvarid)
+c$$$      IF (ierr .NE. NF_NOERR) THEN
+c$$$         PRINT*, "condsurf: Le champ <FLIC> est absent"
+c$$$         CALL abort
+c$$$      ENDIF
+c$$$#ifdef NC_DOUBLE
+c$$$      ierr = NF_GET_VARA_DOUBLE(nid,nvarid,debut,epais,pctsrf(1,is_lic))
+c$$$#else
+c$$$      ierr = NF_GET_VARA_REAL(nid,nvarid,debut,epais,pctsrf(1,is_lic))
+c$$$#endif
+c$$$      IF (ierr .NE. 0) THEN
+c$$$         PRINT*, "condsurf: Lecture echouee pour <FLIC>"
+c$$$         CALL abort
+c$$$      ENDIF
 c
       ELSE ! test sur newlmt