Ignore:
Timestamp:
Sep 20, 2001, 12:44:47 PM (23 years ago)
Author:
lmdzadmin
Message:

Prise en compte du runoff et changement de signe de l'evap sortie
d'orchidee
LF

File:
1 edited

Legend:

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

    r274 r275  
    262262      allocate(run_off(knon), stat = error)
    263263      if (error /= 0) then
    264         abort_message='Pb allocation runoff'
     264        abort_message='Pb allocation run_off'
    265265        call abort_gcm(modname,abort_message,1)
    266266      endif
     
    270270      abort_message='voir ci-dessus'
    271271      call abort_gcm(modname,abort_message,1)
    272       deallocate(coastalflow, stat = error)
    273       allocate(coastalflow(knon), stat = error)
    274       if (error /= 0) then
    275         abort_message='Pb allocation coastalflow'
    276         call abort_gcm(modname,abort_message,1)
    277       endif
    278       deallocate(riverflow, stat = error)
    279       allocate(riverflow(knon), stat = error)
    280       if (error /= 0) then
    281         abort_message='Pb allocation riverflow'
    282         call abort_gcm(modname,abort_message,1)
    283       endif
    284       deallocate(run_off, stat = error)
    285       allocate(run_off(knon), stat = error)
    286       if (error /= 0) then
    287         abort_message='Pb allocation run_off'
    288         call abort_gcm(modname,abort_message,1)
    289       endif
    290     endif
    291     run_off = 0.
     272    endif
     273    coastalflow = 0.
     274    riverflow = 0.
    292275!
    293276! Calcul age de la neige
     
    970953  alb_new(:) = albedo_out(:,1)
    971954
    972 ! LF essai sensible
     955! Convention orchidee: positif vers le haut
    973956  fluxsens = -1. * fluxsens
    974957  fluxlat  = -1. * fluxlat
     958  evap     = -1. * evap
    975959
    976960  if (debut) lrestart_read = .false.
     
    12531237       &                          + tauy(ig)        / FLOAT(nexca)
    12541238      cpl_rriv(ig,cpl_index) = cpl_rriv(ig,cpl_index) &
    1255        &                          + 0.     / FLOAT(nexca)/dtime
     1239       &                          + riverflow      / FLOAT(nexca)/dtime
    12561240      cpl_rcoa(ig,cpl_index) = cpl_rcoa(ig,cpl_index) &
    1257        &                          + 0.     / FLOAT(nexca)/dtime
     1241       &                          + coastalflow    / FLOAT(nexca)/dtime
    12581242    enddo
    12591243  endif
     
    20121996  if (check) write(*,*)'Entree ', modname,' surface = ',nisurf
    20131997
    2014   if (size(run_off) /= knon .AND. nisurf == is_ter) then
     1998  if (size(coastalflow) /= knon .AND. nisurf == is_ter) then
    20151999    write(*,*)'Bizarre, le nombre de points continentaux'
    20162000    write(*,*)'a change entre deux appels. J''arrete ...'
Note: See TracChangeset for help on using the changeset viewer.