Changeset 3698 for trunk/LMDZ.PLUTO/libf/phypluto/callcorrk_pluto_mod.F90

Timestamp:

Mar 26, 2025, 11:05:58 AM (5 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

File:

• trunk/LMDZ.PLUTO/libf/phypluto/callcorrk_pluto_mod.F90 (modified) (3 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