Changeset 2693 for LMDZ5/trunk/libf


Ignore:
Timestamp:
Nov 1, 2016, 10:08:00 AM (8 years ago)
Author:
oboucher
Message:

Cosmetic changes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90

    r2644 r2693  
    2020
    2121  USE infotrac_phy
     22  USE YOMCST
    2223
    2324  IMPLICIT NONE
     25
    2426  include "clesphys.h"
    25   include "YOMCST.h"
    26 
    2727
    2828  ! Input arguments
     
    104104     !--convert to ug m-3 unit for consistency with offline fields
    105105     !
    106      do i=1,nbtr
    107         select case(trim(solsym(i)))
    108            case ("ASBCM")
     106     DO i=1,nbtr
     107        SELECT CASE(trim(solsym(i)))
     108           CASE ("ASBCM")
    109109              id_ASBCM = i
    110            case ("ASPOMM")
     110           CASE ("ASPOMM")
    111111              id_ASPOMM = i
    112            case ("ASSO4M")
     112           CASE ("ASSO4M")
    113113              id_ASSO4M = i
    114            case ("ASMSAM")
     114           CASE ("ASMSAM")
    115115              id_ASMSAM = i
    116            case ("CSSO4M")
     116           CASE ("CSSO4M")
    117117              id_CSSO4M = i
    118            case ("CSMSAM")
     118           CASE ("CSMSAM")
    119119              id_CSMSAM = i
    120            case ("SSSSM")
     120           CASE ("SSSSM")
    121121              id_SSSSM = i
    122            case ("CSSSM")
     122           CASE ("CSSSM")
    123123              id_CSSSM = i
    124            case ("ASSSM")
     124           CASE ("ASSSM")
    125125              id_ASSSM = i
    126            case ("CIDUSTM")
     126           CASE ("CIDUSTM")
    127127              id_CIDUSTM = i
    128            case ("AIBCM")
     128           CASE ("AIBCM")
    129129              id_AIBCM = i
    130            case ("AIPOMM")
     130           CASE ("AIPOMM")
    131131              id_AIPOMM = i
    132            case ("ASNO3M")
     132           CASE ("ASNO3M")
    133133              id_ASNO3M = i
    134            case ("CSNO3M")
     134           CASE ("CSNO3M")
    135135              id_CSNO3M = i
    136            case ("CINO3M")
     136           CASE ("CINO3M")
    137137              id_CINO3M = i
    138            end select
    139      enddo
    140 
     138           END SELECT
     139     ENDDO
    141140
    142141     bcsol(:,:)        =   tr_seri(:,:,id_ASBCM)                         *zrho(:,:)*1.e9  ! ASBCM
     
    171170     !
    172171     ! Read and interpolate sulfate
    173      IF ( flag_aerosol .EQ. 1 .OR. &
    174           flag_aerosol .EQ. 6 ) THEN
     172     IF ( flag_aerosol .EQ. 1 .OR. flag_aerosol .EQ. 6 ) THEN
    175173
    176174        CALL readaerosol_interp(id_ASSO4M_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, sulfacc, sulfacc_pi,loadso4)
     
    178176        sulfacc(:,:) = 0. ; sulfacc_pi(:,:) = 0.
    179177        loadso4=0.
    180      END IF
     178     ENDIF
    181179
    182180     ! Read and interpolate bcsol and bcins
    183      IF ( flag_aerosol .EQ. 2 .OR. &
    184           flag_aerosol .EQ. 6 ) THEN
     181     IF ( flag_aerosol .EQ. 2 .OR. flag_aerosol .EQ. 6 ) THEN
    185182
    186183        ! Get bc aerosol distribution
     
    192189        bcins(:,:) = 0. ; bcins_pi(:,:) = 0.
    193190        loadbc=0.
    194      END IF
    195 
     191     ENDIF
    196192
    197193     ! Read and interpolate pomsol and pomins
    198      IF ( flag_aerosol .EQ. 3 .OR. &
    199           flag_aerosol .EQ. 6 ) THEN
     194     IF ( flag_aerosol .EQ. 3 .OR. flag_aerosol .EQ. 6 ) THEN
    200195
    201196        CALL readaerosol_interp(id_ASPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomsol, pomsol_pi, load_tmp3)
     
    206201        pomins(:,:) = 0. ; pomins_pi(:,:) = 0.
    207202        loadoa=0.
    208      END IF
    209 
     203     ENDIF
    210204
    211205     ! Read and interpolate csssm, ssssm, assssm
    212      IF (flag_aerosol .EQ. 4 .OR. &
    213           flag_aerosol .EQ. 6 ) THEN
     206     IF (flag_aerosol .EQ. 4 .OR. flag_aerosol .EQ. 6 ) THEN
    214207
    215208        CALL readaerosol_interp(id_SSSSM_phy ,itap, pdtphys,rjourvrai, &
     
    228221
    229222     ! Read and interpolate cidustm
    230      IF (flag_aerosol .EQ. 5 .OR.  &
    231           flag_aerosol .EQ. 6 ) THEN
     223     IF (flag_aerosol .EQ. 5 .OR. flag_aerosol .EQ. 6 ) THEN
    232224
    233225        CALL readaerosol_interp(id_CIDUSTM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, cidust, cidust_pi, loaddust)
     
    299291     DO i = 1, klon
    300292        pdel(i,k) = paprs(i,k) - paprs (i,k+1)
    301      END DO
    302   END DO
     293     ENDDO
     294  ENDDO
    303295
    304296!--new aerosol properties
Note: See TracChangeset for help on using the changeset viewer.