Changeset 3970 for trunk/LMDZ.PLUTO/libf

Timestamp:

Nov 21, 2025, 2:27:46 PM (3 weeks ago)

Author:

tbertrand

Message:

PLUTO PCM:
fixing a bug in newstart.F when setting phisfi using newstart
TB

File:

• trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/newstart.F (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/newstart.F

@@ -601 +601 @@
         ithfi(:,:)=inertiedat(:,:)
         ! copy topo
-        CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,phisold_newgrid,phisfi)
+        ! CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,phisold_newgrid,phisfi)
         ! Initialize n2frac
         n2fracfi(1:ngridmx)=1.
@@ -1458 +1458 @@
                 ! Loop where min
                 actualmin=1000.
-                do j=1,i
+                do j=1,i-1
                    if ( (array_angle(j).lt.angle+0.52).and.
      &                  (array_angle(j).gt.angle-0.52).and.
@@ -3150 +3150 @@
 
           ! get min max lon
-          print*, latfi*180/pi
-          print*, '***************************'
-          print*, '***************************'
-          print*, '***************************'
-          print*, '***************************'
-          print*, '***************************'
-          print*, lonfi*180/pi
-          print*, 'iip1=',iip1
-          do ig=2,ngridmx-1
+          ! do ig=2,ngridmx-1
+          do ig=2+iip1,ngridmx-1-iip1
              IF  (lonfi(ig)*180./pi.eq.-180.) THEN
-                 print*, lonfi(ig),lonfi(ig+iip1-2)
-                 phisfi(ig)=(phisfi(ig+1)+phisfi(ig+iip1))/2
+                 phisfi(ig)=(phisfi(ig-1)+phisfi(ig-iip1+2))/2.
              ENDIF
           enddo