Changeset 3473

Timestamp:

Oct 24, 2024, 1:10:47 PM (6 weeks ago)

Author:

afalco

Message:

Pluto PCM: fixed glacier simple initializations bug with gfortran.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/spreadglacier_simple.F90 (modified) (9 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/spreadglacier_simple.F90

@@ -1 +1 @@
       subroutine spreadglacier_simple(ngrid,nq,pqsurf,timstep)
-          
+
       use geometry_mod, only: latitude, longitude, cell_area
       use radinc_h, only : naerkind
@@ -12 +12 @@
 !     -------
 !     Spreading the glacier in a circular crater using a speed of flow
-!  
+!
 !     Inputs
-!     ------ 
+!     ------
 !     ngrid                 Number of vertical columns
-!     pqsurf(ngrid,nq)      N2 ice tracer on surface kg/m2 
-!     
+!     pqsurf(ngrid,nq)      N2 ice tracer on surface kg/m2
+!
 !     Outputs
 !     -------
-!     pqsurf(ngrid,nq)      new value for N2 ice tracer on surface kg/m2 
-!     
+!     pqsurf(ngrid,nq)      new value for N2 ice tracer on surface kg/m2
+!
 !     Authors
-!     ------- 
+!     -------
 !     Tanguy Bertrand (2015)
-!     
+!
 !==================================================================
 
@@ -35 +35 @@
       INTEGER ngrid, nq, ig, ig2
       REAL :: pqsurf(ngrid,nq)
-      REAL timstep 
+      REAL timstep
 !-----------------------------------------------------------------------
 !     Local variables
@@ -48 +48 @@
       LOGICAL firstcall
       DATA firstcall/.true./
-      INTEGER indglac(ngrid) 
-      INTEGER nglac
+      INTEGER,SAVE,ALLOCATABLE :: indglac(:)
+      INTEGER, SAVE :: nglac
 
 !---------------- INPUT ------------------------------------------------
-! Defined for  Tombaugh crater : 
+! Defined for  Tombaugh crater :
 
       dist=300.  ! reference radius km
@@ -62 +62 @@
 !--------------- FIRSTCALL : define crater index -----------------------
       IF (firstcall) then
+         ALLOCATE(indglac(ngrid))
+         nglac=0
          indglac(:)=0
          DO ig=1,ngrid
@@ -68 +70 @@
                indglac(nglac)=ig
             ENDIF
-         ENDDO 
-         !write(*,*) 'inglac=',indglac
-         write(*,*) 'nglac=',nglac
- 
-         firstcall=.false. 
+         ENDDO
+         firstcall=.false.
 
       ENDIF
-
 !--------------- redistribution -------------------------------------
       DO ig=1,nglac
@@ -83 +81 @@
          ! and close to the reservoir
          DO ig2=1,nglac
-           IF ((pqsurf(indglac(ig2),igcm_n2).lt.pqsurf(indglac(ig),igcm_n2)*0.8)) THEN    
+           IF ((pqsurf(indglac(ig2),igcm_n2).lt.pqsurf(indglac(ig),igcm_n2)*0.8)) THEN
              !write(*,*) 'loop=',latitude(indglac(ig2))*180/3.14,longitude(indglac(ig2))*180/3.14
              !.and.(phisfi(indglac(ig2)).le.phisfi(indglac(ig)))) THEN
@@ -98 +96 @@
            ENDIF
          ENDDO
-        ENDIF 
+        ENDIF
       ENDDO
 
@@ -117 +115 @@
       c = 2 * atan2( sqrt(a), sqrt(1-a) )
       dist_pluto2 = rad * c
-      return 
+      return
       end function dist_pluto2