Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

    r138 r139  
    286286    endif   
    287287!
     288! Remplissage des pourcentages de surface
     289!
     290    pctsrf_new(:,nisurf) = pctsrf(:,nisurf)
     291
    288292  else if (nisurf == is_oce) then
    289293
     
    424428       enddo
    425429       alb_new = alb_neig*zfra + 0.6 * (1.0-zfra)
     430!
     431! Remplissage des pourcentages de surface
     432!
     433    pctsrf_new(:,nisurf) = pctsrf(:,nisurf)
    426434
    427435  else
     
    429437    abort_message = 'Index surface non valable'
    430438    call abort_gcm(modname,abort_message,1)
     439  endif
     440
     441  if (check) then
     442    write(*,*)'In ',modname
     443    do ii = 1, nbsrf
     444      write(*,*) 'surface, pctsrf_new',ii,pctsrf_new(:,ii)
     445    enddo
    431446  endif
    432447
     
    10051020      endif
    10061021    enddo
     1022    if (check) write(*,*)'In ',modname
     1023    if (check) write(*,*)' surface, pctsrf_new',is_oce,pctsrf_new(:,is_oce)
     1024    if (check) write(*,*)' surface, pctsrf_new',is_sic,pctsrf_new(:,is_sic)
    10071025  else if (nisurf == is_sic) then
    10081026      call cpl2gath(read_sit, tsurf_new, klon, knon,iim,jjm, knindex)
Note: See TracChangeset for help on using the changeset viewer.