Ignore:
Timestamp:
Oct 6, 2006, 12:38:30 PM (18 years ago)
Author:
Laurent Fairhead
Message:

Corrections bugs divers YM
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/V3_test/libf/phylmd/physiq.F

    r717 r720  
    28802880
    28812881c
     2882cjq - introduce the aerosol direct and first indirect radiative forcings
     2883cjq - Johannes Quaas, 27/11/2003 (quaas@lmd.jussieu.fr)
     2884      IF (ok_ade.OR.ok_aie) THEN
     2885         ! Get sulfate aerosol distribution
     2886         CALL readsulfate(rjourvrai, debut, sulfate)
     2887         CALL readsulfate_preind(rjourvrai, debut, sulfate_pi)
     2888
     2889         ! Calculate aerosol optical properties (Olivier Boucher)
     2890         CALL aeropt(pplay, paprs, t_seri, sulfate, rhcl,
     2891     .        tau_ae, piz_ae, cg_ae, aerindex)
     2892cym
     2893      ELSE
     2894        tau_ae(:,:,:)=0.0
     2895        piz_ae(:,:,:)=0.0
     2896        cg_ae(:,:,:)=0.0
     2897cym     
     2898      ENDIF
     2899
     2900c
    28822901cIM calcul nuages par le simulateur ISCCP
    28832902c
     
    29943013     . EXP((Lheat *qsat2m(i))/(RCPD*zt2m(i)))
    29953014      ENDDO
    2996 c
    2997 cjq - introduce the aerosol direct and first indirect radiative forcings
    2998 cjq - Johannes Quaas, 27/11/2003 (quaas@lmd.jussieu.fr)
    2999       IF (ok_ade.OR.ok_aie) THEN
    3000          ! Get sulfate aerosol distribution
    3001          CALL readsulfate(rjourvrai, debut, sulfate)
    3002          CALL readsulfate_preind(rjourvrai, debut, sulfate_pi)
    3003 
    3004          ! Calculate aerosol optical properties (Olivier Boucher)
    3005          CALL aeropt(pplay, paprs, t_seri, sulfate, rhcl,
    3006      .        tau_ae, piz_ae, cg_ae, aerindex)
    3007 cym
    3008       ELSE
    3009         tau_ae(:,:,:)=0.0
    3010         piz_ae(:,:,:)=0.0
    3011         cg_ae(:,:,:)=0.0
    3012 cym     
    3013       ENDIF
     3015
    30143016
    30153017#ifdef INCA
Note: See TracChangeset for help on using the changeset viewer.