Changeset 139 for LMDZ.3.3/branches/rel-LF/libf/phylmd

Timestamp:

Sep 15, 2000, 3:15:42 PM (24 years ago)

Author:

lmdzadmin

Message:

pctsrf_new n'est pas rempli pour terre et glacier

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90 (modified) (4 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

@@ -286 +286 @@
     endif    
 !
+! Remplissage des pourcentages de surface
+!
+    pctsrf_new(:,nisurf) = pctsrf(:,nisurf)
+
   else if (nisurf == is_oce) then
 
@@ -424 +428 @@
        enddo
        alb_new = alb_neig*zfra + 0.6 * (1.0-zfra)
+!
+! Remplissage des pourcentages de surface
+!
+    pctsrf_new(:,nisurf) = pctsrf(:,nisurf)
 
   else
@@ -429 +437 @@
     abort_message = 'Index surface non valable'
     call abort_gcm(modname,abort_message,1)
+  endif
+
+  if (check) then
+    write(*,*)'In ',modname
+    do ii = 1, nbsrf
+      write(*,*) 'surface, pctsrf_new',ii,pctsrf_new(:,ii)
+    enddo
   endif
 
@@ -1005 +1020 @@
       endif
     enddo
+    if (check) write(*,*)'In ',modname
+    if (check) write(*,*)' surface, pctsrf_new',is_oce,pctsrf_new(:,is_oce)
+    if (check) write(*,*)' surface, pctsrf_new',is_sic,pctsrf_new(:,is_sic)
   else if (nisurf == is_sic) then
       call cpl2gath(read_sit, tsurf_new, klon, knon,iim,jjm, knindex)