Changeset 3698 for trunk/LMDZ.PLUTO/libf

Timestamp:

Mar 26, 2025, 11:05:58 AM (10 months ago)

Author:

emillour

Message:

Pluto PCM:
Some fixes to enable runnnig with dynamico:

• add "strictboundcorrk" flag in callcorrk_pluto to enable running even if outside of kmatrix temperatures (when strictboundcorrk=.true.)
• add premature exiting of writediagsoil if not with lon-lat grid
• while at it, turned surfprop.F90 into a module
• in physiq, enforce the possibility to output subsurface-related field in most cases, not only when "fast=.true."
• adapt reference xml files: subsurface quantities need to be defined on a dedicated grid, otherwise XIOS will generate misleading garbage values. Updated files are put in "deftank/dynamico" for now.

EM

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callcorrk_pluto_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (3 diffs)
• surfprop.F90 (modified) (2 diffs)
• writediagsoil.F90 (modified) (2 diffs)

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

@@ -28 +28 @@
                               strobel,vmrch4_proffix,specOLR,vmrch4fix,&
                               haze_radproffix,callmufi
+      use callkeys_mod, only: strictboundcorrk
       use optcv_pluto_mod, only: optcv_pluto
       use optci_pluto_mod, only: optci_pluto
@@ -210 +211 @@
       real,save :: levdat(Nfine),vmrdat(Nfine)
       real :: vmrch4(ngrid,nlayer)              ! vmr ch4 from vmrch4_proffix
+      character(len=100) :: message
+      character(len=15),parameter :: subname="callcorrk_pluto"
 
 !=======================================================================
@@ -634 +637 @@
           if(tlevrad(k).lt.tgasmin)then
             print*,'Minimum temperature is outside the radiative'
-            print*,'transfer kmatrix bounds, exiting.'
+            print*,'transfer kmatrix bounds'
             print*,'t(',k,')=',tlevrad(k),' < ',tgasmin
-            ! call abort
+            if (strictboundcorrk) then
+              message="Minimum temperature outside of kmatrix bounds"
+              call abort_physic(subname,message,1)
+            else
+              print*,'***********************************************'
+              print*,'we allow model to continue with tlevrad<tgasmin'
+              print*,'  ... we assume we know what you are doing ... '
+              print*,'  ... but do not let this happen too often ... '
+              print*,'***********************************************'
+            endif
           elseif(tlevrad(k).gt.tgasmax)then
             print*,'Maximum temperature is outside the radiative'
-            print*,'transfer kmatrix bounds, exiting.'
+            print*,'transfer kmatrix bounds'
             print*,'t(',k,')=',tlevrad(k),' > ',tgasmax
-            ! call abort
+            if (strictboundcorrk) then
+              message="Maximum temperature outside of kmatrix bounds"
+              call abort_physic(subname,message,1)
+            else
+              print*,'***********************************************'
+              print*,'we allow model to continue with tlevrad>tgasmax'
+              print*,'  ... we assume we know what you are doing ... '
+              print*,'  ... but do not let this happen too often ... '
+              print*,'***********************************************'
+            endif
           endif
          enddo

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

@@ -70 +70 @@
       use callcorrk_mod, only: callcorrk
       use callcorrk_pluto_mod, only: callcorrk_pluto
+      use surfprop_mod, only: surfprop
       use vdifc_mod, only: vdifc
       use vdifc_pluto_mod, only: vdifc_pluto
@@ -1958 +1959 @@
  !     Info about Ls, declin...
       IF (fast) THEN
-         if (is_master) write(*,*),'Ls=',zls*180./pi,' dec=',declin*180./pi
-         if (is_master) write(*,*),'zday=',zday,' ps=',globave
+         if (is_master) write(*,*)'Ls=',zls*180./pi,' dec=',declin*180./pi
+         if (is_master) write(*,*)'zday=',zday,' ps=',globave
          IF (lastcall) then
-            if (is_master) write(*,*),'lastcall'
+            if (is_master) write(*,*)'lastcall'
          ENDIF
       ELSE
-         if (is_master) write(*,*),'Ls=',zls*180./pi,'decli=',declin*180./pi,'zday=',zday
+         if (is_master) write(*,*)'Ls=',zls*180./pi,'decli=',declin*180./pi,'zday=',zday
       ENDIF
 
@@ -2139 +2140 @@
       if (fast) then
          call write_output("fluxrad","fluxrad","W m-2",fluxrad)
-         call write_output("fluxgrd","fluxgrd","W m-2",fluxgrd)
          ! call write_output("dplanck","dplanck","W.s m-2 K-1",dplanck)
+      endif
+
+      if (callsoil) then
          ! "soil" variables
-         call write_output("capcal","capcal","W.s m-2 K-1",capcal)
+         call write_output("capcal","Surface Heat Capacity","W.s m-2 K-1",capcal)
          call write_output("tsoil","tsoil","K",tsoil)
          call write_output("therm_inertia","therm_inertia","S.I.",therm_inertia)

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

@@ +1 @@
+module surfprop_mod
+
+implicit none
+
+contains
+
       subroutine surfprop(ngrid,nq,pfra,qsurf,tsurface,       &
                      capcal,adjust,dist,fluold,ptimestep,zls, &
@@ -507 +513 @@
       end subroutine surfprop
 
+end module surfprop_mod

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

@@ -23 +23 @@
 use mod_phys_lmdz_para, only : is_mpi_root, is_master, gather
 use mod_grid_phy_lmdz, only : klon_glo, Grid1Dto2D_glo, &
-                              nbp_lon, nbp_lat
+                              nbp_lon, nbp_lat, &
+                              grid_type, unstructured
 
 implicit none
@@ -78 +79 @@
 real px2(ngrid)
 #endif
+
+! 0. This routine shoul only be used in lon-lat case
+if (grid_type==unstructured) then
+  return
+endif
 
 ! 1. Initialization step