Changeset 1064


Ignore:
Timestamp:
Dec 17, 2008, 1:27:20 PM (16 years ago)
Author:
Laurent Fairhead
Message:

Correction: initialisation rh2m, qsat2m en dehors boucle nbsrf
IM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90

    r1009 r1064  
    1 !
    2 ! $Header$
    31!
    42MODULE pbl_surface_mod
     
    152150
    153151    IF (type_ocean /= 'slab  ' .AND. type_ocean /= 'force ' .AND. type_ocean /= 'couple') THEN
    154       WRITE(lunout,*)' *** Warning ***'
    155       WRITE(lunout,*)'Option couplage pour l''ocean = ', type_ocean
    156       abort_message='option pour l''ocean non valable'
    157       CALL abort_gcm(modname,abort_message,1)
     152       WRITE(lunout,*)' *** Warning ***'
     153       WRITE(lunout,*)'Option couplage pour l''ocean = ', type_ocean
     154       abort_message='option pour l''ocean non valable'
     155       CALL abort_gcm(modname,abort_message,1)
    158156    ENDIF
    159157
     
    255253    INCLUDE "temps.h"
    256254    INCLUDE "control.h"
    257 
    258255! Input variables
    259256!****************************************************************************************
     
    545542    ytsoil = 999999.
    546543
     544    rh2m(:)        = 0.
     545    qsat2m(:)      = 0.
    547546!****************************************************************************************
    548547! 3) - Calculate pressure thickness of each layer
     
    642641!
    643642!****************************************************************************************
    644     
     643   
    645644    loop_nbsrf: DO nsrf = 1, nbsrf
    646645
     
    988987       trmb2(:,nsrf)  = 0.        ! inhibition
    989988       trmb3(:,nsrf)  = 0.        ! Point Omega
    990        rh2m(:)        = 0.
    991        qsat2m(:)      = 0.
    992989
    993990#undef T2m     
Note: See TracChangeset for help on using the changeset viewer.