Changeset 3243

Timestamp:

Feb 28, 2024, 6:55:38 PM (9 months ago)

Author:

afalco

Message:

Pluto PCM:
Get rid of real_area (should use cell_area probably)
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 (modified) (20 diffs)

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

@@ -16 +16 @@
       use write_field_phy, only: Writefield_phy
 !!
+      use ieee_arithmetic
       use ioipsl_getin_p_mod, only: getin_p
       use radinc_h, only : L_NSPECTI,L_NSPECTV,naerkind, corrkdir, banddir
@@ -29 +30 @@
       use comsaison_h, only: mu0, fract, dist_star, declin, right_ascen
       use comsoil_h, only: nsoilmx, layer, mlayer, inertiedat
-      use geometry_mod, only: latitude, longitude, cell_area, real_area
+      use geometry_mod, only: latitude, longitude, cell_area
       USE comgeomfi_h, only: totarea, totarea_planet
       USE tracer_h, only: noms, mmol, radius, rho_q, qext, &
@@ -55 +56 @@
                               n2cond,nearn2cond,noseason_day,conservn2, &
                               kbo,triton,paleo,paleoyears, &
-                              fast, carbox, methane, UseVdifcPlutold, &
+                              carbox, methane, UseVdifcPlutold, &
                               aerohaze,haze_proffix,source_haze, &
                               season, sedimentation,generic_condensation, &
@@ -286 +287 @@
       REAL flusurf(ngrid,nq)     ! flux cond/sub kg.m-2.s-1
       REAL flusurfold(ngrid,nq)  ! old flux cond/sub kg.m-2.s-1
+      REAL zplev(ngrid,nlayer+1),zplay(ngrid,nlayer)
 
 
@@ -842 +844 @@
 ! II. Compute radiative tendencies
 !---------------------------------
+!     Saving qsurf to compute paleo flux condensation/sublimation
+      DO iq=1, nq
+         DO ig=1,ngrid
+             IF (qsurf(ig,iq).lt.0.) then
+                qsurf(ig,iq)=0.
+             ENDIF
+             qsurf1(ig,iq)=qsurf(ig,iq)
+         ENDDO
+      ENDDO
+
 
       ! Compute local stellar zenith angles
@@ -1653 +1665 @@
       ps(1:ngrid) = pplev(1:ngrid,1) + pdpsrf(1:ngrid)*ptimestep
 
+      ! pressure density !pluto specific
+      IF (.not.fast) then !
+         do ig=1,ngrid
+            do l=1,nlayer
+               zplev(ig,l)=pplev(ig,l)/pplev(ig,1)*ps(ig)
+               zplay(ig,l)=pplay(ig,l)/pplev(ig,1)*ps(ig)
+               rho(ig,l) = zplay(ig,l)/(r*zt(ig,l))
+            enddo
+            zplev(ig,nlayer+1)=pplev(ig,nlayer+1)/pplev(ig,1)*ps(ig)
+         enddo
+      ENDIF
+
 
       ! Surface and soil temperature information
@@ -1795 +1819 @@
          ENDIF
        endif
- 
+
        if (carbox) then
          IF (fast) then
@@ -1814 +1838 @@
          ENDIF
        endif
- 
+
        zrho_haze(:,:)=0.
        zdqrho_photprec(:,:)=0.
@@ -1825 +1849 @@
          ENDDO
        ENDIF
- 
+
        IF (fasthaze) then
         DO ig=1,ngrid
@@ -1832 +1856 @@
         ENDDO
        ENDIF
- 
+
  !     Info about Ls, declin...
        IF (fast) THEN
@@ -1843 +1867 @@
          if (is_master) write(*,*),'Ls=',zls*180./pi,'decli=',declin*180./pi,'zday=',zday
        ENDIF
- 
+
        lecttsoil=0 ! default value for lecttsoil
        call getin_p("lecttsoil",lecttsoil)
@@ -1854 +1878 @@
           CLOSE(13)
        ENDIF
- 
- 
+
+
       if (is_master) print*,'--> Ls =',zls*180./pi
 
@@ -1893 +1917 @@
 
                ! special case if we sublime the entire reservoir
-                IF (-qsurfyear(ig,iq).gt.qsurf(ig,iq)) THEN
-                  masslost(iq)=masslost(iq)+real_area(ig)*   &
-                        (-qsurfyear(ig,iq)-qsurf(ig,iq))
-                  qsurfyear(ig,iq)=-qsurf(ig,iq)
-                ENDIF
-
-                IF (qsurfyear(ig,iq).gt.0.) THEN
-                  massacc(iq)=massacc(iq)+real_area(ig)*qsurfyear(ig,iq)
-                ENDIF
+               !! AF: TODO : fix following lines (real_area), using line below:
+            ! call planetwide_sumval((-qsurfyear(:,iq)-qsurf(:,iq))*cell_area(:),masslost)
+
+               !  IF (-qsurfyear(ig,iq).gt.qsurf(ig,iq)) THEN
+               !    masslost(iq)=masslost(iq)+real_area(ig)*   &
+               !          (-qsurfyear(ig,iq)-qsurf(ig,iq))
+               !    qsurfyear(ig,iq)=-qsurf(ig,iq)
+               !  ENDIF
+
+               !  IF (qsurfyear(ig,iq).gt.0.) THEN
+               !    massacc(iq)=massacc(iq)+real_area(ig)*qsurfyear(ig,iq)
+               !  ENDIF
+
 
                ENDDO
@@ -2068 +2096 @@
 
       !! Pluto outputs
-      call writediagfi(ngrid,"rice_ch4","ch4 ice mass mean radius","m",3,rice_ch4)
+      ! call writediagfi(ngrid,"rice_ch4","ch4 ice mass mean radius","m",3,rice_ch4)
       call writediagfi(ngrid,"dist_star","dist_star","AU",0,dist_star)
 
@@ -2075 +2103 @@
          !AF: TODO which outputs?
       else
+         if (check_physics_outputs) then
+            ! Check the validity of updated fields at the end of the physics step
+            call check_physics_fields("HERE physiq:", zt, zu, zv, pplev, zq)
+         endif
+
        call writediagfi(ngrid,"temp","temperature","K",3,zt)
        call writediagfi(ngrid,"teta","potential temperature","K",3,zh)
@@ -2229 +2262 @@
                call writediagfi(ngrid,"zrho_ch4","zrho_ch4","kg.m-3",3,zrho_ch4(:,:))
             endif
-    
+
              ! Tendancies
              call writediagfi(ngrid,"zdqch4cloud","ch4 cloud","T s-1",&
@@ -2241 +2274 @@
              call writediagfi(ngrid,"zdqadj_ch4","zdqadj ch4","",&
                             3,zdqadj(:,:,igcm_ch4_gas))
-             call writediagfi(ngrid,"zdqsed_ch4","zdqsed ch4","",&
-                            3,zdqsed(:,:,igcm_ch4_gas))
-             call writediagfi(ngrid,"zdqssed_ch4","zdqssed ch4","",&
-                            2,zdqssed(:,igcm_ch4_gas))
-            call writediagfi(ngrid,"zdtch4cloud","ch4 cloud","T s-1",&
+             if (sedimentation) then
+               call writediagfi(ngrid,"zdqsed_ch4","zdqsed ch4","",&
+                              3,zdqsed(:,:,igcm_ch4_gas))
+               call writediagfi(ngrid,"zdqssed_ch4","zdqssed ch4","",&
+                              2,zdqssed(:,igcm_ch4_gas))
+             endif
+             if (metcloud) then
+               call writediagfi(ngrid,"zdtch4cloud","ch4 cloud","T s-1",&
                            3,zdtch4cloud)
-    
-         endif
-    
+             endif
+
+         endif
+
          if (carbox) then
             call writediagfi(ngrid,'co_iceflux','co_iceflux',&
@@ -2258 +2295 @@
             endif
          endif
-    
+
          if (haze) then
    !         call writediagfi(ngrid,"zrho_haze","zrho_haze","kg.m-3",3,zrho_haze(:,:))
@@ -2270 +2307 @@
                      3,zdqhaze(:,:,igcm_prec_haze))
             if (igcm_haze.ne.0) then
-            call writediagfi(ngrid,"zdqhaze_haze","zdqhaze_haze","",&
-                     3,zdqhaze(:,:,igcm_haze))
-            call writediagfi(ngrid,"zdqssed_haze","zdqssed haze",&
-                  "kg/m2/s",2,zdqssed(:,igcm_haze))
+               call writediagfi(ngrid,"zdqhaze_haze","zdqhaze_haze","",&
+                        3,zdqhaze(:,:,igcm_haze))
+               if (sedimentation) then
+                  call writediagfi(ngrid,"zdqssed_haze","zdqssed haze",&
+                     "kg/m2/s",2,zdqssed(:,igcm_haze))
+               endif
             endif
             call writediagfi(ngrid,"zdqphot_ch4","zdqphot_ch4","",&
@@ -2282 +2321 @@
    !     &                   2,zdqhaze_col(:))
          endif
-   
+
          if (aerohaze) then
             call writediagfi(ngrid,"tau_col",&