Ignore:
Timestamp:
Jul 22, 2024, 9:29:09 PM (4 months ago)
Author:
abarral
Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/phylmd/readaerosol_interp.F90

    r2840 r5099  
    11! $Id$
    2 !
     2
    33SUBROUTINE readaerosol_interp(id_aero, itap, pdtphys, r_day, first, pplay, paprs, t_seri, mass_out, pi_mass_out, load_src)
    4 !
     4
    55! This routine will return the mass concentration at actual day(mass_out) and
    66! the pre-industrial values(pi_mass_out) for aerosol corresponding to "id_aero".
    77! The mass concentrations for all aerosols are saved in this routine but each
    88! call to this routine only treats the aerosol "id_aero".
    9 !
     9
    1010! 1) Read in data for the whole year, only at first time step
    1111! 2) Interpolate to the actual day, only at new day
     
    2929  INCLUDE "clesphys.h"
    3030
    31 !
    3231! Input:
    3332!****************************************************************************************
     
    4039  REAL, DIMENSION(klon,klev+1),INTENT(IN):: paprs  ! pression between model layers
    4140  REAL, DIMENSION(klon,klev), INTENT(IN) :: t_seri ! air temperature
    42 !     
     41
    4342! Output:     
    4443!****************************************************************************************
     
    4645  REAL, INTENT(OUT) :: pi_mass_out(klon,klev) ! Mass of preindustrial aerosol (monthly mean data,from file) [ug AIBCM/m3]
    4746  REAL, INTENT(OUT) :: load_src(klon) ! Load of aerosol (monthly mean data,from file) [kg/m3]
    48 !     
     47
    4948! Local Variables:
    5049!****************************************************************************************
     
    9796!****************************************************************************************
    9897! Initialization
    99 !
     98
    10099!****************************************************************************************
    101100
     
    167166! 1) Read in data : corresponding to the actual year and preindustrial data.
    168167!    Only for the first day of the year.
    169 !
     168
    170169!****************************************************************************************
    171170  IF ( (first .OR. iday==0) .AND. lnewday ) THEN
     
    270269
    271270!    Calendar initialisation
    272 !
     271
    273272     DO i = 2, 13
    274273       month_len(i) = REAL(ioget_mon_len(year_cur, i-1))
     
    291290! - 2) Interpolate to the actual day.
    292291! - 3) Interpolate to the model vertical grid.
    293 !
     292
    294293!****************************************************************************************
    295294
     
    297296!****************************************************************************************
    298297! 2) Interpolate to the actual day
    299 !
     298
    300299!****************************************************************************************
    301300    ! Find which months and days to use for time interpolation
     
    385384!****************************************************************************************
    386385! 3) Interpolate to the model vertical grid (target grid)
    387 !
     386
    388387!****************************************************************************************
    389388
     
    527526!****************************************************************************************
    528527! 4) Test for negative mass values
    529 !
     528
    530529!****************************************************************************************
    531530     IF (MINVAL(var_day(:,:,id_aero)) < 0.) THEN
     
    571570!****************************************************************************************
    572571! Copy output from saved variables
    573 !
     572
    574573!****************************************************************************************
    575574
Note: See TracChangeset for help on using the changeset viewer.