Changeset 3630 for trunk/LMDZ.COMMON/libf/dyn3d

Timestamp:

Feb 19, 2025, 2:26:21 PM (5 months ago)

Author:

afalco

Message:

Pluto: larger size for array for High Resolution.
AF

File:

• trunk/LMDZ.COMMON/libf/dyn3d/interp_horiz.F (modified) (10 diffs)

trunk/LMDZ.COMMON/libf/dyn3d/interp_horiz.F

@@ -3 +3 @@
 c
       subroutine interp_horiz (varo,varn,imo,jmo,imn,jmn,lm,
-     &  rlonuo,rlatvo,rlonun,rlatvn)  
+     &  rlonuo,rlatvo,rlonun,rlatvn)
 
 c===========================================================
@@ -14 +14 @@
 c===========================================================
 
-      IMPLICIT NONE 
+      IMPLICIT NONE
 
 c   Declarations:
@@ -21 +21 @@
 c  ARGUMENTS
 c  """""""""
-        
+
        integer imo, jmo ! dimensions ancienne grille (input)
        integer imn,jmn  ! dimensions nouvelle grille (input)
@@ -27 +27 @@
        real rlonuo(imo+1)     !  Latitude et
        real rlatvo(jmo)       !  longitude des
-       real rlonun(imn+1)     !  bord des 
+       real rlonun(imn+1)     !  bord des
        real rlatvn(jmn)     !  cases "scalaires" (input)
 
@@ -42 +42 @@
 c    Info sur les ktotal intersection entre les cases new/old grille
        integer kllm, k, ktotal
-       parameter (kllm = 400*200*10)
+       parameter (kllm = 1400*200*10)
        integer iik(kllm), jjk(kllm),jk(kllm),ik(kllm)
        real intersec(kllm)
@@ -53 +53 @@
        data aire_ok /.true./
 
-       
+
 
 
@@ -75 +75 @@
         end do
        end do
-      end do 
-       
+      end do
+
 c Interpolation horizontale
 c -------------------------
@@ -83 +83 @@
 c
       PRINT *, 'ktotal 1 = ', ktotal
-     
+
       do k=1,ktotal
         do l=1,lm
-         varn(iik(k),jjk(k),l) = varn(iik(k),jjk(k),l) 
+         varn(iik(k),jjk(k),l) = varn(iik(k),jjk(k),l)
      &        + varo(ik(k), jk(k),l)*intersec(k)/airen(iik(k),jjk(k))
         end do
@@ -99 +99 @@
              totn = totn + varn(ii,1,l)
              tots = tots + varn (ii,jmn+1,l)
-           end do 
+           end do
            do ii =1, imn+1
              varn(ii,1,l) = totn/REAL(imn+1)
              varn(ii,jmn+1,l) = tots/REAL(imn+1)
-           end do 
+           end do
        end do
-           
+
 
 c---------------------------------------------------------------
-c  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST 
+c  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST
 !!       if (.not.(firsttest)) goto 99
 !!       firsttest = .false.
@@ -116 +116 @@
 !!           airetest(ii,jj) =0.
 !!         end do
-!!       end do 
+!!       end do
 !!       PRINT *, 'ktotal = ', ktotal
 !!       PRINT *, 'jmn+1 =', jmn+1, 'imn+1', imn+1
-!! 
+!!
 !!       do k=1,ktotal
-!!          airetest(iik(k),jjk(k))= airetest(iik(k),jjk(k)) +intersec(k) 
+!!          airetest(iik(k),jjk(k))= airetest(iik(k),jjk(k)) +intersec(k)
 !!       end DO
-!! 
-!! 
+!!
+!!
 !!       PRINT *, 'fin boucle'
 !!       do jj =1 , jmn+1