Changeset 3813

Timestamp:

Jun 24, 2025, 5:29:54 PM (5 days ago)

Author:

afalco

Message:

Pluto: non-oro gravity waves: imported fix from LMDZ terrestrial model for poles. Allowed to output tendencies.
AF

Location:

trunk/LMDZ.PLUTO/libf

Files:

• dynphy_lonlat/phypluto/grid_noro1.F (modified) (3 diffs)
• phypluto/physiq_mod.F90 (modified) (1 diff)

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

@@ -308 +308 @@
       zmaxthe=0.
       zminthe=0.
-c     print 100,' '
-c100  format(1X,A1,'II JJ',4X,'H',8X,'SD',8X,'SI',3X,'GA',3X,'TH') 
-       print*,'OK1'
-       DO ii = 1, imar
-       DO jj = 1, jmar
-c       print*,'ok0'
-         IF (weight(ii,jj) .NE. 0.0) THEN
-c  Orography moyenne:
-c         print*,'ok1'
-           zmea (ii,jj)=zmea (ii,jj)/weight(ii,jj)
-           zxtzx(ii,jj)=zxtzx(ii,jj)/weight(ii,jj)
-           zytzy(ii,jj)=zytzy(ii,jj)/weight(ii,jj)
-           zxtzy(ii,jj)=zxtzy(ii,jj)/weight(ii,jj)
-           ztz(ii,jj)  =ztz(ii,jj)/weight(ii,jj)
-c         print*,'ok2'
-c  Deviation standard:
-           zstd(ii,jj)=sqrt(amax1(0.,ztz(ii,jj)-zmea(ii,jj)**2))
+      DO ii = 1, imar
+      DO jj = 1, jmar
+        IF (weight(ii,jj) .NE. 0.0) THEN
+          zmea (ii,jj)=zmea (ii,jj)/weight(ii,jj)
+          zxtzx(ii,jj)=zxtzx(ii,jj)/weight(ii,jj)
+          zytzy(ii,jj)=zytzy(ii,jj)/weight(ii,jj)
+          zxtzy(ii,jj)=zxtzy(ii,jj)/weight(ii,jj)
+          ztz(ii,jj)  =ztz(ii,jj)/weight(ii,jj)
+          zstd(ii,jj)=sqrt(amax1(0.,ztz(ii,jj)-zmea(ii,jj)**2))
+        ENDIF
+      ENDDO
+      ENDDO
+      !--- CORRECT VALUES OF HORIZONTAL SLOPE NEAR THE POLES:
+      zxtzx(:,   1)=zxtzx(:,     2)
+      zxtzx(:,jmar)=zxtzx(:,jmar-1)
+      zxtzy(:,   1)=zxtzy(:,     2)
+      zxtzy(:,jmar)=zxtzy(:,jmar-1)
+      zytzy(:,   1)=zytzy(:,     2)
+      zytzy(:,jmar)=zytzy(:,jmar-1)
+
+      DO ii = 1, imar
+      DO jj = 1, jmar
 c  Coefficients K, L et M:
            xk=(zxtzx(ii,jj)+zytzy(ii,jj))/2.
@@ -335 +340 @@
            if(xq.le.xw) xq=xw
            if(abs(xm).le.xw) xm=xw*sign(1.,xm)
-c          print*,'ok3'
 c pente: 
            zsig(ii,jj)=sqrt(xq)
@@ -344 +348 @@
            zthe(ii,jj)=57.29577951*atan2(xm,xl)/2.
 
-c          print 101,ii,jj,
-c    *           zmea(ii,jj),zstd(ii,jj),zsig(ii,jj),zgam(ii,jj),
-c    *           zthe(ii,jj)
-c101  format(1x,2(1x,i2),2(1x,f7.1),1x,f7.4,2x,f4.2,1x,f5.1)     
-c          print*,'ok4'
-         ELSE
-c           PRINT*, 'probleme,ii,jj=', ii,jj
-c          print*,'ok1b'
-         ENDIF
       zmaxmea=amax1(zmea(ii,jj),zmaxmea)
-c         print*,'oka'
       zmaxstd=amax1(zstd(ii,jj),zmaxstd)
-c         print*,'okb'
       zmaxsig=amax1(zsig(ii,jj),zmaxsig)
-c         print*,'okc'
       zmaxgam=amax1(zgam(ii,jj),zmaxgam)
-c         print*,'okd'
       zmaxthe=amax1(zthe(ii,jj),zmaxthe)
-c         print*,'oke'
       zminthe=amin1(zthe(ii,jj),zminthe)
-c      print*,'ok5'
        ENDDO
        ENDDO

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

@@ -2253 +2253 @@
       call write_output('dtaui_17','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,17))
 
+      if (calllott) then
+        call write_output("zdugw","tendancy U lott","m s-2",zdugw)
+        call write_output("zdvgw","tendancy V lott","m s-2",zdvgw)
+        call write_output("zdtgw","tendancy T lott","K s-1",zdtgw)
+      endif
+      
       ! Output tracers.
       if (tracer) then