Ignore:
Timestamp:
Feb 19, 2025, 2:26:21 PM (5 months ago)
Author:
afalco
Message:

Pluto: larger size for array for High Resolution.
AF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.COMMON/libf/dyn3d/interp_horiz.F

    r1 r3630  
    33c
    44      subroutine interp_horiz (varo,varn,imo,jmo,imn,jmn,lm,
    5      &  rlonuo,rlatvo,rlonun,rlatvn) 
     5     &  rlonuo,rlatvo,rlonun,rlatvn)
    66
    77c===========================================================
     
    1414c===========================================================
    1515
    16       IMPLICIT NONE 
     16      IMPLICIT NONE
    1717
    1818c   Declarations:
     
    2121c  ARGUMENTS
    2222c  """""""""
    23        
     23
    2424       integer imo, jmo ! dimensions ancienne grille (input)
    2525       integer imn,jmn  ! dimensions nouvelle grille (input)
     
    2727       real rlonuo(imo+1)     !  Latitude et
    2828       real rlatvo(jmo)       !  longitude des
    29        real rlonun(imn+1)     !  bord des 
     29       real rlonun(imn+1)     !  bord des
    3030       real rlatvn(jmn)     !  cases "scalaires" (input)
    3131
     
    4242c    Info sur les ktotal intersection entre les cases new/old grille
    4343       integer kllm, k, ktotal
    44        parameter (kllm = 400*200*10)
     44       parameter (kllm = 1400*200*10)
    4545       integer iik(kllm), jjk(kllm),jk(kllm),ik(kllm)
    4646       real intersec(kllm)
     
    5353       data aire_ok /.true./
    5454
    55        
     55
    5656
    5757
     
    7575        end do
    7676       end do
    77       end do 
    78        
     77      end do
     78
    7979c Interpolation horizontale
    8080c -------------------------
     
    8383c
    8484      PRINT *, 'ktotal 1 = ', ktotal
    85      
     85
    8686      do k=1,ktotal
    8787        do l=1,lm
    88          varn(iik(k),jjk(k),l) = varn(iik(k),jjk(k),l) 
     88         varn(iik(k),jjk(k),l) = varn(iik(k),jjk(k),l)
    8989     &        + varo(ik(k), jk(k),l)*intersec(k)/airen(iik(k),jjk(k))
    9090        end do
     
    9999             totn = totn + varn(ii,1,l)
    100100             tots = tots + varn (ii,jmn+1,l)
    101            end do 
     101           end do
    102102           do ii =1, imn+1
    103103             varn(ii,1,l) = totn/REAL(imn+1)
    104104             varn(ii,jmn+1,l) = tots/REAL(imn+1)
    105            end do 
     105           end do
    106106       end do
    107            
     107
    108108
    109109c---------------------------------------------------------------
    110 c  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST 
     110c  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST
    111111!!       if (.not.(firsttest)) goto 99
    112112!!       firsttest = .false.
     
    116116!!           airetest(ii,jj) =0.
    117117!!         end do
    118 !!       end do 
     118!!       end do
    119119!!       PRINT *, 'ktotal = ', ktotal
    120120!!       PRINT *, 'jmn+1 =', jmn+1, 'imn+1', imn+1
    121 !! 
     121!!
    122122!!       do k=1,ktotal
    123 !!          airetest(iik(k),jjk(k))= airetest(iik(k),jjk(k)) +intersec(k) 
     123!!          airetest(iik(k),jjk(k))= airetest(iik(k),jjk(k)) +intersec(k)
    124124!!       end DO
    125 !! 
    126 !! 
     125!!
     126!!
    127127!!       PRINT *, 'fin boucle'
    128128!!       do jj =1 , jmn+1
Note: See TracChangeset for help on using the changeset viewer.