Changeset 1151 for LMDZ4


Ignore:
Timestamp:
Apr 28, 2009, 6:33:49 PM (16 years ago)
Author:
jghattas
Message:
  • readaerosol : ENDIF mal place. Bug existant depuis la merge avec la branche LMDZ4_V3_patches. Le bug existe sur cette branche.
  • etat0_netcdf : manque des variables dans call conf_phys.
Location:
LMDZ4/branches/LMDZ4-dev/libf
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4-dev/libf/dyn3d/etat0_netcdf.F

    r1140 r1151  
    126126      logical              :: ok_journe, ok_mensuel, ok_instan, ok_hf
    127127      logical              :: ok_LES
    128       LOGICAL              :: ok_ade, ok_aie, aerosol_couple
     128      LOGICAL              :: ok_ade, ok_aie, aerosol_couple, new_aod
     129      INTEGER              :: flag_aerosol
    129130      REAL                 :: bl95_b0, bl95_b1
    130131      real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
    … …  
    169170     &                 iflag_ratqs,ratqsbas,ratqshaut,                  &
    170171     &                 ok_ade, ok_aie, aerosol_couple,                  &
     172     &                 flag_aerosol, new_aod,                           &
    171173     &                 bl95_b0, bl95_b1,                                &
    172174     &                 iflag_thermals,nsplit_thermals,tau_thermals,     &

  • LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/etat0_netcdf.F

    r1140 r1151  
    127127      logical              :: ok_journe, ok_mensuel, ok_instan, ok_hf
    128128      logical              :: ok_LES
    129       LOGICAL              :: ok_ade, ok_aie, aerosol_couple
     129      LOGICAL              :: ok_ade, ok_aie, aerosol_couple, new_aod
     130      INTEGER              :: flag_aerosol
    130131      REAL                 :: bl95_b0, bl95_b1
    131132      real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
    … …  
    169170     &                 iflag_ratqs,ratqsbas,ratqshaut,                  &
    170171     &                 ok_ade, ok_aie, aerosol_couple,                  &
     172     &                 flag_aerosol, new_aod,                           &
    171173     &                 bl95_b0, bl95_b1,                                &
    172174     &                 iflag_thermals,nsplit_thermals,tau_thermals,     &

  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol.F90

    r1150 r1151  
    143143                    WRITE(lunout,*) 'get_aero  iyr=', iyr,'   ',cyear
    144144                    CALL get_aero_fromfile(cyear, var_1, name_aero(id_aero))
    145                  ELSE
    146                     ! Read in data:
     145
     146                 ELSE ! Read data from 2 decades
    147147                    ! a) from actual 10-yr-period
    148148                    IF (iyr.LT.1900) THEN
    … …  
    159159                    WRITE(cyear,'(I4)') iyr1
    160160                    WRITE(lunout,*) 'get_aero  iyr=', iyr,'   ',cyear
    161                    
    162161                    CALL get_aero_fromfile(cyear, var_1, name_aero(id_aero))
    163162                   
    164                  ENDIF
    165                  ! If to read two decades:
    166                  IF (.NOT.lonlyone) THEN
    167                    
    168                     ! b) from the next following one
    169                     WRITE(cyear,'(I4)') iyr2
    170                     WRITE(lunout,*) 'get_aero  iyr=', iyr,'   ',cyear
    171                    
    172                     CALL get_aero_fromfile(cyear, var_2, name_aero(id_aero))
    173                    
    174                     ! Interpolate linarily to the actual year:
    175                     DO it=1,12
    176                        DO k=1,klev
    177                           DO j=1,jjm
    178                              DO i=1,iim
    179                                 var_1(i,j,k,it) = &
    180                                      var_1(i,j,k,it) - FLOAT(iyr-iyr1)/FLOAT(iyr2-iyr1) * &
    181                                      (var_1(i,j,k,it) - var_2(i,j,k,it))
     163                    ! If to read two decades:
     164                    IF (.NOT.lonlyone) THEN
     165                      
     166                       ! b) from the next following one
     167                       WRITE(cyear,'(I4)') iyr2
     168                       WRITE(lunout,*) 'get_aero  iyr=', iyr,'   ',cyear
     169                       
     170                       CALL get_aero_fromfile(cyear, var_2, name_aero(id_aero))
     171                       
     172                       ! Interpolate linarily to the actual year:
     173                       DO it=1,12
     174                          DO k=1,klev
     175                             DO j=1,jjm
     176                                DO i=1,iim
     177                                   var_1(i,j,k,it) = &
     178                                        var_1(i,j,k,it) - FLOAT(iyr-iyr1)/FLOAT(iyr2-iyr1) * &
     179                                        (var_1(i,j,k,it) - var_2(i,j,k,it))
     180                                ENDDO
    182181                             ENDDO
    183182                          ENDDO
    184183                       ENDDO
    185                     ENDDO
    186184                   
    187                  ENDIF     !lonlyone
     185                    ENDIF ! lonlyone
     186                 ENDIF ! iyr .LT. 1850       
    188187              ENDIF ! aer_type
    189188               
Note: See TracChangeset for help on using the changeset viewer.