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/dyn3dmem/friction_loc.f90

    r5136 r5158  
    6767    IF (pole_sud) jje = jj_end
    6868
    69     do j = jjb, jje
    70       do i = 1, iip1
     69    DO j = jjb, jje
     70      DO i = 1, iip1
    7171        u2(i, j) = ucov(i, j, 1) * ucov(i, j, 1) * unscu2(i, j)
    7272      enddo
     
    7878    IF (pole_sud) jje = jj_end - 1
    7979
    80     do j = jjb, jje
    81       do i = 1, iip1
     80    DO j = jjb, jje
     81      DO i = 1, iip1
    8282        v2(i, j) = vcov(i, j, 1) * vcov(i, j, 1) * unscv2(i, j)
    8383      enddo
     
    9090    IF (pole_sud) jje = jj_end - 1
    9191
    92     do j = jjb, jje
    93       do i = 2, iip1
     92    DO j = jjb, jje
     93      DO i = 2, iip1
    9494        modv(i, j) = sqrt(0.5 * (u2(i - 1, j) + u2(i, j) + v2(i, j - 1) + v2(i, j)))
    9595      enddo
     
    103103      vpoln = 0.
    104104
    105       do i = 2, iip1
     105      DO i = 2, iip1
    106106        zco = cos(rlonv(i)) * (rlonu(i) - rlonu(i - 1))
    107107        zsi = sin(rlonv(i)) * (rlonu(i) - rlonu(i - 1))
     
    111111      enddo
    112112      vpn = sqrt(upoln * upoln + vpoln * vpoln) / pi
    113       do i = 1, iip1
     113      DO i = 1, iip1
    114114        ! modv(i,1)=vpn
    115115        modv(i, 1) = modv(i, 2)
     
    121121      upols = 0.
    122122      vpols = 0.
    123       do i = 2, iip1
     123      DO i = 2, iip1
    124124        zco = cos(rlonv(i)) * (rlonu(i) - rlonu(i - 1))
    125125        zsi = sin(rlonv(i)) * (rlonu(i) - rlonu(i - 1))
     
    129129      enddo
    130130      vps = sqrt(upols * upols + vpols * vpols) / pi
    131       do i = 1, iip1
     131      DO i = 1, iip1
    132132        ! modv(i,jjp1)=vps
    133133        modv(i, jjp1) = modv(i, jjm)
     
    143143    IF (pole_sud) jje = jj_end - 1
    144144
    145     do j = jjb, jje
    146       do i = 1, iim
     145    DO j = jjb, jje
     146      DO i = 1, iim
    147147        ucov(i, j, 1) = ucov(i, j, 1) &
    148148                - cfric * pdt * 0.5 * (modv(i + 1, j) + modv(i, j)) * ucov(i, j, 1)
     
    155155    IF (pole_sud) jje = jj_end - 1
    156156
    157     do j = jjb, jje
    158       do i = 1, iip1
     157    DO j = jjb, jje
     158      DO i = 1, iip1
    159159        vcov(i, j, 1) = vcov(i, j, 1) &
    160160                - cfric * pdt * 0.5 * (modv(i, j + 1) + modv(i, j)) * vcov(i, j, 1)
     
    173173
    174174    !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    175     do l = 1, llm
     175    DO l = 1, llm
    176176      ucov(1:iip1, jjb:jje, l) = ucov(1:iip1, jjb:jje, l) * &
    177177              (1. - pdt * kfrict(l))
     
    185185
    186186    !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    187     do l = 1, llm
     187    DO l = 1, llm
    188188      vcov(1:iip1, jjb:jje, l) = vcov(1:iip1, jjb:jje, l) * &
    189189              (1. - pdt * kfrict(l))
Note: See TracChangeset for help on using the changeset viewer.