Ignore:
Timestamp:
Jul 19, 2024, 6:40:44 PM (13 months ago)
Author:
Laurent Fairhead
Message:

Reverting to r4065. Updating fortran standard broke too much stuff. Will do it by smaller chunks
AB, LF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/filtrez/inifgn.F

    r5079 r5084  
    2828      pi = 2.* ASIN(1.)
    2929C
    30       DO i=1,iim
     30      DO 5 i=1,iim
    3131       dlonu(i)=  xprimu( i )
    3232       dlonv(i)=  xprimv( i )
    33       END DO
     33   5  CONTINUE
    3434
    35       DO i=1,iim
     35      DO 12 i=1,iim
    3636      sddv(i)   = SQRT(dlonv(i))
    3737      sddu(i)   = SQRT(dlonu(i))
    3838      unsddu(i) = 1./sddu(i)
    3939      unsddv(i) = 1./sddv(i)
    40       END DO
     40  12  CONTINUE
    4141C
    42       DO j=1,iim
    43       DO i=1,iim
     42      DO 17 j=1,iim
     43      DO 17 i=1,iim
    4444      vec(i,j)     = 0.
    4545      vec1(i,j)    = 0.
    4646      eignfnv(i,j) = 0.
    4747      eignfnu(i,j) = 0.
    48       END DO
    49       END DO
     48  17  CONTINUE
    5049c
    5150c
    5251      eignfnv(1,1)    = -1.
    5352      eignfnv(iim,1)  =  1.
    54       DO i=1,imm1
     53      DO 20 i=1,imm1
    5554      eignfnv(i+1,i+1)= -1.
    5655      eignfnv(i,i+1)  =  1.
    57       END DO
    58       DO j=1,iim
    59       DO i=1,iim
     56  20  CONTINUE
     57      DO 25 j=1,iim
     58      DO 25 i=1,iim
    6059      eignfnv(i,j) = eignfnv(i,j)/(sddu(i)*sddv(j))
    61       END DO
    62       END DO
    63       DO j=1,iim
    64       DO i=1,iim
     60  25  CONTINUE
     61      DO 30 j=1,iim
     62      DO 30 i=1,iim
    6563      eignfnu(i,j) = -eignfnv(j,i)
    66       END DO
    67       END DO
     64  30  CONTINUE
    6865c
    6966#ifdef CRAY
Note: See TracChangeset for help on using the changeset viewer.