Ignore:
Timestamp:
Aug 2, 2024, 9:58:25 PM (3 months ago)
Author:
abarral
Message:

Put dimensions.h and paramet.h into modules

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/incloud_scav.f90

    r5144 r5159  
    1010  USE lmdz_yomcst
    1111
     12USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
     13  USE lmdz_paramet
    1214  IMPLICIT NONE
    1315
    14   INCLUDE "dimensions.h"
    1516  INCLUDE "chem.h"
    16   INCLUDE "paramet.h"
     17
    1718
    1819  !============================= INPUT ===================================
     
    4445
    4546  DO it = 1, nbtr
    46     !
     47
    4748    DO i = 1, klon
    4849      aux_var2(i) = his_dhlsc(i, it)
     
    5455      ENDDO
    5556    ENDDO
    56     !
     57
    5758    IF (lminmax) THEN
    5859      CALL minmaxqfi(aux_var1, qmin, qmax, 'avt inscav')
    5960      !nhl      CALL minmaxqfi(tr_seri(1,1,it),qmin,qmax,'avt inscav')
    6061    ENDIF
    61     !
     62
    6263    !nhl      CALL inscav_spl(pdtphys,it,masse(it),henry(it),kk(it),0.5e-3,
    6364    !nhl     .             prfl,psfl,zrho,zdz,t_seri,tr_seri(1,1,it),
     
    6566    CALL inscav_spl(pdtphys, it, masse(it), henry(it), kk(it), 0.5e-3, &
    6667            prfl, psfl, zrho, zdz, t_seri, aux_var1, aux_var2)
    67     !
     68
    6869    IF (lminmax) THEN
    6970      CALL minmaxqfi(aux_var1, qmin, qmax, 'depot humide lsc')
    7071      !nhl      CALL minmaxqfi(tr_seri(1,1,it),qmin,qmax,'depot humide lsc')
    7172    ENDIF
    72     !
    73     !
     73
     74
    7475    !-scheme for convective in-cloud scavenging
    75     !
     76
    7677    !nhl      CALL inscav_spl(pdtphys,it,masse(it),henry(it),kk(it),1.e-3,
    7778    !nhl     .             pmflxr,pmflxs,zrho,zdz,t_seri,tr_seri(1,1,it),
     
    7980    CALL inscav_spl(pdtphys, it, masse(it), henry(it), kk(it), 1.e-3, &
    8081            pmflxr, pmflxs, zrho, zdz, t_seri, aux_var1, aux_var3)
    81     !
     82
    8283    IF (lminmax) THEN
    8384      CALL minmaxqfi(aux_var1, qmin, qmax, 'depot humide con')
    8485      !nhl      CALL minmaxqfi(tr_seri(1,1,it),qmin,qmax,'depot humide con')
    8586    ENDIF
    86     !
     87
    8788    DO j = 1, klev
    8889      DO i = 1, klon
     
    9596    ENDDO
    9697
    97     !
     98
    9899  ENDDO !--boucle sur it
    99100
Note: See TracChangeset for help on using the changeset viewer.