Ignore:
Timestamp:
Aug 2, 2024, 2:12:03 PM (3 months ago)
Author:
abarral
Message:

Add missing klon on strataer_emiss_mod.F90
Correct various missing explicit declarations
Replace tabs by spaces (tabs are not part of the fortran charset)
Continue cleaning modules
Removed unused arguments and variables

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_plume_6A.f90

    r5119 r5158  
    170170    ! dans une couche l>1
    171171    !------------------------------------------------------------------------------
    172     do ig = 1, ngrid
     172    DO ig = 1, ngrid
    173173      ! Le panache va prendre au debut les caracteristiques de l'air contenu
    174174      ! dans cette couche.
    … …  
    189189    !boucle de calcul de la vitesse verticale dans le thermique
    190190    !==============================================================================
    191     do l = 2, nlay - 1
     191    DO l = 2, nlay - 1
    192192      !==============================================================================
    193193
    … …  
    195195      ! On decide si le thermique est encore actif ou non
    196196      ! AFaire : Il faut sans doute ajouter entr_star a alim_star dans ce test
    197       do ig = 1, ngrid
     197      DO ig = 1, ngrid
    198198        active(ig) = active(ig) &
    199199                &                 .AND. zw2(ig, l)>1.e-10 &
    … …  
    215215      ztemp(:) = zpspsk(:, l) * ztla(:, l - 1)
    216216      CALL thermcell_qsat(ngrid, active, pplev(:, l), ztemp, zqta(:, l - 1), zqsat(:))
    217       do ig = 1, ngrid
     217      DO ig = 1, ngrid
    218218        !       PRINT*,'active',active(ig),ig,l
    219219        IF(active(ig)) THEN
    … …  
    252252            ! afin de faire moins de calcul dans la boucle
    253253            !--------------------------------------------------
    254             do while (zlmelup>zltup)
     254            DO while (zlmelup>zltup)
    255255              lt = lt + 1
    256256              zlt = zlev(ig, lt)
    … …  
    334334            zdz2 = zlev(ig, lt) - zlev(ig, l)
    335335
    336             do while (lmel>zdz2)
     336            DO while (lmel>zdz2)
    337337              lt = lt + 1
    338338              zlt = zlev(ig, lt)
    … …  
    440440      !-------------------------------------------------
    441441
    442       do ig = 1, ngrid
     442      DO ig = 1, ngrid
    443443        IF (active(ig)) THEN
    444444          !          zw2m=max(0.5*(w_est(ig,l)+w_est(ig,l+1)),0.1)
    … …  
    530530
    531531      activetmp(:) = active(:) .AND. f_star(:, l + 1)>1.e-10
    532       do ig = 1, ngrid
     532      DO ig = 1, ngrid
    533533        IF (activetmp(ig)) THEN
    534534          Zsat = .FALSE.
    … …  
    545545      ztemp(:) = zpspsk(:, l) * ztla(:, l)
    546546      CALL thermcell_qsat(ngrid, activetmp, pplev(:, l), ztemp, zqta(:, l), zqsatth(:, l))
    547       do ig = 1, ngrid
     547      DO ig = 1, ngrid
    548548        IF (activetmp(ig)) THEN
    549549          ! on ecrit de maniere conservative (sat ou non)
    … …  
    614614
    615615      nbpb = 0
    616       do ig = 1, ngrid
     616      DO ig = 1, ngrid
    617617        IF (zw2(ig, l + 1)>0. .AND. zw2(ig, l + 1)<1.e-10) THEN
    618618          !               stop 'On tombe sur le cas particulier de thermcell_dry'
    … …  
    658658
    659659    !on recalcule alim_star_tot
    660     do ig = 1, ngrid
     660    DO ig = 1, ngrid
    661661      alim_star_tot(ig) = 0.
    662662    enddo
    663     do ig = 1, ngrid
    664       do l = 1, lalim(ig) - 1
     663    DO ig = 1, ngrid
     664      DO l = 1, lalim(ig) - 1
    665665        alim_star_tot(ig) = alim_star_tot(ig) + alim_star(ig, l)
    666666      enddo
    … …  
    817817    ! Definition de l'alimentation
    818818    !-------------------------------------------------------------------------
    819     do l = 1, nlay - 1
    820       do ig = 1, ngrid
     819    DO l = 1, nlay - 1
     820      DO ig = 1, ngrid
    821821        IF (ztv(ig, l)> ztv(ig, l + 1) .AND. ztv(ig, 1)>=ztv(ig, l)) THEN
    822822          alim_star(ig, l) = MAX((ztv(ig, l) - ztv(ig, l + 1)), 0.)  &
    … …  
    827827      enddo
    828828    enddo
    829     do l = 1, nlay
    830       do ig = 1, ngrid
     829    DO l = 1, nlay
     830      DO ig = 1, ngrid
    831831        IF (alim_star_tot(ig) > 1.e-10) THEN
    832832          alim_star(ig, l) = alim_star(ig, l) / alim_star_tot(ig)
    … …  
    845845    ! dans une couche l>1
    846846    !------------------------------------------------------------------------------
    847     do ig = 1, ngrid
     847    DO ig = 1, ngrid
    848848      ! Le panache va prendre au debut les caracteristiques de l'air contenu
    849849      ! dans cette couche.
    … …  
    864864    !boucle de calcul de la vitesse verticale dans le thermique
    865865    !==============================================================================
    866     do l = 2, nlay - 1
     866    DO l = 2, nlay - 1
    867867      !==============================================================================
    868868
    … …  
    870870      ! On decide si le thermique est encore actif ou non
    871871      ! AFaire : Il faut sans doute ajouter entr_star a alim_star dans ce test
    872       do ig = 1, ngrid
     872      DO ig = 1, ngrid
    873873        active(ig) = active(ig) &
    874874                &                 .AND. zw2(ig, l)>1.e-10 &
    … …  
    891891      CALL thermcell_qsat(ngrid, active, pplev(:, l), ztemp, zqta(:, l - 1), zqsat(:))
    892892
    893       do ig = 1, ngrid
     893      DO ig = 1, ngrid
    894894        !       PRINT*,'active',active(ig),ig,l
    895895        IF(active(ig)) THEN
    … …  
    922922      !-------------------------------------------------
    923923
    924       do ig = 1, ngrid
     924      DO ig = 1, ngrid
    925925        IF (active(ig)) THEN
    926926          zw2m = max(0.5 * (w_est(ig, l) + w_est(ig, l + 1)), 0.1)
    … …  
    959959      !---------------------------------------------------------------------------
    960960      activetmp(:) = active(:) .AND. f_star(:, l + 1)>1.e-10
    961       do ig = 1, ngrid
     961      DO ig = 1, ngrid
    962962        IF (activetmp(ig)) THEN
    963963          Zsat = .FALSE.
    … …  
    975975      CALL thermcell_qsat(ngrid, activetmp, pplev(:, l), ztemp, zqta(:, l), zqsatth(:, l))
    976976
    977       do ig = 1, ngrid
     977      DO ig = 1, ngrid
    978978        IF (activetmp(ig)) THEN
    979979          ! on ecrit de maniere conservative (sat ou non)
    … …  
    10031003
    10041004      nbpb = 0
    1005       do ig = 1, ngrid
     1005      DO ig = 1, ngrid
    10061006        IF (zw2(ig, l + 1)>0. .AND. zw2(ig, l + 1)<1.e-10) THEN
    10071007          !               stop 'On tombe sur le cas particulier de thermcell_dry'
    … …  
    10461046
    10471047    !on recalcule alim_star_tot
    1048     do ig = 1, ngrid
     1048    DO ig = 1, ngrid
    10491049      alim_star_tot(ig) = 0.
    10501050    enddo
    1051     do ig = 1, ngrid
    1052       do l = 1, lalim(ig) - 1
     1051    DO ig = 1, ngrid
     1052      DO l = 1, lalim(ig) - 1
    10531053        alim_star_tot(ig) = alim_star_tot(ig) + alim_star(ig, l)
    10541054      enddo
Note: See TracChangeset for help on using the changeset viewer.